Path integral treatment of coherence effects in charmonium production in nuclear ultra-peripheral collisions

Abstract

We present for the first time a revised study of charmonium production in nuclear ultra-peripheral collisions (UPC) based on a rigorous Green function formalism. This formalism allows for the proper incorporation of the effects of color transparency, as well as the quantum coherence inherent in the higher twist quark shadowing related to the $Q\bar{Q}$ Fock component of the photon. The significance of this effect gradually decreases towards forward and/or backward rapidities. In the LHC kinematic region we additionally incorporate within the same formalism the leading twist gluon shadowing corrections related to higher multi-gluon photon fluctuations. They represent a dominant source of nuclear phenomena in the mid-rapidity region. Model predictions for the rapidity distributions $d\sigma/dy$ are in good agreement with available UPC data on coherent charmonium production at RHIC and the LHC. They can also be verified by future measurements at the LHC, as well as at EIC.

1 Introduction

The study of heavy quarkonium production in heavy-ion collisions can help us to extend our understanding of the QCD dynamics, as well as manifestations of various nuclear effects, such as quantum coherence, color transparency, gluon shadowing and gluon saturation. The production of charmonia in photo-nuclear reactions ($Q^{2}\approx 0$) is recently intensively studied in ultra-peripheral Pb-Pb collisions at the Large Hadron Collider (LHC) and in Au-Au UPC at Relativistic Heavy Ion collider (RHIC). Here, the gluon shadowing (GS) corrections at LHC energies and the reduced effects of quantum coherence in the forward/backward rapidity region are often neglected or incorporated inaccurately in most of the current calculations. In this contribution, we refer on our study of coherent (elastic) photoproduction of charmonia, $\gamma A\to VA$ [V = $J\!/\!\psi(1S),\psi^{\,\prime}(2S)$], in UPC, extended to charmonium electroproduction, $\gamma^{*}A\to VA$, in expected Electron-Ion-Collider (EIC) kinematic region within the light-front (LF) color dipole approach. We focus on the proper treatment of reduced effects of the coherence length using the Green function technique.

2 Green function formalism for the coherent photo-nuclear reaction

The cross section for photo-nuclear reaction in a heavy-ion UPC, derived in the one-photon-exchange approximation and expressed in the rest frame of the target nucleus $A$ [1], is defined as follows:

$$k\frac{d\sigma}{dk}=\int\,d^{2}\tau\int\,d^{2}b\,\,n(k,\vec{b}-\vec{\tau},y)\,\frac{d^{2}\sigma_{A}(s,b)}{d^{2}b}~{}~{}+~{}~{}\Bigl{\{}y\rightarrow-y\Bigr{\}}\,.$$ (1)

Here the rapidity variable $y=\ln\bigl{[}s/(M_{V}\sqrt{s_{N}}~{})\bigr{]}\approx\ln\bigl{[}(2kM_{N}+M_{N}^{2})/(M_{V}\sqrt{s_{N}}~{})\bigr{]}$ with $M_{N}$ and $M_{V}$ being the nucleon and vector meson mass, respectively, $\sqrt{s_{N}}$ is the collision energy and $k$ is the photon energy related to the square of the photon-nucleon center-of-mass (c.m.) energy $W^{2}=2M_{N}k+M_{N}^{2}-Q^{2}\approx s-Q^{2}$, where $s$ is the photon-nucleon c.m. energy squared. The variable $\vec{\tau}$ is the relative impact parameter of a nuclear collision and $\vec{b}$ is the impact parameter of the photon-nucleon collision relative to the center of one of the nuclei. It holds that $\tau>2R_{A}$, for a UPC of identical nuclei with the nuclear radius $R_{A}$. The variable $n(k,\vec{b})$ represents the photon flux induced by a projectile nucleus. The corresponding formula for the coherent production cross section has the following form [2, 3]

