Correct Definition of Color Singlet P-Wave Non-Perturbative Matrix Element of Heavy Quarkonium Production

Recently we have proved factorization of infrared divergences in NRQCD S-wave heavy quarkonium production at high energy colliders at all orders in coupling constant. One of the problem which still exists in the higher order pQCD calculation of color singlet P-wave heavy quarkonium production/anihillation is the appearance of non-canceling infrared divergences due to real soft gluons exchange, although no such infrared divergences are present in the color singlet S-wave heavy quarkonium. In this paper we find that since the non-perturbative matrix element of the color singlet P-wave heavy quarkonium production contains derivative operators, the gauge links are necessary to make it gauge invariant and be consistent with the factorization of such non-canceling infrared divergences at all orders in coupling constant.


I. INTRODUCTION
divergence in this case, 1) infrared divergence due to the Coulomb gluon exchange and 2) infrared divergence due to real soft gluon emission/absorption (eikonal gluon) [5].
For S-wave and P-wave heavy quarkonium the appearance of infrared divergence due to Coulomb gluon exchange (which is also known as 1 v divergence [6][7][8][9] in the limit v → 0 where v is the relative velocity between heavy quark and antiquark) is usually handled by normalization of the bound state wave function similar to QED [10] where the Coulomb potential is used.
However, one may expect non-canceling infrared divergence due to emission/absorption of real soft gluon (eikonal gluon). In the higher order pQCD calculation of the color singlet S-wave heavy quark-antiquark production/anihillation one finds the absence of such infrared divergence [6,11]. Hence the non-perturbative effects can be factored into the nonperturbative wave function at the origin R(0) for the S-wave heavy quarkonium. In terms of heavy quark fields one finds that the definition of the gauge invariant non-relativistic wave function at the origin of the color singlet S-wave heavy quarkonium is given by [12] where χ is the two component Dirac spinor field that creates a heavy quark and ξ is the two component Dirac spinor field that annihilates a heavy quark and a † H is the creation operator of the heavy quarkonium H. In eq. (2) the non-perturbative matrix element < 0|χ † ξa † H a H ξ † χ|0 > is evaluated at the origin. However, unlike higher order pQCD calculation of color singlet S-wave heavy quarkantiquark production/anihillation, the higher order pQCD calculation of color singlet P-wave heavy quark-antiquark production/anihillation contains non-canceling infrared divergence due to real soft gluons emission/absorption [7][8][9]. Hence the non-perturbative effects can not be factored into the derivative of the non-relativistic wave function at the origin [12] |R ′ (0)| 2 = 2π 27 where and the non-perturbative matrix element < 0|χ †∇ ξa † H · a H ξ †∇ χ|0 > is evaluated at the origin. This is also easy to see from eq. (3) because, unlike |R(0)| 2 in eq. (2) which is gauge invariant, the |R ′ (0)| 2 in eq. (3) is not gauge invariant. Since the issue of gauge invariance and non-canceling infrared divergences due to real soft gluons emission/absorption are related [13,14], one expects that the gauge links need to be supplied in eq. (3) to make |R ′ (0)| 2 gauge invariant and consistent with the factorization of such non-canceling infrared divergences due to the real soft gluons emission/absorption. In this paper we will investigate this in detail.
We find that the correct definition of the gauge invariant non-perturbative matrix element of the color singlet P-wave heavy quarkonium production which is consistent with factorization of infrared divergences is given by is the gauge link in the fundamental representation of SU (3) and We find that the long-distance behavior of the color singlet P-wave non-perturbative matrix element < 0|χ † Φ∇Φ † ξa † H · a H ξ † Φ∇Φ † χ|0 > in eq. (5) is independent of the light-like vector l µ at all order in coupling constant in QCD where l µ is used to define the gauge link in eq. (6).
Note that there is no gauge link in |R(0)| 2 in eq. (2) in case of S-wave heavy quarkonium because there are no derivative operators like that in eq. (5) for P-wave heavy quarkonium.
The |R(0)| 2 in eq. (2) is gauge invariant without gauge links which is consistent with the fact that there are no uncanceled infrared divergences due to real soft gluons emission/absorption in case of color singlet S-wave heavy quarkonium.
In this paper we will provide a derivation of eq. (5).
The paper is organized as follows. In section II we briefly review the appearance of infrared divergences in higher order pQCD calculation in color singlet P-wave heavy quarkonium. In section III we discuss the lack of gauge invariance in S-wave NRQCD color octet non-perturbative matrix element which is used to address the issue of non-canceling infrared divergences in the color singlet P-wave heavy quarkonium production in the literature. In section IV we describe the infrared divergence due to real soft gluon (eikonal gluon) emission/absorption in pQCD. In section V we prove the factorization of infrared divergences in color singlet P-wave heavy quarkonium production at all orders in coupling constant at high energy colliders. In section VI we derive the correct definition of the gauge invariant non-perturbative matrix element of the color singlet P-wave heavy quarkonium production as given by eq. (5). Section VII contains conclusions.

