Spin and Polarization in High-Energy Hadron-Hadron and Lepton-Hadron Scattering

The role of spin degrees of freedom in high-energy hadron-hadron and lepton-hadron scattering is reviewed with emphasis on the dominant role of soft, diffractive, non-perturbative effects. Explicit models based on analyticity and Regge-pole theory, including the pomeron trajectory (gluon exchange in the $t$ channel) are discussed. We argue that there is a single, universal pomeron in Nature, manifest as relatively"soft"or"hard", depending on the kinematics considered. Regge trajectories are non-linear, complex functions, predicting a finite number of resonances. Explicit models and fits are presented that may be used in recovering generalized parton distributions from deeply virtual Compton scattering and electoproduction of vector mesons.

The paper is organized as follows. In Section II we discuss elastic N N scattering with spin. Recall that the dominant point of view is that the role of spin decreases as energy increases, and it can be ignored at multi-GeV energies, i.e., beyond those of the ISR, and even more so at the LHC. Spin effect may be important in the region of the diffraction minimum, i.e., around t = −1 GeV 2 . The point is that the origin and mechanism of the dip-bump phenomenon is still disputable. There are many models but no theoretical understanding of the origin of this important phenomenon. Predictions, mostly based on variants of the Regge-eikonal approach, failed in predicting the position and depth of the dip at the LHC. Spin effects are not the favourite but still a viable possibility, alternative to the dominant ones, based on unitarity (Section II C).
The importance and properties of the Regge trajectories are repeatedly stressed throughout this paper, with a dedicated Section II B. Dual models have shown that Regge trajectories are sort of dynamical variables. The trajectories are non-linear, complex function. Their thresholds and asymptotic behaviour should satisfy known constraints. For example, the lowest, 4m π threshold affects the behaviour of the differential cross section in the Coulomb interference region, competing with possible spin effects.
In Section III we focus on the procedure of recovering Generalized Parton Distributions (GPDs) from electron-proton scattering with spin, measured at the JLab 3 . Here again the pomeron exchange is a basic element of the theory. Different from hadron-hadron scattering, dominated by a "soft" pomeron exchange, here the pomeron is "hard", i.e., its intercept is much larger (α(0) ≈ 1.3) than that considered in Section II. In Section IV we introduce a model for Deep Inelastic Scattering (DIS) amplitude with Q 2 -dependent pomeron intercept. That model will be generalized to deeply virtual Compton scattering (DVCS) in Section V, applicable both to "soft", e.g., N N scattering Section II, and "hard" DVCS.
Various aspects related to the origin of nucleon's spin are discussed in Section VII. We share the point of view that there is no "crisis", over-dramatized after the EMS measurements. The proton spin 1/2 can be collected from that of its constituents and their motion.
Finally, let us mention recent findings in spin physics revealed in ultra-relativistic heavy-ion collisions; STAR collaboration 4 discovered a significantly nonzero global polarization of Lambda hyperons produced in non-central Au-Au collisions in the RHIC Beam Energy Scan (BES) Program 4 . Different hydrodynamic models [5][6][7] generally reproduce the magnitude of the measured polarization. In the hydrodynamic models, the Lambda hyperons produced at "particlization" (fluid to particle transition) hypersurface acquire polarization via a thermodynamic spin-vorticity coupling mechanism. This effect is interesting as a possible manifestation of the most vortic fluid ever made.
Interestingly, this is not the first case that Λ polarization raises interest in the high energy community. In the nineteens of the past century, much discussion was related 8 to the origin of polarization in the scattering of unpolarized particles. According to Ref. 9 polarization of inclusively produced lambdas arises due to the spin-orbit interaction in a scalar field in which quarks recombine into hadrons.
The present paper to a large extent reflects the author's vision of the spin effects in the "soft" region and his personal contribution in the field. For a much wider panorama we recommend to the reader the excellent papers 2,10,11 .

