Uncovering the Scaling Laws of Hard Exclusive Hadronic Processes in a Comprehensive Endpoint Model

We show that an endpoint overlap model can explain the scaling laws observed in exclusive hadronic reactions at large momentum transfer. The model assumes one of the valence quarks carries most of the hadron momentum. Hadron form factors and fixed angle scattering are related directly to the quark wave function, which can be directly extracted from experimental data. A universal linear endpoint behavior explains the proton electromagnetic form factor, proton-proton fixed angle scattering, and the t-dependence of proton-proton scattering at large s>>t. Endpoint constituent counting rules relate the number of quarks in a hadron to the power-law behavior. All proton reactions surveyed are consistent with three quarks participating. The model is applicable at laboratory energies and does not need assumptions of asymptotically-high energy regime. A rich phenomenology of lepton-hadron scattering and hadron-hadron scattering processes is found in remarkably simple relationships between diverse processes.


Experimental Regularities
The experimental study of differential cross sections of hard exclusive hadronic reactions at high energy reveals a remarkable pattern: They are described by power laws [1][2][3]. A model explanation exists [4][5][6][7], yet it is not satisfactory [8] at the energies of experimental measurements. We are driven to find a consistent explanation of experimental regularities by re-examining all the facts from a fresh point of view.
"Hard" reactions are those which depend on a single large scale Q 2 > GeV 2 , or several large scales with a fixed ratio. It is remarkable that the proton electromagnetic form factor F 1 (Q 2 ) agrees well with a decreasing power of Q 2 for Q 2 5 GeV 2 [9]. For large momentum transfer, it is remarkable that pp → pp fixed-angle cross section dσ/dt agrees well with a decreasing power of Q 2 ∼ s [10], where √ s is the center of mass energy. There are many other examples.
We have re-evaluated the phenomenology of power-law dependence or "scaling laws" for exclusive reactions. Due to history, the most simple and plausible explanation failed to be developed. The model appears in the literature as the Feynman process, also known as the Drell-Yan model, also known as the endpoint overlap model [11][12][13]. There is actually much to be learned and much that is new when the model is objectively explored.

