New contributions to gluon poles in direct photon production

We consider the direct photon production in two hadron collision with one of the hadrons being transversely polarized. By using the contour gauge for gluon fields, we find that there are new twist-$3$ terms present in the hadron tensor of the considering process in addition to the standard twist-$3$ terms. In this work, we demonstrate that the significance of these new terms are two-fold: first, they are crucial to both the QED and QCD gauge invariance and, second, their contributions to the hadron tensor are at least the same as those from the standard ones. We also study the resuting effects which are responsible for the universality breaking of the corresponding twist-$3$ parton distributions.


I. INTRODUCTION
The problem of the electromagnetic (QED) gauge invariance in the deeply virtual Compton scattering (DVCS) and similar exclusive processes has intensively been discussed during the last few years (for example, see [1][2][3][4][5]). This development explored the similarity with the earlier studied inclusive spin-dependent processes [6]. The gauge invariance of hard process amplitudes is ensured by twist three contributions and by the use of the equations of motion that provide a possibility to exclude the three-particle (quark-gluon) correlators from the amplitude. So that, after combining with the two-particle correlator contributions, one gets the gauge invariant expressions for the physical amplitudes or, in the case of lepton-hadron processes, for the corresponding hadron tensors [6]. This scheme was originally developed in the case of the particular inclusive processes with transverse polarized hadrons, like structure function g 2 in DIS [6] and Single Spin Asymmetry (SSA) [7] due to the soft quark (fermionic poles [8]). Also, the QCD gauge invariance of the so-called gluonic poles contributions [9] has been a subject of studies in [10] where the methods that are used rely on the Wilson exponentials [11][12][13][14].
We have shown in our recent work [15] that, in order to ensure the QED gauge invariance of the transverse polarized Drell-Yan (DY) hadron tensor, it is mandatory to include the contribution from the extra diagram originating from the non-trivial imaginary part of the corresponding twist-3 function B V (x 1 , x 2 ). Before this study [15], the function B V (x 1 , x 2 ) has been argued as real, and the imaginary part of amplitude was ensured by means of a specially introduced "propagator" 1 in the hard part of the hadron tensor [17]. However, we have explained in [15] that the B V -function does, in fact, have an imaginary part, and the existence of this imaginary part can be realized with a help of the contour gauge. Moreover, the B V -function with the complex prescription induces a new contribution to the hadron tensor. As has been stressed, this extra contribution leads to the amplification of corresponding tensor by a factor of 2. This finding of ours has independently been confirmed in [16] by using of a different approach. Besides, from the point of view of phenomenology, the corresponding SSAs in the DY process and the role of gluon pole contributions have previously been discussed in [18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].
In the present paper, we extend our approach used in [15] to the case of the direct photon production in two hadron collision where one hadron is transversely polarized. We derive the hadron tensor for this process and study the effects which lead to the soft breaking of factorization (or the universality breaking) through the QED and QCD gauge invariance. In a similar manner as in [15], the special role is played by the contour gauge for gluon fields. We demonstrate that the prescriptions for the gluonic poles in the twist 3 correlators are dictated by the prescriptions in the corresponding hard parts. We argue that the prescriptions in the gluonic pole contributions differ from each other depending upon the initial or final state interactions of the diagrams under consideration. Moreover, the different prescriptions are needed to ensure the QCD gauge invariance. We treat this situation as a breaking of the universality condition resulting in factorization soft breaking. The extra diagram contributions, which naively do not have an imaginary phase, is discussed in detail.
In the paper, we also show that the new ("non-standard") terms do contribute to the hadron tensor exactly as the "standard" terms known previously. This is exactly similar to the case of Drell-Yan process studied in [15].

