On-line version ISSN 1678-4448
Braz. J. Phys. vol.38 no.3b São Paulo Sept. 2008
An analysis on single and central diffractive heavy flavour production at hadron colliders
M. V. T. Machado
Universidade Federal do Pampa, Centro de Ciencias Exatas e Tecnológicas Campus de Bage, Rua Carlos Barbosa, CEP 96400-970, Bage, RS, Brazil
In this contribution results from a phenomenological analysis for the diffractive hadroproduction of heavy flavors at high energies are reported. Diffractive production of charm, bottom and top are calculated using the Regge factorization, taking into account recent experimental determination of the diffractive parton density functions in Pomeron by the H1 Collaboration at DESY-HERA. In addition, multiple-Pomeron corrections are considered through the rapidity gap survival probability factor. We give numerical predictions for single diffractive as well as double Pomeron exchange (DPE) cross sections, which agree with the available data for diffractive production of charm and beauty. We make estimates which could be compared to future measurements at the LHC.
Keywords: Heavy flavour production; Pomeron physics; Single diffraction; Quantum Chromodynamics
For a long time, diffractive processes in hadron collisions have been described by Regge theory in terms of the exchange of a Pomeron with vacuum quantum numbers . However, the nature of the Pomeron and its reaction mechanisms are not completely known. Currently, a promising way to clarify these questions is using the hard scattering to resolve the quark and gluon content in the Pomeron , regarding that a parton structure is natural in a modern QCD approach to the strongly interacting Pomeron. Form the experimental point of view, observations of diffractive deep inelastic scattering (DDIS) at HERA have increased the knowledge about the QCD Pomeron, where its diffractive distributions of singlet quarks and gluons have been phenomenologically determined .
In hadronic collisions, a single diffractive event is characterized by one of the colliding hadrons emitting a Pomeron that scatters off the other hadron. Hard diffractive events with a large momentum transfer are also set by the absence of hadronic energy in certain angular regions of the final state phase space (large rapidity gaps). For events in which both colliding hadrons remain intact as they each emit a Pomeron we have the so-called central diffractive events. In the latter case, also known as double Pomeron exchange (DPE) processes, both incoming hadrons are quasi-elastically scattered and the final states system in the center region is produced by Pomeron-Pomeron interaction.
Here, we concentrate on the single diffractive processes P + p() → P + J/Ψ [ϒ] + X and the central diffractive reactions, p + p() → p + J/Ψ [ϒ] + p(). The diffractive heavy quarkonium production has drawn attention because their large masses provide a natural scale to guarantee the application of perturbative QCD. There are several mechanisms proposed for the quarkonium production in hadron colliders [4, 5], as the color singlet model, the color octet model and the color evaporation model. An important feature of these perturbative QCD models is that the cross section for quarkonium production is expressed in terms of the product of two gluon densities at large energies. This feature is transferred to the diffractive quarkonium production, which is now sensitive to the gluon content of the Pomeron at small-x and may be particularly useful in studying the different mechanisms for quarkonium production.
In order to do so, we will use the hard diffractive factorization, where the diffractive cross section is the convolution of diffractive parton distribution functions and the corresponding diffractive coefficient functions. At high energies there are important contributions from unitarization effects to the single-Pomeron exchange cross section. These absorptive corrections cause the suppression of any large rapidity gap process, except elastic scattering. In the black disk limit the absorptive corrections may completely terminate those processes. This partially occurs in (anti)proton-proton collisions, where unitarity is nearly saturated at small impact parameters . These multiple-Pomeron contributions depends on the particular hard process and it is called survival probability factor, which are important for the reliability of theory predictions.
The present contribution is organized as follows. In next section, we present the main formulas to compute the inclusive and diffractive cross sections for heavy quarkonium hadroproduction. We also present the parameterization for the diffractive partons distribution in the Pomeron, extracted recently in DESY-HERA, and theoretical estimations for the gap survival probability factor. In the last section we present the numerical results for Tevatron and perform predictions to future measurements at the LHC experiment. The compatibility with data is analyzed and the comparison with other approaches is considered. As an anticipation of the main results, at the Tevatron the single and central diffractive J/Ψ and ϒ production are observable with a single diffractive ratio RSD (Tevatron) that is between 1 % (charmonium) and 0.5 % (bottomonium), with lower values at the LHC. Also, it is predicted that = 1.12± 0.19 %. The central diffractive cross sections are still feasible to be measured, despite the very small diffractive ratios. Therefore, we stress that the theoretical model dependence is either large and more detailed studies are timely.
2. DIFFRACTIVE HADROPRODUCTION OF HEAVY QUARKONIUM
For our purpose we will use the Color Evaporation Model (CEM) . The main reasons for this choice are its simplicity and fast phenomenological implementation, which are the base for its relative success in describing high energy data[4, 5]. In such an approach, the cross section for a process in which partons of two hadrons, h1 and h2, interact to produce a heavy quarkonium state, h1 + h2 H(nJCP) + X, is given by the cross section of open heavy-quark pair production that is summed over all spin and color states. All information on the non-perturbative transition of the pair to the heavy quarkonium H of quantum numbers JPC is contained in the factor FnJPC that a priori depends on all quantum numbers ,
where σ() is the total hidden cross section of open heavy-quark production calculated by integrating over the pair mass from 2mQ to 2mO, with mO is the mass of the associated open meson. The hidden cross section can be obtained from the usual expression for the total cross section to NLO. These hadronic cross sections in pp collisions can be written as
where x1 and x2 are the fractional momenta carried by the colliding partons and are the proton parton densities. The partonic cross sections are known up to NLO accuracy . Here, we assume that the factorization scale, µF, and the renormalization scale, µR, are equal. We also take µ = 2mQ, using the quark masses mc = 1.2 GeV and mb = 4.75 GeV, which provide an adequate description of open heavy-flavour production . The invariant mass is integrated over in the charmonium case and for ϒ production. The factors FnJPC are experimentally determined  to be F11 ≈ 2.5 x 10-2 for J/Ψ and F11 ≈ 4.6 x 10-2 for ϒ. These coefficients are obtained with NLO cross sections for heavy quark production .
The agreement with the total cross section data is fairly good. The low x region is particularly relevant for J/Ψ production at the LHC as well as at Tevatron. For charmonium production, the gg process becomes dominant and information on the gluon distribution is of particular importance. In Fig. 1, the differential cross section is shown as a function of center of mass energy, . The Tevatron data indicates that correction for gluon saturation may be important in bringing the theoretical curve closer to the experimental result. In particular, the typical values of Bjorken-x would be for Tevatron and x ≈ 10-4 at the LHC taken at a relatively low momentum scale . In Fig. 2, the differential cross section is presented as a function of the center of mass energy (solid line). In addition, the measured cross sections for the sum of the three Y states (ϒ. + ϒ' + ϒ") in the dilepton decay channel are shown . The agreement of the CEM model with accelerator data is very good and an extrapolation to the LHC is presented. The results give us considerable confidence in the extrapolation to the LHC energy.
2.1. Diffractive cross section - singlePomeron exchange
For the hard diffractive processes we will consider the Ingelman-Schlein (IS) picture , where the Pomeron structure (quark and gluon content) is probed. In the case of single diffraction, a Pomeron is emitted by one of the colliding hadrons. That hadron is detected, at least in principle, in the final state and the remaining hadron scatters off the emitted Pomeron. The diffractive cross section of a hadron-hadron collision is assumed to factorise into the total Pomeron-hadron cross section and the Pomeron flux factor . The single diffractive event, h1 + h2 → h1,2 + H [nJCP] + X, may then be written as
where the Pomeron kinematical variable xIP is defined as , where is the center-of-mass energy in the Pomeron-hadron j system and the center-of-mass energy in the hadron i-hadron j system. The momentum transfer in the hadron i vertex is denoted by A similar factorization can also be applied to central diffraction, where both colliding hadrons can in principle be detected in the final state. The central quarkonium production, h1 + h2 → h1 + H [nJCP] + h2, is characterized by two quasi-elastic hadrons with rapidity gaps between them and the central heavy quarkonium products. The central diffractive cross section may then be written as,
Here, we assume that one of the hadrons, say hadron hi, emits a Pomeron whose partons interact with partons of the hadron h2. Thus the parton distribution xi fi/h1 (x1, /µ2) is replaced by the convolution between a distribution of partons in the Pomeron, βfa/IP(β,µ2), and the "emission rate" of Pomerons by the hadron, fIP/h(xIP,t). The last quantity, fIP/h(xIP,t), is the Pomeron flux factor and its explicit formulation is described in terms of Regge theory. Therefore, we can rewrite the parton distribution as
Using the substitution given in Eq. (4), the hidden heavy flavour cross section can be obtained from Pomeron-hadron cross sections for single and central diffraction processes,
In the numerical calculations, we will consider the diffractive pdf's recently obtained by the H1 Collaboration at DESY-HERA . The Pomeron structure function has been modeled in terms of a light flavour singlet distribution Σ(z), consisting of u, d and s quarks and anti-quarks and a gluon distribution g(z). The Pomeron carries vacuum quantum numbers, thus it is assumed that the Pomeron quark and antiquark distributions are equal and flavour independent:, where ΣIP is a Pomeron singlet quark distribution and Nf is the number of active flavours. Moreover, for the Pomeron flux factor, introduced in Eq. (4), we take the experimental analysis of the diffractive structure function , where the xIP dependence is parameterized using a flux factor motivated by Regge theory ,
where the Pomeron trajectory is assumed to be linear, .
2.2. Multiple-Pomeron exchange corrections
Here, we consider the suppression of the hard diffractive cross section by multiple-Pomeron scattering effects. This is taken into account through a gap survival probability factor, < |S|2 >, which can be described in terms of screening or absorptive corrections . This suppression factor of a hard process accompanied by a rapidity gap depends not only on the probability of the initial state survive, but is sensitive to the spatial distribution of partons inside the incoming hadrons, and thus on the dynamics of the whole diffractive part of the scattering matrix. The survival factor of a large rapidity gap (LRG) in a hadronic final state is the probability of a given LRG not be filled by debris, which originate from the soft re-scattering of the spectator partons and/or from the gluon radiation emitted by partons taking part in the hard interaction. Let A(s, b) be the amplitude of the particular diffractive process of interest, considered in the impact parameter, b, space. Therefore, the probability that there is no extra inelastic interaction is
where Ω is the opacity (or optical density) of the interaction.
We will consider the theoretical estimates for < |S|2 > from Ref.  (labeled KMR), which considers a two-channel eikonal model and rescattering effects. The survival probability factor is computed for single, central and double diffractive processes at several energies, assuming that the spatial distribution in impact parameter space is driven by the slope B of the pomeron-proton vertex. We will consider the results for single diffractive processes with 2B = 5.5 GeV~2 (slope of the electromagnetic proton form factor) and without N * excitation, which is relevant to a forward proton spectrometer (FPS) measurement. Thus, we have and for = 1.8 TeV (Tevatron) = 14 TeV (LHC)]. There are similar theoretical estimates, as the GLM approach , which also considers a multi-channel eikonal approach.
3. RESULTS AND SUMMARY
For the numerical calculations, the new H1 parameterization for the diffractive pdf's  has been used. The 'H1 2006 DPDF Fit A' is considered, with the cut xIP < 0.1. The singlePomeron results are presented in Figs. 1 and 2 for J/Ψ and ϒ production (dot-dashed curves), respectively. The single diffractive contribution is large, being of order 5-6 % from the inclusive cross section. The results corrected by multiple-Pomeron suppression factor are a factor about 1/10 lower than the single-Pomeron ones. It should be stressed that sizable uncertainties are introduced by changing, for instance, quark masses and/or the renormalization scale. However, our purpose here is to estimate the diffractive ratios σD/σtot, which are less sensitive to a particular choice. For sake of illustration, we have
nb for the energy of = 2 TeV and 159 (1770) nb for LHC energy, = 14 TeV. For ϒ we have pb and 35 (386) pb for Tevatron and LHC, respectively. Numbers between parentheses correspond to values not corrected by survival probability gap factor.
Concerning the central diffractive cross-sections, the predictions give small values. For instance, we have = 18 nb and B = 0.8 pb at = 2 TeV. This gives little room to observe central diffractive T events at the Tevatron but could be promising for the LH . The values reach = 45 nb and B = 3 pb at = 14 TeV. Once again, the results corrected by multiple-Pomeron suppression factor are a factor about 1/10 lower than the usual IS model.
Let us now compute the diffractive ratios. We have defined the single diffractive ratio as and the central diffractive ratio as . The results are summarized in Table I, where the diffractive ratios for heavy quarkonium production are presented for Tevatron and LHC energies. The multiple-Pomeron correction factors are taken from KMR model. The numbers between parentheses represent the single-Pomeron calculation. Based on these results we verify that the J/Ψ and ϒ production in single diffractive process could be observable in Tevatron and LHC, with a diffractive ratio of order of 1 % or less. This is a similar ratio measured in W and Z production at the Tevatron . The predictions for central diffractive scattering are still not very promising, giving small ratios. However, the study of these events is worthwhile since their experimental signals are quite clear. Finally, our theoretical prediction is in good agreement with the experimental measurement of CDF , which found = L45 ± 0.25 %. As the main theoretical uncertainty in determining the diffractive ratio is the value for the survival factor, we consider the theoretical band 0.210. i5 for Tevatron energy. Therefore, our prediction gives = 1.12±0.19 %, which is consistent with the Tevatron determination.
Our calculation can be compared to available literature in diffractive heavy quarkonium production. For instance, in Ref.  the large pT J/Ψ production in hard diffractive process is computed using the color octet fragmentation mechanism and the normalized Pomeron flux . They found a SD J/Ψ cross section at large pT (> 8 GeV) of order 10 pb and a diffractive ratio RSD = 0.65 ± 0.15%. This is closer to our calculation for the SD ratio for Tevatron within the theoretical errors. Afterwards, in Ref.  the color-octet mechanism combined with the two gluon exchange model(in LO approximation in QCD) for the diffractive J/Ψ production is considered. Now, the SD cross section is a(p J/ΨX) = 66 nb. The comparison between these calculations shows the size of the large theoretical uncertainty. A related calculation, the J/Ψ + y diffractive production, appeared in Ref.  based on IS model (with normalized flux ) and factorization formalism of NRQCD for quarkonia production. They found σ(p→[J/Ψ + γ]X) = 3.0 pb (8.5 pb) and diffactive ratio RSD = 0.5 (0.15) % in central region at the Tevatron (LHC). These results are somewhat still compatible with present calculation.
Concerning central diffraction, there are some theoretical studies in literature. In Ref. , the DPE process, p + p() → p + Xj + p(), is calculated using two-gluon exchange model in perturbative QCD. It is found the following DPE cross sections: σ(Xco) = 735 nb and σ(xb0) = 0.88 nb. In the same work, it is found that
= 2 nb/GeV and = 0.5 pb/GeV. Recently, in Ref.  the double-diffractive production of xc and xb mesons has been studied using also the Regge formalism and pQCD (including unitarity corrections). They found dσ/dy\y=0 = 130 (340) nb for xc production and dσ/dy\y=0 = 0.2 (0.6) nb for xb production for Tevatron (LHC). It was verified that the exclusive production is a factor 10 larger than the inclusive rate. Finally, the diffractive X meson production bas been computed  using the Bialas-Landshoff formalism and makes use of the DPEMC Monte-Carlo simulation for small-mass diffractive production. They found = 65(1.6) % and = 22 (1.83) % for Tevatron (LHC). In summary, the theoretical predictions for exclusive meson production are still quite distinct and more detailed studies are deserved.
In summary, we have presented predictions for diffractive heavy quarkonium production at the Tevatron and the LHC. We use Regge factorization, corrected by unitarity corrections modeled by a gap survival probability factor (correction for multiple-Pomeron exchange). The perturbative formalism for meson hadroproduction is based on the CEM model, which is quite successful in describing experimental results for inclusive production. For the Pomeron structure function, recent H1 diffractive parton density functions have been used. The results are directly dependent on the quark singlet and gluon content of the Pomeron. We estimate the multiple interaction corrections taking the theoretical prediction a multichannel model (KMR), where the gap factor decreases on energy. That is, < |S|2 >~ 15 % for Tevatron energies going down to < |S|2 >~ 9 % at LHC energy (for single diffractive process). We found that at the Tevatron single and central diffractive J/Ψ and ϒ production is observable with a single diffractive ratio RSD(Tevatron) that is between 1 % (charmonium) and 0.5 % (bottomonium), with lower values at the LHC. In particular, we predict = 1.12 ± 0.19%, which is in agreement with Tevatron measurements. The central diffractive cross sections for quarkonium production are still feasible to be measured, despite the very small diffractive ratios. In this case, the theoretical model dependence is still very large and detailed studies are deserved.
This work was supported by CNPq, Brazil. The author is grateful to Uri Maor and David dEnterria for useful discussions.
 P. D. B. Collins, An Introduction to Regge theory and high energy physics (Cambridge University Press, Cambridge, England, 1977). [ Links ]
 G. Ingelman and P. E. Schlein, Phys. Lett. B 152, 256 (1985). [ Links ]
 H1 Collab., A. Aktas etat., Eur. Phys. J. C 48, 715 (2006). [ Links ]
 J. P. Lansberg, Int. J. Mod. Phys. A 21, 3857 (2006). [ Links ]
 M. Krämer, Prog. Part. Nucl. Phys. 47, 141 (2001). [ Links ]
 J. F. Amundson, O. J. P. Eboli, E. M. Gregores, and F. Halzen, Phys. Lett. B 390, 323 (1997); [ Links ] C. Brenner Mariotto, M. B. Gay Ducati, and G. Ingelman, Eur. Phys. J. C 23, 527 (2002). [ Links ]
 M. L. Mangano, P. Nason, and G. Ridolfi, Nucl. Phys. B 373, 295 (1992); [ Links ] A. Petrelli, M. Cacciari, M. Greco, F. Maltoni, and M. L. Mangano, Nucl. Phys. B 514, 245 (1998); [ Links ] S. Frixione et at., Adv. Ser. Direct High Energy Phys. 15, 607 (1998). [ Links ]
 G. A. Schuler and R. Vogt, Phys. Lett. B 387, 181 (1996). [ Links ]
 A. B. Kaidalov, V. A. Khoze, A. D. Martin, and M. G. Ryskin, Eur. Phys. J. C 21,521 (2001). [ Links ]
 E. Gotsman, E. Levin, U. Maor, and A. Prygarin, arXiv:hep-ph/0511060. [ Links ]
 A. A. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 87, 241802 (2001). [ Links ]
 M. B. Gay Ducati, M. M. Machado, and M. V. T. Machado, Phys. Rev. D 75, 114013 (2007). [ Links ]
 F. Yuan and K. T. Chao, Phys. Rev. D 57, 5658 (1998). [ Links ]
 K. Goulianos, Phys. Lett. B 358, 379 (1995). [ Links ]
 F. Yuan, J. S. Xu, H. A. Peng, and K. T. Chao, Phys. Rev. D 58, 114016 (1998). [ Links ]
 J. S. Xu and H. A. Peng, Phys. Rev. D 59, 014028 (1999). [ Links ]
 F. Yuan, Phys. Lett. B 510, 155 (2001). [ Links ]
 V. A. Khoze, A. D. Martin, M. G. Ryskin, and W. J. Stirling, Eur. Phys. J. C 35,211 (2004). [ Links ]
 M. Rangel, C. Royon, G. Alves, J. Barreto, and R. Peschanski, Nucl. Phys. B 774, 53 (2007). [ Links ]
(Received on 20 March, 2008)