$$\frac{d^{2}\sigma_{A}^{coh}(s,b)}{d^{2}b}=\frac{1}{4}\,\Bigl{|}\int\limits_{-\infty}^{\infty}\,dz\,\rho_{A}({b},z)\,H_{1}(s,{b},z)\,\Bigr{|}^{2}\,,$$ (2)

where $\rho_{A}(b,z)$ is the nuclear density distribution employed in the realistic Wood-Saxon form. The function $H_{1}(s,b,z)$ can be expressed within a rigorous path integral technique as

$$\!\!\!\!\!\!\!\!H_{1}(s,b,z)\!\!\! = \!\!\!\!\!\int_{0}^{1}\!\!d\alpha\int\!\!d^{2}r_{1}\,d^{2}r_{2}\,\Psi^{*}_{V}(\vec{r}_{2},\alpha)\,G_{Q\bar{Q}}(z^{\prime}\!\!\to\!\infty,\vec{r}_{2};z,\vec{r}_{1})\,\sigma_{Q\bar{Q}}(r_{1},s)\,\Psi_{\gamma(\gamma^{*})}(\vec{r}_{1},\alpha)\,.$$ (3)

Here, $\Psi_{V}(\vec{r},\alpha)$ is the LF wave function for heavy quarkonium, obtained by the Terent’ev boosting prescription [4] and including the correction due to the Melosh spin rotation effect [5], and $\Psi_{\gamma(\gamma^{*})}(\vec{r},\alpha)$ is the LF distribution of the $Q\bar{Q}$ Fock component of the quasi-real (transversely polarized) or virtual photon. Variable $\vec{r}$ is the transverse separation of the $Q\bar{Q}$ fluctuation (dipole) and $\alpha=p_{Q}^{+}/p_{\gamma}^{+}$ is the boost-invariant fraction of the photon momentum carried by a heavy quark (or antiquark). The dipole cross section $\sigma_{Q\bar{Q}}(r,s)$ in Eq. (3) describes the interaction of the $Q\bar{Q}$ dipole with the nucleon target. The Green function $G_{Q\bar{Q}}(z^{\prime}\to\infty,\vec{r}_{2};z,\vec{r}_{1})$, which describes the evolution of a $Q\bar{Q}$ pair between points $\vec{r}_{1},z$ and $\vec{r_{2}},z^{\prime}$, satisfies the two-dimensional Schrödinger equation [6]

$$i\frac{d}{dz_{2}}\,G_{Q\bar{Q}}(z_{2},\vec{r}_{2};z_{1},\vec{r}_{1})=\left[\frac{\eta^{2}-\Delta_{r_{2}}}{2\,k\,\alpha\,(1-\alpha)}+V_{Q\bar{Q}}(z_{2},\vec{r}_{2},\alpha)\right]G_{Q\bar{Q}}(z_{2},\vec{r}_{2};z_{1},\vec{r}_{1})\,,$$ (4)

where $\eta^{2}=m_{Q}^{2}+\alpha(1-\alpha)Q^{2}$, $m_{Q}$ is the heavy quark mass and Laplacian $\Delta_{r_{2}}$ acts on the coordinate $r_{2}$. Important ingredient of the Schrödinger equation is the complex $Q\bar{Q}$ potential in the LF frame, $V_{Q\bar{Q}}(z_{2},\vec{r}_{2},\alpha)$. The imaginary part of the LF potential controls the attenuation of the $Q\bar{Q}$ pair in the medium and the corresponding real part describes the interaction between the $Q$ and $\bar{Q}$. The values of the real part of the LF potential for arbitrary model of $Q-\bar{Q}$ interaction are obtained numerically, as described in Ref. [3], by solving the LF Schrödinger equation for $\Psi_{V}(\vec{r}_{2},\alpha)$.

The kinetic term in Eq. (4) includes the coherence length (CL), $l_{c}=1/q_{L}$, with $q_{L}$ being the longitudinal momentum transfer defined as