The Endpoint Overlap Model
In their evaluation of the endpoint region, Brodsky and Mueller [14] wrote that "its contribution depends sensitively on the hadronic wave functions". The discussion discovered no actual fault in the endpoint contribution. Instead of finding a flaw, the section ends with a weak suggestion to assume validity of a short-distance perturbative model, as "at least plausible", adding "in any case there is currently no comprehensive alternative theory of these processes".
We suggest that the lack of a comprehensive alternative theory came out of a historical failure of the endpoint contribution to be fully appreciated and developed. In the concluding remarks we link this to the history of early development of perturbative QCD, which is an era long past.
In reviewing the current status we noticed several facts: • The predictions of all models depend on the wave functions. For reasons we believe are obsolete, the opportunity to learn about wave functions using data was bypassed in the promotion of short-distance (SD) models [4][5][6][7].
• Great emphasis was imposed early on asymptotic limits. The motivation was not to learn about hadrons, but an attempt to make hadrons irrelevant for the goal of establishing QCD.
• The asymptotic limits of QCD are now understood to be of negligible experimental relevance. The asymptotic limits of QCD predictions have also never actually been established. Instead limits of models have been established. None of the work assuming a model has gone an inch beyond the boundaries of the model itself. In particular, it has not been shown anywhere that the pion form factor of the general theory known as QCD necessarily falls faster than 1/Q 2 , despite considerable effort to force such a conclusion. Careful reading is needed to verify this. For example, Farrar and Jackson [4] claim an asymptotic limit in opening lines, without actually supporting the claim: The contrary information about regions outside the model is buried in the footnote labeled Ref. 13 in the paper.
• Now that QCD is established, every integration regime, including those contradicting the assumptions of SD models, needs to be considered. Interest in the larger theory and hadron structure has eclipsed the goal of exhibiting a model based solely on perturbation theory.
• The main reason for early interest in perturbative models was power-law behavior. Inexperience with more general models created a folk-lore that "soft" non-perturbative wave functions would lead to exponential dependence on a large Q 2 scale. This is false. As we review, power-law behavior is generic from the endpoint region.
• Divisions in the field separated groups into two camps. Relying on perturbation theory appears to be more theoretically ambitious, but it is actually less general than representing dynamics with wave functions. For one thing, an arbitrary order of perturbation theory can be subsumed into equivalent wave functions, but not vice versa.
• Calculations in perturbation theory use the Fock state basis of free field theory. Considerable effort has been dedicated to making the endpoint region of perturbative calculations go away at asymptotically high energies, towards demonstrating perturbative self-consistency. None of that work is relevant to non-perturbative wave functions, which use a different basis of fully-interacting quanta. It is not logically self-consistent to extend the asymptotic pQCD-based suppression of the endpoint region to wave functions extracted from experimental data.
• Despite years of study, very little is known about pions and protons. The proof comes from the dearth of definite information about pions and protons in terms of non-perturbative wave functions. Contrary to the bias of perturbative QCD, it is definitely possible and absolutely productive to use experiments to learn about non-perturbative wave functions.
Recently Chang, Cloèt, Roberts, Schmidt, and Tandy [15,16] have computed the pion form factor with a method described as self-consistent for all space-like Q 2 . The paper highlights the asymptotic SD model's prediction being about 3 times smaller than the experimental form factor. Ref. [15] states the asymptotic estimate is incapable of converging to a realistic value below Q 2 ∼ 1000 GeV 2 . This is typical of asymptotic SD estimates: The estimates require fabulously high momentum transfer to apply. We agree that the assumptions made in setting up SD models are contradicted by the application of the models to existing momentum transfers.
The calculations of [15] are made in an overlap model emphasizing first-principles predictions of the pion wave function. They are similar in spirit to overlap models in the relativistic impulse approximation [35][36][37]. These models are very successful in describing data at low Q 2 . Our approach accepts the validity of these models as integral representations of form factors, while also extending the scope to other processes. We differ from most models by not attempting to know wave functions in advance. Many studies have been restricted to estimating the endpoint region to order of magnitude. We will show that the end-point region of the wave functions is not only determined, but over-determined, by experimental regularities found in power laws. By measuring the wave functions rather than predicting them we find the endpoint overlap model is a consistent comprehensive description. Rather unexpectedly, our approach extracting information from wave functions is quite consistent with the trend predicted by [15]. Section 2 derives the endpoint overlap constituent counting rules for form factors. Section 3 extends the rules to exclusive hadron-hadron reactions. These rules explain why scaling laws should be observed at the limited energies of laboratory experiments, and how scaling laws are correlated with the number of quarks scattering. We cannot explain why this predictive regularity of all exclusive reactions surveyed has been overlooked. Concluding remarks with brief historical commentary are given in Section 4.

Endpoint Power Counting Rules
In this Section we derive the endpoint constituent counting rules, which predict the scaling power of Q 2 in terms of the number of constituents. Comparing experimental data to these rules finds consistency with three (3) quarks scattered in every proton reaction we have surveyed. The discussion will be organized in increasing levels of detail, beginning with the simplest case of the pion electromagnetic form factor. The form factor is defined by

