Unifying approaches: derivation of Balitsky hierarchy from the Lipatov effective action

We consider a derivation of the hierarchy of correlators of ordered exponentials directly from the Lipatov’s effective action (Lipatov in Nucl Phys B 452:369, 1995; Phys Rep 286:131, 1997; Subnucl Ser 49:131, 2013; Int J Mod Phys Conf Ser 39: 1560082, 2015; Int J Mod Phys A 31(28/29):1645011, 2016; EPJ Web Conf 125: 01010, 2016) formulated in terms of interacting ordered exponentials (Bondarenko and Zubkov in Eur Phys J C 78(8), 617 2018; Bondarenko et al. in Eur Phys J C 81(7):61, 2021). The derivation of the Balitsky equation (Balitsky in Nucl Phys B 463:99, 1996; Phys Rev D 60:014020, 1999; At the frontier of particle physics, vol. 2, p. 1237–1342; Nucl Phys B 629:290, 2002; Phys Rev D 72:074027, 2005) from the hierarchy is discussed as well as the way the sub-leading eikonal corrections to the Balitsky equation arise from the transverse field contribution and sub-leading eikonal corrections to the quark propagator. We outline other possible applications of the proposed calculation scheme.


Introduction
The Lipatov's effective action approach [1][2][3][4][5][6] and the Balitsky-Kovchegov (BK)-Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov and Kovner (JIMWLK) approaches , see also examples of the phenomenological applications in , are both intended to describe amplitudes of QCD scattering at high-energies, despite the fact that the two independent approaches are formulated in terms of different degrees of freedom. The Lipatov's action deals with scattering amplitudes in terms of reggeized gluons, it allows to calculate quasi-elastic amplitudes of high-energy scattering processes in the multi-Regge kinematics. There are many applications of this action for the description of high energy processes and calculation of sub-leading unitarizing correca e-mail: sergeyb@ariel.ac.il (corresponding author) tions to the amplitudes and production vertices, see for example . In the BK-JIMWLK approach on the other hand  one treats the scattering problems in terms of correlators of the ordered exponentials written through external fields.
In our previous studies [60][61][62][63] we demonstrated that the Color Glass Condensate (CGC) and corresponding JIMLK approaches are interconnected with the effective action of Lipatov. In our further studies [7,8,63] we made a step towards the unification of the two approaches, where the Lipatov's action was reformulated in terms of interacting ordered exponentials, which are main degrees of freedom of Balitsky hierarchy of equations, [9][10][11][12][13]. Our ultimate goal is to unify the approaches of Lipatov and Balitsky to the description of the scattering at high-energies in the framework valid for both variants of the theory. In particular, we consider the derivation of Balitsky hierarchy for the quark (fundamental) Wilson lines from the action of Lipatov with the use of [7,8,63] results. The similar calculations were done in [97], there the comparison between the approaches was performed on the base of the diagrammatic calculations with the use of the vertices and Feynman rules obtained from the Lipatov's, [1][2][3][4][5][6], action. Nevertheless, in our calculations we consider an another way to derive the hierarchy of the correlators. This derivation we implement as follows. Concerning the Lipatov's effective action, we can derive the action beginning from the QCD action and introducing there the interacting exponentials as an additional part of the action, see discussion and details in [7,8], the exponentials provide a nonlocal rapidity dependence in the description of the scattering. At the next step, the product of the exponentials is eliminated from the action by the introduction of the auxiliary fields, the Reggeons. In the transformed action those exponentials are represented separately as effective covariant currents with Reggeons as their sources [60,61] and this is the original, Lipatov's one, action, see [1][2][3][4][5][6]. The whole rapidity interval of a scattering, consequently, is labeled and divided to the dif-ferent QCD clusters at some values of rapidity variable, the interaction between the clusters are due the Reggeons and this is a classical description of the high energy Balitsky-Fadin-Kuraev-Lipatov [98][99][100] approach to the scattering in the Regge kinematics. Nevertheless, we can perform the same calculations directly on the level of the complimentary variant of the Lipatov's action, [7,8], with only ordered exponentials included considering them as the objects of interests. In this case, applying the light-cone gauge, the gluon's propagator will not affected by the additional terms in the QCD Lagrangian and we obtain a theory of interacting Wilson lines averaged with respect to the gluon and quarks fields for the each separated QCD cluster.
This formulation of the problem has a few advantages in comparison to the standard derivation of the correlators of Wilson lines. First of all, we deal with the usual QCD Lagrangian, the introduction of the different gauges and corresponding calculations of the propagators in this case is a standard and well defined procedure. An important consequence of this is that it allows systematically define and calculate any kind of the many loops corrections to the correlators in correspondence to any gauge chosen in the initial action. An another advantage of the approach is that we obtain the equations for the different correlators directly from the action, without the approximation which concerns the external fields. In turn, the form and structure of the external fields can be borrowed from the complimentary Lipatov's action written in terms of reggeized gluons. The external fields there are reggeized gluons as sources of the effective currents, see details in [60,61] and this is the way how the reggeized gluon fields are entering in the hierarchy of equations of Balitsky. The shock wave approximation to the fields, required for the derivation of the original BK equation, is applied further in this case, the initial equations for the correlator are correct beyond the shock wave approximation as mentioned above.
Another important issue discussed in the paper is the origin and the form of ordered exponentials in the Lipatov's action. In the papers [7,8,101,102] it was demonstrated how exponentials are arising in the effective QCD action. Following to the [9][10][11][12][13], we can consider this additional part of the action as arising from a truncated propagator of the energetic particle in an external field, see Sect. 2. The additional corrections to the ordered exponential in this case, both sub-leading eikonal and sub-leading non-eikonal, arise in the definition of the effective currents. The usual Wilson line is a leading order (LO) value of such object but in general we have to distinguish the two types of the Wilson lines, for the quarks (fundamental) and for the gluons (adjoint). In the action the currents appear after the color averaging, the form of the currents is not important consequently to the LO precision and it does not affects on the further calculations. Nevertheless, if we account the contribution from the truncated propagator as the operator in the effective action, the difference is important. The propagators of quark and gluon in an external field are not the same, the sub-leading and non-eikonal corrections they provide are not the same as well. Therefore, for mostly general action, we need to introduce a two forms of the Wilson lines directly in the Lipatov's action and in it's modified form. Still, in the present article we derive the Balitsky hierarchy only for the quark Wilson lines, the same hierarchy for the adjoint Wilson lines will be different only if we will account sub-leading and non-eikonal corrections to the obtained system of the equations.
The paper is written as follows. In Sect. 2 we discuss the appearance and role of the operators related to the ordered exponentials in the action. In Sect. 3 we consider the effective action of Lipatov with those operators included in the Lagrangian and the derivation of general BK-like equation for the correlators of Wilson lines. In Sect. 4 we discuss the gluon's propagator requested for the calculation of the correlator whereas Sect. 5 is dedicated to the application of the shock wave approximation to the correlator's equation and derivation of the original form of the Balitsky's equation. The next two sections, Sects. 6 and 7, are about the calculations of the some sub-leading eikonal corrections to the BK equation where we derive these corrections for demonstrating the method of our calculations. The last Section summarizes our paper where we also discuss the possible applications of the framework.

