Double vector meson production in γγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \gamma $$\end{document} interactions at hadronic colliders

In this paper we revisit the double vector meson production in γγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \gamma $$\end{document} interactions at heavy ion collisions and present, by the first time, predictions for the ρρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho \rho $$\end{document} and J/ΨJ/Ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\Psi J/\Psi $$\end{document} production in proton–nucleus and proton–proton collisions. In order to obtain realistic predictions for rapidity distributions and total cross sections for the double vector production in ultra peripheral hadronic collisions we take into account the description of γγ→VV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \gamma \rightarrow VV$$\end{document} cross section at low energies as well as its behavior at large energies, associated to the gluonic interaction between the color dipoles. Our results demonstrate that the double ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} production is dominated by the low energy behavior of the γγ→VV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \gamma \rightarrow VV$$\end{document} cross section. In contrast, for the double J/Ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\Psi $$\end{document} production, the contribution associated to the description of the QCD dynamics at high energies contributes significantly, mainly in pp collisions. Predictions for the RHIC, LHC, FCC, and CEPC–SPPC energies are shown.

(photon-photon process) or it can interact directly with the other hadron (photon-hadron process) [43][44][45][46][47][48]. In these processes the total cross section can be factorized in terms of the flux of equivalent photons from the hadron projectile and the photon-photon or photon-target production cross section. In this paper we will focus on two-photon interactions in hadronic collisions. Experimental results on exclusive twophoton production of W + W − and + − pairs in γ γ interactions reported by the CMS and ATLAS Collaborations [9][10][11][12] have demonstrated that it is possible to measure such events with the experimental apparatus already available at the LHC, allowing for novel studies of QCD at very high energies and searches for Beyond Standard Model Physics (see, e.g., Ref. [49][50][51][52][53][54][55][56][57][58][59]). This motivates us to revisit the analysis of double vector meson production in ultra peripheral heavy ion collisions performed some time ago in Refs. [60][61][62], taking into account recent improvements in the description of the γ γ → V V (V = ρ, J/ ) cross section at low [63,64] and at high [65] energies. In this work we will derive, for the first time, realistic predictions for double meson production in γ γ interactions at pp and p A collisions at RHIC and LHC energies as well as for the proposed energies for the Future Circular Collider (FCC) at CERN [66] and in the Circular Electron Positron Collider with a subsequent Super proton-proton Collider (CEPC-SPPC) in China [67]. Our goal is to determine if this process, after the inclusion of the recent theoretical improvements, can be used to study the QCD dynamics at high energies, as originally proposed in Ref. [60]. As we will show, this remains true for double J/ production, mainly in pp collisions.
Let us start our analysis presenting a brief review of the main formulas to describe double vector meson production in γ γ interactions at hadronic colliders. In the Equivalent Photon Approximation (EPA) [43][44][45][46][47][48] the total cross section for this process can be written as where √ s is center-of-mass energy of the h 1 h 2 collision (h i = p, A), ⊗ characterizes a rapidity gap in the final state and W = √ 4ω 1 ω 2 is the invariant mass of the γ γ system. Moreover, N (ω, b) is the equivalent photon spectrum of photons with energy ω at a distance b from the hadron trajectory, defined in the plane transverse to the trajectory. The spectrum can be expressed in terms of the charge form factor F as follows: where γ L is the Lorentz factor. The factor S 2 abs (b) is the absorption factor, given in what follows by where R h i is the radius of the hadron h i (i = 1, 2). The presence of this factor in Eq. (1) excludes the overlap between the colliding hadrons and allows one to take into account only ultra peripheral collisions. Remembering that the photon energies ω 1 and ω 2 are related to W and the rapidity (Y = 1 2 (y V 1 + y V 2 )) of the outgoing double meson system by the total cross section can be expressed by (for details see e.g. Ref. [68]) It is important to emphasize that in EPA we disregard the photon virtualities, which is a good approximation, mainly for ions, since the typical virtualities are <1/R h . Moreover, the highest energy of the photons is of the order of the inverse Lorentz contracted radius of the hadron ≈ γ L /R h , with the spectra decreasing exponentially at larger energies. Consequently, for the same Lorentz factor, we have Finally, due to the Z 2 dependence of the photon spectra, for a fixed W the following hierarchy is valid for processes induced by γ γ interactions: In order to estimate this cross section we must describe the γ γ → V 1 V 2 interaction in a large energy range. In the following we will assume that where the LE term is associated to the description of the cross section at low energies W 10 GeV, while the HE term describes the region of larger values of W . Double vector meson production at low energies has been discussed and improved in Refs. [63,64]. As in Ref. [63], we will evaluate double ρ production in γ γ interactions directly from the experimental measurements using a fit to the world data, which describes the experimental data in the region of few GeV. In particular, we will take into account the huge enhancement close to the threshold observed in the data, which is not yet well understood. As demonstrated in Ref. [63] this contribution determines the behavior of double ρ production in A A collisions. In the case of double J/ production, as in Ref. [64], we will consider the contribution associated to the box diagrams, calculated in the heavy quark non-relativistic approximation. In Ref. [64] the authors have also estimated the contribution associated to the two-gluon exchange, which implies a γ γ → J/ J/ cross section independent of the energy. In our analysis, we will not include this contribution in the low energy term, since it is the leading order term in the dipole-dipole interaction present in our formalism to treat the high energy term discussed in the following. Finally, it is important to emphasize that one of the main conclusions from Ref. [64] is that in PbPb collisions the box mechanism significantly dominates over the two-gluon exchange one.
The description of double vector meson production in γ γ interactions at high energies has attracted the attention of several theoretical groups in the last years, with the cross section being estimated in different theoretical frameworks [60,61,[63][64][65][69][70][71][72][73][74][75][76][77][78][79][80], as, for instance, the solution of the BFKL equation and impact factors at leading and next-toleading orders. In particular, in Ref. [65] we have estimated the total γ (Q 2 1 ) + γ (Q 2 2 ) → V 1 + V 2 cross sections for V i = ρ, φ, J/ψ, and ϒ in the color dipole formalism considering the improved treatment of the dipole-dipole cross section proposed in Ref. [81]. We have also taken advantage of the progress in the knowledge of vector meson wave functions. Moreover, in Ref. [65] we have taken into account the non-linear effects in the QCD dynamics, which are expected to be present at large energies. An important aspect of the analysis presented in Ref. [65] is that the main ingredients are constrained by LEP and HERA data. In particular, assum- ing the values for the slope parameter B V 1 V 2 proposed in Ref. [78], it is possible to obtain parameter-free predictions for the ρρ and J/ J/ production cross sections at high energies.
In the case of ultra peripheral hadronic collisions, double meson production is induced by the interaction of real photons and can be represented by the diagram in Fig. 1 in the color dipole formalism. In this approach the γ γ → V 1 V 2 interaction at hadronic colliders can be seen as a succession in time of four factorizable subprocesses (see Fig. 1): (1) the photons are emitted by the incident hadrons, (2) the photons fluctuate into quark-antiquark pairs (the dipoles), (3) these color dipoles interact and, (4) the pairs convert into the vector meson final states. In particular, the γ γ → V 1 V 2 cross section can be expressed as follows: where we have approximated the t-dependence of the differential cross section by an exponential with B V 1 V 2 being the slope parameter. The imaginary part of the amplitude at zero momentum transfer A(W 2 , t = 0) reads where γ and V i are the light-cone wave functions of the photon and vector meson, respectively, and T the transverse polarization. The variable r 1 defines the relative transverse separation of the pair (dipole) and z 1 (1 − z 1 ) is the longitudinal momentum fraction of the quark (antiquark). Similar definitions hold for r 2 and z 2 . The variable Y is the rapidity and will be defined later. The basic blocks are the photon wave function, γ , the meson wave function, V , and the dipole-dipole cross section, σ d d . In contrast to the photon wave function, which is well known in the literature (see e.g. [82]), the description of the vector meson wave functions is still a subject of debate. The simplest approach 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 [82][83][84][85][86]. As in Ref. [65] we will assume that the overlap between the photon and the vector meson wave function, for the transversely polarized case, is given by (for details see Ref. [82]) whereê f is the effective charge of the vector meson, m f is the quark mass, N c = 3, 2 = z(1 − z)Q 2 + m 2 f and φ T (r, z) defines the scalar part of the vector meson wave function. In the following we will consider the Gauss-LC model for φ T (r, z), which is then given by The parameters N T and R T are determined by the normalization condition of the wave function and by the decay width (see Ref. [65] for details). The other main input to calculate the γ γ → V 1 V 2 cross section is the dipole-dipole cross section, σ d d . At lowest order, the dipole-dipole interaction can be described by the two-gluon exchange between the dipoles, with the resulting cross section being energy independent (see, e.g., Ref. [87]). The inclusion of the leading corrections associated to terms ∝ log(1/x) (as described by the BFKL equation) leads to a power-law energy behavior of the cross section, which violates the unitarity at high energies. Unitarity corrections were introduced in Ref. [88], considering the color dipole picture and independent multiple scatterings between the dipoles. These corrections were also addressed in Refs. [89,90] in the context of the Color Glass Condensate (CGC) formalism [91][92][93][94][95].
In the eikonal approximation the dipole-dipole cross section can be expressed as follows: where N (r 1 , r 2 , b, Y ) is the scattering amplitude of the two dipoles with transverse sizes r 1 and r 2 , relative impact parameter b and rapidity separation Y . The interaction of two dipoles of similar sizes is still an open question (see, e.g., Ref. [96]). In a first approximation, it is useful to express N in terms of the solution of the Balitsky-Kovchegov (BK) equation (obtained disregarding the b dependence), which has been derived considering an asymmetric frame where the projectile has a simple structure and the evolution occurs in the target wave function [97][98][99]. A shortcoming of this approach is that, although the unitarity of the S-matrix (N ≤ 1) is respected by the solution of the BK equation, the associated dipole-dipole cross section can still rise indefinitely with the energy, even after the black disk limit (N = 1) has been reached at central impact parameters, due to the non-locality of the evolution. In Ref. [81] we have proposed a more elaborated model for the impact parameter dependence in order to obtain more realistic predictions for the dipole-dipole cross section. Basically, we assumed that only the range b < R, where R = Max(r 1 , r 2 ), contributes to the dipole-dipole cross section, i.e. we assumed that N is negligibly small when the dipoles have no overlap with each other (b > R). Therefore the dipole-dipole cross section can be expressed as follows [81]: (12) where N (r, Y ) is the forward scattering amplitude, which can be obtained as a solution of the BK equation disregarding the impact parameter dependence or from phenomenological models that describe the HERA data. The explicit form of σ dd reads where Y i = ln(1/x i ) and As in Refs. [65,81] we will consider in our calculations the IIM-S model [100,101] for the forward scattering amplitude, which is based on the solutions of the BK equation at small and large dipoles, and is given by where a and b are determined by continuity conditions at r Q s (x) = 2, γ s = 0.6194, κ = 9.9, λ = 0.2545, Q 2 0 = 1.0 GeV 2 , x 0 = 0.2131 × 10 −4 and N 0 = 0.7. As demonstrated in Ref. [81], using this model we can describe the LEP data for the total γ γ cross sections and photon structure functions.
In the following we present our predictions for the rapidity distributions and total cross sections for ρ-ρ and J/ -J/ production through γ γ interactions in pp, pPb, and PbPb collisions. In order to estimate the contribution of the gluonic part associated with the dipole-dipole interaction, we will compare the full predictions, obtained considering the low and high energy contributions (denoted Low Energy + IIM-S hereafter) with those without the high energy contribution (Low Energy hereafter). In order to estimate the equivalent photon spectra for A = Pb, we will consider the a monopole form factor F(q 2 ) = 2 /( 2 + q 2 ), with = 0.088 GeV adjusted to reproduce the root mean square (rms) radius of the nucleus. Moreover, we will assume that R A = 1.2 A 1 3 fm. In the proton case, we will consider that F(q 2 ) = 1/[1 + q 2 /(0.71 GeV 2 )] 2 and R p = 0.7 fm. Finally, as in [78] we will assume B ρρ = 10 GeV −2 and B ψψ = 0.44 GeV −2 . Our results for the rapidity distributions are presented in Figs. 2, 3 and 4. In the case of PbPb collisions, presented in Fig. 2, we see that the Low Energy and Low Energy + IIM-S predictions are almost identical for double ρ production for all considered energies, which indicates that the gluonic contribution for this process is very small, in agreement with the conclusion obtained in Ref. [63]. This result is can be related to the energy behavior of the γ γ → ρρ cross section at high energies, which presents a mild growth with W , expected in a process dominated by large dipoles, and to the fact that in A A collisions we are probing values of W ≤ 160 GeV for √ s = 5.5 TeV. In contrast, for double J/ production, we observe that the gluonic contribution increases with the energy, which is associated to the steep energy behavior of the γ γ → J/ J/ cross section. For √ s = 0.5 TeV the analysis of this process can be useful to probe the box mechanism. On the other hand, for √ s = 5.5 TeV the gluonic contribution implies an enhancement by a factor 2 of the rapidity distribution at Y = 0. When pPb collisions are considered, we obtain the asymmetric rapidity distributions presented in Fig. 3, which is expected since the nuclear equivalent photon spectra is enhanced by a factor Z 2 . Moreover, in this case the energy range probed in the γ γ interactions increases for W ≤ 160 GeV for √ s = 8.8 TeV. As a consequence, we observe that the difference between the Low Energy and Low Energy + IIM-S predictions starts to occur in the case of double ρ production and becomes appreciable for double J/ production. Finally, in Fig. 4 we present our results for pp collisions for different values of √ s. In this case the double vector meson production is induced by γ γ interactions with W ≤ 4500 GeV for √ s = 14 TeV. This large range of energies probed in the interaction implies that the gluonic contribution becomes very important for the description of double vector production, increasing the rapidity distributions at Y = 0 by a factor 2 in the case of ρρ production and by a factor 40 for double In agreement with our previous discussion, we can see that this contribution is small for double ρ production and appreciable for double J/ production, mainly in pp collisions. This result indicates that the analysis of double J/ production in ultra peripheral hadronic collisions can be useful to study the QCD dynamics at high energies, as originally suggested in Ref. [60]. In Tables 1 and 2 we present our predictions for total cross sections for double vector meson production in pp, pPb, and PbPb collisions for the energies of RHIC and LHC as well as for the conceptual design energies of the FCC [66] and CEPC-SPPC [67]. It is important to emphasize that our Low Energy + IIM-S predictions can be considered as a lower bound for the gluonic contribution, since other models for the dipole-dipole cross section or for the description of the QCD dynamics imply larger values for the γ γ → J/ J/ cross section (for a detailed discussion see Ref. [65]). Consequently, we believe that the analysis of this process is feasible in hadronic colliders. Additionally, considering the results from Ref. [65], which indicate that γ γ → V 1 V 2 cross sections for the ρ J/ , φ J/ , ρϒ, J ϒ, and ϒϒ production increase strongly with the energy, we can also expect that these final states could be analyzed in the future. As discussed in detail in Refs. [60][61][62]65], the study of these different final states is important to understand the transition between the soft and hard regimes of the QCD dynamics, since different dipole sizes are probed in each process.
Finally, let us summarize our main conclusions. In recent years, a series of studies have discussed in detail the computation of the total cross section and the exclusive production of different final states in γ γ interactions considering very distinct theoretical approaches. One of the basic motivations for these efforts is the possibility to study the behavior of QCD dynamics at high energies. The ideal laboratory for these studies is the scattering of two off-shell photons at high energy in e + e − colliders, which could be performed in the International Linear Collider (ILC). However, as the schedule for the construction and operation of this collider is still an open question, the analysis of alternative ways to study the γ γ interactions is an important theme. The study of double vector meson production in γ γ interactions in ultra peripheral heavy ion collisions as a probe of the QCD dynamics was proposed in Ref. [60] and developed in Refs. [61,62]. However, these studies focused only on the high energy regime and disregarded the low energy mechanisms for double vector production. As emphasized in Refs. [63,64], the contribution of these mechanisms is important in A A collisions, since the maximum center of mass energies probed in the γ γ interactions is not large and the main contribution of the equivalent photon spectrum comes from photons with low energy. However, these studies have disregarded the effects of the QCD dynamics discussed in Refs. [61,62] and recently updated in Ref. [65]. In this paper we have combined these two approaches and derived predictions for the γ γ → V V cross section which are valid in the full kinematical range. We have obtained realistic predictions for the total cross sections in hadronic collisions and estimated the relative contribution of the low and high energy regimes. In particular, the results for pp and pPb have been derived by the first time. Our results demonstrated that double ρ production is dominated by low energy mechanisms. On the other hand, the gluonic contribution for double J/ production strongly increases with the energy, the study of this process becomes feasible in hadronic collisions (mainly in pp collisions) and it may be useful to constrain the QCD dynamics at high energies, as proposed originally in Ref. [60].