$$q_{L}(z)=\frac{M_{Q\bar{Q}}^{2}(z)+Q^{2}}{2\,k}\Rightarrow\frac{\eta^{2}-\Delta_{r_{2}}}{2\,k\,\alpha(1-\alpha)}\,,$$ (5)

where $M^{2}_{Q\bar{Q}}=(m_{Q}^{2}+q_{T}^{2})/\alpha(1-\alpha)$ is the $Q\bar{Q}$ effective mass and $q_{T}^{2}\Rightarrow-\Delta_{r}$ is the quark transverse momentum. The condition $l_{c}\geq R_{A}$ controls the onset of the quantum coherence effect related to $Q\bar{Q}$ Fock component of the photon and leading to higher twist quark shadowing. Another effect, known as the color transparency (CT), is responsible for the final state absorption of produced quarkonium [2, 3]. The onset of CT is controlled by the transverse size evolution of a $Q\bar{Q}$ pair propagating through the medium and is also incorporated in the Green function formalism.

In the electroproduction of heavy quarkonia at EIC energies, we evaluate the total cross section for quarkonium production on nuclear target as follows,

$$\sigma_{\gamma^{*}A\to J\!/\!\psi A}(s,b)=\int d^{2}b~{}~{}\frac{d^{2}\sigma_{A}^{coh}(s,b)}{d^{2}b}\,.$$ (6)

Then, we express the nucleus-to-nucleon ratio (nuclear transparency) as

$$R_{A}^{coh}(J\!/\!\psi)=\frac{\sigma_{\gamma^{*}A\to J\!/\!\psi A}}{A~{}\sigma_{\gamma^{*}N\to J\!/\!\psi N}}$$ (7)

where $\sigma_{\gamma^{*}N\to VN}$ is the total cross section for $J/\Psi$ production on a nucleon target. More details about the formalism can be found in Ref. [3].

We also include in Eq. (3) a small correction due to the real part of the $\gamma^{*}N\to J\!/\!\psi N$ amplitude by performing the following replacement [7, 8, 9],

$$\sigma_{Q\bar{Q}}(r,s)\Rightarrow\sigma_{Q\bar{Q}}(r,s)\,\left(1-i\,\frac{\pi}{2}\,\Lambda\right),\qquad\qquad\Lambda=\frac{\partial\,\ln\,{\sigma_{Q\bar{Q}}(r,s)}}{\partial\,\ln s}\,.$$ (8)

Another phenomenon that affects heavy quarkonium production in photo-nuclear reactions in UPC is gluon shadowing. The effect of GS has to be additionally included as a shadowing correction corresponding to higher Fock components of the photon containing gluons. It was incorporated as a reduction of $\sigma_{Q\bar{Q}}(r,s)$ in nuclear reactions with respect to processes on the nucleon [14],

$$\sigma_{Q\bar{Q}}(r,x)\Rightarrow\sigma_{Q\bar{Q}}(r,x)\cdot R_{G}(x,b)\,.$$ (9)

The correction factor $R_{G}(x,b)$ was calculated within the Green function formalism as a function of the nuclear impact parameter $b$ and the Bjorken variable $x$.

3 Results

We calculated the coherent charmonium photoproduction in UPC at LHC energies according to Eq. (1), as well as the electroproduction at EIC energies using Eqs. (6) and (7).

Figure 1: Rapidity distributions of coherent cross section for photoproduction of $J/\psi$ (left) and $\psi^{\prime}$ (right) in UPC at LHC energies $\sqrt{s_{N}}=5.02\,\,\mbox{TeV}$ (top panel) and 5.5 TeV (bottom panel). The data for $d\sigma/dy$ from ALICE [10], LHCb [11, 12] and CMS [13] Collaborations are shown for comparison. Figure is adopted from Ref. [3].