The Pion Form Factor
where J µ em is the electromagnetic current operator. In a gauge-invariant local field theory the photon interaction involves one (1) struck parton. The minimum number of constituents in the pion is two. The dynamical question of elastic scattering is how scattering one constituent can scatter the entire hadron. This is answered by the quantum mechanical overlap of wave functions. Figure 1 conveys the qualitative picture. The endpoint region is dominated by some transverse hadronic scale "Λ" for which the slow parton obeys x 2 Λ/Q, while the struck parton obeys 1 − x 1 Λ/Q. For the entire region the endpoint contribution is such that the transverse momentum integrations do not contribute any power of Q 2 to the form factor. This happens to be the feature causing the short-distance model to be impossible to justify in this region.
The form factor is The electromagnetic interaction matrix element T contains the quark charges, a gamma matrix, and momentum-conserving delta functions. Wave functions with several independent spin structures are possible, with discussion postponed to Section 2.1.1. We define the delta functions to include factors representing momentum conservation, hence k 1 + k 2 = P, etc. in a frame we now specify. Let P, P be the 4-momenta of the pions, with P = P + q, q 2 = −Q 2 . Choose a Lorentz frame with Cartesian labels (E, p X , p Y , p Z ) where E is the energy. Thus q = (0, Q, 0, 0); Let k µ 1 (k µ 1 ) be the momenta of the struck parton before (after) scattering. Due to the change of direction of the fast momenta, the meaning of symbols " k i " must be adapted to be orthogonal to each hadron's direction. We introduce a basis of transverse vectors adapted to the particular hadronic momenta :ŷ =(0, 0, 1, 0) =ŷ ; P ·ŷ = P ŷ = 0; With these coordinates, the components of the quark momenta are There are only three free parameters, and the quanta are not strictly constrained to the perturbative mass shell. By hypothesis, amplitudes are concentrated in the kinematic region shown.
Integrating over a fourth (minus, or virtuality parameter) concentrated on the region is equivalent. Momentum conservation of the un-struck spectator is k 2 = k 2 . Momentum conservation of the struck quark k 1 = k 1 + q yields four constraints: Solving gives See Figure 2. We emphasize these are exact kinematic relations of the model, regardless of the value of m π . P P' Figure 2: Endpoint kinematics in the pion case. Pion momenta are shown as dashed arrows, while quark momenta are solid arrows. Two isosceles triangles representing energy and momentum conservation must close. The transverse and longitudinal momenta of one spectator covers the difference. By inspection, k x = xQ/2, where x is the momentum fraction of the slow quark.
Evaluating F π gives A soft non-perturbative wave function means that Φ π (x, k x ) is a rapidly falling function with a scale k x Λ, where Λ ∼ 300 MeV. To see how power-law behavior emerges, consider an exponential function multiplied by a function of x: As Q → ∞ the exponential is driven toward the endpoint x 2 Λ/Q. Integrals dominated by their endpoints 1 have an asymptotic series expansion [38] developed 1 The same expansion applies to a wide class of integrals with saddle points in x approaching the endpoint with integration by parts. The first term is given by The last line dropped the exponentially small term from the upper limit. Using repeated integration by parts gives Here φ (n) = ∂ n φ/∂x (n) (0). It is generally expected that light cone wave functions have non-zero x-derivatives at the endpoints. Data from parton distributions supports this. The perception that soft wave function overlaps are incompatible with power-law dependence at large Q 2 is in error. The result does not depend strongly on the exponential dependence or factored form of Eq. (7). For example a Gaussian dependence of Φ ∼ e −k 2 x /Λ 2 → e −x 2 Q 2 /Λ 2 uses ∂/∂x 2 integration by parts. Pursuing that, the region x 2 Q 2 Λ 2 produces just the same power expansion as the region xQ Λ with different constants. It is also always possible to write This simply definesφ(x, k x ) as a function with an exponential dependence removed, namely a function that varies more slowly. Integration by parts proceeds as before, withφ (n) replacing φ (n) . The result is a series in powers of 1/Q which is qualitatively unchanged so long asφ is slowly varying.

Electromagnetic Gauge Invariance
Phenomenological quark models often violate current conservation. Projecting a γ µ vertex with g µν − q µ q ν /q 2 hides the problem without actually curing it. The endpoint overlap model with massless quarks satisfies electromagnetic gauge invariance automatically, passing the test < P |q µ J µ em |P >= 0. This is a very detailed topic, explaining why we postponed details of the Dirac algebra in Section 2.1.
Let be the difference of quark momenta in one pion. The most general wave function for By permuting γ 5 through the wave function, the terms with one gamma ("chirally even") have antiparallel helicity and are antisymmetric in spin quantum numbers, and vice versa for chirallyodd. The A, B terms are proportional to the large number P. These track the Fermion fields which are the largest under a Lorentz boost, making them leading order in energy. Thus A is leading in power of Q, with zero orbital angular momentum while B is also leading in Q, while representing one unit along the momentum axis. (Power counting in SD models is different. Only the A wave function is large as → 0. ) The C, D terms are sub-leading on every basis.
The form factor's integrand contains a trace over the Dirac indices of wave functions from two pions, three quark propagators, and the perturbative photon vertex (γ µ ). Chiral selection rules give zero unless even-even and odd-odd chirality wave functions are composed. There are six non-zero combinations denoted AA , BB , CC , DD , AD , BC , plus their complex conjugates. Our demonstration of gauge invariance comes from contracting q µ with the vertex and computing the terms one by one.
The final result is zero when evaluated using the kinematic conditions of Eq. 4. The result is not obvious from manipulating perturbative Ward identities, and the full calculation is quite extensive. The zero-result requires massless quark propagators, namely chirality conservation, as we assume. Apparently a feature of the global chiral symmetry of the model protects the gauge symmetry of the model. We would like to understand this better, while the issue goes beyond the scope of this paper.