II. INFRARED DIVERGENCE IN HIGHER ORDER PQCD CALCULATION IN COLOR SINGLET P-WAVE HEAVY QUARKONIUM
In the higher order pQCD calculation of heavy quark-antiquark annihilation in the color singlet P-wave state, the non-canceling infrared divergences were found for the color singlet P-wave heavy quarkonium decay to hadrons [7][8][9]. Similarly in the higher order pQCD calculation of hadroproduction of heavy quark-antiquark in the color singlet P-wave state, the non-canceling infrared divergences were found [5]. A satisfactory explanation of these non-canceling infrared divergences in the case of color singlet P-wave heavy quarkonium production/decay has not been given in the literature, partly because of the lack of understanding of the correct definition of the gauge invariant non-perturbative matrix element of the color singlet P-wave heavy quarkonium which is required to cancel these infrared divergences (factorization).
In addition to this it should be noted that the diagrammatic pQCD calculation is not always sufficient to predict the correct definition of non-perturbative matrix element in QCD, no matter how many orders of pQCD calculation is performed. On the other hand the nonperturbative QCD method (the path integral formulation of QCD) can correctly predict the correct definition of the non-perturbative matrix element in QCD. We will use path integral formulation of QCD in this paper.
In case of h c decays to hadrons the partonic process in which the non-canceling infrared divergence is found is given by [15] 1 P 1 → ggg.
Similarly in case of χ cJ decay to hadrons the partonic process in which the non-canceling infrared divergence is found is given by [7][8][9] 3 P J → qqg.
Note that the collinear divergences are canceled by the virtual corrections to the two gluon decay process [5]. Hence the non-canceling infrared divergence is due to the emission/absorption of real soft gluon (eikonal gluon).
Similar to the hadronic decay of color singlet P-wave heavy quarkonium, the non-canceling infrared divergences are found in the hadroproduction of color singlet P-wave heavy quarkonium. For example, the non-canceling infrared divergence appears in the partonic process [5] qq → 3 P J g.
Note that unlike decay process in eq. (9) there is no collinear divergence in the production process in eq. (10) and hence there is no virtual correction present either. The non-canceling infrared divergence is due to emission/absorption of real soft gluon (eikonal gluon).
It can be noted that in the original calculation [7][8][9] the binding energy was used as infrared cutoff which can be mapped into the corresponding infrared divergence analysis in the dimensional regularization scheme, see for example [16].
Hence whether it is color singlet P-wave heavy quarkonium production or decay, the non-canceling infrared divergences occur due to real soft gluons emission/absorption. In this paper we will derive the correct definition of the gauge invariant color singlet P-wave non-perturbative matrix element which cancels these infrared divergences at all order in coupling constant in QCD.
We emphasize here that we are not calculating the finite part of the cross section of the color singlet P-wave heavy quarkonium production in the partonic level processes like that in eq. (10). What we are investigating in this paper is to study the exact behavior of the infrared divergences due to real soft gluons emission/absorption without modifying the finite value of the cross section. Our main aim is to derive the correct definition of the gauge invariant color singlet P-wave non-perturbative matrix element which cancels these infrared divergences at all order in coupling constant.
The color singlet P-wave non-perturbative matrix element which was used in the earlier studies is given by eq. (3) which is not gauge invariant. Hence it is no surprise that the noncanceling infrared divergences were present in these earlier studies of color singlet P-wave heavy quarkonium production/decay. It is expected that if the gauge links are supplied in the definition of the non-perturbative matrix element to make it gauge invariant then it can accommodate these non-canceling infrared divergences [13,14]. Hence when we use the gauge invariant definition of the color singlet P-wave non-perturbative matrix element as given by eq. (5) then we do not expect any uncanceled infrared divergences due to real soft gluons emission/absorption in the color singlet P-wave heavy quarkonium production.
It is useful to mention here that the heavy quarkonium is also an useful signature [17] to study production and detection of quark-gluon plasma at RHIC and LHC. In high energy heavy-ion collisions at RHIC and LHC the study of heavy quarkonium production is more complicated than that in pp collisions due to the presence of QCD medium in nonequilibrium.

