Generalized Parton Distributions: confining potential effects within AdS/QCD

Generalized Parton Distributions are investigated within a holographic approach where the string modes in the fifth dimension describe the nucleon in a bottom-up or AdS/QCD framework. The aim is to bring the AdS/QCD results in the realm of phenomenology in order to extract consequences and previsions. Two main aspects are studied: i) the role of the confining potential needed for breaking conformal invariance and introducing confinement (both: classic Soft-Wall and recent Infra-Red potentials are investigated); ii) the extension of the predicted GPDs to the entire range of off-forward kinematics by means of Double-Distributions. Higher Fock states are included describing the nucleon as a superposition of three valence quarks and quark-antiquark pairs and gluons.


I. INTRODUCTION
Generalized Parton Distribution functions (GPDs) are a source of fundamental information encoding essential aspects of the nucleon structure [1][2][3][4][5][6] as basic ingredients in the description of hard exclusive processes [2]. They are generalization of the well known parton distribution functions and, at the same time, as correlation functions they incorporate quite non trivial aspects of hadrons in the non-perturbative regime like: electromagnetic form factors, spin and angular momentum of the constituents and their spatial distribution [1,4]. Their functional structure is usually written as function of the longitudinal momentum fraction of the active quark (x), the momentum transferred in the longitudinal direction (ξ or skewedness) and the invariant momentum (square) t = −∆ 2 . The Fourier transform of GPDs (at ξ = 0) in the transverse direction encodes information on the partonic distributions in the transverse plane and it translates in a quantitative information (because of the probabilistic interpretation as density functions) on the separation of the struck quark and the center of momentum of the nucleon [3]. The detailed map of quarks and gluons in the nucleon interior is often called "nucleon tomography" since the traditional information from elastic and deepinelastic scattering provide static coordinates or momentum space pictures, separately, while GPDs provide pictures of dynamical correlations in both coordinate and momentum spaces [7].
Amplitudes of different hard exclusive processes (like deeply virtual Compton scattering (e.g. [8]), and virtual vector meson production (e.g. [9]) in the new generation of CLAS experiments at Jefferson Lab.) contain GPDs as essential components. On the other hand the experimental results already collected have shed a fundamental light on their role in different processes and kinematical regimes (for example the H1 [10] and ZEUS [11] at HERA, HERMES at DESY [12], Hall A and Hall B at Jefferson Lab. [13], COMPASS at CERN [14]).
GPDs are non-perturbative objects and their evaluation lies in the realm of non-perturbative QCD. The successes are, till now, strongly limited [15]. An alternative approach is the Holographic Light Front technique. Its fundamentals are in the correspondence between string theory developed in a higher dimensional ant-de-Sitter (AdS) space and conformal field theory (CFT) in Minkowski physical space-time [16][17][18]. Several consequent models have been constructed and they can be divided in top-down and bottom-up approaches. Starting from some brane configuration in string theory, one can, indeed, try to reproduce basic features of QCD following top-down paths (e.g. ref. [19]). On the way up one starts from low-energy properties of QCD (like chiral symmetry breaking and quark confinement) to infer elements for a gravity frame with asymptotically anti-de Sitter (AdS) space, the models are therefore indicated as AdS/QCD (e.g. ref. [20,21] and references therein). In particular within the bottom-up approach two successful models have been constructed: i) the Hard-Wall model which uses a sharp cut-off in the extra dimension to confine the (dual) hadron field [22,23]. The model is simple, analytic and appealing, but it does not reproduce the linear Regge behavior of the meson masses; ii) in the Soft-Wall model [24] a (quadratic) dilation field is added in the meson sector in order to successful reproduce the Regge behavior, however chiral symmetry breaking cannot be consistently realized. In particular it has been shown [25] that the spontaneous chiral symmetry breaking in vacuum and its restoration at finite temperature, can be realized only within a careful choice of the dilaton profile (see also ref. [26]).
Consequently several authors are investigating how to improve the SW description to incorporate the largest number of QCD properties [27][28][29][30][31].
An example particularly interesting in the present perspective is the Infra-Red improved soft-wall AdS/QCD model proposed in ref. [28]: it is constructed for baryons, taking into account a specific baryonic property of the spectrum, namely the parity-doublet pattern of the excited baryons. It shows consistent properties also in the meson sector [29]. This simplified model is taken, in the present paper, as a prototype to investigate GPDs and illustrating, at the same time, a procedure valid to study Generalized Parton Distributions and other observables in a generically modified confining potential.
Within the AdS/QCD approach Deep Inelastic Scattering (DIS) has been first addressed by Polchinski and Strassler in refs. [33], and GPDs have been investigated by many authors both within the Hard-Wall [34] and Soft-Wall [35,36] models. Because of the nature of the AdS-QCD analogy in the region of DIS, the results are restricted to the forward limit (ξ = 0) (cfr. Section III).
In the present work an attempt for a step forward is investigated and in two directions: i) generalizing the study of GPDs for confining potentials more complex than the simple Soft-Wall model; ii) extending the GPDs results to the off-forward region, ξ > 0, by means of a technique called Double-Distributions [37].
In Section II the procedure to evaluate the nucleon holographic wave function in the modified confining potential is discussed and the numerical results illustrated. In Section III the relation between sum rules and the ξ = 0 components of the GPDs is investigated and generalized to include, within a unified framework: i) the effects of the modified confining potential; ii) the contributions of higher Fock states. Numerical results for both helicity-independent and -dependent GPDs are discussed in Section IV and compared with a Light-Front approach. Section V is devoted to the application of Double-Distribution techniques [37] to the AdS/QCD predictions for the Soft-Wall model. It is shown how AdS/QCD can become predictive in the whole kinematical range (x, ξ > 0, t). Conclusions and perspectives in Section VI.