Measuring the Wave Function
Experimentally very little is known about the actual large Q 2 dependence of F π (Q 2 ). In the space-like region of momentum transfer Q 2 F π (Q 2 ) rises and appears to approach a plateau in the region of Q 2 ∼ f ew GeV 2 [39][40][41]. That represents a serious problem for SD models, because Q 2 is far too small for the asymptotic regime to set in. The attempt to interpret the flaw as a "bonus," namely a confirmation of the model's predictions despite contradicting the model's requirements, does not stand up to careful examination. π π γ * , ρ... It is a serious matter that existing data is well fit by vector meson dominance, associated with the pole at time-like Q 2 of the ρ meson [42][43][44]. The importance of the ρ is a problem for SD models. The models [4,45] are based on expanding the wave function of the pion about zero quark separation. That expansion contains no terms capable of representing the ρ meson's singularity: See Figure 3. The asymptotic contribution of the short distance model [4] is MeV is a test of the model. The normalization is too small compared to the experimental data. More ambitious work exploiting self-consistent Bethe-Salpeter wave functions is more convincing [46][47][48]. Yet none of the debate about asymptotic limits matters for our approach, which is based on extracting information from what is observable at laboratory energies.

An Informative Sum Rule
We mentioned in the introduction that folklore about exponentially-dependent form factors is in error. Analysis produces a sum rule for the endpoint-overlap model using Mellin transforms. Represent Positivity of Φ gives one constraint. The integral over N is done over a contour in a strip where Φ N is analytic. Then The right hand side is the Mellin transform of F π with respect to Q: The formula pinpoints what knowledge of F determines about Φ. The large Q dependence is determined by singularities at N > 0. These come from any N k > 0 singularities of the k dependence, and any N x < 1 singularities of the x dependence. Thus: • If Φ falls exponentially with increasing k then Φ N (x) is analytic for N > 0. All large Q singularities are found from the x dependence. If F π ∼ 1/Q 2 , then Φ(k ∼ 0, x) ∼ x 1 , modulo logarithmic factors. Term by term the x dependence determines the powers and terms in the asymptotic series in Q 2 , while the k dependence contributes to the coefficients.
• Suppose Φ has a power-law tail at large k. If Φ ∼ k −2 modulo logarithms, then Φ N (x) will have a pole at N = 2, producing Q −2 dependence. Term by term the k dependence determines the powers and terms in the asymptotic series in Q 2 , while the x dependence contributes to the coefficients, provided Φ is not as large as x 1 • The combination of both x 1 and k −2 can produce a double pole, which translates to possible The last two options have been extensively explored by the short-distance models. By calculating the large k dependence in the first step, they select in advance one of the regions the analysis shows are possible in general. The perturbative calculations also agree with the general analysis by producing Q −2 /log(Q 2 ) behavior from wave functions with corresponding power law behavior. After many years and a considerable investment in manpower, the same calculations have been found not to apply at finite Q 2 , and also to produce contributions that are too small to explain experimental data [8]. That indicates the endpoint region as being dominant at large Q 2 .
That leaves the first option. Experimentally Q 2 F π (Q 2 ) ∼ 0.4 GeV 2 , which leads to with the leading singularities determined by Φ(x) ∼ x 1 . This contradicts the perturbative SD model wave function (sometimes inappropriately called the "asymptotic" wave function). There are no compelling arguments to calculate the x-dependence of a wave function with perturbation theory. Ref. [15] fits a numerical calculation to the equivalent of Φ(x) ∼ x 1.6 , which seems to be reasonably consistent.