II. GETTING STARTED: CASE OF DRELL-YAN PROCESS
For pedagogical reasons, we remind briefly our findings for the Drell-Yan (DY) process where one of hadrons possesses the transverse polarization, see [15] for all details. As usual, the Drell-Yan process with the transversely polarized nucleon is defined as where the virtual photon producing the lepton pair (l 1 + l 2 = q) has a large mass squared (q 2 = Q 2 ) while the transverse momenta are small and integrated out. The dominant light-cone directions for the DY process ( Fig. 1) are defined as p 1 ≈ n * Q/(x B √ 2) and p 2 ≈ n Q/(y B √ 2) with the dimensionless light-cone vectors n * µ = (n * + , 0 − , 0 ⊥ ) and n µ = (0 + , n − , 0 ⊥ ).
Since we deal with a large Q 2 in the process under consideration, it is possible to apply the factorization theorem to get the corresponding hadron tensor factorized in the form of convolution: Usually both the hard and soft parts in Eqn. (2) are independent of each other, UV-and IR-renormalizable. Moreover, various parton distributions which parametrize the soft part have to manifest the universality property. Based on DY-process, it is convenient to study the role of twist 3 by exploring of different kinds of asymmetries, for instance, the left-right asymmetry. This left-right asymmetry means the transverse momenta of the leptons are correlated with the direction S × e z where S µ implies the transverse polarization vector of the nucleon and e z is a beam direction [33]. Generally speaking, any single spin asymmetries (SSAs) can be presented in the following symbolical form (at this moment, the exact expression for SSA is irrelevant): where L µν is an unpolarized leptonic tensor and, consequently, has only the real part; H µν stands for the hadronic tensor which is also real. Since one of hadrons is transversely polarized, the corresponding matrix element which forms the soft part of hadron tensor reads [15] where the light-cone vectorñ is a dimensionful analog of vector n. Therefore, in order to provide for the condition, H µν ∈ e, the complex i in the r.h.s. of (4) has to be compensated either (a) by the complexness of the hard part (this is the standard contribution): or (b) by the complexness of the soft part In Eqns. (5) and (6), F ∼ and O(ψ, ψ, A) are the shorthands for the Fourier transformation and the corresponding quark-gluon operator, respectively (Eqn. (4)). In general, the hard parts in Eqns. (5) and (6) differ from each other. For instance, the hard part of the diagram presented in Fig. 1(a) contains the quark propagator in contrast to the hard part of the diagram in Fig. 1 However, in the previous studies (for example, [10,17,26,27]), B V (x 1 , x 2 )-function has been assumed to be a purely real function: with the function T (x 1 , x 2 ) ∈ e which parametrizes the corresponding projection of ψ G αβ ψ . Therefore, the scenario (b) (Eqn. (6)) will never be realized if the B V (x 1 , x 2 )-function has the form as in Eqn. (7). As a result, the QED gauge invariance of the DY hadron tensor is in question. Indeed, having analyzed the hard subprocess in the context of the QED gauge invariance, we can conclude that only the sum of two diagrams in Fig. 1 (a) and (b) ensures the QED gauge invariance of the hadron tensor. As shown in [15], the contribution of H µν is generated by the diagram presented in Fig. 1(b).
Thus, we may infer that the real function B V (x 1 , x 2 ) in the form as in Eqn. (7) leads to the violation of the QED gauge invariance of the hadron tensor.
To solve this discrepancy, it is instructive to remind the reason which leads to the representation in Eqn. (7). The conclusion that B V is a real function has come from the ambiguity in the solutions of the differential equation (provided that A + = 0): Indeed, the formal solutions of Eqn. (8) have the following forms: From the first glance, Eqns. (9) and (10) seem to be equivalent each other. However, as we will see below, this is not true. Now, if we insert Eqns. (9) and (10) into Eqn. (4) and use the following parametrization 2 The corresponding δ-functions appeared in the hadron tensor and expressed the momentum conservation low should be also referred to the hard parts. This statement was argued in [34] in the context of the so-called factorization links.
we get the following representations: respectively. In Eqns. (12) and (13), the corresponding prescriptions ±i arise from the integral representation for the theta-function: Further, if we suppose that the representations (12) and (13) are equivalent to each other, we can calculate the plus and minus combinations of (12) and (13) resulting in The ambiguity in the solutions of (8) ( (12) and (13)) ultimately gives us the standard representation (7) with the real function B V provided the asymmetric boundary condition for gluons is given by . In fact, the representations (12) and (13) are not equivalent ones [15]. To see that it is necessary to remember that all axial-type gauges are the particular cases of the most general contour gauge (see Appendix A and the text below for details). Using the contour gauge conception, one can easily check that the representation (12) belongs to the gauge [x, −∞] = 1, while the representation (13) belongs to the gauge [+∞, x] = 1. Therefore, one has no any reason to believe that (12) and (13) are equivalent ones (for details, see (A.12), (A.13)),i.e.
Roughly speaking, it resembles the trivial situation where two different vectors have the same projection on the certain direction. In this context, the well-known axil gauge A + = 0 can be treated as some type of "projection" which corresponds to two different "vectors" represented by two different contour gauges (see Appendix A). We now consider the QED gauge invariant hadron tensor for the DY-process with the transverse polarization. In order to understand which contour gauge we need to deal with, before imposing the condition A + = 0, it is necessary to take into account the contributions of p 1 , S T |ψ γ + A + ψ|S T , p 1 in the standard hadron tensor (see Fig. 1(a)). Based on the analysis of the γ-structure of this diagram as shown in [15], the Feynman causal prescription in the quark propagator (see, eqns. (2), (6) and (8) of [15]) uniquely leads to the Wilson line in the quark-gluon correlator: Eqn. (18) suggests that we have to use the contour gauge [−∞ − , 0 − ] = 1. Therefore, the contour gauge defined by [−∞ − , 0 − ] = 1 destroys the ambiguity, and the function B V (x 1 , x 2 ) has to be described by the following representation (see, (A.12) and [15] for details) As a result, the diagram represented in Fig. 1(b) (or the hadron tensor (6)) does contribute to the hadron tensor, and this together with the first diagram represented in Fig.1(a) form the gauge invariant (GI) hadron tensor (see, Eqn. (34) of [15]): From Fig. 1, we can also realize that the representation of B V (x 1 , x 2 ) with the complex prescription +i in the gluonic pole, see (A.12), corresponds to the initial state interaction (ISI) with respect to the hard subprocess.
We want to stress that in the case of DY-process we deal with the initial state interaction only as opposed to the case of the direct photon production which is studied below. For the DY-process, the initial state interaction generates −2 + k − 2 + i (see, the diagram in Fig. 1(a)) in the quark propagator which, in turn, leads to (i) the contour gauge [−∞ − , 0 − ] = 1 and, then, to (ii) the function B V + with the certain complex prescription (A.12). The latter ensures the QED gauge invariance for the hadron tensor. Schematically, the mentioned logical chain can be presented as We can see that the prescription in the quark propagator of the hard part gives information on the contour gauge for gluons from the soft part. In other words the hard and soft parts are not fully independent of each other. Despite of this, the DY hadron tensor has formally been factorized with the mathematical convolution, and the parton distributions, such as the twist-3 function B V (x 1 , x 2 ), still satisfy the universality condition. In contrast to the DY-process, as we will see in the next section, the direct photon production tensor is built with the functions B V (x 1 , x 2 ) that will not manifest the universality.
We will refer to a soft breaking of factorization when the factorization procedure results in the mathematical convolution between the finite hard and soft parts, but there is no universality for the soft functions or the hard and soft parts are not totally independent.
To conclude this section, let us note that the situation with the gauge invariance for the DY-process (21) is very similar to what has been discussed for the vector meson electroproduction in [39] where the final state interaction also predetermined the correct prescription for the spurious gluon pole for the validity of QCD factorization.