S-matrix element of quark's propagation in an external gluon field
In order to understand an appearance of the Wilson line operators in the high energy scattering amplitudes let us consider a propagation of quark in the external field, we do not consider here the same problem with gluons included. The task, of course, is well known, see [103][104][105][106][107][108][109][110][111][112][113][114][115] references, additional examples can be found in [9][10][11][12][13][116][117][118] as well. Introducing J i and J f as sources of quark's creation and absorption, we define the S-matrix element for the general case of quark's propagation in an external field: with the processes of creation and absorption of the quark in y and x correspondingly and S(x, y) as quark's propagator in the external field. For example, for the quark asymptotically free at x, y → ± ∞ we have: Integration by parts provides in turn: with the limits p 2 , q 2 → m 2 assumed to be taken at the end. Further we consider the processes with the quark propagation in an external field, the form of the propagator is complicated in comparison to the propagators of scalar or vector particles. Therefore, in order to make the calculations easier, we will follow to [119] and introduce the two component spinors: where correspondingly with the chiral basis for the gamma matrices used: Similarly we can consider a right-handed spinor: with In any case, writing the Eq. (3) quarks propagator in the two component form we will stay with only G L or G R component in the reduced Eq. (3) expression for the case of the propagation of a quark with chirality preservation. 1 So, let's rewrite the spinor in terms of R and defining we write for the T-matrix element which corresponds to the Eq. (3) S-matrix element: here G R 0 is a free propagator at the absence of the external field. In the case, more interesting for us, when the quark is created at some y and asymptotically free at x, we will obtain instead Eq. (11): with properly redefined J (q) from initial Eq. (1). Therefore, we need to find the expression for G R (x, y) component of S(x, y), we have for the four component propagator: with v as an external gluon field and a coupling constant included in it's definition. Correspondingly we obtain: Resolving this system we have for the G R propagator: where we made G R /m → G R substitution for the sake of simplicity, further we also will use the anti-hermitian representation of gluon field v → ı v. Therefore we rewrite Eq.
as 1 The change of chirality is provided by matrix element of the R,L . . . L ,R with more complicated propagator, we do not consider this case in the paper.
Finally, defining see Eq. (6) definition, we obtain for the propagator: Our next step, consequently, is to write the propagator in the path integral form and isolated from it the requested part which reproduces the Wilson line currents in the effective action of Lipatov. This task is important because whereas the presence of the simple ordered exponentials in the Lipatov's effective action is explained in the [7,8] paper, we face a problem to account the sub-leading and non-eikonal correction there. Now, because the effective currents in the action can be interpreted as a contribution from the truncated propagators, using the precise form of the truncated propagator we consequently can account the sub-leading and non-eikonal corrections to the processes described by initial Lipatov's action.
The different presentation of the quark's propagator in an external field in the form of the path integral can be found in [103][104][105][106][107][110][111][112][113][114][115][116][117], see also [120][121][122][123][124] for examples, partially we will use these results. We begin from the standard representation of the Eq. (19) propagator: Assuming parametrization or the equivalent and extracting the bare propagator from the expression, i.e. leaving in the expression only part of the propagator which depends on the external field, we obtain:Ĝ here the substitution p =ẋ straight was made. This substitution can be understood as a straight line approximation for the description of the high energy quark's propagation in an external field. Namely, the dx μ integral in the Eq. (23) exponential must be defined as an integral on the the straight line quark trajectory, see Sect. 7 further. Considering the propagation of the quark free at x, see Eq. (12), we need an expression with the ( p 2 − m 2 ) −1 pole isolated, see [115] for additional examples. Proceeding with Eq. (23), we restore the regularization of the Eq. (17) operator adding ı ε term in the inverse propagator, and write for the factor which determines the answer for the Eq. (12): After the integration by parts we obtain: Now we can take p 2 → m 2 limit in the exponential in the integral obtaining finally see also [115]. We consider this "phase" operator is a main object of our interests, after the proper normalization it coincides to LO precision with the simplest Wilson line operator. The Eq. (27) expression can be directly used in the effective action of Lipatov, see further. From the expression it is clear, that the form of the same object for the gluon propagation in the external field will coincide with the Eq. (27) only to the LO precision providing different sub-leading and non-eikonal corrections. A different answer we will obtain if we will interesting in the T-matrix element for the quark free at x i and x f , in this case we need to isolate two poles ( p 2 − m 2 ) −1 and (q 2 − m 2 ) −1 , see Eq. (11) expression. We again begin from Eq. (23) and expanding P-exponential rewrite it as following: with correspondingly defined φ function. At the next step we perform variables change in Eqs. (3) and (28): Therefore we will obtain for Eq. (28) expression: with properly redefined argument of φ function. Redefining the variables in the integral we correspondingly can write: and proceeding similarly to done in Eq. (24) we again will obtain an expression for the T-matrix element, see [103][104][105][106][107][120][121][122].