Counting Rule 1
Summarizing this section, under the universal assumption that Φ(x, k x ) falls rapidly to constrain x 1/Q and φ(x) ∼ x A , the cost of overlapping to retain one spectator is a "slow quark probability" factor of 1/Q A+1 .

The Proton Electromagnetic Form Factor F 1 : The Probability of Two Slow Quarks
While our objective focuses on power-counting, we believe value is added by including considerable details in the calculation. Here we compute the proton electromagnetic form factor, F 1 , assuming end point domination. There have been several earlier calculations of this form factor in different models [49][50][51][52][53][54][55][56][57][58][59]. The basic diagram for calculating proton electromagnetic form factors is shown in Fig. 4. The momenta P, P , Q are the same as in Eq.(2), and quark momenta use the same notation as Eq. 3. Let Y be the proton wave function to three quarks, and let the electromagnetic vertex be (10) where N, N are Dirac spinor functions. The momentum space wave functions with leading power of P is written as [60,61] Here α, β, γ are Dirac indices, V, A, T are scalar functions of the quark momenta(k i ), N c is the number of colors, C the charge conjugation operator, σ µν = i 2 [γ µ , γ ν ], and f N is a normalization. P P' k 1 k' 1 q Figure 5: Endpoint kinematics in the proton case. Proton momenta are shown as dashed lines, quark momenta as solid lines. Isosceles triangles representing energy and momentum conservation close as in Figure 2. The momenta of all spectator constituents sum to cover the differences P − k 1 and P − k 1 .
Several combinations appearing in Y Γ µ Y give leading order contributions to F 1 . Consider, for example, Collecting all coefficients proportional to this term written gives Inserting Eq. 13 in Eq. 10 gives where, The delta functions lead to kinematics similar to Eq.(4). Momentum conservation requires k µ 1 = k µ 1 + q µ and k 2 = k 2 . One transverse momentum of the struck quark is unconstrained except by wave functions, while the other transverse momentum is constrained by the relation previously found for the pion: The measure is then replaced by Integrating over the unconstrained variables gives That leaves the integration depending on Q as Once more consider an exponential ansatz where ψ( k T ) ∼ exp − k T 2 /Λ 2 . That leads to Evaluated at k 1x = (x 2 + x 3 )Q/2, both x 2 and x 3 range over intervals of size Λ/Q. For wave functions that are uncorrelated products, ψ = ψ(x 1 )ψ(x 2 ), the probability of finding two slow quarks is precisely the product of two slow quark probabilities. For quark wave functions going like Experimental data finds that Q 4 F 1 ∼ constant for Q 2 ≥ GeV 2 . The data indicates that A ∼ 1, namely that quark wave functions must go like x 1 i near x i ∼ 0. We emphasize that this result does not require extremely large Q 2 (large logarithms of Q 2 ). The estimates are based on comparison with the transverse size of the proton, Q 2 >> Λ 2 for Λ 2 ∼ 0.1 GeV 2 .

A Typical Perturbative Question
The approach based on short-distance perturbation theory typically asks how these relations behave when "soft gluons" are added.
The notion of adding soft gluons is tied to a basis of Fock state wave functions used in perturbation theory. Since the interacting theory is not the free theory, Feynman diagrams represent the interactions with gluons of all momenta. That does not represent our model: By construction, we are concerned with the full wave functions of the interacting theory. Thus the calculation using the fully interacting wave functions is self-consistent without adding gluons.
The question of soft gluons is vital for the internal consistency of the SD model and its (conceptually different) estimate of the endpoint contribution. Assuming one desires a perturbative description, experience indicates that any finite number of soft gluon internal diagrams will not revise the leading power of Q 2 . Considerable effort has gone into arranging calculations that would be simultaneously compatible with the assumptions of short distance. That has led to statements that Sudakov effects suppress the endpoint region. These statements refer to the short-distance model, not ours. If it is true that the endpoint region of the SD model is negligible, it has no bearing on the endpoint region of all possible models expressed in different quantummechanical bases. None of it is our concern once the focus is on extracting non-perturbative information from experiments.
Nevertheless it is interesting to check that soft-gluons do not change the leading power behavior. We feel that a specific calculation is more convincing than an estimate, and present one in the Appendix.