A. Kinematics
In this section, we study the two hadron collisions, where one of the hadrons possesses the transverse polarization, which produce the direct photon in the final state in: The gluonic poles are being manifested in this process in the similar way as for the Drell-Yan process [26]. We perform our calculations within a collinear factorization, and it is convenient (see, e.g., [34]) to fix the dominant light-cone directions as The hadron momenta p 1 and p 2 have the plus and minus dominant light-cone components, respectively. Accordingly, the quark and gluon momenta k 1 and lie along the plus dominant direction while the gluon momentum k 2 -along the minus direction. The final on-shell photon and quark(anti-quark) momenta can be presented as The Mandelstam variables for the process and subprocess are defined as The amplitude of process (22) involves the contributions from • the leading (LO) diagrams: two diagrams with a radiation of the photon before (A LO 1 ) and after (A LO 2 ) the quark-gluon vertex with the gluon going to the lower blob, see the right side of Fig. 2; • the next-to-leading order (NLO) diagrams: eight diagrams constructed from the LO diagrams by insertion of all possible radiations of the additional gluon which together with the quark goes to the upper blob, see the left side of Fig. 2.
We denote the sum of diagrams as that generates the hadron tensor related to the corresponding asymmetry: In this work, we mainly dwell on the discussion of the hadron tensor rather than the asymmetry itself. Diagrammatically (Fig. 2), the hadron tensor can be presented in the form of an interference between the LO and NLO diagrams: In Fig. 2, the upper blob determines the matrix element of the twist-3 quark-gluon operator while the lower blob -the matrix element of the twist-2 gluon operator related to the unpolarized gluon distribution.