II. FROM THE SOFT-WALL TO THE INFRA-RED IMPROVED MODEL
The AdS/QCD framework relates a gravitationally interacting theory in the ant-de-Sitter space AdS d+1 with a conformal gauge theory in d-dimensions defined at the boundary. The needed breaking of conformal invariance (QCD is not a conformally invariant theory) of that correspondence for the baryonic case is obtained introducing, in addition to the dilaton term ϕ(z), an effective interaction ρ(z) in the action of the Dirac field (propagating in AdS d+1 ) [20,21]: Maximal symmetry is restored for ϕ(z) = ρ(z) = 0. One has (g) = R z d+1 while e M A is the inverse vielbein, e M A = z R δ M A . D M is the covariant derivative and the Dirac matrices anti-commute [Γ A , Γ B ] = 2η AB . A Dirac-like wave equation can be derived from Eq.(1) and the dynamical effect due to the dilaton field reabsorbed rescaling the spinor Ψ → e ϕ(z)/2 Ψ. For that reason the term e ϕ(z) is sufficient to break maximal symmetry for mesons but not for the baryon sector. The additional interaction term ρ(z) provides the needed breaking (and confining) contributions to generate the correct baryon spectrum [20,32]. The absence of dynamical effects of the dilaton background field has a particular disappointing side effect in the lack of guidance from gravity to solve the equations. A solution is given by a Light-Front holographic mapping where the LF wave equation can be identified with the equation of motion. In the case of d = 4, Γ A = (γ µ , iγ 5 ) and V (z) = R z ρ(z), the holographic variable z can be identified with the transverse impact variable ζ of the n − 1 spectator system with respect the active parton in a n-parton bound state (z = ζ). In 2 × 2 chiral spinor representation one obtains two coupled differential equations (e.g. ref. [20]) where ν can be identified with the light-front angular momentum, i.e. the relative angular momentum between the active parton and the spectator cluster. Eqs. (2) and (3) are easily reduced to the equivalent system of second order differential equations: A. Linear Soft-Wall potential For a quadratic interaction (and z = ζ within the holographic model), ρ(ζ) ∼ ζ 2 , V (ζ) = α 2 ζ (the so called Soft-Wall linear potential) and Eqs.(4) (5) become:  [28,29] (see Eq. (11) and Section II B for comments). Also a hard-wall potential at ζ0 ∼ 1/ΛQCD is sketched (dotted line).
L l n (x) are the associated Laguerre polynomials and one identifies common eigenvalues M 2 = 4α 2 (n + ν + 1). The linear confining potential generates a mass gap of the order of α. ν is related to the h.o. angular momentum by l + = ν, l − = ν + 1. In the following α 2 = (0.41) 2 GeV 2 will be selected, a value which interpolates among different choices in the literature (cfr. ref. [20] and references therein) and it gives a good fit to the Form factors [36]. A critical analysis of the influence of the α's value on the results of the present approach will be given in Section IV B 1.

