## Vector meson production in ultra-peripheral collisions at the LHC

R. Fiore a ♣ , L. Jenkovszky bglyph[star] , V. Libov c glyph[diamondmath] , and Magno V. T. Machado d ♠

a Dipartimento di Fisica, Universit` della Calabria and a

Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Cosenza

I-87036 Arcavacata di Rende, Cosenza, Italy Kiev, 03680 Ukraine c

b Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine

Deutsches Elektronen-Synchrotron, Hamburg, Germany

d HEP Phenomenology Group, CEP 91501-970, Porto Alegre, RS, Brazil

## Abstract

By using a Regge-pole model for vector meson production (VMP), successfully describing the HERA data, we analyse the correlation between VMP cross sections in photon-induced reactions at HERA and those in ultra-peripheral collisions at the Large Hadron Collider (LHC). Predictions for future experiments on production of J/ψ and ψ (2 S ) are presented.

♣

e-mail address: fiore@fis.unical.it

glyph[star] e-mail address:

jenk@bitp.kiev.ua

glyph[diamondmath] e-mail address: vladyslav.libov@desy.de

♠ e-mail address:

magnus@if.ufrgs.br

## 1 Introduction

Following the shut-down of the HERA collider at DESY (Hamburg), exclusive diffractive production of vector mesons in ultra-peripheral collisions of protons and nuclei became among the priorities of the present and future studies at the CERN LHC [1, 2, 3], triggering a large number of theoretical investigations [4, 5, 6, 7, 8]. For a relevant review paper see, e.g., Ref. [9]. First results on vector meson production (VMP), in particular on J/ψ , are already published [1, 2].

In this study of VMP at the CERN LHC we scrutinize possible changes in the energy dependence of the cross sections when moving from HERA to the LHC, in particular we are interested in the change from light vector mesons to heavy vector mesons ( φ, J/ψ, Υ, etc. ).

## 2 VMP at HERA

At HERA, VMP was studied in details both by the H1 and ZEUS collaborations. Most of the events were chosen in the kinematical region corresponding to diffractive scattering, which means that the processes can be described by a Pomeron exchange (see Fig. 1). Pomeron dominance is especially clean in J/ Ψ production, where, by the Zweig (OZI) rule, any exchange of secondary trajectories, made of quarks, is forbidden, thus leaving alone the uncontaminated Pomeron exchange. This does not mean that the dynamics is simple, but we have the opportunity to study in this class of reactions the nature of the Pomeron, a complicated and controversial object. The main problem is the twin nature of the Pomeron:

it seems to be 'soft' or 'hard' depending on the virtuality of the incident photon and/or the mass of the produced vector meson.

Figure 1: Diagrams of DVCS (a) and VMP (b); (c) DVCS (VMP) amplitude in a Reggefactorized form.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-389c8e61-a6a2-46fd-b8c0-9199b05acb33.jpg)

In most of the papers on the subject the existence of two Pomerons is assumed: a hard (or QCD) Pomeron, resulting from perturbative quantum chromodynamic calculations, and a soft one, somewhat misleadingly called 'non-perturbative'. We instead believe that there is only one Pomeron in Nature, but it has two components, whose relative weight is regulated by relevant ˜ Q 2 -dependent factors in front of the two, where the measure of the 'hardness', ˜ Q 2 = Q 2 + M 2 V , is the sum of the squared photon virtuality Q 2 and the squared mass M 2 V of the produced vector meson.

A specific model realizing this idea was constructed and tested against the experimental data recently (see Ref. [10] and earlier references therein). The relevant VMP amplitude reads

$$( 1 )$$

$$\ell ^ { 2 } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \\ A & ( s, t, Q ^ { 2 }, M _ { v } ^ { 2 } ) = \widetilde { A _ { s } } e ^ { - i \frac { \pi } { 2 } \alpha _ { s } ( t ) } \left ( \frac { s } { s _ { 0 } } \right ) ^ { \alpha _ { s } ( t ) } e ^ { b _ { s } t - n _ { s } \ln \left ( 1 + \widetilde { Q } ^ { 2 } _ { s } \right ) } \Delta _ { Q ^ { 2 } _ { s } } \right ) \\ & + \widetilde { A _ { h } } e ^ { - i \frac { \pi } { 2 } \alpha _ { h } ( t ) } \left ( \frac { s } { s _ { 0 } } \right ) ^ { \alpha _ { h } ( t ) } e ^ { b _ { h } t - ( n _ { h } + 1 ) \ln \left ( 1 + \widetilde { Q } ^ { 2 } _ { s } \right ) + \ln \left ( \frac { \widetilde { Q } ^ { 2 } _ { s } \right ) } { Q ^ { 2 } _ { h } } \right ) }, \\ \ell ^ { 2 } \cdot \text{and} \, \alpha, ( t ) \, \text{are the off and hard Domargon troisiocoftorio} \, \alpha = W ^ { 2 } \, \text{W} \,.$$

where α s ( ) t and α h ( ) t are the soft and hard Pomeron trajectories, s = W 2 , W being the energy of the VMP. Let us stress that the Pomeron is unique in all reactions, but its components (and parameters) vary. Examples with detailed fits can be found in the papers of Ref. [10].

The integrated (called also total) VMP cross section is given by

$$\sigma _ { \gamma p \rightarrow V p } ( \widetilde { Q } ^ { 2 }, W ) = \int _ { t _ { m } } ^ { t _ { t h r } } \frac { d \sigma ( \widetilde { Q } ^ { 2 }, W, t ) } { d t },$$

where the upper limit is t thr = 0 GeV 2 and the lower limit is the kinematical one, t m = -s/ 2. The total cross section can be simply calculated, without integration for an exponential diffraction cone, according to the formula

$$\sigma _ { e l } ( s ) = \frac { 1 } { B ( s ) } \frac { d \sigma } { d t } \Big | _ { t = 0 },$$

where B is the forward slope.

Since our primary goal is the comparison between the energy dependence of VMP at HERA and the LHC, we start with very simple ans¨ atz for the γp → V p cross section, postponing the use of the advanced model given by Eq. (1) to a future study.

## 3 VMP at the LHC

## 3.1 Distribution in rapidity

Vector meson production (VMP) cross section in hadronic collisions (see the relative Feynman diagram in Fig. 2) can be written in a factorized form, according to Refs. [5, 9] (e.g. Eqs. (1) and (9) in [5]a)). The distribution in the rapidity Y of the production of a vector meson V in the reaction h 1 + h 2 → h V h , 1 2 (where h may be a hadron, e.g. proton, or a nucleus, pPb, PbPb,...) is calculated according to a standard prescription based on the factorization of the photon flux and photon-proton cross section (see below).

Figure 2: Feynman diagram of vector meson production in a hadronic collision.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-77232ca7-d631-4ce6-a6ac-86979040a510.jpg)

Generally speaking, the γp differential cross section depends on three variables: the total energy of the γp system, W , the squared momentum transfer at the proton vertex, t , and ˜ Q 2 = Q 2 + M 2 V , where, as previously said, Q 2 is the photon virtuality and M V is the mass of the produced vector meson. Since, by definition, in ultra-peripheral collisions we have b &gt;&gt; R 1 + R 2 , where b is the impact parameter, i.e. the closest distance between the centres of the colliding particles/nuclei with radii R i (i = 1,2), photons are nearly real, Q 2 glyph[similarequal] 0, and M 2 V remains the only measure of 'hardness'. Notice that this might not be true for peripheral collisions, where b ∼ R 1 + R 2 , and in the Pomeron or Odderon exchange instead of the photon. Finally, the t -dependence (shape of the diffraction cone) is known to be nearly exponential. It can be either integrated, or kept explicit. Extending this parametrization to include a t -dependent exponential is easy (see below). Concerning total cross section σ γp → V p ( ˜ Q , t, W 2 ), it is well known from HERA.

We use a simple parametrization of the σ γp → V p ( W ) cross section, suggested by Donnachie and Landshoff [12]: σ γp → V p ( W ) ∼ W , δ δ ≈ 0 8 (more involved models, e.g. . of Ref. [10, 11] will be considered elsewhere).

The differential cross section as function of the rapidity Y reads

$$\frac { d \sigma ( h _ { 1 } + h _ { 2 } \to h _ { 1 } + V + h _ { 2 } ) } { d Y } = \omega _ { + } \frac { d N _ { \gamma / h _ { 1 } } ( \omega _ { + } ) } { d \omega } \sigma _ { \gamma h _ { 2 } \to V h _ { 2 } } ( \omega _ { + } ) \\ + \omega _ { - } \frac { d N _ { \gamma / h _ { 2 } } ( \omega _ { - } ) } { d \omega } \sigma _ { \gamma h _ { 1 } \to V h _ { 1 } } ( \omega _ { - } ),$$

$$+ \, \omega _ { - } \frac { d N _ { \gamma / h _ { 2 } } ( \omega _ { - } ) } { d \omega } \sigma _ { \gamma h _ { 1 } \to V h _ { 1 } } ( \omega _ { - } ),$$

where dN γ/h ( ω ) dω is the 'equivalent' photon flux [9], defined as (see, for instance, Ref. [5] and references therein)

$$\frac { d N _ { \gamma / h } ( \omega ) } { d \omega } = \frac { \alpha _ { e m } } { 2 \pi \omega } [ 1 + ( 1 - \frac { 2 \omega } { \sqrt { s } } ) ^ { 2 } ] ( \ln \Omega - \frac { 1 1 } { 6 } + \frac { 3 } { \Omega } - \frac { 3 } { 2 \Omega ^ { 2 } } + \frac { 1 } { 3 \Omega ^ { 3 } } ),$$

and σ γp → V p ( ω ) is the total cross section of the vector meson photoproduction subprocess (the same as e.g. at HERA, see Ref. [11]). Here ω is the photon energy, ω = W / 2 γp 2 √ s , where √ s denotes now the c.m.s. energy of the proton-proton system; ω min = M / 2 V (4 γ L m , p ) where γ L = √ s/ (2 m p ) is the Lorentz factor (Lorentz boost of a single beam); e.g., for pp at the LHC for √ s = 7 TeV, γ L = 3731. Furthermore, Ω = 1 + Q /Q 2 0 2 min , with Q 2 min = ( ω/γ L ) 2 and Q 2 0 = 0 71 . GeV 2 , Y ∼ ln(2 ω/M V ). The differential cross section in W can be calculated as dσ/dW = ( dσ/dY ) /W .

For definiteness we fix: a) the colliding particles are protons; b) the produced vector meson V is J/ψ , and 3) the collision energy √ s = 7 TeV. We comprise the constants in A = α em / (2 π , ) c = Q γ 2 0 2 L . (Note that the shape of the distribution in Y is very sensitive to the value and the sign of the constant c ). The i = ± signs of ω correspond to the first or second term in Eq. (4), respectively, ω ± ∼ e ± Y .

## 3.2 Corrections for rapidity gap survival probabilities

The above results may be modified by initial and final state interactions, alternatively called as rescattering corrections. The calculation of these corrections is by far not unambiguous, the result depending both on the input and on the unitarization procedure chosen. The better (more realistic) the input, the smaller the unitarity (rapidity gap survival probability) corrections. Since this is a complicated and controversial issue per se deserving special studies beyond the scope of the present paper, to be coherent with the 'common trend', here we use only familiar results from the literature: the standard prescription is to multiply the scattering amplitude (or the cross section) by a factor (smaller than one), depending on the energy and eventually other kinematical variables [6]. In this work we used the constant value of 0.8.

## 4 Results of the fitting procedure

In this Section theoretical predictions are presented and compared to experimental data. IN our alculations we use two models for the exclusive J/ Ψ photoproduction cross section, σ γp → J/ Ψ : the simple power-law, σ γp → J/ Ψ ∼ W δ , with δ = 0 8, and the so-called Reggeometric . model. The latter, in its first version, was suggested and applied to deeply virtual Compton scattering (DVCS) in Ref. [11]. Apart from W and t , it contains also the dependence on virtuality, Q 2 . The model was fitted to the HERA data on DVCS, but it can be applied also to the VMP by refitting its parameters. In this paper we use the first (simpler) version of the Reggeometric model of VMP and DVCS, suggested in Ref. [10].

Figure 3: Differential cross section of exclusive J/ Ψ production at the LHC as function of the rapidity Y , plot (a), and of the photon-proton c.m.s. energy W , plot (b). The red curves and the blue curves represent, respectively, the parameterization of the photon-proton cross section according to the simple power-law and the Reggeometric model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-f973b232-2c50-4b94-ae63-c845f9b845c5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-634c3e47-b8b4-4a2a-b608-b2716478d262.jpg)

## 4.1 Fitting the Reggeometric model to LHCb data

The first version of the Reggeometric model [10] a) applies to photoproduction ( Q 2 = 0). The total photoproduction cross section, Eq. (11) in Ref. [10] a), is given by

$$\sigma _ { \gamma p \to J / \Psi p } = A _ { 0 } ^ { 2 } \frac { ( W / W _ { 0 } ) ^ { 4 ( \alpha _ { 0 } - 1 ) } } { ( 1 + \tilde { Q } ^ { 2 } / Q _ { 0 } ^ { 2 } ) ^ { 2 n } [ 4 \alpha ^ { \prime } \ln ( W / W _ { 0 } ) + 4 \left ( \frac { a } { Q ^ { 2 } } + \frac { b } { 2 m _ { N } ^ { 2 } } \right ) ] },$$

where ˜ Q 2 = Q 2 + M 2 V . The parameters, fitted [10] a) to the J/ψ photoproduction, quoted in Table II of Ref. [10] a), are: A 0 = 29 8 . ± 2 8 . , Q 2 0 = 2 1 . ± 0 4 . , n = 1 37 . ± 0 14 . , α 0 = 1 20 . ± 0 02 . , α ′ = 0 17 . ± 0 05 . , a = 1 01 . ± 0 11 . , b = 0 44 . ± 0 08 . , W 0 = 1 and the relevant dimensions here again are implied. Note that compared to the original formula, s was replaced by W 2 , since W is used in this paper to denote the photon-proton centre-of-mass energy (in contrast, √ s in this paper is the proton-proton centre-of-mass energy). √

Fig. 3(a) shows the predicted differential cross section of J/ Ψproduction at LHC ( s = 7 TeV) as a distribution in rapidity Y . On the whole, the power law and the Reggeometric model yield similar distributions, the latter is flatter though. Fig. 3(b) shows the W -dependence of J/ψ photoproduction differential cross section. Again, the two models give generally similar behaviours.

The LHCb collaboration at the LHC has recently [1, 2] measured the production cross section of J/ Ψas a function of the rapidity Y. Fig. 4(a) shows a comparison of our calculations with these data. As it can be seen, the data are somewhat steeper than both the curves. LHCb also extracted the basic photon-proton photoproduction cross section as a function of W from the data. The result is compared to our predictions in Fig. 4(b). Also the ZEUS and H1 data are shown.

Figure 4: The differential cross section of J/ Ψ production at LHC as a function of the rapidity Y together with the LHCb data is shown in plot (a). Also shown, in plot (b), is the J/ Ψ photoproduction ( γp → J/ Ψ ) cross section as a function of the photon-proton c.m.s. p energy compared with the LHCb, ZEUS and H1 data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-c8c16fee-5b67-4d18-8647-2467a84ef4b6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-88c849f8-9f00-4e18-8555-fd8fd1f31bfe.jpg)

## 4.2 Fitting the power law to the LHCb data

As discussed above, the LHCb rapidity data have a steeper shape than our predictions. Here we investigate the role of the power δ in the rapidity cross section. In particular, we perform a least-square fit to these data. The free parameters are the power δ and the overall normalization. The best fit gives δ =0.37 and is shown in Fig. 5. Indeed, the prediction obtained with this value gives a much better description of the rapidity distribution. It is interesting make a comparison with other parameterizations, in particular by using the logarithmic parameterization, which leads to a less steep line than that obtained in the power law case . This is shown in Fig. 6. The logarithmic parameterization provides a good description as well. Fig. 7 shows, for completeness, the predictions obtained from other groups, together with the LHCb rapidity data. The agreement with the data is reasonable.

Our predictions for the J/ Ψ photoproduction, also compared to the H1 and ZEUS data, are shown in Fig. 8. As it can be seen, the Reggeometric model provides the best description. The power law with the standard δ = 0 8 agrees reasonably as well. . However, the one with δ = 0 37, as well as the logarithmic parameterization fail to describe the data. .

Ψ(2 S ) differential cross section as a function of rapidity Y is compared with the LHCb data in Fig. 9. The power law model was used with parameters as obtained from the fit. Good description is observed.

Figure 5: Differential cross section of J/ψ production at the LHC as function of the rapidity Y , according to the power law together with LHCb data. Both calculations use the simple power law model, with δ = 0 8 (red line) and . δ = 0 37 (green line). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-618b905f-d054-40c7-a0ac-cdf2998fe23d.jpg)

Figure 6: Differential cross section of J/ψ production at the LHC as a function of rapidity Y , according to the various models for the photon-proton cross section used: power law with δ = 0 8 and . δ = 0 37 (red and green lines, respectively), logarithmic parameterization (pink . line) and the Reggeometric model (blue line). Also the LHCb data are shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-4f890cdd-d4be-4ce5-bfd5-f139ea973990.jpg)

Figure 7: Differential cross section of J/ψ production at LHC as a function of rapidity Y together with the LHCb data [2].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-75d4a940-6df6-4b48-aaf8-905e27e33cb6.jpg)

Figure 8: J/ψ photoproduction ( γp → J/ Ψ ) cross section as a function of the photonp proton c.m.s. energy compared with the LHCb, ZEUS and H1 data. Other details are as in Fig. 6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-78ce2bc1-62dd-4142-b42a-4156837237e4.jpg)

Figure 9: Ψ(2 S ) differential cross section as a function of rapidity Y compared with the LHCb data. The power law model was used with parameters as shown in the figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0530v1-figure-17a4fa8a-2e69-4cb4-b002-b38225062eab.jpg)

## 5 Conclusions and Outlooks

In this paper predictions for exclusive J/ Ψmeson production in ultra-peripheral collisions at the LHC are presented and compared with the recent experimental data collected by the LHCb. The HERA experimental data are also shown and compared with the H1 and ZEUS measurements. The simple power-law and a more advanced geometric model, the so called Reggeometric model, were used to describe the photon-proton cross section. The LHCb for the cross section as a function of the rapidity are steeper than predictions. A better description can be obtained by tuning the power, however this makes the predictions inconsistent with the HERA data.

The present study will be extended by

- · the inclusion of the t dependence of the differential cross section - both exponential (corresponding to linear Regge trajectories) and with deviations from the exponential (non-linear Regge trajectories);
- · the inclusion of the dependence of the σ γp → V p cross section on Q 2 , negligible in γ -but important in the case of Reggeon (Pomeron, Odderon,...) exchanges;
- · studies of inelastic processes, i.e. those in which additional particles are produced due to either gluon radiation and/or (c,d) proton dissociation;
- · more advanced studies of corrections due to the rapidity gap survival probability;
- · generalizations to include nuclear collitions.

## Acknowledgements

Wethank A. Salii for discussions. L.J. is grateful to the Dipartimento di Fisica dell'Universit'a della Calabria and the Istituto Nazionale di Fisica Nucleare - Gruppo Collegato di Cosenza, where part of this work was done, for their hospitality and support. He was supported also by the grant 'Matter under extreme conditions' of the National Academy of Sciences of Ukraine, Dept. of Astronomy and Physics and by the DOMUS Curatorium of the Hungarian Academy of Sciences.

## References

- [1] LHCb Collab., R. Aaij, J. Phys. G40 (2013) 045001, arXiv:1301.7084.
- [2] LHCb Collab., R. Aaij, J. Phys. G41 (2014) 055002, arXiv:1401.3288.
- [3] ATLAS Collab., CERN-PH-EP-2014-095; arXiv:1407.5535.
- [4] A. Sch¨fer, L. Mankiewicz, and O. Nachtmann, Phys. Lett. a B272 (1991) 419.
- [5] a) V.P. Goncalves and M.M. Machado, Phys. Rev. D85 (2012) 054019, arXiv:1112.3500;
- b) V.P. Goncalves and M.M. Machado, Eur. Phys. J. C72 (2012) 2231, arXiv:1207.5273;
- c) V.P. Gonsalves and W.K. Suter, Eur. Phys. J. A47 (2011) 117, arXiv: 1004.1952; d)
- V.P. Goncalves and M.M. Machado, Phys. Rev. C84 (2011) 011902, arXiv:1106.3036;
- e) V.P. Goncalves, Nucl. Phys. A902 (2013) 32, arXiv:1211.1207; f) M.B. Gay Ducati,

- M.T. Griep, and M.V.T. Machado, Phys. Rev. D88 (2013) 017504, arXiv:1305.4611; g) V.P. Goncalves and M.M. Machado, Phys. Rev. C80 (2009) 054901, arXiv:0907.4123. e) G. Sampaio dos Santos, M.V.T. Machado, Light vector meson photoproduction in hadron-hadron and nucleus-nucleus collisions at the energies available at the CERN Large Hadron Collider , arXiv:1407.4148; G. Sampaio dos Santos and M.V.T. Machado, azXiv: 1407.4148 [hep-ph].
- [6] a) V.A. Khoze, A.D. Martin, and M.G. Ryskin, Eur. Phys. J. C24 (2002) 459, hepph/0201301; b) S.P. Jones, A.D. Martin, M.G. Ryskin, and T. Teubner, Probes of the small x gluon via exclusive J/psi and Upsilon , arXiv:1307.7099; c) S.P. Jones, A.D. Martin, M.G. Ryskin, and T. Teubner, J. Phys. G41 (2014) 055009, arXiv:1312.6795.
- [7] a) A. Bzdak, L. Motyka, L. Szymanowski, and J.-R. Cudell, Phys. Rev. D75 (2007) 094023, hep-ph/0702134; b) L. Motyka and G. Watt, Phys. Rev. D78 (2008) 014023, arXiv:0805.2113.
- [8] a) W. Sch¨fer and A. Szczurek, Phys. Rev. a D76 (2007) 094014, arXiv: 0705.2887; b) W. Sch¨ afer, G. Slipek, and A. Szczurek, Phys. Lett. B688 (2010) 185, arXiv:0805.0717.
- [9] G. Baur et al. , Phys. Rept. 364 (2002) 359, hep-ph/011221; K. Hencken et al. , Phys. Rept. 458 (2008) 1, arXiv:0706.3356.
- [10] a) S. Fazio, R. Fiore, L. Jenkovszky, and A. Salii, Acta Phys. Polonica B44 (2013) 1333, arXiv:1304.1891; b) S. Fazio, R. Fiore, L. Jenkovszky, and A. Salii, arXiv:1312.5683, Phys. Rev. D , to be published.
- [11] M. Capua et al. , Phys. Lett. B645 (1997) 161, hep-ph/0605319.
- [12] A. Donnachie and P. Landshoff, Phys. Lett. B296 (1992) 227.
- [13] R. Fiore, L.L. Jenkovszky, and F. Paccanoni, EPJ C10 (1999) 461, hep-ph/9812458.
Roberto Fiore Laszlo Jenkovszky Vladyslav Libov Magno Machado
Vector meson production in ultra-peripheral collisions at the LHC
By using a Regge-pole model for vector meson production (VMP), successfully
describing the HERA data, we analyse the correlation between VMP cross sections
in photon-induced reactions at HERA and those in ultra-peripheral collisions at
the Large Hadron Collider (LHC). Predictions for future experiments on
production of $J/\psi$ and $\psi(2S)$ are presented.