B. Factorization procedure
Since the collinear factorization is our main tool, let us outline the main stages of factorization. The factorization procedure contains the following steps: • the decomposition of loop integration momenta around the corresponding dominant direction: k i = x i p + (k i · p)n + k T within the certain light cone basis formed by the vectors p and n (in our case, n * and n); • the replacement: that introduces the fractions with the appropriated spectral properties; • the decomposition of the corresponding propagator products around the dominant direction: • the use of the collinear Ward identity, if it requests by the needed approximation: • performing of the Fierz decomposition for ψ α (z)ψ β (0) in the corresponding space up to the needed projections.
Notice that, for our purposes, it is enough to be limited by the first order of decomposition in the third item.
As a result of this procedure, we should reach the factorized form for the considered subject, see (2).

C. Hadron tensor: QED gauge invariance
At the first item, we want to discuss the QED gauge invariance of the hadron tensor. To check the QED gauge invariance, it is sufficient to consider the typical Feynman diagrams H1, H3 and H5 represented in Fig.  3. For the definiteness, we pay our attention on the anti-quark contribution. All our results can trivially be extended to the quark contribution as well.
Before factorization, the H1-diagram in Fig.3 leads to the following expression: where C 2 implies the corresponding colour factor. Here and in what follows the coupling constants are not shown explicitly. In Eqn. (28) the gluon (unpolarized) twist-2 parameterizing function is defined as while the quark-gluon twist-3 parameterizing function is given bȳ We now carry out the standard factorization procedure and, after some algebra, obtain the following expression: In the diagrams H1, H3 and H5, we deal with the final state interaction (FSI) with respect to the hard part. Therefore, the B V -function has the representation as in (A.13) (cf. [15]), i.e.
In the similar way as in the preceding section, we have to restore the path in the Wilson line with A + -fields in the quark-gluon correlators which appear in the diagrams H1, H3 and H5 in Fig.3. After straight-forward calculations, we derive the Wilson line in the form [+∞ − , z − ] that suggests us to use the representation (A.13) for our B V -function. The contribution of H5-diagram in Fig.3 is equal to zero due to the fact that the photon momentum has the dominant minus light-cone component and, therefore, the γ-structure gives (γ − ) 2 = 0.
We now calculate the hadron tensor term associated with the H3-diagram in Fig.3, it reads where the B V -function is also given by the representation (33) or (A.13).
After some γ-algebra, we can check that the contribution of (34) is equal to the contribution of (31) but with an opposite sign. Therefore, the sum of all contribution gives us zero (we remind that while the second diagram contribution is equal to zero itself) i.e.
In fact, the identity (35) reflects the QED gauge invariance for the hadron tensor. We emphasize that the QED gauge invariance takes place owing to the same complex prescriptions in the definitions of B V -function (Eqn. 33). In turn, the same complex prescriptions emanate from the final state interaction presented in H1, H3 and H5 of Fig.3. Finally, using the similar logical chain as for the DY-process, we can write that We would like to emphasize that the concrete sign of the gluonic pole prescription is not so crucial for the QED gauge invariance because here we deal with only one type of interaction which is the final state interaction (this distinguishes the case of QCD gauge invariance which is considered below). It is more important to have the same prescriptions in all gluonic poles. In other words, from the point of view of contour gauge, we may use the same wrong +i prescription in the representation of B V -function for diagrams H1, H3 and H5 depicted in Fig.3. But this wrong prescription still leads to the QED gauge invariance. However, after calculation of the imaginary parts for the corresponding asymmetry, the wrong prescription definitely plays a negative role.