B. The IR-improved Soft-Wall model and its solutions
The Infra-Red improved Soft-Wall AdS/QCD model proposed in ref. [28] (in the following: IR) exhibits a confining potential of the form shown in Fig.1. The numerical values of the parameters are as in Table I. The potential (11) belongs to the class of potentials obeying V (ζ → 0) = α 2 ζ, and V (ζ ≫ µ −1 g ) = α 2 ζ, i.e. la class of potentials matching the linear wall both in the IR and UV regimes [27].
Therefore the potential V IR must reduce to the linear soft-wall potential in the limiting case λ B = 0, and one has: parameters used in Fig.1 and in the following.
The IR potential has been constructed to reproduce, with good accuracy, both the meson and the baryon masses. In particular it gives consistent predictions for the mass spectra of scalar, pseudoscalar, vector and axialvector mesons, and both confinement and chiral symmetry breaking are well characterized [29]. In the case of baryons the parameters λ A and λ B are fixed by fitting the masses of the first low-lying baryons with even parity (including nucleon). The predicted masses for odd-parity baryons and high excited states of even-parity baryons are consistently reproduced [28] by using the same values of the parameters.
Let us introduce the form (11) in the equations (4), (5), one gets A convenient technique to solve Eqs. (13) is an expansion on the basis of φ ± nl± of Eqs. (8,9), in this way one can keep all the already established properties of the solutions (8,9) within a linear combination of them. Consequently where ν = 3 and l + = ν and l − = ν + 1 for the lowest three quark Fock state of the nucleon [38]. The natural parameter to be chosen to minimize M 2 IR looking for the ground state wave function (the nucleon) is the harmonic oscillator constant which has to be diversifyed in two components α → α ± in order to respect the essential (15) Of course the restricted Hilbert space used in solving the minimization will result in an upper bound for M 2 IR . However, as it will be more clear in the next Section, the convergence is rapid and one has to expect only few percent deviations.

C. Numerical results
The minimization procedure is performed in the two components φ ± ν,IR varying the parameter α − and reaching the minimum value for (M − IR ) 2 = 2.61 GeV 2 for α − = 2.65 fm −1 , with the corresponding (M + IR ) 2 = 2.61 GeV 2 for α + = 2.35 fm −1 , and involving 17 oscillator quanta (n max = 16). The harmonic oscillator angular momentum quantum numbers l − = ν + 1 and l + = ν are fixed by the twist operator ν = 3 for the lowest number of active quarks [38] (of course nmax n=0 (a ± n,l± ) 2 = 1). In table II the actual vales of the coefficients a ± n,l± . One can appreciate the rapid convergence. The basis is in fact the maximum numerical basis supported by the Matlab code used for the minimization, however it is evident that remaining within n max = 10 is a quite good approximation. The numerical calculations of the next Sections will make use of the restricted basis n max = 10.

III. GPDS AND SUM RULES AT ξ = 0
In order to introduce the explicit calculations of the GPDs, let us concentrate first on the chiral even (helicity conserving) distribution H q (x, ξ, Q 2 , t) for partons of q-flavor at the scale where one is assuming valid the calculation for the related amplitudes. For example, the amplitude for deeply virtual Compton scattering where a virtual photon of momentum q µ is exchanged by a lepton to a nucleon of momentum P µ and a real photon of momentum q ′µ is produced (together with a recoiling nucleon P ′µ ).The space-like virtuality is therefore Q 2 = −q µ q µ and it identifies the scale of the process. The invariant momentum square is t = −∆ 2 = (P ′µ − P µ ) 2 and the skewedness ξ encodes the change of the longitudinal nucleon momentum (2ξ = ∆ + /P + , with 2P µ = (P µ + P ′µ )). In the following the common notation of simply three variables (x, ξ, t) instead of (x, ξ, Q 2 , t) is assumed 1 . In addition only the limit ξ = 0 will be discussed and therefore one can remain in the 0 ≤ x ≤ 1 region.