Perturbative calculation of Wilson lines correlators
Our next step is a consideration of the variant of the Lipatov's effective action written in the form of interacting Wilson lines: here S 0 as the QCD gluon's action on which the averaging of the operators is performed, C(R) is an eigenvalue of the Casimir operator for the corresponding gluons representation and This version of the Lipatov's action was derived in [7,8,63], it's form is different from the initial Lipatov's action but, as was demonstrated in [7,8], it is fully equivalent to it. The use of this variant of the action is dictated by the presence of the Wilson lines operators in this variant of the action, they are main degrees of freedom in the BK-JIMWLK approaches as well. Therefore, it allows relatively simple comparison between the approaches, in all of them we will operate with the same type of the degrees of freedom. Now, introducing auxiliary currents and applying the light-cone gauge, we obtain: In the following we consider the Regge like kinematics of the interactions, i.e. shock wave approximation of the scattering. Therefore, we define the currents as following: Taking two derivatives of log Z [J ] with respect to the currents we have: further, in sake of simplicity, we will not write the integrals on x − , y − variables precisely and will account the integrals only in the calculations of the final contributions to BK equation. The Eq. (38) correlator is defined through the usual Wilson lines: We note, that the form of the effective action, Eq. (33), allows to rewrite the T ± (v ± ) operators in terms of Eq. (27) phase factor, it can be included directly in the Eq. (33) action. In this case, the action allows systematically account sub-leading and sub-eikonal correction to the Eq. (38) correlator. We have to mention also, that the [9][10][11][12][13] approach operates with the effective current in the action given by the quark-antiquark (quark) loop, i.e. by the physical currents which include full quark's propagators, consequently in the expression for the amplitude the corresponding impact-factor is arising. In contrast to that, the effective currents in the Lipatov's effective action are untended to reconstruct the correct BFKL correlators of reggeized gluons with proper rapidity dependence, see [60,61]. Therefore, whereas in the Balitsky's approach we need to know the form of the whole quark's propagator, in the form Eq. (20) for example, in the Lipatov's approach we need a "phase" operator in the form of Eq. (27), with the properly defined normalization of the path integral on ξ variable. This issue we also discuss in the Sect. 7. Now, taking the gluon field as the classical solution (Reggeon field) of the equations of motion of the effective action of Lipatov, see [60,61], plus the quantum fluctuation around we can expand the V operator in respect to these fluctuations: Here we defined with D as covariant derivative operator. Finally we can define the correlator of interest writing precisely the LO terms only in the expression: The same we can write in the matrix form: This system of equations is similar (almost) to the Balitsky, [9][10][11][12][13], hierarchy of the correlators. There is a complete system of equations for the Wilson line correlators written in terms of the correlators itself and correlators of gluon and quark 2 fields. Nevertheless, we do not assume nor any special forms of gluon's correlators in Eq. (46) neither any kind of large N c approximations, in general the system can not be be reduced directly to the standard BK equation. For example, the O T x T c O x ik operator can not be written through Eq. (38) Wilson lines which are present in the l.h.s. of the equation in the general case. Therefore, further approximations must be done in the framework if we want to reproduce the standard Balitsky equation. But, before to proceed with the shock wave approximation, 3 in the next Section we discuss the form of the gluon's propagator which appears in the r.h.s. of Eq. (46).