In Figure 1, top panels, we present comparison of our model predictions for the rapidity distributions $d\sigma/dy$ of coherent $J\!/\!\psi(1S)$ (left panel) and $\psi^{\,\prime}(2S)$ (right panel) photoproduction in UPC at c.m. collision energy $\sqrt{s_{N}}=5.02\,\,\mbox{TeV}$ with available data. Our results were obtained for charmonium wave functions generated by two $c$-$\bar{c}$ potential models, POW [15, 16] (dotted and dot-dashed lines) and BT [17] (solid and dashed lines). For the dipole cross sections $\sigma_{Q\bar{Q}}$ we employed the GBW [18] (solid and dot-dashed lines) and KST [2] (dashed and dotted lines) parametrizations. Our predictions based on the Green function formalism manifest a rather good agreement with data. Bottom panels of Figure 1 show the prediction for rapidity distributions of coherent charmonium cross section at larger $\sqrt{s_{N}}=5.5\,\,\mbox{TeV}$, which is planned to be measured at the LHC. Here, solid lines represent calculations within a rigorous path integral technique. The dotted lines are based on the Eikonal approximation ($l_{c}\gg R_{A}$ in Eq. (3)) for production cross section and without GS corrections. The difference between solid and dotted lines at midrapidities illustrates the net GS effect. At forward/backward rapidities, the GS vanishes and the difference between the two lines is caused by the reduced effects of quantum coherence when $l_{c}\mathrel{\hbox to0.0pt{\lower 4.0pt\hbox{\hskip 1.0pt$\sim$}\hss}\raise 1.0pt\hbox{$<$}}R_{A}$.

Figure 2: Ratio $R_{A}^{coh}$ for $J/\psi$ coherent production on the gold target as a function of c.m. energy $W$ at several fixed values of photon virtuality $Q^{2}=0$, 5, 20 and 50 GeV² (left panel). Right panel shows the ratio $R_{A}^{coh}$ as a function of $Q^{2}+M_{J/\Psi}^{2}$ at several fixed values of $W=20$, 30, 50, 100 and 200  GeV. Figure is adopted from Ref. [3].

In Figure 2, we present the ratio $R_{A}^{\rm coh}$ from Eq. (7) as a function of the c.m. energy $W$ (left panel) and as a function of $Q^{2}+M_{J/\Psi}^{2}$ (right panel) for electroproduction of $J\!/\!\psi(1S)$ on a gold target. The values of $R_{A}^{coh}$ were obtained using the charmonium LF wave function generated by the POW $c$-$\bar{c}$ interaction potential and the GBW model for $\sigma_{Q\bar{Q}}$ fitted to data for photon virtualities $Q^{2}\leq 50$ GeV² [18]. In the left panel of Figure 2, our results show significant leading twist corrections, rising with the photon energy, as differences between solid (the Green function formalism with GS) and dashed (Eikonal approximation without GS) lines at large $W$. At small photon energies, up to the position of maximal $R_{A}^{\rm coh}$ values, the difference is caused by the effect of reduced coherence length, as $l_{c}\mathrel{\hbox to0.0pt{\lower 4.0pt\hbox{\hskip 1.0pt$\sim$}\hss}\raise 1.0pt\hbox{$<$}}R_{A}$. Right panel of Figure 2 represents the $Q^{2}+M_{J/\Psi}^{2}$-behavior of $R_{A}^{\rm coh}$ at different values of photon energy $W$. At large $W=100$ and $200$ GeV, CL is long enough, $l_{c}>R_{A}$, to neglect its variation with $Q^{2}$ and the rise of $R_{A}^{\rm coh}$ with $Q^{2}$ is a net manifestation of CT. The effect of reduced CL is visible at lower energies, $W<50$ GeV, as decrease of $R_{A}^{\rm coh}$ with $Q^{2}$.

Our predictions for the onset of reduced coherence effects and gluon shadowing can be tested by future experiments at EIC. The proper treatment of the shadowing and absorption effects in the nuclear medium within a rigorous path integral technique may be crucial for the future conclusive evidence of expected gluon saturation effects at large energies.