NON-PERTURBATIVE MATRIX ELEMENT
It can be mentioned here that in the NRQCD approach of heavy quarkonium study the S-wave color octet non-perturbative matrix element is used to factorize the above mentioned non-canceling infrared divergences which arise due to real soft gluons emission/absorption (eikonal gluons) in the pQCD partonic level calculation in the color singlet P-wave heavy quarkonium production [12]. The NRQCD definition of the color singlet P-wave non-perturbative matrix element used for this purpose is given by [12] < which is gauge invariant where D is the covariant derivative.
Note that in QCD the definition of the color singlet P-wave non-perturbative matrix element of the heavy quarkonium is given by eq. (3) which contains the ordinary derivative ∇ whereas in NRQCD the definition of the color singlet P-wave non-perturbative matrix element of the heavy quarkonium is given by eq. (12) which contains covariant derivative D. The main motivation to use covariant derivative D in NRQCD instead of ordinary derivative ∇ was to make the the definition in eq. (12) gauge invariant. However, it has to be remembered that it is the relative momentum q of the heavy quark-antiquark pair which appears in the definition of the color singlet P-wave non-prturbative matrix element [7][8][9] which can be obtained if ordinary derivative ∇ operates on the heavy quark field which can be seen from the definition in eq. (3) in QCD.
Of course the appearance of ordinary derivative ∇ makes the definition of the color singlet P-wave non-perturbative matrix element of heavy quarkonium in eq. solve infrared divergences problem in color singlet mechanism, we do not bring color octet mechanism to solve infrared divergences problem in color singlet mechanism.
In any case let us briefly discuss the issue of lack of gauge invariance of the S-wave color octet non-perturbative matrix element of the heavy quarkonium production in NRQCD.
In order for the non-canceling infrared divergences due to real soft gluons (eikonal gluons) emission/absorption in the color singlet P-wave heavy quarkonium production to be factorized into the definition of the NRQCD S-wave color octet non-perturbative matrix element as given by eq. (11), the two non-perturbative matrix elements in NRQCD in eqs. (11) and (12) are related by the equation (see eq. (6.5) of [12]) which has the problem of gauge invariance in NRQCD because the left hand side of eq. (13) is not gauge invariant whereas the right hand side of eq. (13) is gauge invariant where Λ is the ultraviolet cutoff of NRQCD. This is because the definition of the NRQCD S-wave color octet non-perturbative matrix element as given by eq. (11) in NRQCD in [12] is not gauge invariant. Note that the gauge links are supplied in the NRQCD S-wave color octet non-perturbative matrix element in [1,13,14] to make it gauge invariant and consistent with factorization of infrared divergences.
It can be mentioned here that in NRQCD, the (the long distance matrix elements) LDMEs are organized as powers of the relative velocity (v) of the heavy quarks in the heavy-quarkonium states. At leading order in v 2 , |R ′ (0)| 2 in eq. (3), which is usually calculated in potential models in Coulomb gauge, is gauge invariant and so is the evolution equation eq. (13).