Propagator of gluon field
In this Section we remind the framework of [61,64] and calculate the gluon field propagator in the shock wave which we will use in the derivation of BK equation further. First of all, we consider an expansion of the gluon field around classical solution provided by Lipatov's effective action in the light-cone gauge: v a here B ± (x ± , x ⊥ ) are the Reggeon fields, see [60][61][62][63][64] papers. Expanding the QCD Lagrangian with respect to the fluctuations, we obtain the following expression: Here we defined (M i ) ac μ ν ∝ g i and note that The bare gluon propagator G 0 ν μ is defined as and correspondingly the full gluon propagator we write as following: The Eq. (52) also can be written in the form of the following perturbative series: To the first order of the approximation we can take v b cl i = 0. Therefore, the correlator (propagator) we need has the following form: Here Now we reproduce the Balitsky-Belitsky result, see [125], performing the full re-summation of Eq. (54) expression using the v cl + = B + classical field. Firstly consider the first order contribution of the classical field to Eq. (54) propagator: The second order contribution has the following form: In the Eikonal approximation we have and correspondingly obtain: Writing the integral dp 1 + 2π we will account further the θ( p 1 − ) contribution only for the sake of simplicity. Now, performing integrating on p 1 , we obtain finally for Eq. (63): The third order contribution is the following one: Using again Eq. (64) we obtain: The value of p, k integrals in the expressions above are defined by poles that allows to leading order to take at s → ∞: We note, that this approximation is equivalent to the shock wave representation for the classical field when it takes as v + (x + , x ⊥ ) ∝ δ(x + ), namely it will lead to the same result if we will substitute the classical fields in the form of the shock wave in the expression and will perform a resummation of the perturbative terms. After the similar calculations of all order terms in Eq. (54) expression, we will obtain for the propagator of interests: with for the adjoint representation of the gluon field. This is a gluon propagator in an external field obtained in [125], further we identify to LO precision v cl + = B + . Following to [125], we also note that due the δ( p − − k − ) in the expression, the integration on p + and k + variables provides θ(x + ) θ(−y + ) and θ(−x + ) θ(y + ) in the answer, we will use this result in the derivation of BK equation.

