Diffractive incoherent vector meson production off protons: a quark model approach to gluon fluctuation effects

The European Physical Journal C, Apr 2019

Fluctuations play an important role in diffractive production of vector mesons. It was in particular recently suggested, based on the Impact-Parameter dependent Saturation model (IPSat), that geometrical fluctuations triggered by the motion of the constituent quarks within the protons could explain incoherent diffractive processes observed at HERA. We propose a variant of the IPSat model which includes spatial and symmetry correlations between constituent quarks, thereby reducing the number of parameters needed to describe diffractive vector meson production to a single one, the size of the gluon cloud around each valence quark. The application to \(J/\varPsi \), \(\rho \) and \(\phi \) diffractive electron and photon production cross sections reveal the important role of geometrical fluctuations in incoherent channels, while other sources of fluctuations are needed to fully account for electroproduction of light mesons, as well as photo production of \(J/\varPsi \) mesons at small momentum transfer.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Diffractive incoherent vector meson production off protons: a quark model approach to gluon fluctuation effects

The European Physical Journal C April 2019, 79:327 | Cite as Diffractive incoherent vector meson production off protons: a quark model approach to gluon fluctuation effects AuthorsAuthors and affiliations Marco Claudio TrainiJean-Paul Blaizot Open Access Regular Article - Theoretical Physics First Online: 11 April 2019 2 Shares Abstract Fluctuations play an important role in diffractive production of vector mesons. It was in particular recently suggested, based on the Impact-Parameter dependent Saturation model (IPSat), that geometrical fluctuations triggered by the motion of the constituent quarks within the protons could explain incoherent diffractive processes observed at HERA. We propose a variant of the IPSat model which includes spatial and symmetry correlations between constituent quarks, thereby reducing the number of parameters needed to describe diffractive vector meson production to a single one, the size of the gluon cloud around each valence quark. The application to \(J/\varPsi \), \(\rho \) and \(\phi \) diffractive electron and photon production cross sections reveal the important role of geometrical fluctuations in incoherent channels, while other sources of fluctuations are needed to fully account for electroproduction of light mesons, as well as photo production of \(J/\varPsi \) mesons at small momentum transfer. 1 Introduction Fluctuations play an essential role in the diffractive production of vector mesons. It was recently suggested that these fluctuations could be dominated by those, event by event, of the constituent quark positions inside the proton, and that these could be constrained by the incoherent diffractive photoproduction of \(J/\varPsi \) mesons off protons [1]. Such fluctuations, of essentially geometrical origin, are commonly referred to as “geometrical fluctuations”. They are the analog of the fluctuations linked to the positions of the nucleons in high energy nucleus-nucleus collisions [2]. As we shall see, a crucial ingredient entering the calculation of the diffractive processes is the cross section of a small color dipole crossing the proton at a given impact parameter. The interaction of the dipole with the proton is directly sensitive to the total density of gluons that it “sees” on its path through the proton. Although we have experimental information about the total (integrated over the impact parameter) density of gluons in a proton, the dependence on the impact parameter is much less constrained. A simple dipole model that includes the physics of saturation and takes explicitly into account the impact parameter dependence of gluon distributions is the Impact-Parameter dependent Saturation model (IPSat) [3, 4, 5, 6, 7]. In the IPSat model the impact parameter dependence of the amplitude is simple to implement and it can be easily generalized from Deep Inelastic Scattering (DIS) off protons to DIS off nuclei [8, 9, 10, 11]. Other excellent probes of the high energy saturation regime are the exclusive diffractive processes in the electron-proton collisions: exclusive vector meson production and deeply virtual Compton Scattering (DVCS) are the prominent examples. Our main interest, in the present work, is the physics of exclusive diffractive meson production, since an interesting new piece of information can be extracted from such reactions, namely how much the spatial gluon distribution fluctuates, event-by-event, within a proton. Experimentally one can access this information via exclusive incoherent diffractive meson production, i.e. events connected with a dissociated proton [12]. Including the analysis of coherent diffractive processes where the proton remains intact, both the impact parameter dependence and the fluctuations of the gluon distribution in the proton can be constrained [1]. Different final states depend in different ways on the impact parameter, where intrinsically non-perturbative physics may become relevant. Thus the fluctuations of the shape of the gluon distribution may be influenced by non-perturbative physics and the aim of the present work is a detailed study of some of such non-perturbative effects. To this end we present a self-consistent approach where the spatial quark and gluons distributions are consistently calculated. The number of parameters drastically reduces and the predictive power of the IPSat model increases since it is based on calculated properties of the quark wave functions. The paper is organized as follows. In Sect. 2 we review the approach used by Mäntysaari and Schenke in Ref. [1] to calculate the vector meson production cross sections. In particular, we stress the role of geometrical fluctuations in the description of incoherent photoproduction of \(J/\varPsi \) mesons. In Sect. 3.1 we present and discuss the quark correlations which are relevant in the description of diffractive processes. These are obtained in a specific quark model that allows for a simple determination of the quark wave functions of the nucleons. Fluctuations in the density of gluons are introduced, as in Ref. [1], by attaching a gluon cloud around each valence quark. In Sect. 4 we use DGLAP evolution equations to various degrees of precision in order to relate quarks, gluons and sea quark degrees of freedom at the initial non-perturbative scale to their values at the large experimental scale. In Sect. 5 we present our main results for the coherent and incoherent \(J/\varPsi \) photoproduction, while \(J/\varPsi \), \(\rho \) and \(\phi \) electron photoproduction is discussed in Sect. 6. Finally, conclusions are drawn in Sect. 7. 2 Diffractive deep inelastic scattering in the dipole picture In deep inelastic lepton–proton scattering, the exclusive production of vector mesons (V) proceeds via the exchange of pomerons in the case of a diffractive process where no color is exchanged between the proton and the produced system. The absence of colored strings leads to a rapidity gap (a region in rapidity with no produced particles) which characterizes experimentally the diffractive events. If the scattered proton remains intact, the process is called coherent, while for incoherent processes the final proton breaks up (see Ref. [13] for an introduction to diffractive processes and their description within perturbative QCD). Explicitly, using the notation of Ref. [1], we write the coherent diffractive cross section as $$\begin{aligned} {d \sigma _{T,L}^{\gamma ^* p \rightarrow V p} \over dt} = {(1+\beta ^2) \over 16 \pi } \left| \langle {{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta )\rangle \right| ^2\,, \nonumber \\ \end{aligned}$$ (1) where \({{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta )\) is the scattering amplitude, \({x_{\mathbb {P}}}=(P-P')\cdot q/(P\cdot q)\) the fraction of the longitudinal momentum of the proton transferred to the pomeron (\({\mathbb {P}}\)), and the momentum transfer (square) is \(t = -(P'-P)^2\) with P and \(P'\) the initial and final proton four-momenta. The virtual photon-proton scattering is characterized by a total center-of-mass-energy squared \(W^2 = (P+q)^2\), (\(Q^2 = -q^2\)). Finally \({\varDelta }=(P'-P)_\perp \) is the transverse momentum transfer.1 The amplitude \({{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta )\) for diffractive vector meson production assumes the form [4, 5] $$\begin{aligned}&{{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta ) = i \int d^2 \mathbf{r} \int d^2 \mathbf{b} \int {dz \over 4 \pi } \nonumber \\&\quad \times \left( \varPsi ^* \varPsi _V \right) _{T,L}(Q^2,\mathbf{r},z) \, e^{-i[\mathbf{b}-(1-z) \mathbf{r}] \cdot {\varDelta }} \nonumber \\&\quad \times {{d \sigma _{q {\bar{q}}} \over d^2 \mathbf{b}}}(\mathbf{b},\mathbf{r}, x_{x_{\mathbb {P}}}), \end{aligned}$$ (2) where the subscripts T and L refer to transverse and longitudinal polarization of the exchanged virtual photon. This expression is based on the dipole picture: the photon fluctuates into a quark-antiquark pair, a color dipole, with transverse size \(\mathbf{r}\), while z is the fraction of the photon’s light-cone momentum carried by a quark. This picture holds in a frame where the dipole lifetime is much longer than the interaction time with the target proton. The \(\gamma ^* p\) scattering then proceeds through three steps: (1) the incoming virtual photon fluctuates into a quark–antiquark pair; the splitting of the photon is described by the virtual photon wave function \(\varPsi \), which can be calculated in perturbative QED (see e.g. Ref. [14]). (2) The q–\({\bar{q}}\) pair scatters on the proton, with a cross section \(\sigma _{q {\bar{q}}}\) to be discussed below. This cross section is Fourier transformed into momentum space, with the transverse momentum transfer \(\varDelta \) conjugate to \(\mathbf{b}-(1-z) \mathbf{r}\) (distance, in the transverse plane, from the center of the proton to the center-of-mass of the dipole [5]). (3) The scattered dipole recombines to form a final state, in the present case the vector meson with wave function \(\varPsi _V\) (cf. Appendix A). The factor \((1+\beta ^2)\) in Eq. (1), is described in Appendix B together with other phenomenological corrections. In Eq. (1) the amplitude is averaged over the proton ground state, as indicated by the angular brackets. When breakup processes are included, the square of the average amplitude leaves the place to a sum over intermediate states. Ignoring in that sum the contribution of the ground state, which yields the coherent part of the cross section, we are left with the incoherent cross section. This takes the form [1] $$\begin{aligned} {d \sigma _{T,L}^{\gamma ^* p \rightarrow V p'} \over dt}= & {} {(1+\beta ^2) \over 16 \pi }\left[ \left\langle \left| {{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta )\right| ^2 \right\rangle \right. \nonumber \\&- \left. \left| \left\langle {{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}(x_{x_{\mathbb {P}}},Q^2, \varDelta ) \right\rangle \right| ^2\right] \end{aligned}$$ (3) and involves the variance of the amplitude. Note that, as written in Eq. (2), the amplitude \({{\mathscr {A}}}\) is averaged over the dipole size and the impact parameter. Open image in new window Fig. 1 Coherent photoproduction (\(Q^2 = 0\)) cross section within the kinematical conditions of the HERA experiments: \(x_{\mathbb {P}}\approx 9.6 \times 10^{-4}\) for \(\langle W\rangle = 100\) GeV. The parameter of Eq. (5) is chosen to be \(B _G= 4\) \(\hbox {GeV}^{-2}\) . The H1 data are from [16, 17] 2.1 Coherent production In this paper, we shall rely on the the IPSat model [3, 4, 5, 6, 7], which has been very successful in describing a wide range of data from HERA. In this model the dipole cross section is given by (see e.g. [15]) $$\begin{aligned}&{d \sigma _{q {\bar{q}}} \over d^2 \mathbf{b}} \nonumber \\&\quad = 2 \left[ 1-\exp \left( - {\pi ^2 \over 2 N_c} \mathbf{r}^2 \alpha _S(\mu ^2)\,{x_{\mathbb {P}}} g({x_{\mathbb {P}}},\mu ^2)\,T(\mathbf{b}) \right) \right] \nonumber \\ \end{aligned}$$ (4) where the proton (transverse) spatial profile function \(T_G(\mathbf{b})\) is assumed to be Gaussian in a first approximation, viz. $$\begin{aligned} T(\mathbf{b}) = T_G(\mathbf{b}) = {1 \over 2 \pi B_G} e^{-{\mathbf{b}^2 / (2 B_G)}}. \end{aligned}$$ (5) The scale \(\mu \) in the gluon distribution function \({x_{\mathbb {P}}} g({x_{\mathbb {P}}},\mu ^2)\) is related to the size \(\mathbf{r}\) of the dipole $$\begin{aligned} \mu ^2 = \mu ^2(\mathbf{r}^2) = \mu _0^2 + {4 \over \mathbf{r}^2}, \end{aligned}$$ (6) and the gluon distribution is parameterized as $$\begin{aligned} x g(x,\mu _0^2) = A_g \, x^{-\lambda _g}\,(1-x)^{5.6}. \end{aligned}$$ (7) Loosely speaking, what the IPSat model does in Eq. (4), is to take the integrated gluon distribution (7), and redistribute the gluons in transverse plane according to the phenomenological profile \(T(\mathbf{b})\) given in Eq. (5). As an illustration of the results obtained within such an approach, we display in Fig. 1 the cross section for the coherent \(J/\varPsi \) diffractive photoproduction (\(Q^2 = 0\), real photons, and therefore transverse response only). The results shown in Fig. 1 reproduce those shown in Fig. 6 of Ref. [1], for \(B=4\) \(\hbox {GeV}^{-2}\). The scale \(\mu _0^2\) entering the initial condition for the DGLAP evolution of the gluon distribution \({x_{\mathbb {P}}} g({x_{\mathbb {P}}},\mu ^2)\) [3], is taken from Ref. [6] (\(m_c = 1.4\) GeV is used for the charm quark mass). 2.2 Incoherent \(J/\varPsi \) diffractive production The incoherent component of the diffractive cross section for vector meson production involves the fluctuation of the amplitude (see Eq. (3)). Following the authors of Ref. [1], we assume that these fluctuations have a geometrical origin, i.e., they are dominated by the fluctuations, event by event, of the locations of the constituent quarks in the transverse plane. We then consider the density \(T(\mathbf{b})\) in Eq. (4) as resulting from the sum of the contributions of the individual quarks, i.e., $$\begin{aligned} T(\mathbf{b}) \rightarrow {1 \over N_q} \sum _{i=1}^{N_q}T_q(\mathbf{b}-\mathbf{b}_i), \end{aligned}$$ (8) with $$\begin{aligned} T_q(\mathbf{b}) = {1 \over 2 \pi B_q}e^{-\mathbf{b}^2/(2 B_q)} \end{aligned}$$ (9) with parameter \(B_q\). That is, we assume that each constituent quark is surrounded by a cloud of gluons, assumed also to be Gaussian, and represented by \(T_q\) in Eq. (9). In practice one starts sampling the constituent quarks’ positions in the transverse plane (\(\mathbf{b}_i\), \(i=1,2,3\)), from a Gaussian distribution with width parameter \(B_{qc}\), neglecting any possible correlations between the quarks [1]. For fixed \(N_q\) (\(N_q = 3\)) the degree of fluctuations is controlled by the relative sizes of the parameters \(B_{qc}\) and \(B_q\). In Fig. 3 an example of a “lumpy” proton configuration is shown: it corresponds to a relatively broad distribution of constituent quarks, \(B_{qc}= 3.3\) \(\hbox {GeV}^{-2} = (0.3585\,\mathrm{fm})^2\), and a small size gluon cloud around each valence quarkx, \(B_q = 0.7\) \(\hbox {GeV}^{-2} = (0.1651\,\mathrm{fm})^2\). In contrast, Fig. 2 shows a “smooth” proton that has little fluctuations: this corresponds to a compact distribution of constituent quarks, \(B_{qc}= 1.0\) \(\hbox {GeV}^{-2} = (0.1973\,\mathrm{fm})^2\), with a broad distribution of gluons around each constituent quark, \(B_q = 3.0\) \(\hbox {GeV}^{-2} = (0.3418\,\mathrm{fm})^2\)). The parameters are chosen in such a way that the two-dimensional gluon root mean square radius of the proton is kept at the fixed value $$\begin{aligned} \sqrt{\langle \mathbf{b}^2\rangle }= & {} \sqrt{2 B} = \sqrt{2(B_{qc}+B_q)} \nonumber \\= & {} 2 \sqrt{2}\; \mathrm{GeV^{-1}} \approx 0.55\; \mathrm{fm}. \end{aligned}$$ (10) Open image in new window Fig. 2 Example of a smooth density profile obtained from Eq. (8) with the parameters indicated in the panel Open image in new window Fig. 3 Examples of a “lumpy” density profile obtained from Eq. (8) with the parameters indicated in the panel Open image in new window Fig. 4 Coherent and incoherent \(J/\varPsi \) photoproduction cross section at \(\langle W\rangle =100\) GeV. The full line show the results for the incoherent scattering and the lumpy configuration parameters of Fig. 3 (\(B_{qc} = 3.3\) \(\hbox {GeV}^{-2}\), \(B_q=0.7\) \(\hbox {GeV}^{-2}\)). The same parameters used for the coherent scattering give the result shown by the dashed line. For comparison also the coherent results obtained without geometric fluctuations (\(B=B_G=4\) \(\hbox {GeV}^{-2}\) (dot-dashed line)) are shown. Coherent H1 data from Refs. [16, 17] (circles), incoherent data (triangles) from H1 and ZEUS experiments of Refs. [17, 18]. The single square refers to the total cross section of the H1 experiment at large momentum transfer [19] The configurations obtained via the sampling procedure just described represent the basic ingredients for a complete calculation of the coherent and incoherent diffractive vector meson production. The number of configurations considered in the present study for the evaluation of Eqs. (1) and (3), is \(N_{\mathrm{conf}}=10{,}000\). We have checked that the results of our simulations are stable when \(N_{\mathrm{conf}}\) is increased beyond this value. We show in Fig. 4 the results obtained for the photoproduction cross sections, in the kinematical conditions of the HERA experiments, and for the “lumpy” configurations of Fig. 3. As can be seen the coherent as well as incoherent data of the HERA experiments are well reproduced (see the captions of Fig. 4 for more details). For comparison also the coherent results obtained without geometric fluctuations and an average Gaussian profile (with \(B = B_G = 4\) \(\hbox {GeV}^{-2}\)) are shown (cf. Eq. (5) and Fig. 1). Open image in new window Fig. 5 Coherent and incoherent \(J/\varPsi \) photoproduction cross section at \(\langle W\rangle =100\) GeV. The full line shows the result for the incoherent scattering using the “smooth” configuration parameters (see legend and Fig. 2). Data as in Fig. 4 Finally, we consider the respective influence of the fluctuations on coherent and incoherent cross sections. If one smoothens the strength of the fluctuations by choosing as Gaussian parameters the values of Fig. 2 (i.e. \(B_{qc} = 1.0\) \(\hbox {GeV}^{-2}\) and \(B_{q} = 3.0\) \(\hbox {GeV}^{-2}\)) and calculates again coherent and incoherent cross sections for diffractive photon-production at HERA kinematical conditions, the results of Fig. 5 are obtained. The incoherent cross section is largely underestimated, while the calculated coherent cross section reproduces the HERA data. This just confirms the conclusion of Ref. [1] regarding the sensitivity of the incoherent scattering to the strength of the (geometrical) gluon fluctuations. 3 A quark model based approach to diffractive scattering The description of incoherent diffractive vector meson production that has been discussed in the previous section relies on simple Gaussian approximations for the quark distribution as well as the gluon distribution around each constituent quark. They have revealed the large sensitivity of the process to the fluctuations in these distributions. However the calculation, which essentially duplicates that of Ref. [1], completely neglects correlations between the constituent quarks. Such correlations could however affect the gluon fluctuations. Our goal in the next sections is to develop a simple treatment of these correlations, based on a quark model for the nucleon wave function (QMBA). As an outcome of this approach, we shall see that the number of free parameters to describe the diffractive scattering is drastically reduced and the predictions are more directly related to the quark and parton dynamics. The correlations among (constituent) quarks are induced by their mutual interaction, in particular by the One-Gluon-Exchange. In the non-relativistic limit, this yields the so-called hyperfine interaction [23] $$\begin{aligned} V_{hyp}= & {} {2 \over 3}{\alpha _S \over m^2} \sum _{i < j} {8 \pi \over 3} \Bigg \{ \mathbf{S}_i \cdot \mathbf{S}_j \, \delta (\mathbf{r}_{ij}) \nonumber \\&+\, {1 \over r_{ij}^2} \left[ 3 \,(\mathbf{S}_i \cdot {\hat{r}}_{ij}) \, (\mathbf{S}_j \cdot {\hat{r}}_{ij}) - \mathbf{S}_i \cdot \mathbf{S}_j\right] \Bigg \}\,. \end{aligned}$$ (11) This interaction introduces a spin dependence in the quark wave function. In particular the contact term of Eq. (11) (which is the most relevant) is repulsive in \(S = 1\) states (uu pairs in protons and dd pairs in neutrons) and attractive in \(S = 0\) (ud pairs). It contains also a tensor component expressed in terms of the quark spin \(\mathbf{S}_i\) and the relative coordinates \(\mathbf{r}_{ij}\). The \(N - \varDelta \) mass difference (fixed at about 300 MeV) also fixes the value of \(\alpha _S\) (see e.g. Ref. [24]). 3.1 The Isgur and Karl model and SU(6) breaking The presence of the hyperfine interaction naturally breaks SU(6) symmetry and leads to a description of the proton as a superposition of different SU(6) configurations (multiplets 56, 70). A specific realization is given by the model introduced by Isgur and Karl, where, by diagonalizing the Hamiltonian in a harmonic oscillator (h.o.) basis up to \(2 \hbar \omega _0\) states, one finds the following nucleon wave function [25, 26] $$\begin{aligned} |N\rangle= & {} a_S |56, N\, ^2S_{1/2}\rangle _S + a'_{S} |56',N\, ^2S'_{1/2}\rangle _M \nonumber \\&+ \, a_M |70, N\, ^2S_{1/2}\rangle _M + a_D |70, N\, ^4D_{1/2}\rangle _M\,. \end{aligned}$$ (12) The first state in Eq. (12) is in the \(0\hbar \omega _0\)-shell, while the remaining ones are all \(2\hbar \omega _0\) states. The explicit values of the parameters obtained by Isgur and Karl are $$\begin{aligned} a_S = 0.931,\; a'_S = - 0.274,\; a_M = - 0.233,\; a_D = -0.067\,. \end{aligned}$$ Neglecting the breaking of the SU(6) symmetry would give \(a'_S = a_M = a_D = 0\) and the spatial distributions of the u and d valence quarks cannot reproduce the charge distribution in the neutron. 3.2 The two-harmonic-oscillator (2 h.o.) model In the present study we simplify the picture with in mind the description of the scattering properties. We will describe the SU(6)-breaking effects induced by the hyperfine interaction within a harmonic-oscillator model and introduce two different force constants between u and d quarks. For the nucleons, the procedure is as follows (see Ref. [27]): the nucleons p and n are constructed from the two types of constituent quarks, u and d, which are considered to be distinct and not to be permuted. The internal quark wave functions are written as p(uud) and n(ddu), in each case taking the first two quarks to be identical. Given spin-dependent forces, the third (unlike) quark will have a different interaction with the first two (like) quarks than these two will have with each other. The justification for using these wave functions has been discussed in detail by Franklin [28] many years ago, and applied by Capstick and Isgur [29] to construct a relativized quark model for baryons. The two-body potential takes the form $$\begin{aligned} V&= {1 \over 2} K \mathbf{r}\,^2_{12} + {1 \over 2} K' \left( \mathbf{r}\,^2_{13} + \mathbf{r}\,^2_{23} \right) \nonumber \\&= {1 \over 2}\left( 2 K + K'\right) \varvec{\rho }\,^2 + {3 \over 2}K' \varvec{\lambda }\,^2\,, \end{aligned}$$ (13) where \(\varvec{\rho }= \left( \mathbf{r}_1 - \mathbf{r}_2 \right) /\sqrt{2}\) and \(\varvec{\lambda }= \left( \mathbf{r}_1 + \mathbf{r}_2 - 2 \mathbf{r}_3 \right) /\sqrt{6}\) are Jacobi coordinates. Two h.o. constants can be defined: $$\begin{aligned} \alpha ^2= & {} m \omega _{0 \rho } = \left( m (2K + K')\right) ^{1/2}\,,\nonumber \\ \beta ^2= & {} m \omega _{0 \lambda } = \left( 3mK'\right) ^{1/2}\,. \end{aligned}$$ (14) The three-quark wave function is then written as $$\begin{aligned} \varPsi _{3q}= & {} \left[ \varPhi _{\mathrm{color}}\right] _A \times {1 \over \sqrt{2}} \, \left( \phi _N \times \chi _{M A} + \phi _N \times \chi _{M S}\right) \nonumber \\&\times {\alpha ^{3/2}\,\beta ^{3/2} \over \pi ^{3/4}}e^{-(\alpha ^2 \varvec{\rho }\,^2 + \beta ^2 \varvec{\lambda }\,^2)/2 }\,, \end{aligned}$$ (15) where \(\phi _p=|uud\rangle \), \(\phi _n = |ddu\rangle , \) and $$\begin{aligned} \chi _{MA} = (\uparrow \downarrow - \downarrow \uparrow ) \uparrow /\sqrt{2},\qquad \chi _{MS} = (\uparrow \downarrow + \downarrow \uparrow ) \uparrow /\sqrt{2}, \end{aligned}$$ are the spin components. A physically sensible way to fix the parameters (14) is to relate them to the charge r.m.s radius of the proton and neutron: $$\begin{aligned} \langle r^2 \rangle _p= & {} \left( 0.862 \pm 0.012\right) ^2\,\mathrm{fm}^2 = {1 \over \alpha ^2}\,, \nonumber \\ \langle r^2 \rangle _n= & {} - \left( 11.94 \pm 0.18\right) \cdot 10^{-2}\,\mathrm{fm}^2 \nonumber \\= & {} -\, {1 \over 2}\left( {1 \over \alpha ^2} - {1 \over \beta ^2}\right) \,, \end{aligned}$$ (16) and consequently \(\alpha ^2 \approx 1.35\) \(\hbox {fm}^{-2}\) and \(\beta ^2 \approx 1.99\) \(\hbox {fm}^{-2}\) (corresponding to \(K'/K \approx 5.3\)). The neutron charge distribution can be reproduced by breaking the SU(6) symmetry (\(\alpha \ne \beta \)) and vanishes in the SU(6)-symmetric limit of a single harmonic oscillator potential (\(\alpha = \beta \)). 3.3 The density profile function and sampling procedure The 2 h.o. average transverse profile function is obtained from the spherical density \(\rho _{2ho}(r)\), $$\begin{aligned} T_{2ho}(\mathbf{b}) = {1 \over {{\mathscr {N}}}_u+{{\mathscr {N}}}_d}\int _{-\infty }^{+\infty } dz \, \rho _{2ho}\left( r=\sqrt{z^2+\mathbf{b}^2}\right) ,\nonumber \\ \end{aligned}$$ (17) where (from the wave function (15)) $$\begin{aligned} \rho _{2ho}(r) = {{\mathscr {N}}}_u {\kappa _u^3 \over \pi ^{3/2}}\,e^{-\kappa _u^2 r^2}+{{\mathscr {N}}}_d {\kappa _d^3 \over \pi ^{3/2}}\,e^{-\kappa _d^2 r^2}, \end{aligned}$$ (18) and \(\int d \mathbf{r}\, \rho _{2ho}(r) = {{\mathscr {N}}}_u+{{\mathscr {N}}}_d = 3\). Consequently $$\begin{aligned} T_{2ho}(\mathbf{b})= & {} {{{\mathscr {N}}}_u \over {{\mathscr {N}}}_u+{{\mathscr {N}}}_d}{1\over 2 \pi B_u}\,e^{-\mathbf{b}^2/(2 B_u)} + \nonumber \\&+ {{{\mathscr {N}}}_d \over {{\mathscr {N}}}_u+{{\mathscr {N}}}_d}{1 \over 2 \pi B_d}\,e^{-\mathbf{b}^2/(2B_d)}, \end{aligned}$$ (19) with $$\begin{aligned}&\int d \mathbf{b}\, T_{2ho}(\mathbf{b}) = 1, \nonumber \\&{1 \over 2 B_u} = \kappa _u^2 = {3 \over 2} 4{\alpha ^2 \beta ^2 \over 3 \alpha ^2 + \beta ^2} \approx 2.67\;\mathrm{fm^{-2}};\;\;\;\nonumber \\&\quad B_u \approx 4.8\;\mathrm{GeV^{-2}},\nonumber \\&{1 \over 2 B_d} = \kappa _d^2 = {3 \over 2}{\beta ^2} \approx 2.99\;\mathrm{fm^{-2}}; \;\;\; B_d \approx 4.3\;\mathrm{GeV^{-2}}. \nonumber \\ \end{aligned}$$ (20) Equation (19) explicitly summarizes the effects on the profile function of the correlations between quarks that are due to the SU(6)-breaking component of the One-Gluon-Exchange and the spin-isospin symmetries of the proton wave function. The SU(6)-symmetric limit of a single harmonic oscillator wave function is recovered for $$\begin{aligned} \alpha ^2=\beta ^2 \rightarrow {1 \over 2B_u} = {1 \over 2 B_d} = {3 \over 2}\alpha ^2= {1 \over 2 B_0}, \end{aligned}$$ (21) in which case $$\begin{aligned} T_{2ho}(\mathbf{b}) \rightarrow T_{ho}(\mathbf{b}) = {1 \over 2 \pi B_0}\,e^{-\mathbf{b}^2/(2 B_0)}, \end{aligned}$$ (22) i.e. a Gaussian approximation with \(B_0 = 6.34\) \(\hbox {GeV}^{-2}\). The sum in Eq. (19) can be sampled by a random selection of the single term of the sum, followed by sampling the distribution of that term [30]. In this way the correlated positions of the quarks relative to the origin, \(\mathbf{b}_i\) (\(i=1,2,3\)) are sampled from the 2 h.o. distribution (19). We should emphasize here that the sampling of the one-body density takes into account the correlations among the constituent quarks only to the extent that these modify the one-body density. In principle, since the full wave-function is known, it should be possible to calculate more fully the effect of these correlations, but this is beyond the scope of the present paper. Gluon densities are obtained by adding, as was done earlier, around each constituent quark in the transverse plane as described by the profile (19), a Gaussian gluon distribution with parameter \(B_q\). Examples of gluon transverse density profiles obtained in this way are shown in Fig. 6. The results are analogous to those shown earlier in Fig. 3, but in the present case the only free parameter is the width \(B_q\) of the gluon cloud around each valence quark, the positions of the quarks being determined by the simplified 2 h.o. wave function and the electromagnetic sizes of proton and neutron (cf. Eqs. (16)), with no additional free parameter. Open image in new window Fig. 6 Examples of density profiles from Eq. (19) with \(B_q=0.7\) \(\hbox {GeV}^{-2}\) characterizing the Gaussian shape of each gluon distribution around the constituent quark 4 From quarks to partons The description of the quark states, within an appropriate quantum mechanical approach, as detailed in the previous section, allows us to connect quarks and partons in a consistent way avoiding a new set of parameters entering the gluon distribution (cf. Eq. (7)). In the present section we recall how to connect partons and quarks within a framework which makes use of QCD perturbative evolution. 4.1 Valence quarks and partons A simple description which connects the parton distributions to the momentum density of the constituents has been developed in the past (e.g. [31, 32, 33, 34, 35, 36]). Within that approach the valence quark distribution for the bare nucleon is written as $$\begin{aligned} q_V(x)\large |_{\mathrm{bare}}= & {} {1 \over (1-x)^2}\,\int d^3k\, n(k)\,\delta \left( {x \over 1-x} - {k_+ \over M}\right) \nonumber \\= & {} 2 \pi {M \over (1-x)^2}\,\int _{k_m(x)}^\infty dk\,k\,n(k), \end{aligned}$$ (23) where $$\begin{aligned} k_m(x) = {M \over 2} \left| {x \over 1-x} - \left( {m \over M} \right) ^2 \,{1-x \over x} \right| \,; \end{aligned}$$ (24) \(k_+ = k_0 - k_z\) is the light-cone quark momentum fraction, n(k) the quark momentum density distribution predicted by the specific QM wave functions, and M and m are the nucleon and constituent quark masses, respectively. One can check that \(\int dx \,q_V (x, \mu _0^2) = \int d^3 k \,n(k) = {{\mathscr {N}}}_u + {{\mathscr {N}}}_d = 3\), i.e. the particle sum rule is preserved and the valence quark distributions (23) are defined within the correct support \(0< x < 1\). In detail, the 2 h.o. quark momentum distribution within the proton reads $$\begin{aligned} n(k)= & {} {{\mathscr {N}}}_u\,{1 \over \pi ^{3/2}}\,{1 \over \gamma _u^3}\, e^{-{k^2 \over \gamma _u^2}} + {{\mathscr {N}}}_d \,{1 \over \pi ^{3/2}}\,{1 \over \gamma _d^3}\, e^{-{k^2 \over \gamma _d^2}} = \nonumber \\\equiv & {} n_u(k)+n_d(k) \end{aligned}$$ (25) where: \({{\mathscr {N}}}_u = 2\) and \({{\mathscr {N}}}_d=1\) are the numbers of u and d constituent quarks, while \({1 \over \gamma _u^2} = {3 \over 2}\,{4 \over 3\alpha ^2+\beta ^2}\) and \({1 \over \gamma _d^2} = {3 \over 2}\,{1 \over \beta ^2}\) are the combinations of parameters relevant for u and d momentum densities. One has: $$\begin{aligned} u_V(x)\large |_{\mathrm{bare}}= & {} {{\mathscr {N}}}_u \, {1 \over \sqrt{\pi }} {1 \over \gamma _u}\,{M \over (1-x)^2}\,e^{-{k^2_m(x) \over \gamma _u^2}}\,, \end{aligned}$$ (26) $$\begin{aligned} d_V(x)\large |_{\mathrm{bare}}= & {} {{\mathscr {N}}}_d\, {1 \over \sqrt{\pi }} {1 \over \gamma _d}\,{M \over (1-x)^2}\,e^{-{k^2_m(x) \over \gamma _d^2}} \,; \end{aligned}$$ (27) $$\begin{aligned} q_V(x)\large |_{\mathrm{bare}}= & {} \left[ u_V(x)+d_V(x)\right] _{\mathrm{bare}}\,. \end{aligned}$$ (28) Of course the distributions (26) and (27) refer to an extremely low energy scale where the total amount of momentum is carried by the three valence quarks, with no gluon or sea contributions (bare nucleon). In the next Section a concrete way to include the cloud degrees of freedom will be presented. 4.2 From the meson cloud to sea quark and gluon distributions The quark model can be integrated with its virtual meson cloud incorporating \(q {\bar{q}}\) pairs into the valence-quark picture of the parton distributions described in the previous sub-section, dressing the bare nucleon to a physical nucleon (see e.g. Ref. [37] and references therein). The physical nucleon state is built by expanding it [in the infinite momentum frame (IMF) and in the one-meson approximation] in a series involving bare nucleons and two-particle, meson–baryon (MB) virtual states. The description of deep inelastic scattering (Sullivan process) implies that the virtual photon can hit either the bare proton p or one of the constituents of the higher Fock states. In the IMF, where the constituent of the target can be assumed as free during the interaction, the contribution of those higher Fock states to the quark distribution of the physical proton can be written $$\begin{aligned} \delta q_p(x)= & {} \sum _{BM} \left[ \int _x^1 \frac{dy}{y}\,f_{MB/p}(y)\,q_M\left( \frac{x}{y}\right) \right. \nonumber \\&+ \left. \int _x^1 \frac{dy}{y}\,f_{BM/p}(y)\,q_B\left( \frac{x}{y}\right) \right] \,. \end{aligned}$$ (29) The splitting functions \(f_{BM/p}(y)\) and \(f_{MB/p}(y)\) are the probability of the Fock state containing a virtual baryon (B) with longitudinal momentum y and a meson (M) with longitudinal momentum fraction \(1-y\). The quark distributions in a physical proton are then given by $$\begin{aligned} q(x,\mu _0^2) = Z\,q_p^{\mathrm{bare}}(x) + \delta q_p(x)\,, \end{aligned}$$ (30) where \(q_p^{\mathrm{bare}}\) is given by Eqs. (26) and (27) and \(\delta q_p\) is from Eq. (29). $$\begin{aligned} Z = 1 - \sum _{MB}\,\int _0^1 dy\,f_{BM/p}(y)\,, \end{aligned}$$ (31) Open image in new window Fig. 7 The parton distributions at the scale \(\mu _0^2\) of the physical nucleon. The total non-perturbative sea due to Meson–Baryon fluctuations contains both strange and non-strange components (dashed line). The valence distribution (continuous line) is consistently renormalized (cf. Eqs. (30), (31)). For comparison also the valence distribution of the bare nucleon is shown (dot-dashed line) is the renormalization constant and is equal to the probability to find the bare nucleon in the physical nucleon. In Fig. 7 the results are shown comparing the physical and bare parton distributions. The valence distribution is renormalized by the inclusion of the non-perturbative sea, the total sea distribution includes \(\pi ,\rho , \omega , K, K^*\) Meson–Baryon fluctuations, therefore the total sea is $$\begin{aligned} Sea(x,\mu _0^2)= & {} 2\,{\bar{u}}(x,\mu _0^2)+\, 2\,{\bar{d}}(x,\mu _0^2) + s(x,\mu _0^2) \\&+ {\bar{s}}(x,\mu _0^2), \end{aligned}$$ and strange and non-strange components are considered. The final results for the parton distributions at high resolution scale \(\mu ^2 = \mu _0^2 + 4/r^2\) (cf. Eq. (6)) are then obtained by evolving the initial distribution calculated at the scale \(\mu _0^2\), by means of the DGLAP equations. More details can be found in Ref. [37]. 5 \(J/\varPsi \) photoproduction within the quark model based approach (QMBA) Before showing the complete set of results for the \(J/\varPsi \) diffractive photoproduction in the coherent and incoherent channels, it is perhaps useful to summarize the approach that we have presented in Sects. 3 and 4. 1. We have proposed a generalization of the usual color-dipole picture (IPSat). The aim is to connect the diffractive scattering to proton properties like size, wave function symmetries, avoiding, as far as possible, ad hoc parametrization like in Eqs. (5)–(7). We have constructed a proton wave function in which the SU(6) breaking is simply introduced by means of a two harmonic oscillator potential between constituent quarks whose parameters are fixed by means of the experimental radii of neutron and proton. From that model the parton distributions are calculated at low resolution scale \(\mu _0^2\) including a sea component by means of a well established formalism for the light-cone (perturbative) Meson–Baryon fluctuations. The procedure implies many parameters for the coupling constant, but they are taken from the most recent literature without any specific changes for the description of the diffractive scattering. Standard DGLAP evolution is applied to generate gluon distributions at the scale of the process, \(\mu ^2\). No further parameters are needed.   2. The description of the coherent photoproduction of \(J/\varPsi \) does not need further ingredients and its calculation represent an absolute prediction directly related to a low-energy proton model. To describe incoherent diffraction an additional parameter (\(B_q\)) is needed to relate the fluctuations in the gluon density to the motion of the constituents quarks in the transverse plane, in analogy with Eqs. (8) and (9) as discussed in Sect. 5.2. The parameter \(B_q\) which controls the size of the gluon cloud around each valence quark, is the only adjustable parameter of the approach.   Open image in new window Fig. 8 Coherent and Incoherent photoproduction (\(Q^2 = 0\)) cross sections within the kinematical conditions of the HERA experiments (\(x_{\mathbb {P}}\approx 9.6 \times 10^{-4}\) for \(\langle W\rangle = 100\) GeV). The dashed lines represent the QMBA-2ho predictions for the gluon distributions evolved at LO; the continuous lines represent results with gluons at NNLO (NLO and NNLO predictions cannot be distinguished in the Figure, as emphasized in the text). The incoherent scattering calculations within QMBA-2ho are made with \(B_q = 0.7\) \(\hbox {GeV}^{-2}\). Data as in Fig. 4 5.1 DGLAP evolution The predictions of the cross section for coherent and incoherent diffractive \(J/\varPsi \) photoproduction are compared with HERA data in Fig. 8. The leading order (LO) predictions (dashed lines) refer to a calculation where the gluon distribution that enters the dipole cross section is obtained by evolving the initial parton distribution using DGLAP equation at Leading Order. The next or next-to-next to leading order (NNLO) evolution equations are used for the calculation leading to the full lines (the difference between NLO and NNLO results could not be appreciated in the Figures). Strictly speaking, only the LO calculation is fully consistent with the form of the dipole cross section that we use. In the present calculation, the gluon distributions are predicted at high resolution scale starting from a low resolution physical picture of the nucleon. Their final values depend on the order of the evolution, which is reflected in the (weak) dependence of the diffractive cross sections on the order of the evolution, as shown in Fig. 8. The fact that, as seen in Fig. 8, the experimental results appear to be better reproduced by the higher order evolution may reflect the better determination of the gluon distribution, although the slight inconsistency mentioned above prevents us to draw a too firm conclusion at this stage. We may, minimally, regard the difference between the two sets of calculations as reflecting the intrinsic uncertainties of the theoretical model predictions here discussed. In the following we will present results obtained at NNLO, mainly because they appear to be numerically more stable that the LO ones. 5.2 Fluctuations in incoherent scattering Incoherent diffractive photoproduction can only be described by including gluon fluctuation effects, as we have emphasized earlier. The procedure to include gluon fluctuations extends that used with the Gaussian approximation for the profile functions (see Eqs. (8) and (9), and also Fig. 4). In the case of the QMBA profile, the substitution analogous to Eq. (8) reads $$\begin{aligned}&T_{2ho}(\mathbf{b}) \rightarrow {1 \over N_q} \sum _{i=1}^{N_q}T_{q}(\mathbf{b}-\mathbf{b}_i), \end{aligned}$$ (32) $$\begin{aligned}&T_q(\mathbf{b}) = {1 \over 2 \pi B_q}e^{-\mathbf{b}^2/(2 B_q)}. \end{aligned}$$ (33) From a practical point of view one starts by considering a sampling of the constituent quarks’ positions (\(\mathbf{b}_i\), \(i=1,2,3\), in the transverse plane), from the distribution \(T_{2ho}(\mathbf{b})\) of Eqs. (32). This distribution includes part of the correlations between the quark positions coded in the 2 h.o. wave function. The gluon density around each constituent quark is assumed to be Gaussian in the transverse plane, and is described by the function \(T_q\) (33). For fixed \(N_q\) (\(N_q = 3\)) the degree of fluctuations is controlled by the parameter \(B_q\). In Fig. 6 we have already shown an example of lumpy proton configuration assuming \(B_q = 0.7\) \(\hbox {GeV}^{-2} = (0.1651\,\mathrm{fm})^2\) as suggested by the study of a Gaussian profile (Fig. 3 and related discussion). Open image in new window Fig. 9 Coherent and Incoherent photoproduction (\(Q^2 = 0\)) cross sections within the kinematical conditions of the HERA experiments: upper panel: \(x_{\mathbb {P}}\approx 9.6 \times 10^{-4}\) for \(\langle W\rangle = 100\) GeV; lower panel: \(x_{\mathbb {P}}\approx 1.7 \times 10^{-3}\) for \(\langle W\rangle = 75\) GeV. The dot-dashed lines represent the QMBA-2ho predictions for coherent scattering and the related NNLO gluons. Incoherent scattering calculations within QMBA-2ho are represented by the full lines and dashed lines according to the fluctuation parameter \(B_q\) (see legend). Data as in Fig. 4 Figure 9 shows the relevant results for incoherent scattering comparing the calculations with the coherent component. The relevance of the fluctuations is confirmed, and also the value of the parameter \(B_q\). The Fig. 9 shows in fact that the values 0.7 GeV \(^{-2} \le B_q \le 1.0\) \(\hbox {GeV}^{-2}\) remain the favorite range. A consideration which is now independent from other parameters, specifically the parameter \(B_{qc}\) needed to sample the position of the three quarks within the Gaussian approximation used in Ref. [1]. In fact, in our quark model based approach, the quark positions are sampled directly from the proton profile (19) deduced from the quark model wave function. Of course the inclusion of a gluon distribution surrounding each valence quark does modify the global transverse gluon profile as already discussed in the case of a Gaussian transverse density (see Eq. (10) and the related discussion). In the present case the rôle played by the parameter \(B_{qc}\) of Sect. 5.2 is assigned to the two parameters \(B_u\) and \(B_d\) of Eq. (19). In order to keep the transverse gluon root mean square radius fixed one has to replace (cfr. Eqs. (10), (20)). $$\begin{aligned} B_u\rightarrow & {} B_u - B_q; \nonumber \\ B_d\rightarrow & {} B_d - B_q; \end{aligned}$$ (34) when \(B_q > 0\). In that way the gluon rms radius $$\begin{aligned} \sqrt{\langle \mathbf{b}^2\rangle } = \sqrt{{2 \over 3} 2 (B_u +B_q) + {1 \over 3} 2 (B_d+B_q)} \approx 0.60\;\mathrm{fm}\nonumber \\ \end{aligned}$$ (35) will remain fixed varying \(B_q\). 5.3 Quark correlations In order to illustrate the specific rôle of the SU(6)-breaking symmetry and the related quark correlations, one can compare the results of the present QMBA 2h.o. correlated model with the limiting case of a single h.o. wave function which belongs to the 56-th multiplet. The parameters of the two models are chosen in a consistent way, namely fixing the charge radius of the proton at the experimental value, cfr. Sect. 3.2. Obviously the single harmonic oscillator model will predict a vanishing charge radius of the neutron as discussed in the same Section, just because of the lack of SU(6) configuration-mixing in the neutron wave function. The resulting transverse gluon profile function has been discussed in Sect. 3.3. In particular, forcing the single harmonic oscillator model to reproduce the proton charge radius, will result in a rather large value of the transverse gluon root mean square radius. The corresponding coherent scattering cross section is, therefore, too low as it is evident from Fig. 10. The introduction of fluctuations does not alter the conclusion. On the other hand, the incoherent scattering cross section calculated within the same single harmonic oscillator, SU(6)-symmetric, potential is able to follow the data behavior when a lumpy configuration is chosen (\(B_q = 0.7\) GeV \(^{-2}\), full line in Fig. 10). We could conclude that the incoherent scattering is so strongly dominated by the fluctuations that the rôle of quark correlations is unimportant. However such a conclusion needs to be qualified with the following considerations: 1. from Fig. 9: if one uses a wavefunction which includes the proper correlation effects, fluctuations are essential to reproduce the incoherent cross section and, at the same time, the results are moderately sensitive to the free parameter \(B_q\);   2. from Fig. 10: if one uses a wavefunction poorly correlated (e.g. a single harmonic oscillator), the effects of fluctuations are strongly sensitive to the parameter \(B_q\).   The present calculations reveal therefore an interplay between the effects of correlations and those of fluctuations, the latter remaining however the crucial ingredient. Open image in new window Fig. 10 Coherent and Incoherent photoproduction (\(Q^2 = 0\)) cross section within the kinematical conditions of the HERA experiments (\(x_{\mathbb {P}}\approx 9.6 \times 10^{-4}\) for \(\langle W\rangle = 100\) GeV). The dot-dashed line represents the single-ho predictions for coherent scattering. Incoherent scattering calculations within single-ho are represented by the full and dashed lines according to the fluctuation parameter \(B_q\) (see legend). Data as in Fig. 4 5.4 Small |t| behavior The region at very small |t| deserves a specific comment. Indeed this is the region where our predictions for incoherent scattering appear to deviate significantly from the data. When |t| becomes small, the relevant fluctuations acquire a typical wavelength of the order of the size of the proton, and are not described by the geometrical fluctuations that we calculate. This can be seen from a simple analysis of Eqs. (2) and (4). In the limit where \(\varDelta \rightarrow 0\), the integral over the impact parameter in Eq. (2) becomes unconstrained, and the amplitude becomes proportional to the total dipole cross section, that is to the integral of Eq. (4) over impact parameter. The reason why this kills the fluctuations can be easily understood by recalling how fluctuations are generated through the sampling of the valence quarks configurations, namely Eqs. (32), (33): \(T_{2ho}(\mathbf{b})\) fluctuates because its value at a given \(\mathbf{b}\) depends on whether there are valence quarks in the vicinity of \(\mathbf{b}\), which is controlled by the function \(T_q(\mathbf{b}-\mathbf{b}_i)\). When \(|\mathbf{b}|\) is constrained to be small, i.e., \(|\mathbf{b}|\le R\) with R the nucleon size, the final value of \(T(\mathbf{b})\) is sensitive to the location of the individual quarks and its value fluctuates. But when \(|\mathbf{b}|\) is allowed to vary over distance larger than the nucleon size, which is the situation when \(\varDelta \rightarrow 0\), the value of \(T(\mathbf{b})\) becomes insensitive to the precise location of the quarks. Thus the geometrical fluctuations of the kind discussed in the present paper are effective only at not too low momentum transfer. In the region of small momentum transfer, extra sources of fluctuations are called for. This issue has been discussed in Ref. [1]. There, the authors have argued that fluctuations of the gluon density around each valence quark, that they express in terms of the fluctuations of the saturation momentum, can account for the missing ingredient, and can be tuned to reproduce the data in the small |t| region. Note that such fluctuations of the gluon density of the proton could also be understood in terms of the fluctuations of the size of the dipole going through the proton (see e.g. [38] for a recent discussion of such issues). We have already indicated that such fluctuations are explicitly left out in the present calculation. Note also that the fluctuations of the dipole size appear to be the relevant ones at \(t=0\) in the approach based on cross section fluctuations, as discussed for instance recently in Ref. [39]. The discussion of electroproduction in the next section will provide other indications on the importance of these fluctuations. Open image in new window Fig. 11 Upper panel: Coherent differential cross section (\(\sigma _L + \sigma _T\)) for electroproduction of \(J/\varPsi \) as a function of the total momentum transfer (square) t. Predictions within the quark-model-based approach (2 h.o.) and related NNLO gluon profile are compared with data from HERA (H1) [16]. For \(Q^2=22.4\) \(\hbox {GeV}^2\) also the predictions of the Gaussian profile function (\(B_G=4.0\) \(\hbox {GeV}^{-2}\)), and the related (LO) gluon fit are shown as an example. In the lower panel: the incoherent components of the cross section at different \(Q^2\) and within the QMBA (2 h.o.). The gluon fluctuations are included by using the same gluon distribution around each single quark used for \(J/\varPsi \) photoproduction (\(B_q = 0.7\) \(\hbox {GeV}^{-2}\)) 6 \(J/\varPsi \), \(\rho \) and \(\phi \) electroproduction within the QMBA In the present section we complete the presentation of the QMBA approach to the kinematical conditions of electroproduction, i.e. for \(Q^2 > 0\). Diffractive data exist for \(J/\varPsi \), and lighter vector mesons like \(\rho \) and \(\phi \), and whenever possible, we compare our results to the available data. When appropriate, we also compare with the predictions based on the simple Gaussian profile function introduced in Sect. 2.1. The Boosted wave functions used for the vector mesons are described in Appendix A.2. 6.1 \(J/\varPsi \) electroproduction A systematic comparison of our results with the HERA data is presented in the upper panel of Fig. 11 for the \(J/\varPsi \) coherent electroproduction. The data are rather well reproduced within the QMBA description of the transverse gluon shape (cf. Eq. (19)) with no ad hoc parameters. The slopes of the curves are essentially determined by the geometrical size of the nucleon, fixed by the two parameters of the 2ho wave function (cf. Sect. 3), while the gluon distribution entering the dipole cross section (4) keeps the form determined at \(Q^2 = 0\), i.e. for \(J/\varPsi \) photoproduction (cf. Sect. 4.2). The agreement deteriorates somewhat at large \(Q^2\). For \(Q^2 = 22.4\) \(\hbox {GeV}^2\) the Gaussian-model appears to perform slightly better. The lower panel of Fig. 11 provides predictions for the incoherent electroproduction, for which there are no available data. We have considered two kinematical conditions, \(W=100\) GeV and \(W=75\) GeV, and values of \(Q^2\) that are identical to those of the data for the coherent scattering (upper panel). The geometrical fluctuations are calculated following the procedure discussed for the \(J/\varPsi \) photoproduction in the previous section. We recall that at low transfer, these predictions should not be trusted, for the reasons discussed in Sect. 5.4. 6.2 Lighter meson electroproduction6.2.1 \(\rho \) production A large amount of data exist both for coherent and incoherent diffractive electroproduction of \(\rho \)-mesons. An example is shown in Fig. 12 where the H1 and ZEUS data, within a large range of \(Q^2\) values, are compared with our predictions. The present 2ho QMBA approach and the Gaussian approximation (\(B_G = 4.0\) \(\hbox {GeV}^{-2}\)) is confronted to both data sets (lower and upper panels). Slopes and \(Q^2\) dependence are rather well reproduced except for the largest values of \(Q^2\). As was observed already in the case of the \(J/\varPsi \), the Gaussian profile function seems to leads to a better agreement at large \(Q^2\). We note however that as \(Q^2\) increases, the photon wave functions probably become inaccurate (the \(Q^2\) dependence of the whole cross section is entirely due to the \(Q^2\)-dependence of the photon wave function ( Appendix A.1); also, we use a conservative value for the mass \(m_f\) that enters the meson wave function, as it can be seen from table (1) in Appendix A). Finally, we use here the LO evolution in \(\mu (r)\) (cf. Eq. (6)). All these factors could play a role and further studies would be needed to pin down precisely their respective effects. The analysis of incoherent scattering represents a novelty in the study of diffractive vector meson electroproduction2 and in Fig. 13 we show our main results. Once again fluctuations are the crucial ingredient in order to have a non vanishing incoherent cross section. In the present case, however, the slope can be reproduced with a rather larger fluctuation parameter, i.e. \(B_q = 1.5\) \(\hbox {GeV}^{-2}= (0.2417\) fm\()^2\). Keeping instead the value \(B_q=0.7\) \(\hbox {GeV}^{-2}= (0.1651\,\mathrm{fm})^2\) used for the \(J/\varPsi \) photo and electro-production, would lead to too much fluctuations. For the largest value of \(Q^2\) the QMBA results overestimate the data values while the Gaussian approximation predictions (\(B_{q_c} = (4-1.5)\) \(\hbox {GeV}^{-2}\) and \(B_q=1.5\) \(\hbox {GeV}^{-2}\)) are in better agreement. (We recall that \(B_{q_c} + B_q = 4\) \(\hbox {GeV}^{-2}\), in the Gaussian approximation cfr. Sect. 2.2). Open image in new window Fig. 12 Differential cross sections for coherent electroproduction of \(\rho \) mesons as a function of the total momentum transfer (square) t. The two set of data refer to the HERA experiments: H1 (upper panel) [40] and ZEUS (lower panel) [41]. The data are compared with present calculations within the QMBA (2 h.o.) and the Gaussian (\(B_G\)) approximation profiles Open image in new window Fig. 13 Differential cross sections for incoherent electroproduction of \(\rho \) mesons as a function of the total momentum transfer (square) t. The data refer to the HERA H1 experiment [40] and are compared with present calculations within the QMBA (2 h.o.) (full lines) and the Gaussian approximation (\(B_G\)) profiles (dashed lines). In both cases the fluctuation parameter \(B_q = 1.5\) \(\hbox {GeV}^{-2}\) and LO gluon distributions are used 6.2.2 \(\phi \) production Results for coherent and incoherent \(\phi \) elettroproduction are shown in Figs. 14 and 15. In particular in Fig. 14 the coherent diffractive cross section is shown as evaluated within both the QMBA and the Gaussian profile. As in the case of the \(\rho \) meson, both models fail in reproducing the largest \(Q^2\) data which seem to follow a different slope. Open image in new window Fig. 14 Differential cross sections for coherent electroproduction of \(\phi \) mesons as a function of the total momentum transfer (square) t. The data refer to the HERA H1experiment [40] and are compared with present calculations within the QMBA (2 h.o.) and the Gaussian (\(B_G\)) approximation profiles Open image in new window Fig. 15 Differential cross section for incoherent electroproduction of \(\phi \) mesons as a function of the total momentum transfer (square) t. The data refer to the HERA H1 experiment [40] and are compared with present calculations within the QMBA (2 h.o.) (full lines) and the Gaussian (\(B_G\)) approximation profiles (dashed lines). In the upper panel the “lumpy” density profile selected to describe fluctuations in the \(J/\varPsi \) photo and electro production (\(B_q = 0.7\) \(\hbox {GeV}^{-2}\)) is used and the predictions overestimate the data at large |t|; the lower panel results are obtained with a slightly “smoother” density profile (\(B_q = 1.5\) \(\hbox {GeV}^{-2}\)) and are in good agreement with data. See text In order to emphasize the rôle of the fluctuation parameter \(B_q\) we show in the upper panel of Fig. 15, the incoherent diffractive cross section evaluated within the QMBA and the Gaussian approximation fixing \(B_q\) at the value used for the \(J/\varPsi \) production, i.e. \(B_q=0.7\) \(\hbox {GeV}^{-2}\). As in the case of \(\rho \) meson production, for such value of the parameter \(B_q\) one gets too much fluctuations. However, if one chooses the larger value already fixed in the case of diffractive production of the \(\rho \), namely \(B_q=1.5\) \(\hbox {GeV}^{-2}\), one obtains the results shown in the lower panel of Fig. 15 which are in better agreement with the experimental data. Aside from the issues already pointed out in our discussion of the \(\rho \) meson electroproduction, it appears that large dipoles are playing an important role for light mesons. Now, if this is the case, there are features of our calculation that are not well treated. In particular, the fluctuations of the dipole size may induce additional fluctuations that can in fact contribute to smear out the effect of geometrical fluctuations. Such a smearing is achieved here by increasing the size of the parameter \(B_q\). Such fluctuations of the dipole size appear to be unimportant for heavy quarks, and at not too small values of t, and the \(J/\varPsi \) meson production is dominated by contributions of dipoles of small sizes. In this context, it would clearly be very interesting to have data on electroproduction of \(J/\varPsi \) mesons at various \(Q^2\), to test for instance the predictions in Fig. 11 and in view of new electron-ion collider (e.g. Ref. [42]). 7 Conclusions and perspectives In the first part of the present work we have calculated the diffractive photoproduction of \(J/\varPsi \), for both coherent and incoherent channels, including quark correlations in the evaluation of the gluon transverse density profile. The description of the gluon density in the transverse plane has been achieved through a generalization of the IPSat model. This is based on an explicit quark model for the wave function of the valence quarks, with each constituent quark being surrounded by a gluon cloud. Both spatial correlations, induced by a simple two-harmonic-oscillator potential model, and SU(6)-breaking symmetry correlations are included in the wave function, and these appear to play a rôle in the explicit calculation of the cross sections in the two channels. Since the parameters of the quark model are fixed on low energy properties of the proton and the neutron, no adjustable parameters are needed to calculate coherent diffractive production, and a single parameter needs to be selected to describe incoherent diffractive scattering: the width of the gluon distribution around each valence quark. The integrated gluon density is explicitly calculated from the parton distribution deduced (at low resolution scale) from the quark model and evolved to the experimental, high energy, scale using DGLAP equations. A subtle interplay between quark correlations and geometric fluctuations has been observed. The second part of our work has been devoted to enlarge the domain of our study to diffractive vector meson production at \(Q^2 > 0\), i.e. the electroproduction of \(J/\varPsi \) and lighter mesons like \(\rho \) and \(\phi \). Two new ingredients enter the calculations: (1) the \(Q^2\) dependence of the cross sections; (2) the lighter mass of the mesons together with possible new non-perturbative effects. 1. The \(Q^2\) dependence of the cross sections is determined by the photon wave function, more precisely by the overlaps \((\varPsi ^*_V\varPsi )_{T,L}(Q^2,\mathbf{r},z)\) of Appendix A. The Gaussian Boosted wave functions show their limits in the descriptions of large dipoles as discussed in Sect. 6.2. The effect can be sizable at small as well as high \(Q^2\) because of the interplay with the fluctuation of the dipole size. Calculations are in progress to model the lighter meson wave functions and dipole size fluctuations within a more elaborate approach better describing non-perturbative aspects (see also Refs. [44, 45]).   2. The smaller masses of lighter mesons introduce non-perturbative contamination. The net result is that the only parameter describing incoherent diffractive production in our approach (i.e. the size of the gluon cloud around each quark) differs from the heavy \(J/\varPsi \) meson from that needed for lighter mesons like \(\rho \) and \(\phi \). This points to the relevance, for lighter systems, of fluctuations of different origin than the geometrical fluctuations discussed in this paper. This is the case in particular of the fluctuations in the dipole size, that we have argued could play also a role in the very small t region.   Finally, we note that the method introduced here can also be translated in the discussion of Deep Virtual Compton Scattering. Footnotes 1. Throughout this paper we use bold face letters to denote vectors in the transverse plane. 2. An initial analysis of incoherent diffractive electroproduction of \(\rho \) mesons has been proposed in Ref. [1]. Notes Acknowledgements M.T. thanks the members of the Institut de Physique Théorique, Université Paris-Saclay, for their warm hospitality during a visiting period when the present study was initiated. He thanks also the Physics Department of Valencia University for support and friendly hospitality. Useful remarks by Heikki Mäntysaari are gratefully acknowledged. Appendix A: Overlap functions \(\left( {\varPsi ^*\varPsi _V}\right) _{T,L}(Q^2,\mathbf{r},z)\) Appendix A.1: Forward photon wave function The forward photon wave function has been calculated perturbatively (e.g Ref. [20]) $$\begin{aligned}&\varPsi _{h {\bar{h}},\lambda = 0}(r,z,Q) = e_f e \sqrt{N_c} \delta _{h,-{\bar{h}}}\,2\,Q\,z(1-z)\,{K_0(\epsilon r) \over 2 \pi }, \nonumber \\&\mathrm{for \, longitudinal} \mathrm{\,photon\,polarization\, (\lambda =0)};\end{aligned}$$ (A.1) $$\begin{aligned}&\varPsi _{h {\bar{h}},\lambda = \pm }(r,z,Q) = \pm e_f e \sqrt{2 N_c}\Bigg \{i e^{\pm i \theta _r}\left[ z \delta _{h,\pm }\delta _{{\bar{h}},\mp } \right. \nonumber \\&\qquad \left. - (1-z) \delta _{h,\mp }\delta _{{\bar{h}},\pm }\right] \partial _r + m_f \delta _{h,\pm }\delta _{{\bar{h}},\mp }\Bigg \} {K_0(\epsilon r) \over 2 \pi },\nonumber \\&\mathrm{for \, transverse} \mathrm{\, photon\,polarization\, (\lambda =\pm )}, \end{aligned}$$ (A.2) where \(r=|\mathbf{r}|\), \(e=\sqrt{4 \pi \alpha _{\mathrm{em}}} \approx \sqrt{4 \pi /137}\), the subscripts h and \({\bar{h}}\) are the helicities of the quark and the antiquark, respectively, \(\theta _r\) is the azimuthal angle between the vector \(\mathbf{r}\) and the x-axis in the transverse plane. \(K_0\) is the modified Bessel function of second kind, \(\epsilon ^2 \equiv z(1-z)Q^2+m_f^2\) and \(N_c = 3\) is the number of colors. The flavor dependence f enters through the values of the quark charge \(e_f\) and mass \(m_f\), and \(\partial _r K_0(\epsilon r) = - \epsilon K_1(\epsilon r)\). Appendix A.2: Forward vector meson wave function from Ref. [5] (see also Ref. [20]) The simplest approach to modeling the vector meson wave function is to assume that the vector meson is predominantly a quark-antiquark state and that the spin and polarization structure is the same as in the photon case. The transversely polarized vector meson wave function (in complete analogy to the transverse polarized photon) is $$\begin{aligned}&\varPsi ^V_{h {\bar{h}},\lambda =\pm }(r,z) = \pm \sqrt{2 N_c} {1 \over z(1-z)} \Big \{i e^{\pm i \theta _r}\left[ z \delta _{h,\pm }\delta _{{\bar{h}},\mp } \right. \nonumber \\&\qquad \left. - (1-z) \delta _{h,\mp }\delta _{{\bar{h}},\pm }\right] \partial _r + m_f \delta _{h,\pm }\delta _{{\bar{h}},\mp }\Big \} \nonumber \\&\qquad \times \phi _T(r,z). \end{aligned}$$ (A.3) The longitudinally polarized wave function is slightly more complicated since the coupling of the quarks to the meson, contrary to the photon case, is not local. One has: $$\begin{aligned}&\varPsi ^V_{h {\bar{h}},\lambda =0}(r,z) \nonumber \\&\quad = \sqrt{N_c} \delta _{h,-{\bar{h}}} \left[ M_V + \delta \, {m^2_f -\nabla _r^2 \over M_V z(1-z)} \right] \phi _L(r,z), \end{aligned}$$ (A.4) where \(\nabla _r^2 \equiv (1/r) \partial _r + \partial ^2_r\) and \(M_V\) is the meson mass. The nonlocal term was first introduced in Refs. [21, 22]. The overlaps read then: $$\begin{aligned}&(\varPsi ^*_V\varPsi )_T(Q^2,\mathbf{r},z) = {{\hat{e}}_f e \over \pi } {N_c \over z (1-z)}\left\{ m^2_f K_0(\epsilon r) \phi _T(r,z) \right. \nonumber \\&\quad - \left. \left[ z^2 + (1-z)^2 \right] \epsilon K_1(\epsilon r) \partial _r \phi _T\right\} \end{aligned}$$ (A.5) $$\begin{aligned}&(\varPsi ^*_V\varPsi )_L (Q^2,\mathbf{r},z) = {{\hat{e}}_f e \over \pi } 2 N_c\,Q\,z (1-z) K_0(\epsilon r) \nonumber \\&\quad \times \left[ M_V \phi _L(r,z) + \delta \, {m^2_f -\nabla _r^2 \over M_V z(1-z)} \phi _L(r,z)\right] , \end{aligned}$$ (A.6) where the effective charge \({\hat{e}}_f =2/3, 1/3\), or \(1/\sqrt{2}\), for \(J/\varPsi \), \(\phi \), or \(\rho \) mesons, respectively. In addition \(\delta = 1\) is the natural choice done. The boosted Gaussian wave functions in configuration space are written (see Refs. [20, 21, 22]) $$\begin{aligned} \phi _{T,L}(r,z)= & {} {{\mathscr {N}}}_{T,L}\,z(1-z) \exp \left[ - {m^2_f {{\mathscr {R}}}^2 \over 8 z (1-z)} \right. \nonumber \\&- \left. {2 z (1-z) r^2 \over {{\mathscr {R}}}^2} + {m^2_f {{\mathscr {R}}}^2 \over 2} \right] , \end{aligned}$$ (A.7) and \({{\mathscr {N}}}_{T,L}\) and \({{\mathscr {R}}}\) are fixed by normalization conditions and the decay width (see Ref. [5] for other details and Table 1 for the values of the parameters). Table 1 Parameters of the “boosted Gaussian” vector meson wave functions Meson \(M_V\)/GeV \(m_f\)/GeV \({{\mathscr {N}}}_T\) \({{\mathscr {N}}}_L\) \({{\mathscr {R}}}^2\)/\(\hbox {GeV}^{-2}\) \(J/\varPsi \) 3.097 1.4 0.578 0.575 2.3 \(\phi \) 1.019 0.14 0.919 0.825 11.2 \(\rho \) 0.776 0.14 0.911 0.853 12.9 Appendix B: Phenomenological corrections The derivation of the amplitude for the exclusive vector meson production (2) (or DVCS amplitude if \(V \rightarrow \gamma \), real photon) relies on the assumption that the S-matrix is purely real and, consequently, the exclusive amplitude \({{\mathscr {A}}}\) of Eq. (2) purely imaginary. The corrections due to the presence of the real part is accounted for by the factor \((1+\beta ^2)\) multiplying the differential cross sections (1), (3). \(\beta \) is the ratio of real and imaginary parts of the amplitude and is calculated by means of $$\begin{aligned} \beta = \tan {\pi \lambda \over 2}; \;\;\;\;\; \mathrm{with}\;\;\;\; \lambda \equiv {\partial \ln (\mathrm{Im }\,{{\mathscr {A}}}_{T,L}^{\gamma ^* p \rightarrow V p}) \over \partial \ln (1/x_{x_{\mathbb {P}}})}. \end{aligned}$$ (B.8) In addition for vector meson production (or DVCS) one should use the off-diagonal (or generalized) gluon distributions [43]. Such a “skewed” effect is accounted for (in the limit of small \(x_{x_{\mathbb {P}}}\)), by multiplying the gluon distribution \(x g(x,\mu ^2)\) by a factor \(R_g\) given by [5] $$\begin{aligned} R_g(\lambda _g)= & {} {2^{\lambda _g+3} \over \sqrt{\pi }}{\varGamma (\lambda _g+5/2) \over \varGamma (\lambda +4)}, \nonumber \\ \mathrm{with}\;\;\;\; \lambda _g\equiv & {} {\partial \ln [x_{x_{\mathbb {P}}} g(x_{x_{\mathbb {P}}},\mu ^2)]\over \partial \ln (1/x_{x_{\mathbb {P}}})}. \end{aligned}$$ (B.9) The phenomenological corrections (in particular the skewedness correction) are numerically relevant. Evaluated without fluctuations and for \(J/\varPsi \) photoproduction, their (average) numerical values are around \(10\%\) for the real part corrections and \(40\%\) for the skewedness correction in the kinematical region \(|t| \le 0.5\) \(\hbox {GeV}^2\) (see also Ref. [1]). References 1. H. Mäntysaari, B. Schenke, Revealing proton shape fluctuations with incoherent diffraction at high energy. Phys. Rev. D 94, 034042 (2016). arXiv:1607.01711 [hep-ph]ADSCrossRefGoogle Scholar 2. M.L. Miller, K. Reygers, S.J. Sanders, P. Steinberg, Glauber modeling in high energy nuclear collisions. Ann. Rev. Nucl. Part. Sci. 57, 205 (2007). arXiv:nucl-ex/0701025 ADSCrossRefGoogle Scholar 3. J. Bartels, K.J. Golec-Biernat, H. Kowalski, A modification of the saturation model: DGLAP evolution. Phys Rev. D 66, 014001 (2002). arXiv:hep-ph/0203258 ADSCrossRefGoogle Scholar 4. H. Kowalski, D. Teaney, Impact parameter dipole saturation model. Phys. Rev. D 68, 114005 (2003). arXiv:hep-ph/0304189 ADSCrossRefGoogle Scholar 5. H. Kowalski, L. Motyka, G. Watt, Exclusive diffractive processes at HERA within the dipole picture. Phys. Rev. D 74, 074016 (2006). arXiv:hep-ph/0606272 ADSCrossRefGoogle Scholar 6. A.H. Rezaeian, M. Siddikov, M. Van de Klundert, R. Venugopalan, Analysis of combined HERA data in the impact-parameter dependent saturation model. Phys. Rev. D 87, 034002 (2013). arXiv:1212.2974 [hep-ph]ADSCrossRefGoogle Scholar 7. J. Cepila, J.G. Contreras, J.D. Tapia Takaki, Energy dependence of dissociative \(J/\varPsi \) photoproduction as a signature of gluon saturation at LHC. Phys. Lett. B 766, 186 (2017). arXiv:1608.07559 [hep-ph]ADSCrossRefGoogle Scholar 8. H. Kowalski, T. Lappi, R. Venugopalan, Nuclear enhancement of universal dynamics of high parton densities. Phys. Rev. Lett. 100, 022303 (2008). arXiv:0705.3047 [hep-ph]ADSCrossRefGoogle Scholar 9. H. Kowalski, T. Lappi, C. Marquet, R. Venugopalan, Nuclear diffractive structure functions at high energies. Phys. Rev. C 78, 045201 (2008). arXiv:0805.4071 [hep-ph]ADSCrossRefGoogle Scholar 10. H. Mäntysaari, B. Schenke, Probing subnucleon scale fluctuations in ultraperipheral heavy ion collisions. Phys. Lett. B 772, 832 (2017). arXiv:1703.09256 [hep-ph]ADSCrossRefGoogle Scholar 11. J. Cepila, J.G. Contreras, M. Krelina, Coherent and incoherent \(J/\varPsi \) photonuclear production in a energy dependent hot-spot model. Phys. Rev. C 97, 024901 (2018). arXiv:1711.01855 [hep-ph]ADSCrossRefGoogle Scholar 12. H. Mäntysaari, B. Schenke, Evidence of strong proton shape fluctuations from incoherent diffraction. Phys. Rev. Lett. 117, 052301 (2016). arXiv:1603.04349 [hep-ph]ADSCrossRefGoogle Scholar 13. V. Barone, E. Predazzi, High-Energy Particle Diffraction, vol. 565, Texts and Monographs in Physics (Springer, Berlin, 2002)CrossRefGoogle Scholar 14. Y.V. Kovchegov, E. Levin, Quantum Chromodynamics at High Energy (Cambridge University Press, Cambridge, 2012)CrossRefGoogle Scholar 15. A.H. Mueller, Small-\(x\) behavior and parton saturation: a QCD model. Nucl. Phys. B 335, 115 (1990)ADSCrossRefGoogle Scholar 16. H1 Collaboration, A. Aktas et al., Elastic \(J/\varPsi \) production at HERA. Eur. Phys. J. C 46, 585 (2006). arXiv:hep-ex/0510016 17. H1 Collaboration, C. Alexa et al., Elastic and proton-dissociative photoproduction of \(J/\varPsi \) mesons at HERA. Eur. Phys. J. C 73, 2466 (2013). arXiv:1304.5162 [hep-ex] 18. ZEUS Collaboration, S. Chekanov et al., Measurement of proton dissociative diffractive photoproduction of vector mesons at large momentum transfer at HERA. Eur. Phys. J. C 26, 389 (2003). arXiv:hep-ex/0205081 ADSCrossRefGoogle Scholar 19. H1 Collaboration, A. Aktas et al., Diffractive photoproduction of \(J/\varPsi \) mesons with large momentum transfer at HERA. Phys. Lett. B 568, 205 (2003). arXiv:hep-ex/0306013 20. J.R. Forshaw, R. Sandapen, F. Shaw, Colour dipoles and \(\rho \),\(\phi \) electroproduction. Phys. Rev. D 69, 094013 (2004). arXiv:hep-ph/0312172 ADSCrossRefGoogle Scholar 21. J. Nemchik, N.N. Nikolaev, B.G. Zakharov, Scanning the BFKL pomeron in elastic production of vector mesons at HERA. Phys. Lett. B 341, 228 (1994). arXiv:hep-ph/9405355 ADSCrossRefGoogle Scholar 22. J. Nemchik, N.N. Nikolaev, E. Predazzi, B.G. Zakharov, Color dipole phenomenology of diffractive electroproduction of light vector mesons at HERA. Z. Phys. C 75, 71 (1997). arXiv:hep-ph/9605231 CrossRefGoogle Scholar 23. A. De Rújula, H. Georgi, S. Glashow, Hadron masses in a gauge theory. Phys. Rev. D 12, 147 (1975)ADSCrossRefGoogle Scholar 24. M.M. Giannini, E. Santopinto, The hypercentral constituent model and its applications to baryon properties. Chin. J. Phys. 53, 020301 (2015). arXiv:1501.03722 [nucl-th] 25. Isgur, G. Karl, \(P\)-wave baryons in the quark model. Phys. Rev. D 18, 4187 (1978)ADSCrossRefGoogle Scholar 26. Isgur, G. Karl, Positive-parity excited baryons in a quark model with hyperfine interactions. Phys. Rev. D 19, 2653 (1979)ADSCrossRefGoogle Scholar 27. L. Conci, M. Traini, Quark Momentum distribution in nucleons. Few Body Syst 8, 123 (1990)ADSCrossRefGoogle Scholar 28. J. Franklin, A model of baryons made of quarks with hidden spin. Phys. Rev. 172, 1807 (1968)ADSCrossRefGoogle Scholar 29. S. Capstick, N. Isgur, Baryons in a relativized quark model with chromodynamics. Phys. Rev. D 34, 2809 (1986)ADSCrossRefGoogle Scholar 30. M.H. Kalos, P.A. Whitlock, Monte Carlo Methods (Wiley-VHC Verlag GmbH & Co KGaA, Weinheim, 2008)CrossRefGoogle Scholar 31. M. Traini, L. Conci, U. Moschella, Deep inelastic parton distributions and the constituent quark model. Nucl. Phys. A 544, 731 (1992)ADSCrossRefGoogle Scholar 32. M. Traini, V. Vento, A. Mair, A. Zambarda, Constituent quarks and parton distributions. Nucl. Phys. A 614, 472 (1997)ADSCrossRefGoogle Scholar 33. A. Mair, M. Traini, Quark models and meson cloud in deep inelastic scattering. Nucl. Phys. A 628, 296 (1998)ADSCrossRefGoogle Scholar 34. S. Scopetta, V. Vento, M. Traini, Towards an unified picture of constituent and current quarks. Phys. Lett. B 421, 64 (1998). arXiv:hep-ph/9708262 ADSCrossRefGoogle Scholar 35. S. Scopetta, V. Vento, M. Traini, Polarized Structure functions in a constituent quark scenario. Phys. Lett. B 442, 28 (1998). arXiv:hep-ph/9804302 ADSCrossRefGoogle Scholar 36. P. Faccioli, M. Traini, V. Vento, Polarized parton distributions and light-front dynamics. Nucl. Phys. A 656, 400 (1999). arXiv:hep-ph/9808201 ADSCrossRefGoogle Scholar 37. M. Traini, Next-to-next-to-leading-order nucleon parton distributions from a light-cone quark model dressed with its virtual meson cloud. Phys. Rev. D 89, 034021 (2014). arXiv:1309.5814 [hep-ph]ADSCrossRefGoogle Scholar 38. J .P. Blaizot, High gluon densities in heavy ion collisions. Rep. Prog. Phys. 80(3), 032301 (2017). arXiv:1607.04448 [hep-ph]ADSMathSciNetCrossRefGoogle Scholar 39. V. Guzey, M. Strikman, M. Zhalov, Nucleon dissociation and incoherent \(J/\psi \) photoproduction on nuclei in ion ultraperipheral collisions at the large hadron collider. Phys. Rev. C 99, 015201 (2019). arXiv:1808.00740 [hep-ph]ADSCrossRefGoogle Scholar 40. H1 Collaboration, F. Aaron et al., Diffractive electroproduction of \(\rho \) and \(\phi \) mesons at HERA. JHEP 1005, 032 (2010). arXiv:0910.5831 [hep-ex] 41. H1 Collaboration, S. Chekanov et al., Exclusive \(\rho ^0\) production in deep inelastic scattering at HERA. PMC Phys. A 1, 6 (2007). arXiv:0708.1478 [hep-ex] 42. D. Bendova, J. Cepila, J.G. Contreras, Dissociative production of vector mesons at electron-ion colliders. Phys. Rev. D 99, 034025 (2019). arXiv:1811.06479 [hep-ph]ADSCrossRefGoogle Scholar 43. A.G. Shuvaev, K.J. Golec-Biernat, A.D. Martin, M.G. Ryskin, The effects of off diagonal parton distributions fixed by diagonal partons at small \(x\) and \(\xi \). Phys. Rev. D 60, 014015 (1999). arXiv:hep-ph/9902410 ADSCrossRefGoogle Scholar 44. M. Ahmady, R. Sandapen, N. Sharma, Diffractive \(\rho \) and \(\phi \) production at HERA using an AdS/QCD holographic light-front meson wavefunction. Phys. Rev. D 94, 074018 (2016). arXiv:1605.07665 [hep-ph]ADSCrossRefGoogle Scholar 45. C.H. Lee, H.-Y. Ryu, I. Zahed, Vector photoproduction using holographic QCD. Phys. Rev. D 98, 056006 (2018). arXiv:1804.09300 [hep-ph]ADSCrossRefGoogle Scholar Copyright information © The Author(s) 2019 Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3 Authors and Affiliations Marco Claudio Traini12Email authorJean-Paul Blaizot31.INFN-TIFPATrentoItaly2.Dipartimento di FisicaUniversità degli Studi di TrentoTrentoItaly3.Institut de Physique ThéoriqueUniversité Paris Saclay, CEAGif-sur-YvetteFrance

This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-019-6826-0.pdf

Marco Claudio Traini, Jean-Paul Blaizot. Diffractive incoherent vector meson production off protons: a quark model approach to gluon fluctuation effects, The European Physical Journal C, 2019, 327, DOI: 10.1140/epjc/s10052-019-6826-0