II. HADRON-HADRON SCATTERING
In this section we collect the basic notions and rules indispensable in describing hadronic reactions. We restrict our discussion of this vast field to elastic nucleon scattering in the nearly forward region, dominated by the diffraction cone, and nominated "soft" or "non-perturbative", implying that perturbative QCD methods here are not applicable. Instead, the basic tools are those based on the analytic S matrix, namely, dispersion relations, unitarity and Regge poles. We relate these means in the context of the subject of the present paper-spin and polarization effects at high energies. Below is a short overview of the relevant tools.
Soft events are described by Regge pole models, sometimes appended by "QCD-inspirations". In most of the papers on the subject, a spin-independent invariant scattering amplitude is used at high-energies. Since spin effects cannot be excluded a priori, below we briefly present an introduction to the relevant formalism.
In the Regge pole model, the trajectories contain most of the basic information on the dynamics. Their form is constrained by unitarity and analyticity, and is intimately connected by the Chew-Frautchi plot to the spectrum of hadrons. This is a central issue in high energy physics, since it is related to the problem of confinement: how do heavy resonances melt producing a quark-gluon soup? In our approach, the real part of Regge trajectories is limited, implying a finite number of resonances in Nature. This issue is related to the Hagedorn spectrum and "ultimate temperature", now interpreted as a the temperature of the phase transition between hadrons and quark-gluon plasma 12 .
Unitarity is both important and technically difficult to be satisfied. Various approaches and approximations to unitarity are presented and compared in this Section. Let us mention that the final result (construction of a viable scattering amplitude) depends both on the input and subsequent unitarization procedure. The better the input ("born term"), the better are chances for its subsequent successful unitarization: unitarity correction to a reasonable input (for example, a dipole pomeron, reproducing itself under unitarization) are small 13 .

A. Elastic Proton-Proton Scattering
Elastic proton-proton scattering is described by five helicity amplitudes 14,15 , that in the Regge limit, t fixed and s → ∞, can be expressed in terms of the Regge-pole contributions as [16][17][18][19][20] where α i ≡ α i (t) are Regge trajectories. In high-energy elastic scattering the dominant trajectory is that of the pomeron. In the simplest case, it is linear, parametrized, for simplicity as α(t) = α(0) + α t. We will use also advanced models of Regge trajectories, to be introduced in the next Section II B. Recently 21 it was used to predict glue and oddballs-resonances made of two or three gluons. The above helicity amplitudes differ only by powers of √ t leaving the Regge propagator (i.e., nature of the Regge pole) intact. This means that the addition of spin-flip amplitudes by themselves will not provide the dip mechanism, although they will slightly modify it. The (single) diffraction minimum followed by maximum was shown 13 to arise even at the Born level by choosing the pomeron to be a dipole (double Regge pole or dipole pomeron (DP) 13 -an interesting and unique alternative to a simple Regge pole). Fits to the the data on pp scattering, including those from the LHC (TOTEM Collaboration) using DP was performed recently 22 . The global fit to all energies, including the dip-bump is good, although it leaves room for further perfection. Inclusion of the spin-flip component may improve the situation.
In most of the relevant papers the dip-bump structure is generated by unitarity corrections, see Section II C, although there is no consensus on the details.
The differential cross sections is given by The helicity amplitudes can be written as comes from the strong interactions, Φ em i (s, t) from the electromagnetic interactions and ϕ(s, t) is their interference. The spin correlation parameters, the analyzing power-A N and the and double-spin parameter A N N , can be extracted from experimental measurements: where ∆σ s and ∆σ d refer to the difference of single-and double-spin-flip cross sections. The expressions for these parameters are

B. Regge Trajectories
Unitarity imposes a severe constraint on the threshold behaviour of the trajectories, as shown on the diagram, Figure 1: while asymptotically the trajectories are constrained by The above asymptotic constraint can be still lowered to a logarithm by imposing wide-angle power behaviour for the amplitude. The trajectory satisfying the above constraints is: where t 0 = 4m 2 π for the pomeron and t 0 = 9m 2 π for the odderon and δ, α 1 , α 2 are parameters fitted to the data with the obvious constraints: α(0) ≈ 1.08 and α (0) ≈ 0.3 (in the case of the pomeron trajectory). Trajectory (9) has square-root asymptotic behaviour, in accordance with the requirements of the analytic S-matrix theory.
The contribution of the lowest threshold in the t channel of the trajectory and the scattering amplitude is shown on the diagram, Figure 1. The presence of this threshold produces a "break" in the diffraction cone and generates the pion cloud of the nucleon, as shown in Figure 2.
For t >> t 0 , |α(t)| → α1 α2 |t|. For t > t 0 (on the upper edge of the cut), mα > 0. The intercept is α(0) = 1 + δ and the slope at t = 0 is Resonances terminate, i.e., they do not appear beyond the maximum value or the finite real part of a given Regge trajectory.
For a recent application of this trajectory to glueball spectroscopy see Ref. 21 .