GLUON) EMISSION/ABSORPTION IN QCD
As mentioned in section II the non-canceling infrared divergences occur in the color singlet P-wave heavy quarkonium production in the processes like that in eq. (10) due to real soft gluons (eikonal gluons) emission/absorption. Hence in this section we will describe the infrared divergence due to real soft gluon (eikonal) emission/absorption in QCD. Before describing infrared divergence due to soft gluon (eikonal gluon) emission/absorption in QCD let us first describe the infrared divergence due to soft photon (eikonal photon) emission/absorption in QED.
Consider the emission of a real photon with four momentum k µ from an electron of four momentum p µ in QED. The Feynman diagram contribution from this real photon emission process in QED is given by [18] 1 where we have used The photon field ǫ µ (k) can be written as where ǫ µ phys (k) is the transversely polarized (physical) photon field and ǫ µ pure (k) = k µ p·k p · ǫ(k) is the (unphysical) longitudinally polarized pure gauge photon field.
Hence we expect that it is the transversely polarized (physical) photon field ǫ µ phys (k) in eq. (16) that contributes to the physical (finite) value of the cross section but the (unphysical) longitudinally polarized pure gauge photon field ǫ µ pure (k) does not contribute to the physical (finite) value of the cross section. This can be seen by using eq. (16) in (14) to find and k ǫ pure (k) 2p · k = 0.
Eq. (17) is the non-eikonal contribution from this soft photon emission Feynman diagram which contributes to the finite value of the physical cross section whereas from eq. (18) we find that the (unphysical) longitudinally polarized pure gauge photon field ǫ µ pure (k) can not contribute to the finite value of the physical cross section. The (unphysical) longitudinally polarized pure gauge photon field ǫ µ pure (k) in eq. (16) accounts for the infrared divergence from this soft photon emission Feynman diagram. This can be seen as follows.
Eq. (19) gives the infrared divergence in the eikonal approximation from this real soft photon (eikonal photon) emission Feynman diagram. From eq. (20) we find that the transversely polarized (physical) photon field ǫ µ phys (k) in eq. (16) does not contribute to the infrared divergence from this real soft photon emission Feynman diagram in the eikonal approximation.
Hence from eqs. (19), (18) and (14) we find that if the charge produces pure gauge field then the study of infrared divergence due to soft photon (eikonal photon) emission from that charge can be simplified in quantum field theory by using pure gauge.
In our calculation the momentum P µ 1 of the heavy quark can be different from the momentum P µ 2 of the heavy antiquark where If the photon interacts with the electron and positron simultaneously then that photon is the virtual photon and we can consider the photon propagator G µν (k) in that situation instead of the polarization vector ǫ µ (k) of the photon which we considered in eq. (16) for the real photon emission from the electron only.
Hence if the photon is attached simultaneously to the electron of momentum p µ and to the positron of momentum q µ then the photon propagator G µν (k) can be split into the physical (transverse) propagator G µν phys (k) and the pure gauge (longitudinal) propagator G µν pure (k) as follows which is equivalent to eq. (16) of our paper where in the infrared (IR) limit (see eq. (2.1) of [18]) and For the virtual photon interacting simultaneously with the electron of momentum p µ and the positron of momentum q µ the eikonal Feynman rule for the infrared (IR) divergences is given by (see, for example, Eq. (29) of [14]).
Hence from eqs. (22), (23), (24) and (25) we find and The above equation is equivalent to eq. (20) of our paper. Hence one finds that the pure gauge part is the source of the infrared divergence.
It is well known that the light-like charge produces pure gauge field at all the time-space points x µ except at the positions perpendicular to the direction of the motion of the charge at the time of closest approach in classical mechanics [19][20][21]. This is also true in quantum field theory which can be shown as follows.
The generating functional for the photon field Q µ (x) in the presence of external source J µ (x) is given by which gives the effective action [1,22] where L ef f (x) is the effective lagrangian density.
The eikonal current density J µ eik (x) of the charge with light-like four velocity l µ can be obtained from the eikonal term of eq. (14) which gives the light-like eikonal current density Using eq. (31) in (29) we find From eq. (14) the non-eikonal contribution gives the non-eikonal current density Using eqs. (34) and (31) in eq. (29) we find that the effective interaction lagrangian density due to the interaction between the (light-like or non-light-like) non-eikonal current and the gauge field generated by the light-like eikonal current is given by which implies that the light-like eikonal line can be replaced by the pure gauge background field to study factorization of infrared divergences in quantum field theory without modifying the finite value of the physical cross section.
Hence we find that since the light-like eikonal current produces pure gauge field, the study of infrared divergence due to soft photon (eikonal photon) emission due to the presence of light-like eikonal line can be simplified in quantum field theory by using pure gauge without modifying the finite value of the cross section. Since the light-like eikonal current produces pure gauge field in classical mechanics and in quantum field theory, the factorization of infrared divergence due to the presence of light-like eikonal line can be studied by using path integral formulation of the background field method of quantum field theory in the presence of pure gauge background field [1,[23][24][25].
In QED the light-like electric charge produces U(1) pure gauge field A µ (x) given by In QCD the light-like color charges produce SU(3) pure gauge field given by [1,23,24] where Φ(x) is the gauge link given by eq. (6) where l µ is the light-like four velocity.

V. PROOF OF FACTORIZATION OF INFRARED DIVERGENCE IN COLOR
SINGLET P-WAVE HEAVY QUARKONIUM PRODUCTION As discussed above the factorization of infrared divergences due to the soft gluons emission/absorption with the light-like eikonal line can be studied by using the path integral formulation of the background field method of QCD [1,23,24] in the presence of SU(3) pure gauge background field A µa (x) as given by eq. (37). We follow the notation of [26][27][28] and denote the quantum gluon field by Q µa and the background field by A µa .
In the path integral formulation of QCD we find [26,29] < 0|Ψ( and in the path integral formulation of the background field method of QCD we find [26][27][28] < 0|ψ( where and is the gauge fixing term and the type I gauge transformation in the background field method of QCD is given by In eqs. (38) and (39) the ψ l is the Dirac field of the light quark of flavor l = 1, 2, 3 = u, d, s of mass m l , the Ψ is the Dirac field of the heavy quark of mass M and∇ is defined by eq. (4).
Changing the integration variable Q → Q − A inside the path integration in eq. (39) we where with factorization of infrared divergences in QCD at all orders in coupling constant is given by which reproduces eq. (5) where the gauge link Φ(x) is given by eq. (6) and∇ is defined in eq. (7).
Since the left hand side of eq. (49) is independent of the light-like vector l µ one finds that the right hand side of eq. (49) is independent of the light-like vector l µ used to define the gauge link Φ(x) in eq. (6). This proves that the long-distance behavior of the color singlet P-wave non-perturbative matrix element < 0|χ † Φ∇Φ † ξa † H · a H ξ † Φ∇Φ † χ|0 > in eq. (50) is independent of the light-like vector l µ at all order in coupling constant in QCD.
This completes the derivation of the correct definition of the non-perturbative matrix element of color singlet P-wave heavy quarkonium production which is gauge invariant and is consistent with the factorization of infrared divergences at all order in coupling constant in QCD.
Note that the relationship between the color singlet P-wave long distance matrix element (LDME) and the derivative of the non-relativistic wave function at origin |R ′ (0)| 2 is very useful because it helps us to reduce the number of non-perturbative input parameters to be determined. Hence it is important to discuss how this relationship is affected by the redefinition in eq. (5). First of all it is found that if one uses the derivative of the nonrelativistic wave function at the origin |R ′ (0)| 2 in eq. (1) then the infrared (IR) divergences do not cancel [5,[7][8][9], i. e., the non-perturbative effects can not be factored into the derivative of the non-relativistic wave function at the origin |R ′ (0)| 2 . In order to cancel the infrared (IR) divergences and to be consistent with the factorization theorem at all orders in coupling constant the definition of the color singlet P-wave non-perturbative matrix element from eq. (5) should be used in eq. (1). Note that since the color singlet P-wave non-perturbative matrix element in eq. (5) is a non-perturbative quantity it can not be studied by using perturbative QCD. Due to these infrared (IR) divergences, factorization and non-perturbative QCD issues there is no simple relation between the color singlet Pwave non-perturbative matrix element in eq. (5) at all orders in coupling constant and the derivative of the non-relativistic wave function at the origin |R ′ (0)| 2 . The color singlet P-wave non-perturbative matrix element in eq. (5) at all orders in coupling constant can be studied by using the non-perturbative QCD. In the phenomenological studies of the infrared sensitive processes in color singlet P-wave heavy quarkonium production, the definition from eq. (5) should be used in eq. (1) to extract the non-perturbative matrix element from the experiments in order to be consistent with the factorization of the infrared (IR) divergences at all orders in coupling constant.

VII. CONCLUSIONS
Recently we have proved factorization of infrared divergences in NRQCD S-wave heavy quarkonium production at high energy colliders at all orders in coupling constant. One of the problem which still exists in the higher order pQCD calculation of color singlet P-wave heavy quarkonium production/anihillation is the appearance of non-canceling infrared divergences due to real soft gluons exchange, although no such infrared divergences are present in the color singlet S-wave heavy quarkonium.
In this paper we have found that since the non-perturbative matrix element of the color singlet P-wave heavy quarkonium production contains derivative operators, the gauge links are necessary to make it gauge invariant and be consistent with the factorization of such non-canceling infrared divergences at all orders in coupling constant. In case of color singlet S-wave heavy quarkonium production the gauge links cancel because of the absence of derivative operators which is consistent with the absence of such non-canceling infrared divergences at all orders in coupling constant.