The Pauli Form Factor F 2
We have investigated the large Q 2 dependence of the Pauli form factor F 2 in the endpoint model. We find 1/Q 5 dependence occurs in more than one way, together with wave functions that go like x(1 − x). This result is surprising and impinges directly on the issue of quark orbital angular momentum [62] and "the shape of the proton" [63,64]. Yet the calculations we have available are complicated, and too detailed to be appropriate to review here. We plan to present them in an future paper [65]. For the purposes of this survey, we can objectively report that it has not been shown that experimental data [34] measuring F 2 at large Q 2 is in conflict with power-counting of the endpoint overlap model. The lack of previous work is itself remarkable because it has been shown that F 2 ∼ 1/Q 5 is incompatible with short distance models.

Endpoint Constituent Counting Rules for Form Factors
We are now in a position to state the leading power endpoint constituent counting rules for form factors. Let there be n I N (n OUT ) quarks in the I N (OUT) state hadron. For now choose n I N = n OUT = n. One constituent is scattered, requiring n − 1 quarks to be slow. The form factor is given by induction, For wave functions with linear dependence as x j → 0, each extra constituent beyond the valence configuration causes a suppression of the form factor by a power of 1/Q 2 . The leading power dependence is These happen to be the same rules as the early "dimensional" counting rules [1], but for entirely different reasons. Hard propagator factors are not the explanation. The explanation lies in the phase space to find quarks available to scatter. The general dominance of phase space over hard scattering is reminiscent of the independent scattering mechanism originally discovered by Landshoff [66,67].
We will be straightforward with what is new in the power law. It has been noticed again and again that the endpoint contribution cropped up and competed with the short-distance model of form factors. Many papers have dealt with the issue as a troublesome instability of short-distance dominance. Yet we are not aware of an explicit, positive statement of the predictive regularity between the number of scattered constituents and the observed power laws. We cannot explain why the universal potential of endpoint overlap models has not been not widely recognized.

