Exclusive dilepton production at forward rapidities in $PbPb$ collisions

The dilepton production in diffractive and exclusive processes at forward rapidities considering ultraperipheral $PbPb$ collisions at the LHC is investigated. Predictions for the $e^+ e^-$, $\mu^+ \mu^-$ and $\tau^+ \tau^-$ cross sections are presented taking into account of realistic cuts that can be implemented by the LHCb Collaboration in a future experimental analysis. Our results indicate that the background associated with the diffractive production can be strongly suppressed and the exclusive processes can be cleanly separated. For the $\tau^+ \tau^-$ production, the semi and purely leptonic decay channels are considered. Our results indicate that a future experimental analysis of the dilepton production at the LHCb is feasible and can be useful to search for BSM physics.

The experimental study of photon -induced interactions in pp/pA/AA collisions became a reality in the last years, motivated by the possibility of improve our understanding of the Standard Model as well as by the opportunity of use these processes as an alternative tool to search for New Physics (for a recent review see e.g. Ref. [1]). Recent results for the exclusive diphoton production in ultraperipheral heavy ion collisions (UPHICs) [2] provided the first observation of the light -by -light (LbL) scattering [3,4], which was one of the most important predictions in the beginning of Quantum Electrodynamics (QED). Another important result was the observation of the exclusive dilepton production by photon -photon interactions in heavy ion collisions at the RHIC [5] and LHC [6][7][8], which has demonstrated the feasibility of study this final state and allowed us to improve the description of the nuclear photon flux and the modeling of the absorption effects, that suppress the strong interactions in these collisions [9,10]. Moreover, both results strongly motivated the studies related to beyond Standard Model (BSM) physics performed in Refs. [11][12][13][14][15], which focused to estimate the impact of the axionlike particle production in the diphoton final state and in the possibility of determine the τ lepton electromagnetic moments in the γγ → τ + τ − interactions present in P bP b collisions at the LHC. In particular, the results obtained in Refs. [14,15] have demonstrated the potentiality of the exclusive dilepton production to constrain BSM physics.
The studies performed in [9,10,14,15] were restricted to the rapidity range covered by the central detectors as e.g. ALICE, ATLAS and CMS. Such selection is justified, since the rapidity distribution is larger for midrapidities. However, for the LHC energies, the exclusive dilepton production for forward rapidities can become of the same order, especially if dileptons with small values of invariant masses and low transverse momenta are measured, as possible in the LHCb detector. One of the goals of this letter is to present, for the first time, predictions for the exclusive e − e + , µ − µ + and τ − τ + production at forward rapidities considering realistic selection cuts that can implemented by the LHCb Collaboration to select the associated events. As we will show below, a future experimental analysis is feasible. Our second goal is to estimate, for the first time, the background associated to dilepton production in double diffractive processes in P bP b, which also are characterized by two intact nuclei and two rapidity gaps in the final state, i.e., empty regions in pseudo-rapidity that separate the intact very forward nuclei from the l + l − state. We will demonstrate that such background can be strongly suppressed by the exclusivity cuts, which implies that future data will be a direct probe of the dilepton production by photon -photon interactions. As a by product, we also will show that the diffractive events can be separated using a different set of cuts (For a similar analysis for the diphoton production see Ref. [16]).
In what follows we will present our predictions for the exclusive dilepton production by photon -photon interactions in ultraperipheral heavy ion collisions, represented in Fig. 1 (a), which are characterized by an impact parameter b greater than the sum of the radius of the incident nuclei. In such collisions, the strong interactions between the incident nucleus are suppressed and the electromagnetic interaction between them become dominant [2,17]. The cross section for the P bP b → P b l + l − P b can be estimated using the equivalent photon approximation [18], which allow us to factorize the elementary γγ → l + l − process and describe the intense electromagnetic fields that accompany the relativistic heavy ions in terms of a spectrum of equivalent photons proportional to the square of the nuclear charge Z. As a consequence, the nuclear cross section is proportional to Z 4 , implying a very large number of events at RHIC and LHC. Moreover, it is assumed that the photons are coherently radiated by the whole nucleus. Such condition imposes that the minimum photon wavelength must be greater than the nuclear radius R. In other words, the photon virtuality must satisfy Q 2 = −q 2 ≤ 1/R 2 , with the photon four -momentum being q µ = (ω, q ⊥ , q z = ω/v), where q ⊥ arXiv:2012.08923v1 [hep-ph] 16 Dec 2020 is the transverse momentum of the photon in a given frame, where the projectile moves with velocity v. Therefore, we have that Q 2 = ω 2 /γ 2 L + q 2 ⊥ . The coherence condition limits the maximum energy of the photon to ω < ω max ≈ γ L /R and the perpendicular component of its momentum to q ⊥ ≤ 1/R. Therefore, the coherence condition sets an upper limit on the transverse momentum of the photon emitted by the hadron, which should satisfy q ⊥ ≤ 1/R, being ≈ 28 MeV for P b beams. One has that the photon virtuality can be neglected and the photons can be considered as being real. Moreover, the dilepton pair will be produced with a very small transverse momentum. In the impact parameter representation, the total cross section can be expressed as follows where √ s is center -of -mass energy of the P bP b collision, ⊗ characterizes a rapidity gap in the final state, W = √ 4ω 1 ω 2 is the invariant mass of the l + l − system and Y its rapidity. The photon energies ω 1 and ω 2 are related to W and to the rapidity Y of the outgoing dilepton system by ω 1 = W 2 e Y and ω 2 = W 2 e −Y . The cross sectionσ is the elementary cross section to produce a pair of leptons, where l = e, µ, τ . Moreover, N (ω i , r i ) is the equivalent photon spectrum with energy ω i at a transverse distance r i from the center of nucleus, defined in the plane transverse to the trajectory, which is determined by the charge form factor of the nucleus. Finally, the factor S 2 abs (b) depends on the impact parameter b of the P bP b collision and is denoted by absorptive factor, which excludes the overlap between the colliding nuclei and allows to take into account only ultraperipheral collisions. As in Refs. [10,16], the nuclear photon flux will be estimated assuming a point -like form factor and that the absorptive factor S 2 abs (b) will be described by the model proposed by Baur and Ferreira -Filho [19] to exclude the strong interactions between the incident nuclei. One has that the predictions for small invariant masses, which is the kinematical range of interest in this letter, are not dependent on these choices, as demonstrated in Ref. [10]. In our calculations we express Eq.(1) in the momentum space in order to impose cuts that depend on the lepton transverse momenta. Moreover, our analysis will be performed using the Forward Physics Monte Carlo (FPMC) [20], which allow us to estimate the associated cross sections and distributions taking into account of the detector acceptances.
One possible background for the exclusive dilepton production in UPHICs is the Drell -Yan (DY) production in double diffractive processes, represented in Fig. 1 (b). As shown in Ref. [21] for pp collisions, such process can dominate in some regions of the phase space. However, its contribution for heavy ion collisions is still an open question. In the diffractive DY process, the dilepton is produced through the qq → l + l − subprocess, with the partons being constituents of the Pomerons emitted by the incident nuclei, which remain intact after the emission. As the Pomeron is a color -singlet object, two rapidity gaps are expected in the final state, as in the case of the exclusive dilepton production. The main difference between these processes is the presence of the Pomeron remnants X and X in the double diffractive case, which are expected to generate additional tracks in the final state, in contrast with the exclusive case, where only the dilepton state is present. Moreover, from the analysis of the diphoton production performed in Ref. [16], we can expect that the transverse momentum distribution of the dileptons in these two channels will be different, with the DY process generating dileptons with larger transverse momentum. Such aspects will be considered in what follows in order to suppress the contribution of the diffractive process in comparison to the exclusive one. As in Ref. [16], we will estimate the cross section for the diffractive DY production using the Resolved Pomeron model [22], which implies that the cross section can be expressed by where S 2 eik is the soft survival factor, which takes into account of the soft interactions which lead to an extra production of particles that destroy the rapidity gaps in the final state [23], and q D i (x i , µ 2 ) andq D i (x i , µ 2 ) are the diffractive quark and antiquark densities of the nucleus i with a momentum fraction x i . Such quantities will be modeled using the approach proposed in [24] (See also Ref. [25]), which have described the diffractive distributions for a nucleus taking into account the nuclear effects associated to the nuclear coherence and the leading twist nuclear shadowing (For details see Ref. [16]). In our analysis we will estimate the associated cross section and distributions using the FPMC generator, where the hard matrix elements are treated by interfacing FPMC with HERWIG v6.5. Moreover, we will assume that S 2 eik = 3.4 × 10 −5 for P bP b collisions at √ s = 5.5 TeV, which was calculated in Ref. [26], using an approach that generalizes the model described in Ref. [27] for coherent double exchange processes in nuclear collisions.
In what follows we will present our results for the exclusive and diffractive dilepton production in P bP b collisions at √ s = 5.5 TeV. In our analysis we will use the FPMC event generator [20] to estimate the processes represented in the Figs. 1 (a) and (b). For the final state selection we use HepMC2 [30] with a plugin called HepPDT[31] which has been designed to be used by any Monte Carlo particle generator or decay package. HepPDT has the function to store particle information such as charge and nominal mass in a table which is accessed by a particle ID number. The particle ID number is defined according to the Particle Data Group's Monte Carlo numbering scheme [28]. The final analysis and distributions are done with ROOT [32]. We will focus on dileptons with invariant mass in the range 1.5 GeV ≤ m ll ≤ 20 GeV, removing the dileptons with invariant mass in the ranges 2.796 GeV ≤ m ll ≤ 3.196 GeV, 3.586 GeV ≤ m ll ≤ 3.786 GeV, and 9.0 GeV ≤ m ll ≤ 10.6 GeV, associated to the J/Ψ, Ψ(2S) and Υ decays, respectively.
Initially, let's consider the dimuon production and present in Fig. 2 our results for the invariant mass (m µµ ), rapidity (y µµ ) and acoplanarity (1 − (∆φ/π)) distributions considering the exclusive and diffractive channels, which were obtained at the generation level, without the inclusion of any selection in the events. We have that the diffractive contribution is strongly suppressed by the soft survival effects, with the associated invariant mass and rapidity distributions being a factor ≈ 10 5 smaller than the photon -induced one in the m µµ and y µµ ranges considered. On the other hand, the results for the acoplanarity distribution indicate that in the exclusive case the muons are dominantly produced in a configuration back -to -back. Such result is expected due to the very small transverse momentum carried by the photons in the initial state. In contrast, for the diffractive case, the transverse momentum of the quarks inside the Pomeron that interact to generate the dimuon can be large, which implies a non -negligible contribution associated to muon pairs with a large acoplanarity. Similar results are obtained for the e + e − and τ + τ − production.
In order to obtain realistic estimates of the dilepton production in P bP b collisions, which can be compared with the future experimental data, we will include in our analysis the experimental cuts that are expected to be feasible in the next run of the LHC by the LHCb detector. Due to the distinct characteristics of the τ decay, we will study this final state using a different strategy to be discussed below. For the e + e − and µ + µ − production, the selection criteria implemented in our analysis of double diffractive and exclusive diphoton processes at forward rapidities are the following: • The initial selection is made on the transverse momentum of the leptons and invariant mass: m ll < 20 GeV and p T (l, l) > 0.4 GeV. The notation p T (l, l) represents that the transverse momentum cut is applied in both leptons separately; • We will impose a cut on the acoplanarity (1 − (∆φ/π) < 0.03) and transverse momentum of the dilepton system (p T (ll) < 2 GeV); • We will select only events where leptons are produced in the rapidity range 2.0 < η(l, l) < 4.5 and no other particle with p T > 0.5 GeV in the range 8.0 < |η| < 5.5 is produced (this corresponds to the HERSCHEL acceptance installed in the LHCb experiment [33]).
The associated predictions are presented in the Table I. One has that the predictions for the dielectron and dimuon production are very similar, which is expected in the invariant mass range considered, and the cut on the transverse momentum of the pair suppress the diffractive contribution. In addition, the HERSCHEL selection also suppresses this contribution and does not impact the exclusive one. Such results are expected, since in the exclusive case there are no additional tracks beyond those associated with the dilepton, in contrast to the diffractive production, where the remnants of the Pomeron are expected to generate new tracks. One important comment is that the cut on rapidity for the kinematical range covered by the LHCb detector reduces the cross section by approximately one order of magnitude. However, the final values are still very large, implying that the experimental analysis is feasible.
The results presented in Fig. 2 indicated that the acoplanarity distribution for exclusive and diffractive processes are distinct. Therefore we can invert the acoplanarity cut used in our previous analysis in order to select the diffractive events instead of the exclusive one. The results for the dimuon production derived using this new set of cuts are presented in the Table II. Similar results are obtained for the dielectron production. One has that the exclusive contribution is strongly suppressed and the measurement of dileptons for these selection criteria is a direct probe of the diffractive channel. Such result opens the possibility to constrain the magnitude of the soft survival factor S 2 eik in the dilepton production, which is one of the main sources of theoretical uncertainty in the treatment of diffractive processes at nuclear collisions. In what follows we will discuss the exclusive and diffractive τ + τ − production in UPHICs. The exclusive case was recently investigated in Refs. [14,15], strongly motivated by the possibility of use this final state to constrain the anomalous magnetic momentum of the τ lepton and search for CP violation effects. These studies focused on the kinematical range covered by the ATLAS and CMS detectors. We will complement these analyzes, by analyzing the possibility of measure this final state in the LHCb detector, which collected data with low pile-up and it is able to reconstruct low momentum tracks. We will restrict our study to the Standard Model subprocesses and postpone for a future publication the investigation of the impact of new physics. Differently from the lighter charged leptons, the very short lifetime of the τ lepton [O(10 −13 s)] implies that it decays into lighter leptons and hadrons before to interact with the detector. As a consequence, the separation of the γγ → τ + τ − events depends on a precise identification of the decay products. The τ decays can be classified as one prong decays (≈ 80%), where one charged particle is present in the final state, and three prong decays (≈ 20%), which are characterized by the production of three charged particles. In our analysis, only the 1-prong decays (semi and purely leptonic) production will be considered. Moreover, in order to suppress the contribution of the diffractive DY process for the purely leptonic channel, we will select events where the leptons in the final state have different flavours. Therefore, the following events will be included in our study: τ τ → (e ± ν e ν τ ) (µ ∓ ν µ ν τ ), (e ± ν e ν τ ) (π ∓ π 0 ν τ ), (µ ± ν µ ν τ ) (π ∓ π 0 ν τ ) . ( The HERWIG 6.5 program [29], which incorporates the TAUOLA MC code [34], is used to model τ decay for our studies. We select events with p T (l, chg) > 1.0 GeV and 2.0 < η l,chg < 4.5 with 0 extra tracks and p T (l + chg) < 2.0 GeV. As considered for lighter leptons, we discard events where particles with p T > 0.5 GeV are in the range 8.0 < |η| < 5.5, corresponding to the HERSCHEL acceptance in the LHCb experiment [33]. Before presenting our predictions, it is important to emphasize that we also have performed the analysis considering the rapidity range covered for central detectors and obtained a very good agreement with the results presented in Ref. [15].
In Fig. 3 we present our results for the correlation between the hard process and reconstruction of the τ + τ − pair considering the 1-prong category at generation level. The reconstruction of rapidity of the pair is done before the exclusivity cuts and establish a benchmark of the effective combinations resulting in the τ + τ − pair products. Our results indicate that the rapidity of the system before decay and the final state charged particle in exclusive and diffractive processes are correlated, with this correlation being larger in the exclusive case. In Table III we present our predictions for the γγ → τ + τ − reaction in P bP b collisions at √ s = 5.5 TeV considering the different decay channels. As already observed for lighter leptons, the exclusive channel dominates the τ + τ − production. Moreover, the semileptonic contribution for the exclusive production is a factor ≈ 2.5 larger than the purely leptonic one. Our results indicate that the total cross section for the exclusive τ + τ − production in the kinematical covered by the LHCb detector will be ≈ 60 nb, which implies that the number of events expected for the next run of the LHC is ≈ 600 per year. Although smaller than the prediction for central detectors, this result indicates that a future experimental analysis is feasible.
Diffractive Exclusive Decay e ± µ ∓ µ + π − µ − π + e + π − e − π + e ± µ ∓ µ + π − µ − π + e + π − e − π + Generation 2958.0 6.7 × 10 8 Channel 152.1 11.1 7.8 11.8 8.6 3.6 × 10 7 5.2 × 10 7 5.2 × 10 7 5.3 × 10 7 5.2 × 10 7 pT (l, chg) > 1.0 GeV 18.7 0.07 0.07 0.03 0.02 5.4 × 10 6 6.2 × 10 6 6.2 × 10 6 6.3 × 10 6 6.6 × 10 6 2.0 ≤ η chg,l ≤ 4. Let's summarize our main results and conclusions. In this letter, we have investigated, for the first time, the dilepton production in diffractive and exclusive production at forward rapidities considering ultraperipheral P bP b collisions at the LHC. We have derived predictions for the e + e − , µ + µ − and τ + τ − cross sections taking into account of realistic cuts that can be implemented by the LHCb experiment in a future data analysis. Our results indicate that the background associated to the diffractive production can be strongly suppressed and the exclusive processes can be cleanly observed. We also demonstrated that selecting events with large acoplanarity, the diffractive events can be studied, allowing to constrain the modeling of the soft survival factor. For the τ + τ − production, the semi and purely leptonic decay channels were considered and the final results indicate that the study of the exclusive production can be performed. Finally, the predictions presented in this indicate that a future experimental analysis of the dilepton production at the LHCb is feasible and can be useful to search for BSM physics.