C. Unitarity
The best known unitarization procedure is that of Regge-eikonal, where the input ("Born term") is a Reggepole amplitude.
It implies two Fourier-Bessel transforms (integrals), one from the Mandelstam variable t (momentum space) to the impact parameter b (coordinate space). While direct integration is usually performed analytically (trivial for linear Regge trajectories and still simple for a single square-root threshold, see Ref. 23 , the inverse integral is much more cumbersome, and can be done only numerically. The technicalities become even more complicated when spin degrees are included (see below).
A simple unitarization procedures is by the use 18,19 of the equation: where y = log(s/s 0 ) and ∆ = 1 − α(0), (α(0) is the intercept of the leading Regge pole trajectory. Its solution is where χ(s, b) ≈ s ∆ is connected with the Born term of the scattering amplitude. The scattering amplitude is The phase χ(s, b) is connected to the interaction potential: If the potential contains a non-spin-flip part and, for example, spin-orbital and spin-spin interactions, the phase is: If however only the spin-flip and spin-non-flip parts are accounted for, the overlap function is Using the representation for the Bessel functions the representation of spin-non-flip and spin flip amplitude becomes

"U -matrix" Unitarization
With an extra coefficient n one gets the so-called U -matrix unitarization form, see 16 and references therein. This way of unitarization, developed mainly in Serpukhov, Dubna and Kiev, is less familiar than the eikonal one, nevertheless it results in a number of interesting predictions, among which is the so-called reflective scattering 24 , resulting in an increasing ratio σ el (s)/σ tot (s) at the LHC. The unorthodox predictions of the U matrix approach may bringing new ideas in the complicated field of "soft physics".
In the impact parameter representation, the properties of the U -matrix were explored in Ref. 16 . In that approach, the hadronic amplitude is given by where χ(s, b) is the same Born amplitude as before.
Comparing Equation (21) with Equation (13), we see that they are of similar rational form, just they differ by the additional coefficient in the denominator. This additional coefficient leads to different analytic properties: the upper bound at which the overlapping function saturates will be twice that compared with the eikonal or the U e representations, and the inelastic overlap function at b = 0 tends to zero at high energies. Furthermore, σ el /σ tot → 1.

Eikonal
To get the standard eikonal representation of the elastic scattering amplitude in the impact parameter representation one uses the non-linear equation where y = log(s/s 0 ) and the subscript "e" implies that the solution N e is of the standard eikonal form The eikonal representation is then where: With account for Equations (17) and (27), one gets respectively for spin non-flip and spin-flip where For Gaussian potentials V 0 and V 1 in the first Born approximation, Φ h 0 and Φ h 1 take the forms In view of these technical complications, of interest is the unique case of a dipole pomeron (DP), reproducing itself under eikonal corrections, as noticed in Ref. 25 , see also the Appendix in 26 . In other words, DP in the impact parameters representation remains stable (intact) under eikonalization. In addition, a DP produces logarithmically rising cross sections at unit pomeron intercept and obeys geometrical scaling, making it attractive as an efficient model of diffraction. Note that higher order Regge poles violate unitarity.
FIG. 2: (a) Typical shape of the high-energy differential (diffractive) cross section, and (b) Its impact-parameter image. A "dictionary" relating the shape of the diffraction cone (in t), left icon to its impact-paramter image (rignt), and highlighting their two silent features: the "break", origin of nucleon's atmosphere and the dip-bump, resulting from absorptions (unitarity!) at small values of the impact parameter, see 13 .

III. NUCLEAR STRUCTURE: FROM DEEPLY VIRTUAL COMPTON SCATTERING (DVCS) TO GENERALIZED PARTON DISTRIBUTIONS (GPDS)
An important part of spin physics is lepton-hadron scattering. Deeply virtual exclusive ep processes provide an important tool in accessing the generalized parton distributions (GPDs) [27][28][29] . The underlying mechanism is Regge-pole exchange in the t-channel. Based on factorization theorems 30,31 , GPDs offer a partonic interpretation of these processes, where unobserved transverse degrees of freedom are integrated out. Thereby, these universal functions, defined in terms of matrix elements of quark and gluon operators or, alternatively, as a non-diagonal overlap of light-cone wave functions 32,33 , encode the non-perturbative aspects of the nucleon 34,35 . In particular, GPDs provide access to the transverse spatial distribution of patrons [36][37][38] , and to the decomposition of the nucleon spin in terms of quark and gluon degrees of freedom 39 .
Phenomenologically, exclusive electroproduction of a real photon, DVCS, depicted in Figure 3 (left), is the golden channel to constrain GPDs as it is theoretically clean and the phase of its amplitude can be measured using the interference with the Bethe-Heitler (BH) amplitude (see Figure 3; cf. Figures 4 and 5).
More details, especially those related to experiments can be found in 40,41 .

A. Deeply Virtual Compton Scattering
The differential photon electroproduction cross section is the sum of the BH amplitude squared, DVCS amplitude squared, and the interference (INT) terms (see Figure 4) Here the +(−) signs correspond to the electron (positron) beam, x B is the common Bjorken scaling variable, φ is the azimuthal angle between lepton and hadron scattering planes, and ϕ = Φ − φ, where Φ(≡ φ S ) is the angle between the , and neglecting the electron mass, the three terms on the r.h.s. of (34) are known in terms of the electromagnetic Pauli form factor F 1 (t) and the Dirac form factor F 2 (t), parameterizing the BH amplitude, and a set of 12 photon helicity dependent CFFs F ab (x B , t, Q 2 ), parameterizing the DVCS amplitude, see Figure 3. These CFFs are labelled by the helicities of the incoming a ∈ {+, 0, −} and outgoing

B. Relating DVCS Observables to GPDs
GPDs, denoted as besides depending on the partonic momentum fraction x and the momentum transfer squared t, depend also on the t-channel longitudinal momentum fraction η, called skewness (often denoted by ξ in the literature), and on the factorization scale µ 2 . The unpolarized parton GPDs are called H and E 42 , where the former (latter) GPD can be loosely associated with a proton helicity (non)conserved distribution. DVCS observables can be evaluated in terms of the helicity CFFs (35). To express them in terms of GPDs, it is appropriate to utilize a conventionally defined GPD-inspired CFF basis, such as the one introduced in 43 : Here, the CFFs H, E, H, and E are associated with twist-two GPDs F ∈ {H, E, H, E} and govern the photon helicity non-flip DVCS amplitude, i.e., at leading twist-two accuracy one has To LO they are calculated from the handbag diagram, depicted in Figure 6, yielding the convolution formula where e i are the fractional quark charges. The variable ξ ∼ x B /(2 − x B ) is a conventionally defined Bjorken-like scaling variable, equated to the longitudinal momentum fraction in the t-channel, and µ 2 ∼ Q 2 being the factorization scale. For twist-two and LO accuracy, and light quarks one adopts the conventions From Regge-pole models, consistent with the data, see Section IV below, the real part of the dominant CFF H in the small-and even moderate-x B region is much smaller than its imaginary part (at least for smaller values of −t). One concludes that the interference term is negligible and one can simplify the t-differential cross section to

IV. MODELLING DVCS
Below we present models of DVCS and vector meson production (VMP), potentially useful in calculating GPDs. Regge pole models provide an adequate framework to describe high-energy, low t scattering phenomena. Being part of the S matrix theory, however, strictly speaking, they are valid only for the scattering of on-mass-shall particles. 6: Factorization of the DVCS amplitude to leading order in perturbative QCD and to leading twist-two accuracy.
Still, the successful application of the Regge pole models in describing the off-mass-shall HERA data opened the way to their use in deep-inelastic scattering, DVCS and vector meson production at HERA. Off-shell extension was realized e.g., by calling the Q 2 -dependent Regge trajectories "effective" ones. A particularly simple and efficient Regge pole model 44 with Q 2 -dependent residues (vertices) is presented below, in Section IV.
In a more advanced Regge-pole model, Section V the Regge trajectories and the residues do not depend on virtuality. Instead, the amplitude contains two (or more) Regge-pole terms, whose relative weight depends on Q 2 , mimicking the multi-pole nature of the so-called QCD pomeron.

Simple Model of DVCS
By Regge factorization, Figure 7, the DVCS amplitude can be written as where A 0 is a normalization factor, V 1 (t, Q 2 ) is the γ * P γ vertex, V 2 (t) is the pP p vertex and α(t) is the exchanged pomeron trajectory, which we assume to be logarithmic: Such a trajectory is nearly linear for small |t|, thus reproducing the forward cone of the differential cross section, while its logarithmic asymptotic provides for the large-angle scaling behaviour, typical of hard collisions at small distances, with power-law fall-off in |t|, compatible with the quark counting rules. Here we are referring to the dominant pomeron contribution eventually appended by a secondary trajectory, e.g., the f -Reggeon. For convenience, and following the arguments based on duality, the t dependence of the pP p vertex is introduced via the α(t) trajectory: V 2 (t) = e bα(t) where b is a parameter. A generalization of this concept will be applied also to the upper, γ * P γ vertex by introducing the trajectory where the value of the parameter α 2 may be different in α(t) and β(z) (a relevant check will be possible when more data are available). Hence the scattering amplitude, with correct signature, becomes where L ≡ ln(−is/s 0 ). The model contains a limited number of parameters. Moreover, most of them can be estimated a priori. The product α 1 α 2 is just the forward slope α of the Reggeon (≈0.2 GeV −2 for the pomeron, but much higher for f and/or for an effective Reggeon). The value of α 1 can be estimated from the wide-angle quark counting rules. For large t (|t| >>1 GeV 2 ) the amplitude goes roughly as ∼ e −α1 ln(−t) = (−t) α1 , where the power α 1 is related to the number of quarks in a collision.
From Equation (44) the slope of the forward cone is which, in the forward limit, t = 0 reduces to Thus, the slope shows shrinkage in s and antishrinkage in Q 2 .
In the Q 2 → 0 limit the Equation (44) becomes where we recognize a typical Regge-behaved photoproduction (or, for Q 2 → m 2 H , on-shell hadronic (H)) amplitude. The related deep inelastic scattering structure function is recovered by setting Q 2 2 = Q 2 1 = Q 2 and t = 0, to get a typical elastic virtual forward Compton scattering amplitude: In the Bjorken limit, when both s and Q 2 are large and t = 0 (with x ≈ Q 2 /s valid for large s), the structure function is given by: where α e is the electromagnetic coupling constant and the normalization is σ t (s) = 4π s A(s, Q 2 ). The resulting structure function has correct (required by gauge invariance) Q 2 → 0 limit and approximate scaling (in x) behavior for large enough s and Q 2 .

V. REGGEOMETRY
In a series of papers (see Refs. 45,46 ), photon virtuality was incorporated in a "geometrical" way, reflecting the observed trend in the decrease of the forward slope as a function of Q 2 . This geometrical approach, combined with the Regge-pole model, was named "Reg-ge-ometry" (a wordplay, pun). A Reggeometric amplitude dominated by a single pomeron shows reasonable agreement with the HERA data on VMP and DVCS, when fitted separately to each reaction.
As further step, to reproduce the observed trend of hardening as Q 2 increases, and following Donnachie and Landshoff 47,48 , a two-term amplitude, characterized by a two-component-"soft" + "hard"-pomeron, was suggested 45 . We stress that the pomeron is unique, but we construct it as a sum of two terms. Then, the amplitude is defined as (s = W 2 is the square of the c.m.s. energy), such that the relative weight of the two terms changes with Q 2 in the right way, i.e., the ratio r = A h /A s increases as the reaction becomes "harder" and v.v. It is interesting to note that this trend is not guaranteed "automatically": both the "scaling" model 44 or the Reggeometric one 45,46 show the opposite tendency, that may not be merely an accident and whose reason should be better understood. This "wrong" trend can and should be corrected, and in fact it was corrected 47,48 by means of additional Q 2 -dependent factors H i ( Q 2 ), i = s, h modifying the Q 2 dependence of the amplitude, in a such way as to provide increasing of the weight of the hard component with increasing Q 2 . To avoid conflict with unitarity, the rise with Q 2 of the hard component is finite (or moderate), and it terminates at some saturation scale, whose value is determined phenomenologically. In other words, the "hard" component, invisible at small Q 2 , gradually takes over as Q 2 increases. An explicit example of these functions is presented below.
Recall that the invariant scattering amplitude is defined as where is the linear pomeron trajectory (the use of more advanced, non-linear trajectories of Section II B is straightforward), a and b are two parameters to be determined from the fitting procedure and m N is the nucleon mass. The coefficient H is a function providing the right behavior of elastic cross sections in Q 2 : where A 0 is a normalization factor, Q 2 0 is a scale for the virtuality and n s is a real positive number. In this model one uses an effective pomeron, which can be "soft" or "hard", depending on the reaction and/or kinematic region defining its "hardness". In other words, the values of the parameters α 0 and α must be fitted to each set of the data. Apart from α 0 and α , the model contains five more sets of free parameters, different in each reaction, as shown in Table I. The exponent in the exponential factor in Equation (51) the differential and integrated elastic cross sections become, and where The results of the fits are presented in Table I. The "mass parameter" for DVCS was set to M = 0 GeV, therefore in this case Q 2 = Q 2 . Each type of reaction was fitted separately. As it can be seen from the right plot of Figure 8, the single-term model fails to fit both the high-and low-|t| regions properly, especially when soft (photoproduction or low Q 2 ) and hard (electroproduction or high Q 2 ) regions are considered. One of the problems of the single-term Reggeometric pomeron model, Equation (51), is that the fitted parameters in this model acquire particular values for each reaction, The VMP results clearly show hardening of the pomeron in the change of α 0 and α when going from light to heavy vector mesons.
It is also interesting to note that the effective pomeron trajectory for DVCS (α 0 = 1.23, α = 0.04, see Table I), contrary to expectations is rather "hard".  Table I.

Two-Component Reggeometric Pomeron
Now we introduce the universal, "soft" and "hard", pomeron model. Using the Reggeometric ansatz of Equation (51), we write the amplitude as a sum of two parts, corresponding to the "soft" and "hard" components of a universal, unique pomeron: Here s 0s and s 0h are squared energy scales, and a i and b i , with i = s, h, are parameters to be determined with the fitting procedure. The two coefficients H s and H h are functions similar to those defined in Refs. 47,48 : where A s and A h are normalization factors, Q 2 s and Q 2 h are scales for the virtuality, n s and n h are real positive numbers. Each component of Equation (56) has its own, "soft" or "hard", Regge (here pomeron) trajectory: The terms "soft" and "hard" may be misleading, alluding to the wide-spead nomenclature by which the "soft" pomeron is associated with the traditional Regge-pole theory, while the "hard" one is the BFKL pomeron derived from QCD, having little to do with the former, by meaning that they are different object. Instead, we insist that there is only one pomeron, although it may be complicated, consisting of more terms etc. In our case the two terms refer to a single object, THE pomeron, while its (two) parts are weighted with virtuality Q 2 .
The "pomeron" amplitude (56) is unique, valid for all diffractive reactions, its "softness" or "hardness" depending on the relative Q 2 -dependent weight of the two components, governed by the relevant factors H s ( Q 2 ) and H h ( Q 2 ).
To reduce the number of free parameters, we simplified the model, by fixing a s,h = 0 and substituting the exponent 2 a s,h (56). The proper variation with Q 2 will be provided by the factors H s ( Q 2 ) and Consequently, the scattering amplitude assumes the form The "Reggeometric" combination 2 a s,h was important for the description of the slope B(Q 2 ) within the single-term pomeron model (see previous Section), but in the case of two terms the Q 2 -dependence of B can be reproduced without this extra combination, since each term in the amplitude (58) has its own Q 2 -dependent factor H s,h (Q 2 ). By using the amplitude (58) and Equation (53), we calculate the differential and elastic cross sections, by setting for simplicity s 0s = s 0h = s 0 , to obtain In these two equations we use the notations Notice that amplitude (58) can be rewritten in the form where the two exponential factors e bst−ns ln 1+ Q 2 can be interpreted as the product of the form factors of upper and lower vertices (see Figure 7c). Interestingly, the amplitude (61) resembles the scattering amplitude of Ref. 44 .

Fitting the Two-Component Pomeron to VMP and DVCS HERA Data
The fitting strategy is based on the minimization of the quantityχ 2 = 1 N N i=1χ 2 i , where N is the number of all reactions involved (i.e., ρ, φ, ω, J/ψ, Υ and γ production);χ i is the mean value of χ 2 for different types of data for selected class of reactions, defined asχ i = 1 Ni Ni k=1χ 2 k,i , whereχ 2 k,i is χ 2 k,i /d.o.f. for i-th class of reactions and k-th type of data, i.e., those relative to σ el (Q 2 ), σ el (W ) and dσ el (t)/dt; N i is number of different type of data for i-th class of reactions. The normalization parameters were fixed at and s 0 was set 1 GeV 2 . DVCS and VMP are similar in the sense that in both reactions a vector particle is produced. However there are differences between the two because of the vanishing rest mass of the produced real photon. The unified description of these two types of related reactions does not work by simply setting M γ = 0. From fits we found M ef f DV CS = 1.8 GeV and a normalization factor f DV CS = 0.091 follows.
Some results of the fits are shown in Figures 9-12 (for σ el (W )) and Figure 13 (for dσ el (t)/dt) for vector meson production, with the values of the fitted parameters given in Table ??. The mean value of the totalχ 2 (see above its definition) is equal to 0.986. The mean values ofχ 2 of the fit for different observables (i.e., σ el (Q 2 ), σ el (W ) or dσ el (t)/dt) and different reactions (VMP or DVCS), together with the numbers of degrees of freedom (number of data points) and the global mean valueχ 2 i , are shown in Table ??. In Table II the parameters of the two-component pomeron model (Equations (59) and (60)) fitted to the combined VMP and DVCS data are quoted, when the pomeron trajectories are fixed to α s (t) = 1.08 + 0.25t and α h (t) = 1.20 + 0.01t.
Next, by using Equation (59) with the values of the parameters from Table ?? and the formula one calculates the forward slopes and compares them with the experimental data on VMP, including those for the Ψ(2S) production. A compilation of all results is presented in Figure 14.
The number of fitted parameters of the two-component pomeron model (Equation (62)) is 12 ( Table ??), with five additional normalization factorsfor six vector particle productions (ρ 0 , φ, ω, J/ψ, Υ and γ). By fixing the pomeron trajectories α s (t) and α h (t) (see Table II), the number of free parameters reduces to 8. More results and discussion on the fitting details can be found in Refs. 45,46 .

VI. BALANCING BETWEEN "SOFT" AND "HARD" DYNAMICS
In this section we illustrate the important and delicate interplay between the "soft" and "hard" components of our unique pomeron. Since the amplitude consists of two parts, according to the definition (50), it can be written as As a consequence, the differential and elastic cross sections contain also an interference term between "soft" and "hard" parts, so that they read and σ el = σ s,el + σ h,el + σ interf,el , according to Equations (59) and (60), respectively. Given Equations (65) and (66), we can define the following ratios for each component: and where i stands for {s, h, interf }. Figure 15 shows the interplay between the components for both σ i,el and R i ( Q 2 , t), as functions of Q 2 , for W = 70 GeV. In Figure 16 both plots show that not only Q 2 is the parameter defining softness or hardness of the process, but such is also the combination of Q 2 and t, similar to the variable z = t − Q 2 introduced in Ref. 44 . On the whole, it can be seen from the plots that the soft component dominates in the region of low Q 2 and low |t|, while the hard component dominates high Q 2 and high |t|.
In other words, the variables t and Q 2 have much in common: both refer to squared momentum transfer and have the meaning of "hardness". The present model provides a valuable laboratory to study this important issue.
Hadron-induced reactions, discussed in Section II differ from those induced by photons at least in two aspects. First, hadrons are on the mass shell and hence the relevant processes are typically "soft". Secondly, the mass of incoming hadrons is positive, while the virtual photon has negative squared "mass". Our attempt to include hadron-hadron scattering into the analysis with our model has the following motivations: (a) by vector meson dominance (VMD) the photon behaves partly as a meson, therefore meson-baryon (and more generally, hadron-hadron) scattering has much in common with photon-induced reactions. Deviations from VMD may be accounted for the proper Q 2 dependence of the amplitude (as we do hope is in our case!); (b) connection between space-and time-like reactions; (c) according to some claims the highest-energy (LHC) proton-proton scattering data indicate the need for a "hard" component in the pomeron (to anticipate, our fits do not found support the need of any noticeable "hard" component in pp scattering).
We do not intend to perform here a high-quality fit to the pp data; that would be impossible without the inclusion of subleading contributions and/or the odderon.
The pp scattering amplitude is written in the form similar to the amplitude (58) for VMP or DVCS, the only difference being that the normalization factor is constant since the pp scattering amplitude does not depend on Q 2 : We fix the parameters of pomeron trajectory at We conclude that, while the data on total cross section are compatible with a small "hard" admixture in the amplitude, the slope parameter with a hard component included seems to manifest a wrong tendency, by slowing down with increasing energy, while the TOTEM measurements show the opposite.
The resulting fits are reasonable, despite the following open problems: • sub-leading Regge contributions must be included in any extension of the model to lower energies (below 30 GeV); • the Q 2 dependence of the scattering amplitude, introduced empirically has to be compared with the results of unitarization and/or QCD evolution.
• as seen from Section VI, the "soft" component of the pomeron dominates in the region of small t and small Q 2 . Hence, any parameter responsible for the "softness" and/or "hardness" of processes, should be a combination of t and Q 2 . A simple solution was suggested in Ref. 44 with the introduction of the variable z = t − Q 2 . The interplay of these two variables remains an important open problem that requires further investigation.
The extension of our formalism to hadronic reactions (pp scattering) shows that the available data can be will described by a single-soft-component.

VII. SPIN OF THE PROTON IN TERMS OF ITS CONSTITUENTS
The origin of spin in the nucleon is a rich subject deserving a dedicated presentation. Here I only briefly sketch some recent developments in the field. A complete and comprehensive treatment of the subject can be found e.g., in Refs. 2,10 .
In the simplest case, if only valence quarks contributed, the nucleon spin would be simple a sum where ∆Σ is the quark spin contribution and L q is the orbital angular momentum (OAM) contribution. Partition between the quark spin and the OAM for a long time were subject of debates generally called "spin crisis". Estimates of shares varied in a wide range. Among the reasons of disagreement and contradictions was the role of gluons, and sea as compared to that of the valence quarks. The phenomenon known as "spin crisis" arose from an EMC experiment at CERN 1 quoting Σ ≈ 0 (with large uncertainty), thus allegedly contradicting the naive quark model. Recent COMPAS, HERMES and Jlab measurements are consistent with a value of Σ ≈ 0.3. Actually, J must obey conservation of the total angular momentum known as Jaffe-Manohar sum rule 49 , An important task in the context of the above some rule is the separation of sea and valence quarks, indistinguishable in DIS observables. In addition, as noted in Ref. 50 , strange and anti-strange sea quarks can contribute differently in to the nucleon spin.
Paper 2 discusses also how DIS data on proton, neutron and deuteron targets can be used to separate the contributions from different quark polarizations assuming validity of SU (3) f .
That paper contains many more ideas, among which is the use of parton-hadron duality of Bloom and Gilman 51another effective tool to be used in studies of spin effects and a possible key to understand better problems related to quark confinement. modified to In the naive parton model, for a fast moving proton with helicity +1/2, a 0 = ∆Σ ≡ (∆u + ∆u) + (∆d + ∆d) + (∆s + ∆s) where the ∆q, ∆q are the first moments of the polarized quark-parton helicity densities. Hence, bearing in mind that ∆q = q + − q − , where ± corresponds to the number densities with spin along or opposite to the proton's momentum, one obtains in the naive parton model and if there is no other source of angular momentum one expects implying, naively, After the famous EMC experiments revealed that only a small fraction of the nucleon spin is due to quark spins, there has been great interest in 'solving the spin puzzle', i.e., in decomposing the nucleon spin into contributions from quark/gluon spin and orbital degrees of freedom.
The EMC experiment gave a 0 ≈ 0 and later experiments confirmed that a 0 1, giving rise to the spin crisis in the (naive) parton model. However, Equation (VII) cannot possibly be true because the right hand side is a fixed number, whereas the left hand side is, beyond the naive level, equal to a 0 (Q 2 ), i.e., a function of Q 2 ! Thus failure of Equation (VII) to hold cannot be used to infer that there is crisis. It is obvious that a correct relation between the spin of a nucleon and the angular momentum of its constituents should include their orbital angular momentum and should also include a contribution from the gluons.

Ji's and J-M's Decompositions
Finally, we briefly mention two sum rules visualising in different albeit complementary ways the spin decomposition in the nucleon.
The Ji decomposition 39 appears to be very useful, as not only the quark spin contributions ∆q but also the quark total angular momenta J q ≡ 1 2 ∆q + L z q (and by subtracting the spin piece also the the quark orbital angular momenta L z q ) entering this decomposition can be accessed experimentally, through generalized parton distributions (GPDs). The terms in (78) are defined as expectation values of the corresponding terms in the angular momentum tensor in a nucleon state with zero momentum. Here i D = i ∂ − g A is the gauge-covariant derivative. Jaffe and Manohar proposed 49 an alternative decomposition of the nucleon spin, which does have a partonic interpretation whose terms are defined as matrix elements of the corresponding terms in the +12 component of the angular momentum tensor In Figure 17 the share of various components to the nucleon spin, according to the two (Ji and Jaffe-Manohar) distribtutions is shown.
What appeared to be spin crisis in the parton model, over three decades years ago, was a consequence of a misinterpretation of the results of the famous European Muon Collaboration experiment on polarized deep inelastic scattering.

VIII. CONCLUSIONS
Interest in spin physics and polarization has moved from hadron-hadron to lepton-hadron processes. This is mainly due to the fact that in proton-proton and proton-antiproton scattering spectacular spin effects were/are not expected at highest energies (ISR, Fermilab, BNL or LHC), except maybe for the dip-bump structure at the diffraction cone. Various aspects of polarization in inelastic hadronic reactions reactions were treated in simple intuitive models in Refs. 9 .
Instead, deeply virtual Compton scattering (DVCS), with spin degrees of freedom is now at the focus of interest as a source of information on generalized parton distributions (GPDs).
In both cases (hadron-hadron and lepton-hadron) diffraction and the pomeron play a crucial role. Therefore much space was dedicated to the theoretical problems of diffractive scattering, properties of the vacuum Regge trajectory (pomeron), and confinement of quarks and gluons.
In the near future, research programs at the Electron-Ion Colliders (EICs) will play an important role in understanding gluon GPDs. Apart from providing constraints on the total quark/gluon contributions to the proton spin, the GPD's will provide important information on the nucleon tomography, for example, the 3D imaging of partons inside the proton. Together with the gravitational form factors extracted from the DVCS, this will deepen our understanding of the nucleon spin structure in return. The EIC may shed light on the quark/gluon orbital angular momentum (OAM) directly through various hard diffractive processes. Pioneer experimental effort to constrain the gravitation form factor from DVCS experiment at JLab has been carried out in Ref. 3 . The quark and gluon helicity contributions to the proton spin with the unique coverage in both x and Q 2 at EIC will provide the most stringent constraints on ∆Σ and ∆G 52,53 .
Last but not least, Regge trajectories Section II B, relating uniquely the spin and mass of particles, are building blocks of the theory, containing the basic information on the dynamics. By crossing symmetry, unitarity and duality they relate various aspects of high-energy strong interaction dynamics.
This research received no external funding