D. Hadron tensor: QCD gauge invariance
We are now in a position to dwell on the QCD gauge invariance of the hadron tensor for the direct photon production. To check this invariance, we have to consider four typical diagrams H1, H5, D1 and H9, depicted in Fig.3, which come from the corresponding ξ-process (see, [40]). Notice that the gluon enters in the quark-gluon correlator as an internal field. For the QCD gauge invariance, we have to assume that all charged particles are on-shell, i.e we deal with the physical gluons only. The substantial differences between this case and a pure perturbative Compton scattering case are discussed in Appendix B.
To write down the Ward identity, we need to replace the gluon transverse polarization T α by the gluon longitudinal momentum L α in the quark-gluon correlator: where a + ( ) stands for the gluon creation operator. The summation over the intermediate states is not shown explicitly. Notice that the parametrization of this correlator through B V -function stays unchanged. Consider now the contribution of the H1-diagram in Fig.3 to the hadron tensor. Before going further, it is instructive to begin with the gluon loop integration corresponding to the mentioned diagram. We have where we do not explicitly write the operators which are irrelevant at the moment (cf. (37)). After factorization, we obtain where we decompose the hard part around the dominant direction and put L = (x 2 − x 1 )p 1 which is actually dictated by the γ-structure and the momentum conservation. Using all these, we get the following expression: where As can be seen, however, that this diagram does not contribute to the Ward identity. Indeed, after calculation of the imaginary part we get the factor (x 2 − x 1 ) in the numerator of (40) which goes to zero owing to δ( Further, calculation of the H5-diagram, presented in Fig.3, gives us while the contribution of the D1-diagram in Fig.3 takes the form Finally, the contribution of the H9-diagram with the three-gluon vertex, see Fig.3, reads We now turn to the contour gauge. First of all, based on our previous discussions notice that even a fleeting glance is sufficient to anticipate the corresponding prescriptions for B V -functions in (42)-(44). The H5-diagram in Fig.3 corresponds to the final state interaction and, therefore, the function B V − should appear here. On the other hand, the D1-and H9-diagrams in Fig.3 correspond to the initial state interaction which leads to the function B V + . Performing the explicit calculations (see also [15]), we can arrive at the conclusion mentioned above by restoring the Wilson lines in the quark-gluon correlators of the mentioned diagrams. The Wilson line, [+∞ − , z − ], enters in the hadron tensor represented by the H5-diagram in Fig.3 while the Wilson line, [z − , −∞ − ], appears in the hadron tensor represented by the D1-and H9-diagrams in Fig.3.
We sum all contributions and get the following final expression: where the function B V − are represented by (A.13) or (33), and the function B V + is given by (A.12) or We are now calculating the imaginary part and, ultimately, we derive the QCD Ward identity in the form We want to stress that the identity (47) takes place provided only the presence of the different complex prescriptions in gluonic poles dictated by the final or initial state interactions: We emphasize on the principle differences (see details in Appendix B) between the considered case and the proof of the QCD gauge invariance for the perturbative Compton scattering amplitude with the physical gluons in the initial and final states. The latter does not need any external condition such as the presence of gluon poles. Thus, the situation which we discuss is again absolutely similar to that one which has been described in [39] for the dijet production. From (48), it is seen that the different diagrams correspond to the different contour gauges and, consequently, to the different functions, B V ± , which parametrize the hadronic matrix element forming the soft part. In this context, we also have a soft breaking of factorization because, first, it spoils the universality principle and, second, the gluonic pole prescriptions in the soft part are traced to the causal prescriptions in the hard part. At the same time, we can use the replacement [39]: and finally get the same function B V + for all diagrams (in other words, we can use the same contour gauge for all diagrams). However, it contains the additional δ(x 1 − x 2 )-term which may lead to the collinear factorization violation in the same manner as in [39]. The full analysis of this case will be implemented in our forthcoming work.

IV. HADRON TENSOR OF THE DIRECT PHOTON PRODUCTION II: NEW CONTRIBUTIONS
In this section we calculate the full expression for the hadron tensor which involves both the standard and new contributions to the gluon pole terms. The full expression for the hadron tensor related to the case we are discussing can be split into two groups: (i) the first type of contributions corresponds to the diagrams H1-H12 depicted in Fig. 3 and, before factorization, takes the following form and (ii) the second type of contributions is given by the diagrams D1-D4 in Fig. 3 can be presented as In eqns. (50) and (51), the corresponding coefficient functions are denoted by H αβ,ρ (k 1 , k 2 , ) and D αβ (k 1 , k 2 ).
The unpolarized twist-2 gluon distribution Φ g (k 2 ) and the twist-3 quark distribution Φ in the standard forms, see (29) and (30). The twist-3 quark distribution which appears in the diagrams D1-D4 presented in Fig. 3 is given by where the sum over the corresponding intermediate states is implied. We now perform the factorization procedure for Eqns. (50) and (51), and we obtain for the first type of contributions and for the second type of contributions. To simplify our calculations without losing generality, we may impose the frame where q 2 ⊥ S. The Mandelstam variable defined for the subprocess,û, is a small variable and can be neglected. It means that the Bjorken fraction y B becomes independent of x B , and one can write y B = −T /S (due toŝ +t +û = 0).
The next stage is to determine the type of the twist-3 function B V (x 1 , x 2 ) which is related to the certain complex prescriptions in gluon poles according to the way described in the preceding sections. The correct definition of B V (x 1 , x 2 )-function should be implemented for each of diagrams. Also, it is instructive to notice that the diagrams H1-H8 would not possess the gluon poles in the case where the function B V assumed to be a real one. However, this is not true in our case.
After computing the corresponding traces and performing simple algebra within the frame we are choosing, it turns out that the only nonzero contributions to the hadron tensor come from the diagrams H1, H7, D4 and H10: and Here, The other diagram contributions disappear owing to the following reasons: (i) the γ-algebra gives (γ − ) 2 = 0; (ii) the common pre-factor T + yS goes to zero, (iii) the diagrams H2 and H5 cancel each other.
Analysing the results for the diagrams H1, H7, D4 and H10 (see Eqns. (55)-(58)), we can see that In other words, as similar to the Drell-Yan process, the new ("non-standard") contributions generated by the diagrams H1, H7 and D4 result again in the factor of 2 compared to the "standard" diagram H10 contribution to the corresponding hadron tensor. This is our principle result.

V. CONCLUSIONS
In this work, we explore both the QED and QCD gauge invariance of the hadron tensor for the direct photon production in two hadron collision where one of hadrons is transversely polarized. We present the important details related to the use of the contour gauge for the Drell-Yan process with essential transverse polarizations. We study the effects which lead to the soft breaking of factorization through the QED and QCD gauge invariance.
We show that the contour gauges for gluon fields play the crucial role for our study. The contour gauge belongs to the class of non-local gauges that depends on the path connecting two points in the correlators. It turns out, in the cases which we consider, the prescriptions for the gluonic poles in the twist 3 correlators are dictated by the prescriptions in the corresponding hard parts in the similar manner as for the DY-process considered in [15].
For the direct photon production, we demonstrate that the prescriptions in the gluonic pole contributions differ from each other depending upon the initial or final state interactions in the related diagrams. We stress that the different prescriptions are needed to ensure the QCD gauge invariance. This situation has been treated as a soft breaking of the universality condition resulting in factorization breaking. Besides, the presence of the complex prescriptions in the gluonic pole contributions allow the extra diagrams to contribute nontrivially to the hadron tensor.
We find that the "non-standard" new terms, which exist in the case of the complex twist-3 B V -function with the corresponding prescriptions, do contribute to the hadron tensor exactly as the "standard" term known previously. This is another important result of our work. We also observe that this is exactly similar to the case of Drell-Yan process studied in [15]. denotes the group defined on the given fiber and π is a transformation of the base R 4 into the fiber bundle P). Each element g(x) of the fiber, with the help of the gluon field A α , defines the gauge-transformed field: The set of these fields for all g(x) forms the orbit of the gauge-equivalent fields. It is well-known that in order to quantize the system of the gauge fields, one has to choose the only element of each orbit. In contrast to the usual way, we first fix an arbitrary point (x 0 , g(x 0 )) 3 in the fiber. Then, we define two directions: one of them in the base, the other in the fiber. The direction in the base R 4 is nothing else than the tangent vector of a curve which goes through the given point x 0 . At the same time, the direction in the fiber can be uniquely determined as the tangent subspace which is related to the parallel transition. After following this procedure, one can uniquely define the point in the fiber bundle. Further, solving the parallel transport equation which is defined on the fiber as one can find the solution in terms of the Wilson line: where the points x 0 and x are connected by the path P. The starting point x 0 is usually fixed, i.e. x 0 is independent on x. Here, g(x 0 ) is chosen to be equal to unity 4 . Note that the fixing of g( The contour gauge condition demands that g(x) is equal to unity for all x belonging to the base, i.e.
Therefore, within the contour gauge, the field A g µ (see, (A.5)) becomes i.e. the gluon field A g µ is a linear functional of the tensor G µν . Let us briefly comment a choice of the boundary conditions for gluons. As is pointed out in [36], the starting point x 0 of the path P(x 0 , x) (see, (A.8)) is not well defined by construction. In turn, this leads to the presence of the so-called residual gauge transformation. To avoid this ambiguity, it is natural to fix the uncertainty in Eqn. (A.8) by A µ (x 0 ) = 0. In other words, the starting point is fixed, and it is independent from the destination of the path P(x 0 , x).
It is easy to see that the representations (9) and (10) (or the representations (12) and (13)) can be derived from (A.8) by fixing of the path P(x 0 , x) as a straight line connecting the point x with ∓ ∞, respectively. In fact, it means that the representations (12) and (13) correspond to two different contour gauges and give us two different representations for two different functions B V − (x 1 , x 2 ) and B V + (x 1 , x 2 ). These two representations are associated with the final and initial state interactions (Eqns. (33) and (46)).
Moreover, in order to get a concrete representation for gluons within the axial gauge, we can explicitly parametrize the straight line between the points x and ±∞ along the "minus" light-cone direction n − : (A.9) Then, using Eqn. (A.8), one gets Thus, we can conclude that the r.h.s. of (12) and (13) correspond to the different functions B V ± (x 1 , x 2 ): It is clear that these two representations can not be equivalent of each other. Also, it is reasonable to assume the zeroth boundary conditions for gluons to be B V A(±∞) = 0 which is in agreement with Refs. [15,36]. Note that the functions B V ± do not possess a certain property under the time-reversal transformation. This property becomes a well-defined one only after calculation of the imaginary part for the hadron tensor.
We thus have no the ambiguity in the solutions of (8). As a result, the B V -function has a non-trivial imaginary part which contributes to H In this Appendix, we discuss an essential difference between the QCD gauge invariance for the "perturbative" Comptom amplitude, q(k)g( )| S |g(k 2 )q(k 1 ) , and the "nonperturbative" analog of Compton amplitude, X(P X )q(k)| S |g(k 2 )A(P ) , where one of gluons is included in the loop integration (see, Figs.4, 5).

QCD gauge invariance of the perturbative Compton amplitude
We first consider the "perturbative" Compton amplitude represented in Fig.4. To proof the QCD gauge invariance, we assume that all gluons are physical ones with the transverse polarizations. The first diagram in Fig.4 gives us To check the gauge invariance with respect to, for example, the physical gluon with momentum , we have to replace the given polarization vector * ⊥ β on the momentum β , provided · * ⊥ = 0, in the amplitude (B.1). We have Using the momentum conservation, k 1 + k 2 = k + , and the corresponding equations of motion, the Eqn. (B.2) reduces to the following form: In the similar way, we can get the contribution of the second diagram presented in Fig.4. It reads Therefore, the sum of these two diagrams gives us the color commutator: In the Feynman gauge 5 , the third diagram in Fig. 4 contributes as (where we again replace the transverse gluon polarization on the gluon momentum: * ⊥ → ) Making use of ⊥ · = * ⊥ · = 0 and The sum of all three diagrams leads to which ensures the QCD gauge invariance. In the case of the "perturbative" Compton amplitude with physical gluons, the momentum conservation, δ (4) (k 1 + k 2 − − k), and the on-shellness of all external particles play important role for checking of the gauge invariance.

QCD gauge invariance of the nonperturbative Compton scattering amplitude
In order to get the diagrams in Fig. 5, we now attach the nonperturbative blobs to the diagrams in Fig. 4. In contrast to the "perturbative" Compton amplitude, the gluon with momentum, , with respect to which we check the gauge invariance, is included in the loop integration. As a result, the momentum conservation for the subprocess has a form: k 1 + k 2 = k.
Let us focus on two diagrams in Fig. 5 which are enough to demonstrate how the color commutator contribution can be formed. Both diagrams are given by X(P X )q(k)| S |g(k 2 )A(P ) =ū(k) ⊥ α (d 4 ξ)(d 4 η)e −ik2ξ+ikη P X | δ 2 S δψ(η)δA α (ξ) |P ψ=...=A=0 . (B.11) The first diagram in Fig. 5 contributes as This expression corresponds to the amplitude before factorization and, therefore, the quark with k 1 + and gluon with are off-shell. Next, we apply the factorization procedure and derive the following expression for the first diagram contribution (here we use the light-cone basis presented above): where we replace the transverse gluon polarization vector by the gluon momentum (x 2 − x 1 )n * . In Eqn. (B.14), we use the Fourier transformations for the quark and gluon fields. a + and b − denote the gluon creation and quark annihilation operators, respectively. In the similar manner, we consider the second diagram contribution of Having applied the factorization procedure, we derive the following expression: Analyzing Eqns. (B.14) and (B.16), one can see that in order to form the color commutator combination we have to insist on the "external" conditions which actually emanate from the presence of gluon poles. Indeed, if the gluon poles are present, the amplitude B(dia.1) is accompanied by the factor of iπδ(x 2 − x 1 ). At the same time, the amplitude B(dia.2) has the factor of −iπδ(x 2 − x 1 ) in according to the FSI and ISI prescriptions, see (48). Thus, to check the QCD gauge invariance, the case of the "perturbative" Compton amplitude with the physical gluons in the initial and final states does not need the external condition for the presence of gluon poles. It is sufficient to use only the momentum conservation, k 1 + k 2 = + k, and the equations of motion for the initial and final quarks. At the same time, in order to demonstrate the QCD gauge invariance for the case represented in Fig. 5 where one of gluon momenta belongs to the loop integration, the existence of the gluon poles with the corresponding FSI and ISI prescriptions must be included.