BK equation
In order to obtain the familiar form of BK equation from Eq. (46) expression, we consider the given terms of the Eq. (46) hierarchy and apply the shock wave approximation there. It means the following, taking any operator of interest in Eq. (46), we have to define the expressions for these operators whereas the classical fields in the operators are distributed around δ(x + ) only. Therefore, we define: and Correspondingly, for the last term in Eq. (46) for example, we obtain: with the regularization of the integral around x + = 0 applied, see also Eq. (70) definition of the propagator which determines the separation of the integration regions. Now, inserting the Eq. (70) expression for the propagator and accounting that we obtain: here we used and considered only θ( p − ) part of the propagator. An another contribution to BK equation we consider corresponds to the self-energy of the Wilson lines. Considering, for example, the fourth term in Eq. (46) we have there: Using Eq. (42) definition of G + (x, w) Green's function, we can write for the case of shock wave: and performing subsequent integration over x + and w + , we will obtain for the Eq. (80) expression: Given the fifth term in Eq. (45), which is symmetrical with respect to the permutation of x + and w + variables, we obtain the following answer for these two terms and for two additional terms with x → y substitution: Therefore, we obtain finally for the correlator of interests: the expression is BK equation in Balitsky's formulation, see [9][10][11][12][13].

Transverse fields sub-leading contributions to the gluon propagator
The simplest corrections to the Balitsky hierarchy which we can account in the present framework it is a correction to the gluon's propagator in the external field. Namely, considering the Eqs. (53)-(55) expressions for the propagator, we note that there are corrections to the propagator arise form the transverse fields contributions, see also Eq. (49). In the Lipatov's approach these fields are not independent, due the non-locality of the amplitudes the only fields of interests are Reggeons, the transverse fields are expressed through them, see [60,61]. Therefore, proceeding with the Eq. (53) expression for the full propagator in the external field, we able to account the LO contribution of the transverse field into the propagator. In turn, the corrections will be expressed in terms of the longitudinal external field preserving the selfconsistency of the hierarchy. Now, first of all, we consider the contribution to the propagators arising from the M i j part of Lagrangian proportional to the transverse fields: this vertex also contributes in G ++ and G i+ propagators. Nevertheless, we put an attention that for the part of G ++ propagator in the external field we have the contribution proportional to after the expression for the bare propagator Eq. (58) is accounted. For the G i+ propagator the similar contributions gives which does not give zero in Eq. (55) expression in general. Nevertheless, we know that the classical transverse field has the following structure: is the second Reggeon field. We see, that to this precision we can write this field in Eq. (88) and correspondingly in Eq. (87) as that also provides zero contribution for the Eq. (87) vertex in the external field part of Eq. (55) propagator. Therefore, we conclude, that the transverse field part of M i j vertex does not contribute to our equation to this precision order An another contribution we have to discuss is contribution from M −i vertex to the corresponding propagators. Whereas the system of Eqs. (54)-(55) equations for G ++ and G i+ full propagators can be resolved and written precisely in the form of infinite perturbative series, we will use the truncated full solution for G ++ propagator. In shortened notations it reads as: This expression is reduced to the previous one by taking M −i = 0, the B − field here is the smallness parameter, it's smallness is related to the asymmetry of the scattering process. Our next observation is that the contribution to the propagator's external field part is coming from p + = 0 momentum, therefore we can neglect all G 0 ++ terms in Eq. (91), see Eq. (54) expression. Therefore, we obtain: The last term in the Eq. (91) expression determines the additional contributions to the BK equation provided by the classical v i gluon fields. Correspondingly, repeating the Eq. (77) calculations, we obtain for the momentum dependent contribution to Eq. (92): here, for the sake of simplicity, we preserved only the first perturbative contribution from the expression of the gluon's propagator. We also use Eq. (88) expression for the transverse field at v cl + = 0 , the color structure of the expression stays the same as in LO BK equation. Performing the simple integrations we will obtain: that in turn gives: We take into account now that for ∂ − v cl c 1 d j (z − 1 , z 1 ⊥ ) term's contribution we obtain: and, therefore, we have: Now, we obtain the following term additionally to Eq. (77) answer: The Eq. (99) expression, in turn, depends on the form of ∂ −1 − operator which appears in Eq. (90). Namely, taking for example we will obtain for Eq. (99) expression with the use of LO value of the transverse field from Eqs. (88)- (90): here the integral's value depends on the parity of the physical source (impact factor) of B − Reggeon field. In the case of the shock wave approximation when B − ∝ δ(w − ) the Eq. (101) is zero of course, there the non-zero contribution will appear only from the correction to the δ(w − ) distribution of the classical external field. We also note, that for the shock wave, the Eq. (88) expression can be rewritten as which provides, for example, additional correction for the ∂ −1 − operator different from Eq. (100).
Performing the similar calculations with respect to Eq. (82) contributions to BK equation, we obtain: Therefore, the NLO sub-leading contributions to the BK equation arise from the classical transverse field contribution have the following form: We see, that this correction is of the eikonal type and its contribution crucially depend on the form of Eq. (100) operator. Taking it symmetrical for the case of the shock wave approximation, the contribution will not be zero only when the corrections to the δ(x ± ) shock wave form of the classical fields will be accounted. In this extend we can consider these corrections as corrections beyond shock wave approximations.