A. Contribution of the valence quarks (ν = 3)
The helicity conserving H q distributions, in the limit t = 0 and ξ = 0 reduce to ordinary parton distributions the unpolarized quark distribution of flavor q and one has where N q fixes the number of valence quarks of flavor q (N u = 2, N d = 1). The integral properties are therefore model independent and strongly constrain the helicity conserving distributions in any model and/or parametrization (the conditions on N q are satisfied within all the models presented). The second moment is related to the momentum sum rule (cfr. Fig.2) and the models discussed differ: the LF-model is based on Light-Front wave functions and obeys M u+d = 1 since the valence contribution is the only component at low momentum scale. The numerical calculations give: M u+d = 0.92 for the SW, while M u+d = 0.91 for the IR when higher Fock states are considered (cfr. Section III B). In addition the first t-dependent moments of the GPDs are related to the nucleon elastic form factors [4], i.e.
where F q 1 (∆ 2 ) and F q 2 (∆ 2 ) are the contribution of quark q to the Dirac and Pauli form factors. The property (19) does not depend on ξ and it holds also in the present approach with ξ = 0 and therefore 0 ≤ x ≤ 1, (cfr refs. [4,7]), and one has: where t = −∆ 2 and isospin symmetry has been assumed.
In terms of the holographic wave functions φ ± derived from AdS/QCD, the Dirac form factors for the nucleons in the present Soft-Wall linear model are given by [20,34,35,38] where κ p/n are the proton and neutron anomalous gyromagnetic factors respectively. The kernels V ± have a simple and analytic integral form [39]: The specific boundary condition V ± (∆ 2 = 0, ζ) = 1 imposes the normalizations: and the results of the SW [35] model are recovered in the limit α ± → α = 0.41 GeV, a ± nl± → a nl± = 1, and therefore : φ ± ν,IR (ζ) → φ ± nl± (ζ) of Eqs. (8,9) The resulting expressions for the GPDs are B. Higher Fock states (ν = 4, ν = 5) The formalism developed in the previous Sections for the solution of the improved IR potential at the lowest twist (ν = 3), can easily accommodate also higher Fock states in the wave functions opening the possibility of studying their effects on the Generalized Parton distributions even in presence of a modified potential. In particular additional gluons (ν = 4) or a quark-antiquark pair (ν = 5) as discussed in ref. [38]. One obtains: with c 3 = 1.25, c 4 = 0.16, and The previous conditions and values are taken from ref. [38] where they are established for the linear SW potential. However the criteria are rather general and directly related to experimental observables, their application also for the IR potential seems quite natural and it represents, in any case, a first sensible approximation. The minimization has to be repeated for ν = 4, 5 in analogy with the numerical analysis of Section II C. All the expressions derived in the previous Section III are generalized in a straightforward way replacing (φ ± ν,IR (ζ)) 2 and φ ± ν,IR (ζ) φ ∓ ν,IR (ζ) with the linear combinations (31). In order to comment in more detail, the generalization of the Eq.(27) is given as an example: One has to notice that the normalization factors N 2 will depend on ν while the harmonic oscillator parameters α ± will not. In fact the baryon masses (fixed by the explicit form of the confining potential) will get their minimum values for the same α ± , as it has been checked numerically. The generalization is straightforward. In the appendix the numerical details.

