in the course of the last few years exclusive photo- and electroproduction of neutral vector mesons ( @xmath5 , @xmath6 , @xmath1 , @xmath7 ) at intermediate and high energies has attained increased theoretical interest , which has been stimulated by corresponding experimental efforts at cebaf  @xcite and hera  @xcite .
the aim of the hera experiments is a more fundamental understanding of pomeron phenomenology in terms of qcd .
vector - meson dominance in combination with pomeron exchange proves to be very successful in describing diffractive photo- and electroproduction of neutral vector mesons , provided the energy is the only large scale in the process  @xcite . in their approach donnachie and landshoff
@xcite assume that the pomeron couples directly to on - shell quarks in the vector meson and the proton , respectively .
a qcd - inspired version of pomeron exchange , in which the pomeron is replaced by two non- perturbative ( abelian ) gluons , goes back to the same authors  @xcite and is usually termed a soft pomeron . in the limit @xmath8
the soft - pomeron exchange becomes very similar to the usual pomeron exchange .
this can be used to interpret the form factor at the quark - pomeron vertex  @xcite .
however , if in addition to the energy another scale , like the photon virtuality or the mass of the vector meson , becomes large perturbative qcd effects tend to enter the game . for large photon virtuality ( and small values of mandelstam @xmath3 ) perturbative qcd becomes applicable due to a factorization theorem which has been proved in ref .
the three building blocks of the factorization formula are the hard photon - parton scattering amplitude , the distribution amplitude of the vector meson and finally a skewed gluon distribution of the proton . as a link between hard inclusive and exclusive reactions which generalizes the ordinary parton distributions ,
skewed parton distributions have recently attracted a great amount of interest ( see , e.g. , ref .
@xcite and references therein ) .
electroproduction with highly virtual photons thus offers a possibility to study such skewed parton distributions . within this perturbative formalism for vector - meson electroproduction with highly virtual photons
the two - gluon aspect of the pomeron emerges automatically in a form which is related to the skewed gluon distribution in the proton .
one of the aims of photoproduction experiments at cebaf is , on the other hand , to understand the transition from non - perturbative to perturbative production mechanisms . in the case of photoproduction
the transverse momentum transfer @xmath9 has to be large in order that perturbative qcd becomes applicable .
the pertinent factorization formula , which in general holds for arbitrary exclusive hadronic processes at large transverse momenta , is the result of a leading - order perturbative analysis .
it represents a hadronic amplitude as a convolution integral of a hard - scattering amplitude with distribution amplitudes ( das )  @xcite .
the process - dependent hard - scattering amplitude is perturbatively calculable and represents the scattering of hadronic constituents in collinear approximation .
the process- independent das contain the non - perturbative bound - state dynamics of the hadronic constituents .
there is no doubt that this hard - scattering approach ( hsa ) gives the correct description of exclusive hadronic processes in the limit of asymptotically large momentum transfer .
but there is increasing evidence that the pure perturbative contribution is not the dominant one at experimentally accessible momentum transfers . a reasonable reproduction of the existing exclusive large-@xmath3 data is only achieved with very endpoint- concentrated hadron das , like those proposed by chernyak and zhitnitsky on the basis of qcd sum rules  @xcite .
such das , however , occur to be problematic for a perturbative calculation since they enhance contributions from the soft end - point regions in the convolution integral  @xcite .
some improvement with regard to the consistency of the perturbative calculation can be achieved by applying a modified factorization scheme in which transverse - momentum effects and sudakov suppressions are included . but
this leads to a substantial reduction of the perturbative contribution  @xcite . for the pion
an end - point concentrated da meanwhile seems even to be ruled out by the recent cleo data  @xcite for the @xmath10-transition form factor , since it provides a much too large perturbative result . on the other hand ,
it has been shown recently for electromagnetic nucleon form factors  @xcite and for wide - angle compton scattering  @xcite off protons that the existing data can be saturated by soft overlap contributions .
furthermore , these contributions can be modeled by a nucleon light- cone wave function which gives rise to a distribution amplitude close to the asymptotic one .
this , however , implies that the perturbative contribution becomes small .
the approach which we are going to apply for hard exclusive photoproduction of @xmath1 mesons is still based on the hsa but includes also a certain amount of soft physics which is modeled by means of diquarks .
diquarks may be considered as an effective way to account for quark - quark correlations in baryon wave functions .
part of the soft overlap contributions are then absorbed into phenomenologically parameterized diquark form factors which occur at the gauge - boson diquark vertices .
this hsa - based diquark model has already been applied successfully to other photon - induced hadronic reactions like magnetic and electric baryon form factors in the space- @xcite and time - like region @xcite , real and virtual compton scattering @xcite , two - photon annihilation into baryon - antibaryon @xcite , and photoproduction of @xmath11-mesons  @xcite .
further applications of the diquark model include the charmonium decay @xmath12 @xcite and the calculation of landshoff contributions in elastic proton - proton scattering @xcite .
the purpose of the present investigation is to extend this foregoing work to @xmath13 .
the reason why we concentrate on the @xmath1 channel is that it is cleaner and simpler than photoproduction of @xmath5 or @xmath6 mesons . within the hsa the valence fock state of the @xmath1 , i.e. the @xmath14-@xmath15 state
, can only be produced via the exchange of two gluons between the @xmath1 and the proton .
quark exchange could only happen via the knockout of an @xmath14-@xmath15 pair from the strange - quark sea in the nucleon which is , however , power suppressed within the hsa .    in the following section we start with a short outline of the diquark model , show how feynman
diagrams are properly grouped according to their propagator singularities ( sec .
[ sec : elamp ] ) , sketch how we deal with the propagator singularities numerically ( sec .  [ sec : propnum ] ) , and present our treatment of hadron - mass effects ( sec .
[ sec : mass ] ) .
section  [ sec : helamp ] summarizes how photoproduction of vector mesons is described in terms of helicity amplitudes .
furthermore , the fixed - angle scaling behavior of these helicity amplitudes is discussed and related to the treatment of mass effects .
analytical expressions for the hard part of the hadron - helicity - conserving amplitudes are also given in this section .
our predictions for photoproduction observables along with predictions from pomeron and soft - pomeron exchange are presented in sec .
[ sec : result ] .
at the end of this section we comment on the prospects of extending our calculation to photoproduction of @xmath7 mesons .
finally , a summary and concluding remarks can be found in sec .
[ sec : sum ] .
as already mentioned before , within the hsa an exclusive scattering amplitude at large momentum transfer is expressed as a convolution integral of a hard - scattering amplitude with hadron distribution amplitudes  @xcite .
for the photoproduction reaction @xmath16 , which we are interested in , this integral takes on the particular form @xmath17 the distribution amplitudes @xmath18 are probability amplitudes for finding the valence fock state in the hadron @xmath19 with the constituents carrying certain fractions of the momentum of their parent hadron and being collinear up to a factorization scale @xmath20 . in our model the valence fock state of an ordinary baryon
is assumed to consist of a quark ( @xmath21 ) and a diquark ( @xmath22 ) .
we fix our notation in such a way that the momentum fraction appearing in the argument of @xmath23 is carried by the quark  with the momentum fraction of the other constituent ( either diquark or antiquark ) it sums up to 1 ( cf . fig .
[ kinem ] ) . for our actual calculations the ( logarithmic ) @xmath20 dependence of the das
is neglected throughout since it is of minor importance in the restricted energy and momentum - transfer range we are interested in .
the hard scattering amplitude @xmath24 is calculated perturbatively .
intrinsic transverse momenta of the constituents are thereby neglected , i.e. , one assumes that the momenta of the constituents are collinear to those of their parent hadron .
@xmath24 consists in our particular case of all possible tree diagrams contributing to the elementary scattering process @xmath25 .
eight of the , altogether , sixteen diagrams which contribute to photoproduction of ( heavy ) quarkonia are depicted in fig .
the subscript @xmath26 of @xmath27 represents a set of possible photon , proton and vector - meson helicities . for our purposes it appears to be more convenient to express the analytical results in terms of massless mandelstam variables @xmath28 , @xmath29 , and @xmath30 than in terms of the usual ( massive ) ones .
a detailed account of our treatment of hadron - mass effects is given in sec .
[ sec : mass ] below .    the diquark model , as applied in refs .
@xcite-@xcite , comprises scalar ( @xmath31 ) as well as axial - vector ( @xmath32 ) diquarks .
the dynamics of @xmath31 and @xmath32 diquarks is determined by their coupling to photons and gluons .
the expressions for the photon - diquark and gluon - diquark vertices correspond ( almost ) to the most general form for the ( parity and time - reversal invariant ) coupling of a spin-1 gauge boson to a spin-0 or a spin-1 particle , respectively  @xcite .
thereby @xmath32 diquarks are allowed to possess an anomalous ( chromo)magnetic moment @xmath33 . in applications of the model feynman graphs
are calculated first with the corresponding feynman rules for point - like diquarks . as a second step , feynman graphs with @xmath34 gauge bosons coupling to the diquark
are combined in a group which we call n - point contribution .
finally the composite nature of diquarks is taken into account by multiplying each n - point contribution with phenomenological vertex functions , the diquark form factors . the particular choice @xmath35 for 3-point contributions and @xmath36 for n - point contributions ( @xmath37 ) ensures that in the limit @xmath38 the scaling behavior of the diquark model goes over in that of the pure quark hsp .
the factor @xmath39 ( @xmath40 for @xmath41 ) provides the correct powers of @xmath42 for asymptotically large @xmath43 . for the running coupling @xmath44 the one - loop result @xmath45
is used with @xmath46 .
in addition @xmath47 is restricted to be smaller than @xmath48 .
@xmath49 and @xmath50 are strength parameters which allow for the possibility of diquark excitation and break - up in intermediate states where diquarks can be far off - shell .
having sketched how diquarks are treated perturbatively , a few words about the choice of the das ( which incorporate the bound - state dynamics of the hadronic constituents ) and the spin - flavor wave functions of the proton and the @xmath1 are in order .
if one assumes zero relative orbital angular momentum between quark and diquark and takes advantage of the collinear configuration ( @xmath51 and @xmath52 ) the valence fock - state wave function of an octet baryon b may be written as @xmath53 the two terms in eq .  ( [ wfcov ] ) represent states consisting of a quark either a scalar or a vector diquark .
the pleasant feature of the covariant wave - function representation , eq .
( [ wfcov ] ) , is that it contains , besides @xmath54 and @xmath55 ( the lorentz index of the vector - diquark polarization vector ) , only baryonic quantities ( momentum @xmath56 , helicity @xmath57 , baryon mass @xmath58 ) .
assuming an su(6)-like spin - flavor dependence , the flavor functions @xmath59 for a proton take on the form @xmath60 } \ , , \;\;\ ; \chi_v^p = \phantom{- } [ u v_{\{u , d\ } } -\sqrt{2 } d v_{\{u , u\ } } ] / \sqrt{3 } \ , .
\label{flavp}\ ] ] in eq .
( [ wfcov ] ) the color part of the quark - diquark wave function has been neglected .
since the quark - diquark state has to be a color singlet the diquark has to be in a color antitriplet state ( like an antiquark ) thus the quark - diquark color wave function is simply @xmath61 .    like the baryon wave function , eq .
( [ wfcov ] ) , also the @xmath62-@xmath63 wave function of a vector meson @xmath64 can be represented in a covariant way @xmath65 with @xmath66 denoting the polarization vector of the vector meson .
the flavor function of the @xmath1 is @xmath67 .    in refs .
@xcite-@xcite a quark - diquark da of the form ( @xmath68 for @xmath31 diquarks ) @xmath69 \ , ,   \label{dap}\ ] ] turned out to be quite appropriate for octet baryons @xmath70 .
the origin of the phenomenological baryon da , eq .
( [ dap ] ) , is a non - relativistic harmonic - oscillator wave function ( cf .
therefore the masses appearing in the exponential have to be considered as constituent masses ( 330 mev for light quarks , 580 mev for light diquarks , strange quarks are 150 mev heavier than light quarks ) .
the oscillator parameter @xmath71 gev@xmath72 is chosen in such a way that the full wave function gives rise to a value of @xmath73 mev for the mean intrinsic transverse momentum of quarks inside a nucleon .
the exponential does not only provide a flavor dependence , but also suppresses the end - point regions @xmath74 in the convolution integral , eq .
( [ convol ] ) .
one thus avoids to pick up sizable contributions from the soft end - point regions and results become less dependent on details of the baryon da . in the actual data fitting the exponential plays only a minor role .
therefore , the constituent masses and the oscillator parameter have not been considered as free parameters but have been fixed in advance .
we stress that the constituent masses do not appear in the hard scattering amplitudes ( cf .
[ sec : mass ] below ) .
a model for the @xmath1 meson da which fulfills qcd sum - rule constraints has been proposed by benayoun and chernyak  @xcite .
since we also want to investigate to which extent our calculations for @xmath1 production may be generalized to @xmath7 production we modify this da by attaching an exponential factor which provides a flavor dependence .
the effect of this factor is negligible for @xmath1s , but it becomes crucial for @xmath7s .
the form of the exponential factor is inspired by a simple relativistic treatment of heavy mesons  @xcite and provides a flavor dependence in accordance with heavy - quark effective theory .
for longitudinally polarized @xmath1s and @xmath7s the da reads ( @xmath75 ) @xmath76\ , .
\label{dabc}\ ] ] the da for transversally polarized @xmath1s and @xmath7s contains an additional polynomial in @xmath77 which enhances the maximum at @xmath78 and makes it narrower ( cf .
@xcite ) @xmath79 \nonumber \\ & & \times ( 5225.2\ , z_1 ^ 6 ( 1 - z_1)^6 + 0.39\ , z_1 ( 1 - z_1 ) ) \ ,
. \label{dabcl}\end{aligned}\ ] ] for the oscillator parameter @xmath80 we take a value of @xmath81  gev@xmath82 in accordance with estimates for the radius of the @xmath7 meson  @xcite .    the usual normalization condition , @xmath83 , fixes the constants @xmath84 and @xmath85 in eqs .
( [ dap ] ) , ( [ dabc ] ) , and ( [ dabcl ] ) .
the quantity @xmath86 showing up in eq .
( [ wfcovv ] ) is related to the experimentally determinable decay constant of the meson @xmath64 . in case of the @xmath1
we get from the radiative decay width of the @xmath87 decay @xmath88  mev .
the analogous constants @xmath89 and @xmath90 for the @xmath21-@xmath22 das of baryon @xmath70 are free parameters of the diquark model .
they are , in principle , determined by the probability to find the @xmath21-@xmath22 state ( @xmath91 ) in the baryon @xmath70 and by the transverse - momentum dependence of the corresponding wave function .
our present study is performed with the set of diquark - model parameters found in ref .
@xcite by fitting elastic electron - nucleon scattering data .
the numerical values of the parameters are : @xmath92 for further details of the diquark model we also refer to the paper of jakob et al .
@xcite .
the introduction of diquark form factors already requires a decomposition of the hard - scattering amplitude @xmath93 into 3- and 4-point contributions @xmath94 } } \widehat{t}_{\{\lambda\}}^{(s)}(x_1,y_1,z_1;\hat{s},\hat{t } ) \nonumber \\   & + & e_s \widehat{t}_{\{\lambda\}}^{(q , v)}(x_1,y_1,z_1;\hat{s},\hat{t } )   \ , .
\label{tn}\end{aligned}\ ] ] @xmath95 and @xmath96}}$ ] are the electric charges of the quark @xmath21 and the s diquark @xmath97}$ ] in units of @xmath98 , respectively .
the superscripts indicate whether the photon couples to the heavy quark @xmath62 in the meson , or the light quark @xmath21 in the proton with the diquark @xmath99 acting as spectator ( 3-point contributions ) , or whether the photon couples to the @xmath31 or @xmath32 diquark in the proton ( 4-point contributions ) .
note that the only @xmath32- diquark contribution is @xmath100 .
the reason is that the charge - flavor factor in front of @xmath101 vanishes and @xmath102 represents already a mass correction of second order , which is neglected in our calculation .
for the numerical treatment of the convolution integral , eq .
( [ convol ] ) , it is advantageous to further split the @xmath103-point contributions into two parts which differ by their propagator singularities : @xmath104   f_d^{(3)}(-x_2 y_2 \hat{t } ) \ , , \nonumber \\
\widehat{t}_{\{\lambda\}}^{(q , d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) & = & \left [ \frac{f_{\{\lambda\}}^{(q , d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) }             { ( \hat{g}_1 ^ 2 +
i \epsilon ) ( \hat{q}_2 ^ 2 +
i \epsilon^\prime ) }    +    \frac{g_{\{\lambda\}}^{(q , d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) }              { ( \hat{g}_1 ^ 2 + i \epsilon ) ( \hat{q}_3 ^ 2 +
i \epsilon^\prime ) } \right ]   f_d^{(3)}(-x_2 y_2 \hat{t } ) \ , , \nonumber \\
\widehat{t}_{\{\lambda\}}^{(d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) & = & \left [ \frac{f_{\{\lambda\}}^{(d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) }             { ( \hat{g}_2 ^ 2 + i \epsilon ) ( \hat{q}_5 ^ 2 +
i \epsilon^\prime ) }    +    \frac{g_{\{\lambda\}}^{(d)}(x_1,y_1,z_1;\hat{s},\hat{t } ) }              { ( \hat{g}_2 ^ 2 +
i \epsilon ) ( \hat{q}_4 ^ 2 +
i \epsilon^\prime ) } \right ]   f_d^{(4)}(-x_2 y_2 \hat{t } ) \ , . \label{tfg}\end{aligned}\ ] ]
the functions @xmath105 and @xmath106 correspond to gauge - invariant subgroups of the 16 feynman diagrams which contribute to the hard - scattering amplitude @xmath93 ( cf .
[ feyn ] ) .
the diagrams entering the @xmath106s are related to those entering the @xmath105s by interchange of the two gluons .
the @xmath107 and @xmath108 denote just those quark and gluon propagators which can go on- shell when integrating over @xmath54 , @xmath109 , and @xmath77 .
explicitly , the propagator denominators read : @xmath110 as already indicated in eq .
( [ tfg ] ) propagator singularities are treated by means of the usual @xmath111 prescription .
the resulting imaginary part in the amplitudes is considered as a non - trivial prediction of perturbative qcd which is unaffected by long - distance effects  @xcite .
as we will see below , it gives rise to non - vanishing predictions for polarization observables which require the flip of a hadronic helicity .
for analytical expressions of the functions @xmath105 and @xmath106 and their symmetry properties we refer to sec .
[ sec : helamp ] .      the numerical difficulties in performing the convolution integral , eq .
( [ convol ] ) , are mainly caused by the occurrence of propagator singularities in the range of the @xmath112-integration .
these give rise to principal value integrals @xmath113 as explained in some detail in ref .
@xcite such integrals can still be treated by means of a rather fast and stable fixed - point gaussian quadrature after carefully separating the principal value parts and doing the corresponding integrations analytically .
here we only want to give a short account of our integration procedure .
we first perform the @xmath54-integration for which the various contributions to the hard scattering amplitude ( cf .
( [ tfg ] ) ) lead to integrals of the general form @xmath114 in order to simplify notations we have neglected helicity labels and the dependence on the mandelstam variables @xmath28 and @xmath29 . furthermore , the distribution amplitudes have been absorbed into the function @xmath115 . by applying a partial fractioning , exploiting eq .
( [ pv ] ) , and separating the principal value part eq .
( [ convol2a ] ) may be rewritten to give @xmath116 @xmath117 and @xmath118 are the zeroes of the propagator denominators @xmath119 and @xmath120 , respectively , and @xmath121 stands for either @xmath122 or @xmath123 , depending for which of the two the limit @xmath124 has been taken .
as long as the two zeroes @xmath125 and @xmath126 ( which still depend on @xmath109 and @xmath127 ) do not coincide in the @xmath128-integration region  this is the case for the contributions of @xmath129 , @xmath130 , @xmath131 , and @xmath132
the @xmath112-integration can be performed very efficiently with the help of eq .
( [ convol3a ] ) . after analytic integration of the two principal - value integrals in eq .
( [ convol3a ] ) ( which is a simple matter ) all the integrands of the remaining integrations are regular functions of the integration variables and can be treated by fixed - point gaussian quadrature .    for @xmath133 and @xmath134
, however , the situation becomes a little bit more complicated . after having performed the @xmath54-integration in eq .
( [ convol3a ] ) one observes that the two propagator singularities @xmath125 and @xmath126 coincide along a trajectory in the @xmath128-integration region .
this leads to a quadratic zero @xmath135 of @xmath136 .
a closer inspection of eq .  ( [ convol3a ] ) reveals that this quadratic zero in the denominator of @xmath137 does not cause any further problem when performing the @xmath128-integration for the real part of @xmath137 .
the reason is that it is completely compensated by a corresponding quadratic zero in the numerator of @xmath138 . in @xmath139 , however , the zero of @xmath140 is only partly compensated by the numerator and one still encounters a single pole . at the pole
@xmath139 becomes proportional to @xmath141 = \{1/(y_1- y_1^{(k_1,k_2 ) } + i \hat{\epsilon } ) + 1/(y_1-y_1^{(k_1,k_2)}- i \hat{\epsilon})\}/2 $ ] . as a consequence @xmath142 becomes a pure principal - value integral ( cf .
( [ pv ] ) ) @xmath143 which can be treated numerically analogous to the principal - value integral in @xmath54 .
one only has to separate the principal - value part and do the corresponding integration analytically . for this purpose the analytical expression of the residue of @xmath144 at the position of the pole @xmath145
is needed .
the remaining integrations are then again amenable to gaussian quadrature .
thus it is possible , by carefully separating the singular contributions , exploiting @xmath146 functions , rewriting principal - value integrals as ordinary integrals plus analytically solvable principal - value integrals , to do all the numerical integrations by means of fixed - point gaussian quadrature . whereas the authors of ref .
@xcite , who have used a different integration procedure , were only able to obtain reliable results for the imaginary part of the convolution integral , eq .
( [ convol ] ) , our numerical results for the imaginary * and * the real part occur to be very stable .
going from an @xmath147 grid of @xmath148 points to a @xmath149 grid changes the results for the amplitudes ( over the whole angular range ) by less than 0.1% .
our numerical calculations were performed on an alphaserver 1000 4/266 . for the larger grid size
the evaluation of @xmath150 ( taking into account the 8 hadronic non- and single - flip amplitudes ) took less than 7 seconds per energy point and angle .
our calculation of the hard - scattering amplitude involves an expansion in powers of @xmath151 which is performed at fixed scattering angle . we keep only the leading order and next - to - leading order terms in this expansion .
hadron masses , however , are fully taken into account in flux and phase - space factors .
the reasons why we want to include hadron - mass effects in the hard scattering amplitude are twofold . on the one hand , we want to apply our model already at momentum transfers of only a few gev , where hadron masses may still play a role . on the other hand ,
we also want to make predictions for spin observables which require the flip of hadronic helicities . within usual perturbative qcd
such observables would vanish since a spin-1 gauge boson which couples to a nearly massless ( current ) quark can not flip the helicity of the quark . in the original diquark model  @xcite violations of hadronic helicity conservation are generated by massive v- diquarks on the account of introducing a mass parameter for the v - diquark . in this paper
we adopt another strategy which allows for a more consistent treatment of mass effects without introducing new mass parameters for the hadronic constituents . like the authors of ref .  @xcite we assume for every hadronic constituent that its 4-momentum is proportional to the 4-momentum of its parent hadron .
this requires to assign to every constituent of the hadron @xmath19 an effective mass @xmath152 , where @xmath153 is the fraction of the four- momentum of the hadron @xmath19 carried by the constituent . keeping in mind that the momentum fractions are weighted by the hadron das in the convolution integral , eq .
( [ convol ] ) , this means that quarks and diquarks acquire ( on the average ) a mass which is rather the mass of a constituent than a current particle .
this may be interpreted in such a way that at intermediate momentum transfers a quark or diquark surrounded by a cloud of @xmath154 pairs and gluons is not resolved but rather acts like a single effective particle with a corresponding mass .
this prescription for the constituent masses applies to the ( on - shell ) particles which correspond to the external legs of the feynman diagrams .
it has already been tacitly exploited when casting the hadronic spin wave functions into covariant form ( cf .
( [ wfcov ] ) and ( [ wfcovv ] ) ) . taking the quark - diquark da of eq .
( [ dap ] ) the average values of the effective masses of quarks and diquarks in a proton become @xmath155 analogously the effective mass of the heavy ( anti)quark in the meson becomes on the average half of the meson mass .
the choice of the masses which occur in the propagators of internal particles is not that obvious , since , e.g. , a quark in the incoming proton has mass @xmath156 whereas the same quark in the outgoing proton has mass @xmath157 .
so what should we take for the mass of this quark between two vertices ? in order to explain our recipe we , e.g. , consider the momentum of the heavy quark between the two gluon vertices in the third diagram of fig .
it is @xmath158 with @xmath159 , @xmath160 and @xmath161 representing the 4- momenta of the hadrons as indicated in fig .
[ kinem ] .
the quark mass in the propagator related to @xmath162 is now obtained by replacing the 4-momenta in eq .
( [ propmom ] ) by the masses of the corresponding particles , i.e. @xmath163 taking again into account that the momentum fractions are weighted by the hadron das , eq .
( [ dap ] ) and eq .
( [ dabc ] ) , the average heavy quark mass in the propagator becomes @xmath164 . extending this recipe to all the quark and diquark propagators and using the covariant wave - function representations , eqs .
( [ wfcov ] ) and ( [ wfcovv ] ) , for the external particles , leads to analytical expressions for the feynman diagrams which contain ( apart of the momentum fractions @xmath165 , @xmath166 , and @xmath167 ) only hadronic momenta , masses , spinors , polarization vectors , and helicities .
the occurring traces and spinor products can be evaluated in the usual way .
the final step is then to reexpress the ( massive ) mandelstam variables @xmath14 , @xmath3 , and @xmath168 in terms of massless ones @xmath169 and to make an expansion in terms of the ( small ) parameters @xmath151 , keeping the scattering angle @xmath170 fixed . in order to facilitate the treatment of propagator singularities
we retain only the leading and next - to - leading order term in this expansion . taking again as an example the propagator related to the momentum @xmath162 ( cf .
( [ propmom ] ) ) the propagator denominator becomes @xmath171 hence the quark propagator is given by @xmath172 i.e. , by neglecting mass corrections of @xmath173 we end up with propagator denominators which still appear like those of massless particles ( cf .
( [ prop ] ) ) .
apart from the simplifications in the treatment of propagator singularities and in the analytical expressions for the hard part of the amplitudes our approximate treatment of mass effects has a few other nice features . by including mass corrections of @xmath174 we obtain non - vanishing results for hadron - helicity - flip amplitudes and thus for polarization observables which require the flip of a hadronic helicity ( see the discussion in the next section ) . like the leading - order hadron - helicity - conserving amplitudes
also the hadron - helicity - flip amplitudes are @xmath175 gauge invariant with respect to the photon and @xmath176 gauge invariant with respect to the gluon . finally , crossing relations for the hadronic helicity amplitudes
are fulfilled up to @xmath174 , as has been shown for compton scattering off baryons and its crossed process @xmath177  @xcite .
for exclusive photoproduction of vector mesons @xmath178 one finds , altogether , 24 complex helicity amplitudes . by virtue of parity invariance only 12 of these helicity amplitudes
are independent . following the notation of ref .
@xcite we denote them by @xmath179 with @xmath180 . for our normalization of the helicity
amplitudes the unpolarized differential cross section takes on the form @xmath181 as examples for single spin observables , for which our model provides nontrivial predictions , we will consider the beam and target asymmetry .
the corresponding expressions are  @xcite @xmath182 for the photon and @xmath183 for the proton , respectively .    within the hard scattering approach
the energy dependence of the helicity amplitudes at fixed cm angle and large @xmath14 is roughly @xmath184 depending on whether the helicity of the hadronic constituents is conserved or flipped by one or two units .
this scaling behavior is modified by logarithms due to the running coupling @xmath44 and , eventually , the evolution of the hadron distribution amplitudes .
further deviations are due to the diquark form factors . for small momentum transfers
diquarks appear nearly point- like and the decay behavior of the helicity amplitudes is weakened by one power of @xmath14 . the energy dependence in eqs .
( [ scaling ] ) is only approached at large enough momentum transfer , where the diquark form factors become fully operational .    for the helicity
amplitudes our treatment of mass effects ( cf .
[ sec : mass ] ) leads to an expansion in powers of @xmath185 .
as compared to the hadron - helicity - conserving amplitudes @xmath186 and @xmath187 , the leading - order term of the single - flip amplitudes @xmath188 , @xmath189 , @xmath190 , and @xmath191 is suppressed by @xmath174 .
the leading - order term of the double - flip amplitudes @xmath192 and @xmath193 is even suppressed by @xmath185 and thus has the same order as the first ( non - vanishing ) mass - correction term of the hadron - helicity - conserving amplitudes . to simplify matters we only consider non- and single - flip amplitudes to leading - order and neglect double - flip amplitudes at all .    in the following we quote the analytical results for the functions @xmath105 and @xmath106 which have been obtained with the help of the symbolic computer program mathematica , in particular with the two program packages ` feynarts`@xcite and ` feyncalc`@xcite . in case of the hadron - helicity conserving amplitude @xmath187 ,
i.e. for the helicity combination @xmath194 , @xmath195 , @xmath196 , @xmath197 , the functions @xmath105 read :    [ fs1 ] @xmath198 \nonumber \\   & & \left . + x_1 y_2 \hat{t}^2 \left [ 4 z_1 ^ 2 ( y_1 - z_1 ) + y_1 ( x_2 + z_2 ) + z_1 ( 2 x_1 z_1 - 3 y_1 ) \right ] \right\}\ , , \\
%   f_{\left\ { 0 ,   - \frac{1 } { 2 } ; + 1 , - \frac{1 } { 2 } \right\}}^{(q , s ) } & = &    - \frac{6 m_p^2   } { \kappa \hat{t } } f_{\left\ { 0 ,   - \frac{1 } { 2 } ; + 1 , - \frac{1 } { 2 } \right\}}^{(q , v ) } = - 4 v   \frac{1 } { x_2 y_1 \hat{t } } \left\ { x_2 z_1 \hat{u } + y_2 z_2 \hat{s } \right\}\ , ,   \\
% f_{\left\ { 0 ,   - \frac{1 } { 2 } ; + 1 , - \frac{1 } { 2 } \right\}}^{(s ) } & = &   - 4 v \frac{1 } { x_1 y_1   \hat{t } } \left\ { x_1 z_1 \hat{u } + y_1 z_2 \hat{s } \right\}\ , , \end{aligned}\ ] ]    with @xmath199 where @xmath200 denotes the color factor and @xmath55 the fine - structure constant , respectively .
the functions @xmath106 are related to the @xmath105s in a simple way :    [ gs1 ] @xmath201    vector 4-point contributions @xmath202 and @xmath203 are , in general ( for all helicity combinations ) , suppressed by @xmath173 or even stronger .
relations analogous to also valid between the @xmath105s and @xmath106s entering the other @xmath204 amplitudes @xmath188 , @xmath186 , and @xmath189 .
the functions @xmath105 contributing to @xmath186  the other hadron - helicity conserving amplitude ( helicity combination @xmath205 , @xmath206 , @xmath196 , @xmath207 )  are obtained from those entering @xmath187 by interchange of the mandelstam variables @xmath208 and the momentum fractions @xmath209 :    [ fs2 ] @xmath210    the analytical expressions for the hadron - helicity - flip amplitudes @xmath188 , @xmath189 , @xmath190 , and @xmath191 are of similar length as the expressions quoted above and can be obtained from the authors on request . like in eqs .
( [ fs1 ] ) the vector 3-point contributions to @xmath190 and @xmath211 are proportional to the corresponding scalar ones .
it is also interesting to observe that due to the ( anti)symmetry properties of the feynman graphs under interchange of the momentum fractions @xmath77 and @xmath212 ( cf .
( [ gs1 ] ) ) , only diagrams in which the photon couples to the @xmath14 quark contribute to the @xmath204 amplitudes .
actually we have computed all the functions @xmath105 and @xmath106 , and relations like ( [ gs1 ] ) and ( [ fs2 ] ) served as a check of our analytical expressions . furthermore , we have verified @xmath175 gauge invariance with respect to the photon and @xmath176 gauge invariance with respect to the gluon .
the proof of gauge invariance is facilitated by the observation that not only the sum of all 16 tree diagrams provides a gauge invariant expression , but rather each of the functions @xmath105 and @xmath106 is gauge invariant by itself .
finally , we have recalculated a few of the diagrams by hand to confirm the outcome of mathematica .
fig .  [ result ] shows the diquark model predictions ( solid line ) for @xmath213 $ ] at @xmath214  gev together with the outcome of the pomeron- exchange mechanism ( dotted line )  @xcite and results of a two - gluon - exchange model ( dash - dotted and dashed line )  @xcite .
a direct comparison of our predictions with experiment is not possible yet since data are only available at low momentum transfers ( @xmath215  gev@xmath216 ) where the perturbative photoproduction mechanism is certainly not the dominant one .
the experimental situation , however , will hopefully improve as soon as the data analysis of the cebaf-93 - 031 experiment will be completed . nevertheless , fig .
[ result ] shows that the diquark model provides a reasonable extrapolation of the low-@xmath3 data .
thereby one should keep in mind that the cross section which decays exponentially in forward direction is expected to flatten around @xmath217 ( @xmath3 and @xmath168 large ) .
the few available data for photoproduction of pseudoscalar mesons and also photoproduction of @xmath5 and @xmath6 mesons which extend to larger values of @xmath3 exhibit just such a behavior  @xcite . as can be seen from fig .
[ result ] the forward cross section for @xmath218 production is reasonably well reproduced by simple pomeron phenomenology ( dotted line ) . as a mechanism for photoproduction of vector mesons donnachie and landshoff
@xcite proposed that the photon fluctuates into a quark - antiquark pair which subsequently recombines to a vector meson .
the resulting quark loop is connected to a quark in the proton via the exchange of a single pomeron which couples to quarks like an isoscalar photon .
the qcd - inspired version of the pomeron exchange of donnachie and landshoff  @xcite , in which the pomeron is replaced by two non - perturbative ( abelian ) gluons , has been applied by laget and mendez - galain  @xcite to photoproduction of @xmath1 mesons ( dash - dotted line ) .
the latter authors consider the two - gluon - exchange model as a link which connects the diffractive with the hard - scattering domain . at @xmath219
gev@xmath216 the two - gluon - exchange cross section exhibits a characteristic node which can be understood as an interference effect between the two feynman diagrams which enter the photoproduction amplitude . however , this picture changes drastically if the two gluons are allowed to couple to different quarks in the proton .
recently laget  @xcite tried to estimate such contributions by modeling quark - quark correlations inside a proton via a simple wave function .
he found that the node in the cross section is completely washed out and that the contributions in which the two gluons couple to different quarks inside the proton start to dominate the photoproduction of @xmath1 with increasing @xmath3 ( dashed line ) .
actually , for @xmath220  gev@xmath216 the result of laget becomes comparable with the diquark - model result .
this is not surprising since we know from the perturbative analysis that the hard - scattering mechanism , i.e. contributions from diagrams without loops in which all the hadronic constituents are connected via gluon exchange , should become dominant for large values of @xmath3 and @xmath168 .
correlations between hadronic constituents are automatically accounted for within such an approach and hence also within our diquark model .    to check whether our model provides a reasonable energy dependence we have also calculated the differential cross section at a much larger photon energy ( @xmath221  gev ) .
the result is depicted in fig  [ result2 ] ( solid line ) along with a two - parameter fit of @xmath222  gev@xmath216 data from zeus  @xcite .
the hatched area indicates the uncertainties in the forward cross section and the slope parameter . keeping again in mind that the cross section is expected to flatten with increasing @xmath3 the magnitude of our prediction occurs still to be within the range of an extrapolation of the low-@xmath3 data . the dashed line in fig .
[ result2 ] represents a leading - order calculation in which mass effects are neglected .
this means that only the two hadron - helicity conserving amplitudes @xmath186 and @xmath187 are taken into account in the leading - order cross section .
a comparison between the full calculation ( solid line ) and the leading - order calculation ( dashed line ) reveals that mass effects are still sizable and amount to @xmath223-@xmath224 of the full differential cross section .
they increase in forward and backward direction and are smallest around @xmath217 .    at this point
we want to comment on another attempt to calculate photoproduction of @xmath1 mesons within the diquark model .
the authors of ref .  @xcite use a different parameterization of the model .
only the hadron - helicity - conserving amplitudes @xmath186 and @xmath187 are taken into account in their calculation . with their numerical method
they encounter difficulties in calculating the real parts of these amplitudes .
but they assert that the real parts are small for their parameterization and therefore neglect them .
their results for differential cross sections are about one order of magnitude smaller than ours .
this is plausible if one looks at their values of @xmath225 and @xmath226 , i.e. the normalization  of the quark - diquark das , and the cutoff masses in the diquark form factors .
both are considerably smaller than ours .
in contrast to the authors of ref .
@xcite we think that the present experimental situation in @xmath1 photoproduction does not allow to discriminate a particular parameterization of the diquark model . however , we want to point out that our parameterization has already been applied successfully to a variety of other exclusive processes  @xcite . on the other hand ,
it is still unclear whether these processes can also be described satisfactorily with the parameterization of carimalo et al .
their parameterization has been determined by means of the proton magnetic form factor .
apart from @xmath1 photoproduction it has only been applied to charmonium decays into the @xmath227 final state  @xcite with the result that the data could not be reproduced .    in the kinematic region , where @xmath3 and @xmath168 are large
, the differential cross section at fixed scattering angle is expected to exhibit a scaling behavior of approximately @xmath228 , which is supposed to indicate that the hard photoproduction mechanism becomes relevant .
[ scalingfig ] shows @xmath229 for photon energies of @xmath230 , @xmath231 , and @xmath232  gev .
5 gev is presently the upper limit of cebaf , 10 gev will be reached by the planned upgrade of cebaf and 20 gev is a value which could be reached at future facilities like elfe or a further extension of cebaf which are presently still under discussion . comparing the diquark - model results for these energies
one observes that one is still far away from the @xmath228 scaling behavior .
the reasons have already been mentioned in sec .
[ sec : helamp ] .
a closer inspection shows that the cross section at @xmath233 scales rather like @xmath234 than @xmath228 for @xmath235  gev .
starting with point - like diquarks and a constant @xmath236 one would expect that the leading - order cross section ( without mass effects ) should scale like @xmath237 . by including the diquark form factors and cross - section contributions from hadron - helicity - flip amplitudes , which are already suppressed by one power of @xmath14 as compared to the leading - order contributions , the scaling behavior changes already to nearly @xmath228 .
finally , the approximate @xmath234 decay is caused by the fact that @xmath236 is taken as a running coupling constant ( with argument @xmath238 ) .
[ effects ] shows how the differential cross section ( at @xmath214  gev ) is influenced by different contributions within our model .
a comparison of the solid ( full calculation ) and the short - dashed ( mass effects neglected ) line reveals , e.g. , the influence of mass effects . for @xmath214  gev and @xmath239 they amount to @xmath240 of the full cross section .
the dash - dotted line represents the pure s - diquark part of the cross section .
it varies between @xmath241 and @xmath242 if one goes from backward to forward angles . since s- and v - diquark contributions add coherently ( and constructively , cf .
( [ fs1 ] ) ) this means that the corresponding amplitudes are of the same order of magnitude .
the differential cross section for the photoproduction of longitudinally polarized @xmath1s is given by the long - dashed curve .
it amounts to @xmath243 of the full cross section at @xmath239 .
this fraction increases in backward direction and decreases in forward direction .
finally , we have also examined the influence of the exponential in the @xmath1 da ( cf .
( [ dabc ] ) ) . neglecting this exponential one ends up with the dotted curve which is nearly indistinguishable from the full calculation .
this also shows that we are not very sensitive to the @xmath244 endpoint regions when performing the convolution integral eq .
( [ convol ] ) .    as examples for
single - polarization observables the beam asymmetry @xmath245 and the target asymmetry @xmath246 are depicted in fig .  [ asymmetry ]
. the beam asymmetry would be already non - zero , if only the leading - order helicity amplitudes @xmath186 and @xmath187 were taken into account ( cf .
( [ gasy ] ) ) .
it would not even be necessary that the helicity amplitudes acquire a non - trivial phase .
on the other hand , a non - trivial target asymmetry demands for both , non - vanishing hadron - helicity - flip amplitudes and non - trivial phases of the amplitudes ( cf .
( [ pasy ] ) ) . as already mentioned at the end of sec .
[ sec : elamp ] phases are generated through propagator singularities , hadron - helicity - flip amplitudes occur in our model as mass corrections .
as it turns out within our model , also the beam asymmetry is considerably affected by the inclusion of the hadron - helicity - flip amplitudes . in a recent paper
@xcite the influence of nucleonic resonance effects on polarization observables in @xmath1-meson photoproduction has been investigated within a constituent - quark model in which the diffractive contribution is produced by t - channel pomeron exchange .
it is interesting to observe that the result of these authors for the target asymmetry at a rather small photon energy of @xmath247  gev is comparable in size to our result at @xmath214  gev .
this is remarkable since two completely different concepts are at work : on one hand nucleonic resonance contributions , on the other hand a purely perturbative mechanism . also for the beam asymmetry some similarities can be found .
the corresponding predictions of both approaches are again comparable around @xmath233 .
a ( somewhat smaller ) negative beam asymmetry has also been found in another model which is based on pomeron exchange , pion exchange , and a direct knockout production from the strange - quark sea in the proton  @xcite .
it would be interesting to see experimentally , whether polarization observables like the target asymmetry @xmath248 , the recoil polarization @xmath249 , or the vector - meson polarization @xmath250 , which vanish in leading - order perturbative qcd , are still sizable at momentum transfers of a few gev .
such observables would be a more sensitive tool than the unpolarized differential cross section to figure out whether the kind of nonperturbative ingredients ( diquarks , effective constituent masses ) which are taken into account in our perturbative approach , are already sufficient , or whether much more non - perturbative physics is still at work at a few gev of momentum transfer .    at the end of this section
we want to comment on the perspective to extend this calculation to photoproduction of @xmath7 .
it is a simple matter , since we have already assumed a flavor dependent da for the vector meson ( cf .
( [ dabc ] ) ) .
the only thing to do is to replace the @xmath1 mass by the @xmath7 mass and the decay constant @xmath251 by @xmath252 .
one , however , has to keep in mind that we have applied an expansion in the parameter @xmath151 . in order
that this parameter is approximately of the same size in @xmath1 and @xmath7 production the photon lab energy has to be increased approximately by a factor @xmath253 when going from @xmath1 to @xmath7 .
a second criterion for the reliability of our approximate treatment of mass effects is that the coefficient of the first mass correction term should become small compared to the leading order .
this coefficient is angular dependent and becomes large at small scattering angle .
that means , if the size of the mass corrections relative to the leading order should be about the same in @xmath1 and @xmath7 production ( @xmath151 comparable ) , one has to consider much larger momentum transfers for @xmath7 than for @xmath1 production
. a more detailed discussion of this point along with some numerical results can be found in ref .
a possible way to apply the perturbative approach to @xmath7 production at a few gev of momentum transfer would be to retain our prescription for the choice of the effective masses , but to refrain from the final expansion in @xmath151 .
this has not been attempted yet since the analytical expressions for the amplitudes become very lengthy and the treatment of propagator singularities in the convolution integral , eq .  ( [ convol ] ) , becomes much more involved .
we have investigated the reaction @xmath254 in the few - gev momentum transfer region where unpolarized differential cross - section data from the cebaf-93 - 031 experiment at jlab are expected to become available very soon .
our theoretical approach is based on perturbative qcd supplemented by the assumption that baryons can be treated as quark - diquark systems .
the same approach has already been successfully applied to other photon - induced hadronic reactions  @xcite .
the motivation for introducing diquarks is , above all , to extend the perturbative approach from large down to moderately large momentum transfers .
there are mainly two reasons why this can be achieved with diquarks .
first , the overall momentum transfer has to be shared between less hadronic constituents than in the pure quark hsa .
hence the gluons which keep the hadronic constituents collinear are on the average harder and thus the running strong coupling constant becomes smaller .
the second reason is that a certain class of non - perturbative effects , namely quark - quark correlations inside a baryon , are effectively accounted for by means of diquarks .
the fact that a heavy quarkonium is produced in the final state entails a considerable reduction in computational effort .
only a small fraction of the feynman diagrams which contribute in general to photoproduction of arbitrary mesons occurs in @xmath1 production . to be more specific , these are just those diagrams in which the photon fluctuates into the heavy quark - antiquark pair which is then connected to the quark and the diquark in the proton via the exchange of two gluons , respectively . both , the scalar and the vector - diquark component of the proton have been taken into account .
we have also improved the diquark model in the sense that hadron - mass effects in the hard - scattering amplitude have been included in a systematic way .
this has been accomplished by assigning to each hadronic constituent an effective mass @xmath255 ( @xmath256 longitudinal momentum fraction of the constituent ) and expanding the scattering amplitude in the small parameter @xmath151 .
only the leading and next - to - leading order terms in this expansion have been kept .
the leading - order term provides for the hadron - helicity - conserving amplitudes and agrees , of course , with the massless result .
the next - to - leading order term represents a power- correction and gives rise to non - vanishing hadron - helicity - flip amplitudes .
it enables us to make predictions for polarization observables which require the flip of a hadronic helicity like , e.g. , the target asymmetry or the recoil polarization .
we want to point out that our treatment of mass effects still preserves gauge invariance with respect to the photon and the gluon .
our numerical studies have been performed with the diquark - model parameters and the quark- diquark das proposed in ref .
the da of the @xmath1 meson has been taken from the literature  @xcite .
it fulfills qcd sum - rule constraints .
hence we have been able to make predictions without introducing new parameters .
we have paid special attention to the correct and numerically robust treatment of propagator singularities . in the absence of large-@xmath3 data
( @xmath257  gev@xmath216 ) our results can only be compared with the available data which go up to @xmath258  gev@xmath216 .
our predictions for the differential cross section at two rather different photon energies ( @xmath214  gev and @xmath259  gev ) which are typical for cebaf and hera , respectively , occur to be still within the range of an extrapolation of the existing low-@xmath3 data .
it is worth mentioning that the predictions of a two - gluon - exchange model , in which the gluons are considered as a remnant of the pomeron  @xcite , become comparable to our result for @xmath257  gev@xmath216 if the gluons are allowed to couple to different quarks in the proton .
this is not surprising because contributions in which all the hadronic constituents are ( approximately ) kept collinear by means of gluon exchange to prevent the hadrons from breaking up become increasingly important as soon as one enters the hard - scattering domain .
just those contributions are assumed to be the dominant ones within our hsa - based diquark model .
we have also investigated within our model whether the ( asymptotic ) @xmath228 scaling behavior of the differential cross section at fixed angle is already realized for @xmath14 values between @xmath231 and @xmath260  gev . after dividing out the strong running coupling ( which occurs to the sixth power in the cross section ) an approximate @xmath228 scaling has been found .
our investigation of hadron - mass effects has revealed that they are still non - negligible in the few- gev momentum transfer range and photon energies of the order of 10  gev .
correspondingly , also polarization observables which require the flip of a hadronic helicity , like the target polarization , occur to be sizable .
finally , we have considered the possibility to extend the present calculation to photoproduction of @xmath7 mesons . in order
that our treatment of mass effects remains reliable for @xmath7s one has to go to much larger energies and momentum transfers where experimental data can not be expected .
the situation with respect to the momentum transfer could perhaps be improved by refraining from the final expansion in @xmath151 .
this , however , will considerably complicate the treatment of propagator singularities and make the analytical expressions for the amplitudes much more complicated
. therefore it has not been attempted as yet .
our diquark - model calculations for @xmath261 will meet the first severe test in the near future when the data analysis of the cebaf 93 - 031 experiment will be completed .
with these data , which go up to @xmath3 values of about 5  gev@xmath216 , one will get a first glimpse how different production mechanisms ( pomeron exchange @xmath262 gluon exchange @xmath262 hard - scattering mechanism ) in high - energy photoproduction evolve with increasing momentum transfer . however , in order to really touch the hard - scattering domain much higher photon energies than those which can be reached at the present stage of cebaf would be desirable .
in addition to unpolarized differential cross sections , polarization observables which require the flip of hadronic helicities , could serve as an indicator for the relative importance of hard and soft photoproduction mechanisms , since hard scattering is closely connected with hadronic helicity conservation .    c. f. berger would like to thank the paul - urban - stipendienstiftung for supporting a visit at the institute of theoretical physics at the university of graz during which part of this paper has been completed .
email address : carola.berger@sunysb.edu email address : wolfgang.schweiger@kfunigraz.ac.at g. audit _
et al . _ ,
cebaf proposal no .
93 - 031 ( 1993 )
. j. a. crittenden , _ exclusive production of neutral vector mesons at the electron - proton collider hera _ , springer tracts in modern physics * 140 * ( springer , berlin heidelberg , 1997 ) .
a. donnachie and p. v. landshoff , phys .
b * 348 * , 213 ( 1995 ) .
a. donnachie and p. v. landshoff , phys .
b * 185 * , 403 ( 1987 ) .
a. donnachie and p. v. landshoff , nucl .
b * 311 * , 509 ( 1989 ) .
j. c. collins , l. frankfurt , and m. strikman , phys .
d * 56 * , 2982 ( 1997 ) .
m. diehl , t. feldmann , r. jakob , and p. kroll , eur .
j. c * 8 * , 409 ( 1999 ) . for a comprehensive review of hard exclusive processes and further references
see , e.g. , s. j. brodsky and g. p. lepage , in _ perturbative quantum chromodynamics _ , edited by a. h. mueller ( world scientific , singapore , 1989 ) . v. l. chernyak and a. r. zhitnitsky , phys
112 * , 173 ( 1984 ) .
n. isgur and c. h. llewellyn smith , nucl .
b * 137 * , 526 ( 1989 ) .
et al . _ ,
c * 66 * , 267 ( 1995 ) .
j. gronberg _
et al . _ ,
d * 57 * , 33 ( 1998 ) .
j. bolz and p. kroll , z. phys .
a * 356 * , 327 ( 1996 ) .
a. v. radyushkin , phys .
d * 58 * , 114008 ( 1998 ) .
r. jakob , p. kroll , m. schrmann , and w. schweiger , z. phys . *
a347 * , 109 ( 1993 ) .
p. kroll , th .
pilsner , m. schrmann , and w. schweiger , phys .
b316 * , 546 ( 1993 ) .
p. kroll , m. schrmann , and p. a. m. guichon , nucl .
a * 598 * , 435 ( 1996 ) . c. f. berger , diploma thesis , technological university graz ( 1997 ) ; p. kroll , m. schrmann , k. passek , and w. schweiger , phys .
d * 55 * , 4315 ( 1997 ) .
r. jakob , phys .
d * 50 * , 5647 ( 1994 ) .
t. huang , nucl .
[ proc . suppl .
] b * 7 * , 320 ( 1989 ) .
m. benayoun and v. l. chernyak , nucl .
b * 329 * , 285 ( 1990 ) .
m. wirbel , b. stech , and m. bauer , z. phys .
c * 29 * , 637 ( 1985 ) .
t. feldmann and p. kroll , phys .
b * 413 * , 410 ( 1997 ) .
g. r. farrar , g. sterman , and h. zhang , phys .
* 62 * , 2229 ( 1989 ) . c. carimalo , n. arteaga - romero , and s. ong , preprint hep - ph/9712485 ( 1997 ) .
m. anselmino , f. caruso , and f. murgia , phys .
d * 42 * , 3218 ( 1990 ) .
m. pichowsky ,  .
avkli , f. tabakin , phys .
c * 53 * , 593 ( 1996 ) .
j. kblbeck , m. bhm , and a. denner , comp .
. commun . * 60 * , 165 ( 1990 ) .
r. mertig , m. bhm , and a. denner , comp .
* 64 * , 345 ( 1991 ) .
laget and r. mendez - galain , nucl .
a * 581 * , 397 ( 1995 ) .
laget , preprint dapnia / sphn-98 - 54 ( 1998 ) .
j. ballam et al . , phys .
d * 5 * , 3150 ( 1973 ) .
barber et al . , z. phys .
c * 12 * , 1 ( 1983 ) .
r. l. anderson et al . , phys .
d * 14 * , 679 ( 1976 ) .
m. derrik et al . , phys .
b * 377 * 259 ( 1996 ) . c. carimalo and s. ong , z. phys .
c * 52 * , 487 ( 1991 ) .
q. zhao , j .-
didelez , m. guidal , and b. saghai , preprint nucl - th/9906013 ( 1999 ) .
a. i. titov , y. oh , s. n. yang , and t. morii , phys .
c * 58 * , 2429 ( 1998 ) . c. f. berger and w. schweiger , preprint hep - ph/9909270 , to appear in proceedings of the int / jlab workshop on _ exclusive and semi - exclusive processes at high momentum transfer _ , newport news , usa , may 1999 .