Sub-leading eikonal corrections to the quark propagator
As mentioned above, we can consider the Wilson line in Eq. (3) Lipatov's action as first order approximation to the Eq. (27) "phase" operator Eq. (27), see for example also [110][111][112][113][114][115][116][117]. This is important observation, because it allows to extend the initial formulation of the effective action beyond the eikonal approximation and additionally allows to account the sub-leading eikonal corrections to the ordinary Wilson line. An another important lesson is that going to account these corrections, we need to write in the action separately two covariant currents with Reggeons as sources: the first one for the fundamental Wilson line and the second one for the adjoint. As discussed, the reason for that is the coincidence of the forms of these Wilson line only to LO precision, in general they are different. We consider further the case of the fundamental Wilson lines, for the adjoint ordered exponential the calculations can be done similarly of course, with only change of the Eq. (27) "phase" operator on the gluon's one. Therefore, going beyond the eikonal approximation, we can use the following expression instead Eq. (39) form of the operator: with Eq. (22) as a parametrization of the trajectory of the quark 4 and the following normalization of the ξ path integral 4 For an arbitrary kinematic and configuration of the gluon fields the equations of motion for the x can be solved firstly that determines a calculation scheme for the path integral in respect to the fluctuations ξ .
in Eq. (105): This normalization, as we mentioned above, corresponds to the definition of the truncation of the quark's propagator required in the Lipatov's formalism. In our case, when we talk about Regge kinematics, we have: and it provides: with the following definitions of the quark's trajectory defined: The variable ξ in the expressions is a fluctuation of the trajectory around the straight line, with F μν as the gluon field strength tensor we have: The Eq.  [116,117,[126][127][128][129].
Further we do not consider these types of the non-eikonal contributions taking this factor equal to one, the calculations of it's contribution to the BK equation we will perform in an additional publication.
In our case, therefore, we have the following expression for the eikonal ordered exponential with sub-leading corrections included: Taking into account the contribution of v + gluon field in Eq. (111) expression and adjusting Eq. (109) definitions we obtain correspondingly: Expanding the Wilson line with respect to the fluctuations around the classical Reggeon, see Eq. (40), we will obtain instead Eq. (41): definition used. Now, correspondingly, we obtain instead Eq. : We note, that due the δ( p − − k − ) function present in the gluon's propagator, the contribution which arises from the action of the ∂ x − + ∂ y − operator on the propagator is zero. The same we obtain for the action of ∂ x i + ∂ y i operator on the terms which represents the Wilson's line self-energy leading order contribution, there δ 2 (x ⊥ − w ⊥ ) function is present. Therefore, we stay with the only following additional contribution to BK equation: which we finally write in the full form as This sub-leading correction is as well of the eikonal type, it describes the quark's spin contribution to BK equation in the processes without helicity flip.