General Features
We are not aware of a previous focused effort to study the contributions of the endpoint overlap model in 2 → 2 hadron-hadron reactions. The power counting for both the SD and endpoint models are expedited by a simple observation. The amplitude of almost all contributions scales like the combination of form factor amplitudes. This observation is quite old [81], and developed for a different purpose, yet it is rather general. If there is another contribution with a qualitatively different behavior, its momentum flow will go by a qualitatively different topology.
In the high energy limit the differential cross section for 2 → 2 reactions with amplitude M is given by The composition of two form factors F A F B with a 1/Q 2 exchange kernel scales like F A F B Q 2 /Q 2 . The numerator factor of Q 2 accounts for the vector vertex factors not contained in F i . There is one significant difference between models, however. Multiple gluon exchanges in the SD model have no strong selection rules from the color singlet nature of hadrons. Scattering a single constituent in the endpoint overlap model requires at least two gluons in a singlet combination. It is well known that the box diagram of two gluon exchanges scales with just the same power of Q 2 as a single gluon exchange, times logarithmic factors that are exactly computable [82]. For our purposes multi-gluon exchanges are indistinguishable, and at most laboratory momentum transfers, probably necessary.
To proceed: For fixed angle kinematics s ∼ t, the endpoint overlap model composing form factors predicts: • For ππ → ππ, • For π p → π p, For pp → pp, • For 2 → 2 scattering of hadrons with n 1 and n 2 valence constituents, Agreement with experiment [5, 10, 74-79, 83, 84] are explained by the endpoint overlap model. Unlike the SD model, no approximations of an asymptotic character are needed. The approximations assume only that Q 2 >> Λ 2 , adds support to the valence state of the pion having two constituents, and the proton having three. As before, scattering constituents beyond the valence components is suppressed by powers of 1/s. The counting is different for the t dependence of amplitudes at fixed s >> GeV 2 . In that case MM * /s 2 ∼ MM * . With |t| << s the leading dependence replaces s → t, and multiplies the results above by t 2 . By far the most important example comes from pp → pp scattering, which displays a stunning experimental dependence falling like t −8 [66,[85][86][87]: exactly the endpointoverlap contribution.

Discussion
We mentioned that a qualitatively different momentum flow could change the counting. The independent scattering model [66] is usually highlighted to explain the t −8 dependence. Landshoff had earlier found the model by not making the same assumptions of the model of Brodsky and Farrar. The independent scattering (IS) model gets its power law partly from the phase space of fast quarks with x ∼ 1/3 to overlap with the wave function in the final state. There are three (3) hard vector exchanges, suppressing the amplitude by corresponding powers of 1/t. In comparison the endpoint overlap contribution uses one (1) hard vector exchange, while obtaining the same powers of 1/t from the probability to find two quarks near the endpoint.
The phenomenology of complex phases and spin dependence are very similar for the IS and endpoint models. When treated in Fock-basis perturbation theory both model have similar Su-P 2 P' 2 =P 2 +q l 1 l' 1 P 1 P' 1 =P 1 -q k 1 k' 1 k 1 -r l 1 +r Figure 6: Proton-proton elastic scattering with a minimal hard-scattering kernel. dakov factors. Such factors may well explain the oscillations seen in pp fixed angle scattering and color transparency. Based on the results of Mueller [88], we conjecture that saddle point interpolation between the endpoint model and the SD model will be the dominant asymptotic amplitude. This is because the Sudakov suppression of one fast quark is less severe than the three fast quarks of the IS model.

Supporting Calculation
As with the form factors, we believe that supporting calculations are at least as important as general arguments.
Consider proton-proton scattering of Fig.6, p(P 1 ) + p(P 2 ) → p(P 1 ) + p(P 2 ) in the limit s = (P 1 + P 2 ) 2 ∼ t = (P 1 − P 1 ) 2 = q 2 = −Q 2 . In the center of mass frame the momenta of the two incoming particles are P 1 =(p, 0, 0, −p), The amplitude for the scattering diagram is given by As with the case of proton form factor, the interaction vertex will have delta functions enforcing the conservation of momentum in the quark interactions, which are implicit in the above expression. Extract the integral over the free momentum r, given by Integration is performed using Feynman parametrization. Simplifying the denominator using Feynman parameters(a i ) leads to a denominator D given by, where l = r − a 1 k 1 − a 3 q + a 4 l 1 ; We may neglect terms of the form k 2 1 , l 2 1 assuming the quarks are nearly light like. Terms of the form k 1 · q, l 1 · q, k 1 · l 1 are of the same order as Q 2 , assuming t ∼ s. Thus the dominant contribution in the ∆ goes like Q 2 .
Terms in the numerator of the form l µ k ν 1 , l µ q ν , l µ l ν 1 . . . vanish upon integration. The other terms can be integrated using the standard substitution In comparison, the other terms in the numerator scale like 1/Q 4 . To leading power we keep the l µ l ν term. The amplitude is given by Here The calculation can be simplified by a Lorentz transformation to a frame where the momenta of the protons becomes equivalent to the momenta of Eq. (2). For example Such a transformation will allow the use of results of the proton form factor calculation. Substituting the wave function from Eq. (11), a single term from the wave function is sufficient to understand the behavior of this integral. We illustrate the term M V going like (/ PC) αβ (γ 5 N) γ V, which is From the calculations following Eq. (14), the integrations become Calculate the cross section using The leading term from simplifying the trace goes like p 4 . Using s ∼ t we find

Concluding Remarks
We have shown that the endpoint overlap model stands as a comprehensive theory of hadronic reactions at large momentum transfer. It explains the observed experimental regularities in all cases we have investigated. The history of endpoint dominance is curious, and possibly explains why the model failed to be completely developed.
In 1970 Drell and Yan [12], and later West [13] (DYW) discussed a partonic model connecting hadronic form factors to deeply inelastic scattering. Using symbol η for the parton momentum fraction since called x, the central region 1 − η ≥ Λ/Q was found to predict a form factor falling too fast to agree with data. From this region F 1 (Q 2 ) ∼ g(Q 2 )/Q 2 , with g(Q 2 ) ∼ exp(−Q 2 /Λ 2 ) is expected. In comparison, the endpoint region 1 − η ≤ Λ/Q was observed to predict F 1 (Q 2 ) ∼ (1/Q 2 ) (p+1)/2 , if the structure function νW 2 ∼ (1 − η) p . The value p = 1 was computed by Drell and Yan [12] in a prototype two constituent calculation. The value p = 1 was also noted as being too small to fit the data needing p ≥ 2. While the two components of the toy model were a pion and a nucleon, the 1970 calculation in our view constitutes the prototype "constituent counting" relation connecting the scaling power with the number of constituents. The paper appears approximately five years before the papers of Brodsky and Farrar [1] and Matveev et al. [2], which found counting rules on a different basis.
Subsequently many workers noticed that a calculation of the endpoint contribution would be revised by a power of 1/Q 2 for each additional constituent added to a given process. Many workers also concluded that the endpoint contribution might be dominating the calculation of their particular process. Yet the endpoint region never saw anything like the degree of development of the short-distance perturbative model. For reasons we cannot explain, we cannot find a reference strongly advocating for the endpoint region, and developing it as a "comprehensive theoretical picture" that explains the observed power law dependence.
There exists a possible explanation coming from the drive of the early era. That time was concerned with testing the Lagrangian of QCD, and without needing to know wave functions. Exclusive reactions were hardly a good testing ground. Imagine trying to test perturbative quantum electrodynamics in a world where the Hydrogen atom bound states had not been solved. In such a Universe, calculations actually depending on unknown wave functions would be "bad." That is, they would be bad for establishing a Lagrangian by perturbation theory. Calculations for Hydrogen-Hydrogen scattering not depending on unknown wave functions would be very difficult to concoct. However, once any scheme self-consistent with perturbation theory was found, it would be "good" whether or not it was incomplete and "wrong." We believe that the lure of a strictly perturbative procedure caused a false perception that "correct physics" could only depend on operator product expansion moments of wave functions. The operator product expansion is an ansatz of great power when it applies, while creating huge gaps when it does not, but this was not obvious right away. Once the attitude was adopted, the opportunity to use data to learn about wave functions was rejected. As a result, the opportunity to actively use data to learn about hadron structure remains a relatively unexplored field. It has taken 30 years of more and more detailed calculations to find that the tiny integration region of the short-distance model can at most be relevant for extremely large Q 2 which need momentum transfers that are tens, hundreds, or thousands of times larger than laboratory scales to be relevant. For emphasis, we do not know of a single calculation that strongly supports the numerical dominance of the short-distance region, relentless advocacy notwithstanding. Indeed the first step towards arranging for short distance dominance has been to banish the endpoint region as perturbatively inconsistent, which (we maintain) is a signal of a concept error.
Through the entire period the endpoint contribution has never gone away. The time has come to accept endpoint contribution, and explore it further.

Appendix: Soft Gluon Effects
Soft gluon effects are an intrinsic difficulty of the SD models. When the effects of perturbation theory produce large corrections they indicate that the first approximations were not dynamically stable. We mentioned that soft gluon effects are not intrinsically present in the non-perturbative quantum mechanical basis we use. Whether or not they are added they do not change the power of Q 2 . We demonstrate this explicitly here by considering a particular two gluon exchange diagram. Consider, as an example, the amplitude shown in Fig. 7. The hard scattering contributions of such diagrams have been analyzed in [56]. We will extract the |Q| dependence of the amplitude in the endpoint region. The momenta of the virtual fermions and gluons are p g 1 = k 3 − k 3 , p f 1 = k 2 + k 3 − k 3 , p g 2 = P − P − k 1 + k 1 , p f 2 = P − P + k 1 . Considering only the endpoint region, it is understood that the momenta transferred to the spectator fermions is soft, hence the gluons and the spectator fermions both have low momenta. One of the terms in the amplitude is Hence we obtain the expected momentum dependence.