IV. GPDS AND CONFINING POTENTIALS: RESULTS AND COMMENTS
Results for H u (x, ξ = 0, t) are shown in Figures 3 and  4. In particular in Fig.3 the results for the valence components H u V (x, ξ = 0, t), i.e. the twist-3 contributions (ν = 3) are shown for both the SW model and the IR improved model. One could imagine that the change in the confining potential encodes just refinements producing only small effects on the observables. This comment is true from the point of view of the baryon spectra, however the modifications induced on the wave functions can show up in a more consistent way in appropriate observables. It is the philosophy of the present work and it is well illustrated in Fig.3: comparing the SW and the IR improved results one can appreciate the effects produced by the tuning of the confining potential (cfr. Fig.1). Analogous effects emerge in the analysis of the response of d-valence quarks (Fig.5). In that case the effects of the IR improved potential seem to be even more evident in the low-x region and for both t = −0.2 GeV 2 and t = −0.5 GeV 2 .
The t-dependence of the H-GPDs can be appreciated comparing the upper and lower panels of Figs.3 and 5 where the responses are shown for two different values of the momentum t = −0.2 GeV 2 and t = −0.5 GeV 2 . In particular in Fig.3 the results of the present AdS/QCD approach are compared with an investigation (cfr. ref. [40]) which makes use of a Light-Front relativistic quark model developed in ref. [41] and based on a q-q-potential with a linear plus a Coulomb-like component: V = − τ r + κ l r . The predictions of the two approaches look rather different. The constrains due to conformal symmetry breaking imposed by the AdS/QCD approach seems to reduce the response considerably (and in the whole x-range) changing, at the same time, their t-dependence in a relevant way. Figures 4 and 5 are devoted to the investigation of the higher Fock states effects. Within the IR-improved potential the ν = 3 and ν = 3, 4, 5 responses are shown and compared. The effects of higher-Fock states is rather weak, but one has to keep in mind the limited validity of the contribution for ξ = 0, the only component here discussed. The role of quark-antiquark and gluon components should show up in a more consistent way in the ξ-dependence of the response [1,42,43]. It would be particularly interesting, in view of the next generation of experiments, to add explicitly such components together with the appropriate perturbative QCD evolution. Work in this direction is in progress.
The comparison with experiments seems also particularly interesting from the point of view of the tdependence of the responses. Often such a dependence is taken following the fall off of the nucleon form factors. Modeling GPDs does not confirm that hypothesis and the results of ref. [40] already questioned such a t-dependence. The results of the AdS/QCD approach show an even stronger t-dependence, a peculiarity which should be explicitly investigated in future experiments.
The integral properties of the helicity non-conserving responses E q are more model dependent: where κ q is the anomalous magnetic moment. Experimentally κ u = 2κ p + κ n = 1.67 and κ d = 2κ n + κ p = Also for the E q distributions the SW and the IR-2 More explicitly, ref. [40] investigates two L-F quark models: i) a Hypercentral potential which includes linear and Coulombian interactions and which is SU (6) symmetric; ii) a model with Goldstone Boson Exchange (GBE) [45] which breaks SU (6). Despite the fact that κ p/n are in principle sensitive to SU (6) breaking effects, the vales of the two models do not differ that much. For details cfr. ref. [40]. improved potentials predict significantly different results as far as their x-dependence is concerning. The comparison with the LF-approach shows also an important difference in t-dependence between the LF and the AdS/QCD approaches. The inclusion of higher Fock states is illustrated in Figs.7 and 8.

The α parameter: a critical analysis
Before discussing some application of the GPDs in AdS/QCD, a critical analysis of the parameter α characterizing the SW potential (cfr. Section II A) could help in fixing the precision one can expect in the present, and analogous, investigations. To this end it is convenient to write explicitly the GPDs within the SW approach, as they result from Eqs. Eq.(26) 3 : (with t = −∆ 2 ). It is evident from Eq.(34), (and analogous expressions can be written for the helicity dependent components) that, for t = 0, the parameter α does not affect the x-dependence of H q , it influences its t-dependence. Such a conclusion has the relevant consequence that the differences one can see in Fig.2 are α-independent. For t < 0 the effects are more complicated correlating in a critical way the x and tdependence 4 and the choice of the α-parameter appears to be critical. The discussion of the SW spectrum for baryons shows that the masses obey the Regge behavior and α ≈ 0.5 GeV is needed to reproduce the nucleon spectrum (cfr. Section II A). In the literature the values α = 0.49 GeV and α = 0.51 are considered the best choices to reproduce, within the holographic AdS/QCD, the nucleon and the ∆ spectra respectively [20]. The freedom in the choice of α is related to the nature of the AdS/QCD approach and the actual value is fixed following physical constraints like the nucleon and the ∆ masses. In the study of the nucleon electromagnetic form factors, α is fixed in order to reproduce their momentum transfer behavior and, to this end, it has been chosen [36] α = 0.4066 ≈ 0.41 GeV. It is physically sensible to remain within this choice in order to study GPDs. However just to give a flavor of the α dependence of the present investigation in Fig.9 the sensitivity of the helicity independent GPDs to the α's values is shown for both the SW and IR improved potential. The values chosen are: i) the choice made in the previous Sections and related to the electromagnetic form factors, α = 0.41 GeV; ii) the value from the best fit of the nucleon masses, α = 0.49 GeV. The variations shown in Fig.9 could represent an upper bound to the absolute theoretical error. However, one cannot consider the range of the results shown in the figure as genuine theoretical error bars; in fact the value α = 0.41 GeV is well constrained to be associated to the electromagnetic interactions as described within AdS/QCD. One has to keep in mind, indeed, that α is the parameter that appears in the dilaton definition used to break conformal invariance in AdS and it affects all fields considered in the model, including the vector massless field which allows the calculation of form factors (and GPDs), cfr. Eq. (24).
The same parameter appears, in the case of the nucleon, in the Soft-Wall potential, in the holographic coordinate V (z) = α 2 z. Consistency is mandatory and the α value has to be fixed by physical constrains connected with the vector massless field dual to the electromagnetic field and the form factors appear a natural choice.
In concluding the present Section IV, a general comment can be added in order to justify the large differences one can see in the IR versus SW potential predictions as well as in the comparison with the LF model. The quite different behavior of the potentials at intermediate vales of z ≈ 0.5 fm (see Fig.1), introduce relevant differences in the high-momentum components of the wave corresponding functions and, consequently, on the H-distributions. In particular, if the behavior of the IR-potential is extrapolated to small distances (< 0.5 fm) to match Coulomb tail like in the case the LF model, the enhancement at small-and intermediate-x values is emphasized as it emerges, for instance, from Figs. (3) and (6). The responses of the LF model (in the region 0 ≤ x ≤ 0.3) are larger than the IR potential ones; the IR responses are, in turn, larger than the SW model distributions: a coherent behavior. On the contrary, because of the sum rule constrains, the responses in the large-x region follow an inverse behavior.

V. MODELING THE ξ-DEPENDENCE WITH DOUBLE DISTRIBUTIONS
In the present Section the results obtained at ξ = 0, are generalized to the whole ξ domain by means of a double distribution approach developed by Radyushkin in ref. [37]. The approach involves a given profile function and the forward parton distribution as evaluated in the previous Sections (or in a generic model). In order to be specific let us concentrate on the chiral even (helicity con-serving) distributions H q (x, ξ, Q 2 , t). One can introduce NonSinglet (valence) and Singlet quark distributions: The analogous distribution for gluons is symmetric in x, with Once again, the Q 2 dependence has been omitted following the common simplified notation, it will be discussed in Section V A when the hadronic scale Q 2 0 will be introduced. Due to the polynomiality property [47] the symmetry characters, (35), (36) and (37), hold also under ξ → −ξ. The Singlet and gluon components mix under evolution, while the NonSinglet distribution evolve independently.
The t-independent part can be parametrized by a two component form [37] with and H q (x, ξ) ≡ H q (x, ξ, t = 0).
The D q contribution in Eq. (39) is defined in the region |x| ≤ ξ and therefore does not contribute in the forward limit. The D-term contributes to the Singlet-quark and gluon distributions and does not contribute to NonSinglet components. Its effect under evolution is restricted at the level of few percent [44] and it will be disregarded in the following.
Following Radyushkin the DD terms entering Eq.(40) are written as where H q (β, 0, 0) = q(β) (cfr. Eq. (16)) and the profile function is parametrized as [48] h(β, α) = Γ(2b + 2) The parameter b fixes the width of the profile function h(β, α) and the strength of the ξ-dependence. In principle it could be used (within the double distribution approach) as a fit parameter in the extraction of GPDs from hard electro-production observables. The favored choice is b NS = b S = 1 (producing a maximum skewedness) and b gluon = 2 [42,48]. In the limiting case b → ∞, h(β, α) → δ(α)h(β) and H q (x, ξ) → H q (x, ξ = 0). The explicit evaluation of H q (x, ξ) in Eq.(39) makes use of the results of the previous Sections within the holographic AdS/QCD approach.
A. Results at low momentum scale: the Soft-Wall model In the present Section the results for the chiral even distributions of the Soft-Wall model valid at ξ = 0 as discussed in Section IV A, are generalized to ξ > 0 by means of the Double Distributions presented in the previous Section. They are defined in the different regions of the (generalized) x-values by the integrals (40) and the combinations (35), (36) and (37) (cfr. also ref. [48]): with The results at low-momentum scale, Q 2 0 , where the Soft-Wall model is supposed to be valid, are shown in Fig.10 for three different values of the skewedness parameter ξ = 0.1, 0.2, and 0.3 and invariant momentum t = −∆ 2 = 0. The value of the low-momentum scale is identified by means of the momentum sum rule. In fact the number of particles are well defined at the initial scale (cfr. Eq.(17)), and the momentum sum rule is not fulfilled by valence quarks only. As matter of fact one has: Differently from a quark model (relativistic or nonrelativistic) based on the presence of only valence quarks at the lowest scale, the holographic approach is intrinsically based on the QCD dynamics. The bound system of valence quarks cannot share momentum among a pure three-quark system. The masses of the quarks are unknown and what is reproduced is the spectrum of the system. The interpretation of the nucleon bound system implies the presence of gluons exchanged among the valence quarks. A natural consequence is an additional gluon distribution filling the gap to the total momentum. A gluon distribution proportional to the valence densities at Q 2 0 (à la Glück, Reya, Vogt [50]) can be a sensible choice such that with A g = 0.091 and dx xg(x, Q 2 0 ) = 0.08. The (small) H g gluon distribution of Eq.(37) as consequence of the density (45), is shown in the lowest panel of Fig.10. The factor 20 is needed to make H g comparable with the results shown in the other panel of the same figure, H S and H N S . The choice (45) is only one of the possible choices one can make. One could assume a different parametrization of the gluon distribution (45), and deduce a different behavior of the (small) H g component of Fig.10. In the previous studies of GPDs within AdS/QCD no mention is made of the fact that the momentum sum rule is not satisfied, i.e. property (44). The main reason to introduce here a conserving momentum sum rule like Eq. (46), is related to the possibility of a detailed investigation of the perturbative QCD evolution properties of the distributions. The simplified assumption made in Eq. (45) is mostly connected to the fact that the perturbative evolution is dominated by the value of the moment carried by the gluon component rather than by the exact form of the distribution.
In order to appreciate the role of the invariant momentum transfer t = −∆ 2 , the results of Fig.10 valid at t = 0 are summarized in Fig.11 and compared with the analogous predictions for t = −0.5 GeV 2 . The SW model gives a non vanishing contribution to quark GPDs in the region |x| < ξ at the lowest scale Q 2 0 without introducing discontinuities at |x| = ξ, the ξ-dependence is rather weak (cfr. Fig.10). One can check, in particular, that H S = H N S at x > ξ, a peculiarity due to the absence of sea contribution at Q 2 0 .

VI. CONCLUSIONS AND PERSPECTIVES
A study of GPDs within a general AdS/QCD framework, has been presented. Two main features have been emphasized and investigated in detail: i) the role of the confining potential in the holographic coordinate as described within the Soft-Wall and within more general potential models; in particular the possibility of introducing high Fock states in the calculation of GPDs.
A method to study effects due to different confining potentials introduced to break conformal symmetry in the AdS/QCD approach to baryons has been proposed in Section II. In several works devoted to the investigation of AdS/QCD wave functions for baryons, often the complementary aspect is stressed: the potentials must manifest iso-spectral properties and therefore their differences have to be adequately mitigated [27]. On the contrary the use of different (almost iso-spectral potentials) in calculating amplitude and responses in Deep Inelastic Scattering, can put in evidence relevant differences that can discriminate among them. The specific observables discussed in relation with Generalized Parton Distributions are a good example. The method implies the use of the Soft Wall solutions as a complete basis to solve more sophisticated potential models. The results are promising: the power of the holographic approach seems to be preserved and observables can be calculated following well established techniques. Higher Fock states can be accommodated showing their relevance in the whole x-region.
ii) the extension of the GPDs results from AdS/QCD methods from the forward (ξ = 0) to the offforward region (ξ > 0). The procedure used (Double -Distributions) enlarge the phenomenological domain of the GPDs predictions opening the concrete evaluation of the Single, Non-Singlet GPDs in the whole (x,ξ,t) domain. The procedure used identifies also the resolution scale of the results. The example developed is restricted to the Soft-Wall, but it is easily generalized to more complex confining potentials.
n a + (ν=4)n a − (ν=4)n a + (ν=5)n a − the new generation of experimental data. In particular the important contributions due additional degrees of freedom like non-perturbative gluon and sea components should be further investigated together with a detailed analysis of the perturbative effects due to QCD evolution. The elegance and the effectiveness of the AdS/QCD approaches has to be integrated in a complete predictive scheme for a large variety of observables in the perspectives of modeling the nucleon structure [46], work in that direction is in progress.

Appendix A: Higher Fock states
In this appendix some details of the method proposed in Section II are illustrated.
The procedures can be generalized in order to accommodate higher Fock states. As discussed in Section III B the values of α ± remain the same also for the wave functions with ν = 4 and ν = 5 (as a numerical check has confirmed). The minimization produces the values of the coefficients shown in table III.