Conclusion
In this paper we demonstrated the interconnection between the Lipatov's QCD effective action and Balitsky description of the correlators of ordered exponentials of classical external fields. The main result of the paper is that we show how the correlators can be obtained directly from the effective action of Lipatov.
There are the following important issues that we addressed in the course of the derivation. The effective currents in the action of Lipatov and correlators of the approach of Balitsky are energetic quark's (or gluon's) propagators in an external field. These propagators, full or truncated, in the leading order approximation can be written in terms of ordered exponentials. Therefore, there are two important issues we address: (a) in which form and (b) at what precision one must account for the propagators. As it was already mentioned in [119], the usual four dimensional spinor form of the quark's propagator is not convenient to use, consequently we consider twodimensional form of the propagator which is easier to use in the processes without flip of the quark helicity. This simplification of the propagator's form makes the calculations simpler, in turn, the form of the propagators clarifies the corrections to the BK equation from both the spin's correction to the propagator and from the deviation of the trajectory of the quarks from the straight line, see Eqs. (105) and (108). Once both corrections are written in the single expression one can develop a perturbative expansion of the propagators with respect to the fluctuations of the trajectory about the straight line. Therefore, taking the effective current in the Lipatov's action in the form of Eq. (105), we introduce these type of the corrections in our approach in a straightforward manner, that extend the area of the application of the Lipatov's effective action in both original, [1][2][3][4][5][6], and new, [7,8], forms. Consequently, there is a new and interesting task to formulate the original version of the action, [1][2][3][4][5][6], in terms of the Reggeon fields with included in it the Eq. (108) operator. This action's generalization will provide the sub-leading and non-eikonal corrections to the amplitudes of the processes with Reggeon fields involved, see also [64]. We also conclude, that considering gluon and quark propagations both together in the action, we will need to add in the effective action, additionally to the Eq. (108) expression, the similar operator arises from the propagation of the gluon in the external field, we plan to consider this task in an additional publication. The processes with the helicity flip is an another example where this extended action can be used. Namely, these processes can be considered similarly to done with the only difference that we will need to use there the propagators corresponding to the flips which can be introduced similarly to discussed in Sect. 2.
It is worth emphasizing that the Lipatov's action does not originally imply the shock wave form of the external fields. The action is written for the Reggeon degrees of freedom, which interaction is averaged with respect to the gluon and quark fields in some rapidity interval. In this case, the whole amplitude of any scattering process can be written as nonlocal interactions of these clusters by the reggeized gluons, the interactions of the Reggeons with the target and projectile at the edges of the interval determine the form of the Reggeon fields. In this extend, we obtained a precise system of equations for the ordered exponentials, similar to the [9-13] hierarchy, directly from the action of Lipatov. The new hierarchy is precise in the sense that it includes the ordered exponentials written in terms of general external fields beyond the shock wave approximation. The Eq. (46) equation can be written fully in terms of O and O T operators, nevertheless it is not clear if we can resolve the equation in terms of the operators in this case. Therefore, in order to obtain the familiar form of the BK equation it is crucial to assume the shock wave form for the external fields and after that one can write the equation in terms of the infinite Wilson line only. This different approach to the problem perhaps allows to introduce the corrections to the equation on some clear base. Namely, in the proposed approach we firstly write the equation and only after that we can expand any object in the equation with respect to any parameters of the problem: it can be finite range of the external field, corrections to the gluon's propagator, non-straight line of the particle's trajectory, all possible corrections together and etc.
The central result of the proposed framework is that it provides a useful instrument for the further perturbative calculations related to the BK equation. First of all, the effective action allows to formulate the problem for the different QCD gauges directly in the action and, in turn, allows to determine corresponding redefinition of the propagators (correlators) of interest which appear in r.h.s. of Eq. (46). Additionally, it clarifies the translation vocabulary between the different approaches. For example, the variation of the action with respect to the auxiliary currents, Eq. (38), correspond to the correlators of the reggeized gluons fields in the complimentary theory written in terms of the Reggeons fields, there the reggeized gluons are sources of the currents, see [64] for example. From that it is clear, that the Eq. (41) expansion of the Wilson line is in correspondence with the trajectory of the reggeized gluon, the correlator of two Wilson lines, Eq. (45), is in correspondence to the BFKL Pomeron, the correlators of the higher orders are in correspondence with the correlators written in terms of the B + reggeized gluon fields. It is interesting to note, as it follows from this relation, that the calculation of the B + B + correlator in the original theory of the effective action, [1][2][3][4][5][6][60][61][62][63][64] with B − = 0, B + = 0 taken at the end and corresponding calculation of B − B − with B − = 0, B + = 0 at the end, will determine the forms of both triple Pomeron vertices in the light-cone gauge, see effective theory for the Pomerons interactions with these vertices in [130][131][132][133][134]. It is important to mention, that because the [1][2][3][4][5][6] effective action can be written in any gauge, the corresponding calculations can be done symmetrically with respect to target and projectile, i.e. symmetrically with respect to Reggeon fields B − and B + or symmetrically with respect to Wilson lines which depend on these two fields. This property can be useful if we consider the processes where the target and projectile are similar. Moreover, concerning the many loops perturbative calculations of the kernel of BK equation, it can be easier, perhaps, to use in the calculations a gauge different from the light-cone. As a result, the new gauge scheme can be introduced directly in the action and the many-loops calculations of the BK kernel further can be performed as defined in the paper, this task we will consider in our further publications.
The proposed framework allows to calculate also the subleading corrections to the BK equation discussed in the paper. First of all, those are corrections to the gluon's propagator which arise from the contribution to it from the classical transverse field. From the Lipatov's formulation of the high energy scattering we know that the action can be written in terms of the physical degrees of freedom which are longitudinal reggeized gluons. Therefore, our local in rapidity action can be fully written in terms of the Reggeon, i.e. the action can reformulated as RFT Regge Field Theory (RFT) action. Practically it means, that the transverse fields are determined in terms of the longitudinal ones, we have only two degrees of freedom after all. The gluon's propagator in this case gains a correction from the transverse field, which, in turn, appears as correction from the longitudinal B − reggeon field, see Eqs. (99)- (101). This correction is definitely sub-leading at the shock wave regime and in the case of a non-symmetrical scattering, i.e. when the B + field is enhanced in comparison to B − due their coupling to the different impact-factors. Nevertheless, the correction can be important when the projectile and target are symmetrical or if we want to account the corrections to BK beyond shock wave approximation, this is an additional interesting issue for the further investigation of the approach. An another correction we discussed is the quark's spin contribution to the BK equation, we consider the processes without helicity flip, similar and more general description of the similar processes can be found in [135][136][137][138][139][140] for example. We note, that in general this process is incorporated in the framework, but of course, when we will need to consider the helicity flip of the quark, then we will need to redefine the quark's propagator in accordance with the Sect. 2 prescription. In summary, we consider the proposed approach a very promising framework for many other calculations in the BFKL and BK physics of high energy scattering.
The authors are indebted to Ian Balitsky for the fruitful discussions and the helpful advises in the course of the preparation of the current paper.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no data to present, the article is a theoretical one.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .