From Heisenberg-Euler Lagrangian to the discovery of Chromomagnetic Gluon Condensation ∗

I reexamine the phenomena of the chromomagnetic gluon condensation in Yang-Mills theory and perform the exact integration of the Heisenberg-Euler Lagrangian in the limit of massless chiral fermions. The extension of the Heisenberg-Euler Lagrangian to the Yang-Mills theory allows to calculate the corresponding effective action, the energy-momentum tensor and demonstrate that the energy density curve crosses the zero energy level of the perturbative vacuum state at nonzero angle and continuously enters to the negative energy density region. At the crossing point and further down the effective coupling constant is small and proves that the true vacuum state of the Yang-Mills theory is below the perturbative vacuum state and is described by the nonzero chromomagnetic gluon condensate. The renormalisation group analyses allows to express the energy momentum tensor, its trace and the vacuum magnetic permeabilities in QED and QCD in terms of effective coupling constant and Callan-Symanzik beta function. These considerations were provided earlier to justify the existence of the chromomagnetic gluon condensation and of the nonzero energy gap between perturbative and non-perturbative vacuum states. We shall present the work of the Niels Bohr Institute theory group and of Curt Flory proving the absence of the imaginary part in a chromomagnetic field. In the vacuum the energy-momentum tensor is proportional to the space-time metric, and it induces a negative contribution to the effective cosmological constant. ∗ Based on lectures at the Leipzig University in occasion of the 80 Years of Heisenberg-Euler Lagrangian 19362016, ITP, Leipzig, November 21, 2016 and 40 Years of Discovery of the Chromomagnetic Gluon Condensation 1977-2017 at the Ludwig-Maximilian University München, Arnold Sommerfeld Colloquium at Center for Theoretical Physics, April 18, 2018. ar X iv :1 91 0. 00 65 4v 1 [ he pth ] 1 O ct 2 01 9


Introduction
In this article we shall analyse the effective action in QED and QCD by using the perturbative loop expansion and renormalisation group equations and discuss the physical consequences which can be derived from their explicit expressions. We shall reexamine the proof of the existence of the chromomagnetic gluon condensation in Yang-Mills (YM) theory and will present the derivation of the new results. The Heisenberg-Euler Lagrangian in QED [1,2,3,4,5,6] is a sum of the one loop diagrams with a vacuum electron-positron pair circulating in the loop and the gluons and quarks in case of QCD [7,8,9,10,11,12,13,14]. The effective action Γ[A] has the following representation: where L is the effective Lagrangian, Γ (n) is a one-particle irreducible (1PI) vertex function, A a µ (x) ≡< 0|A a µ (x)|0 > is the vacuum expectation value of the field operator and W (n) , n = 1, 2, .. represent the terms of the loop expansion.
We shall consider the limit of massless electrons and quarks and demonstrate that the proper time integral in the Heisenberg-Euler Lagrangian can be calculated explicitly by using covariant renormalisation condition [11,13,14] ∂L ∂F | where F = 1 4 G a µν G a µν is the Lorentz and gauge invariant form of the YM field strength tensor G a µν and µ 2 is the renormalisation scale parameter. In the massless limit the QED effective Lagrangian has the exact logarithmic dependence as a function of the invariant F (see Fig.1): where H and E are magnetic and electric fields. This expression should be compared with the one-loop effective Lagrangian in pure SU(N) gauge field theory, which has the form [11,13] (see From (1.2) it follows that the corresponding quark contribution considered in the chiral limit is Figure 1: The graph shows the qualitative behaviour of the QED vacuum energy density ( H 2 ) (1.2), (2.43) and of the effective coupling constantē 2 ( H 2 ) (4.58) as the functions of the background magnetic field. The effective coupling constant is singular at H 0 2 , the "Moscow zero" [83,84].
where N f is the number of quark flavours.
The effective Lagrangian technique allows to calculate the magnetic induction B of the vacuum defined through the derivative of the effective Lagrangian [11]: (1.5) From (1.2), (1.3) and (1.4) it follows that in QED the vacuum responds to the background magnetic field as diamagnet and in QCD as paramagnet with the magnetic permeabilities of the following form [11]: The diamagnetism of the QED vacuum (1.6) means that it repels the magnetic fields by forming induced magnetic field in the direction opposite to that of the applied magnetic field. This phenomenon is similar to the Landau orbital diamagnetism of free electron gas when the counteracting field is formed when the electron trajectories are curved due to the Lorentz force [80]. The paramagnetism of the QCD vacuum (1.7) means that it amplifies the applied chromomagnetic field by generating induced chromomagnetic field in the direction of the applied field. In QCD the large polarisation of the gluon spins is responsible for the amplification of the background field. This phenomenon is similar to the Pauli paramagnetism, an effect associated with the polarisation of the electron spins [79].
The effective Lagrangian approach allows to calculate the quantum-mechanical corrections to the energy momentum tensor by using the formula derived by Schwinger in [5]: . (1.8) In case of the Heisenberg-Euler effective Lagrangian Schwinger presented the expression for the T µν in the fine structure constant α = e 2 /4π expansion: (1.9) with its nonzero trace (1. 10) In massless QED using the one-loop expression (1.2) for T µν one can get The T µν becomes proportional to the space-time metric tensor g µν at the extreme magnetic field H 2 0 = H 2 c exp (6π/α) and therefore induces a positive effective cosmological constant (see Fig.1). To calculate the energy momentum tensor T µν in pure SU (N ) YM theory one should use the expression (1.3) and in the case of QCD, in the limit of chiral fermions, one should also add the quark contribution (1.4) by using the substitution 11N → b = 11N − 2N f : The vacuum energy density T 00 ≡ (F) has therefore the following form [13]: The energy density has its new minimum outside of the perturbative vacuum state < G 2 µν >= 0, at the Lorentz and renormalisation group invariant field strength [13] where b = 11N − 2N f and characterises the dynamical breaking of scaling invariance in YM theory † : Substituting the vacuum field intensity (1.14) into the expression for the energy momentum tensor (1.12) one can get that in the vacuum the tensor T µν is proportional to the space-time metric g µν : (1. 15) In this form the energy momentum tensor represents the relativistically invariant equation of state vac = −P vac , which uniquely characterises the vacuum [15,16] with its negative energy density vac . The vacuum energy momentum tensor (1.15) generates the effective cosmological where the chromomagnetic condensate (1.14) is < 2g 2 F > vac = Λ 4 QCD . The magnetic permeability (1.7) in the vacuum state (1.14) is equal to zero: It is useful to derive the expression of the effective Lagrangian by using the renormalisation group equation [13,14]. The solution of the renormalisation group equation in terms of effective coupling constantḡ(g, t), with the boundary conditionḡ(g, 0) = g, has the following form [13,14]: The derivative (1.18) of the effective Lagrangian has transparent expression in terms of the effective coupling constant and allows to obtain the effective Lagrangian by integration over F in all order of the perturbative expansion: and find out the expressions for the physical quantities beyond the one-loop approximation.
One can calculate different observables of physical interest that will include the effective energy momentum tensor, vacuum energy density, the magnetic permeability, the effective coupling constants and their behaviour as a function of the external fields. In particular, the energy momentum tensor (1.8) will take the following form: where the trace of the energy momentum tensor T µµ is given by the following expression: The last formula provides all-loop expression for the conformal anomaly in gauge field theories ‡ .
As far as the beta functionβ(g) has no zeros, is negative analytical function of the coupling constant and Considering the value of the field strength tensor F 0 at which the vacuum energy density (1.13) vanishes (F 0 ) = 0, the point F 0 shown on Fig.2, one can observe that the effective coupling constant (1.18) at this field strength has the value 96π 2 11N −2N f and tends to zero as N → ∞. The energy density curve (F) (1.13) intersect the horizontal zero energy line at the nonzero angle θ > 0 (see (6.110) and Fig.2). The energy density curve can be continuously extended from the point F 0 deep into the negative energy density region arbitrary close to the value of the vacuum condensate < F > vac by considering a larger values of N and keeping the t'Hooft coupling constant g 2 N small and fixed. This proves that the true vacuum of the Yang-Mills theory is below the perturbative vacuum, the energy density curve (F) is approaching the point of zero energy density from below and that there exists a nonzero energy gap between perturbative and non-perturbative vacuum. The details of the proof are presented in the Section 6.
The article is organised as follows. In the second section we shall use gauge and renormalisation group invariant scheme (1.1) [11,13] to renormalise the massless Heisenberg-Euler Lagrangian and derive the exact one-loop expression for the effective Lagrangian in QED (1.2). In the third section we shall use the renormalisation group equations for the effective Lagrangian to derive all loops results for the vacuum energy density and the traces of the energy momentum tensor. In ‡ If one considers the approximation in whichḡ(t) is field independentḡ(t) ≡ g then (1.22) will reduce to the one given in literature [52,53,54,55].  (4.71). At the intersection point F 0 (6.108) the effective coupling constant is small (6.109) and the intersection angle θ is strictly positive (6.110). The energy density curve can be continuously extended from the point F 0 deep into the negative energy density region arbitrary close to the value of the vacuum condensate < F > vac by considering a larger values of N and keeping the t'Hooft coupling constant g 2 N fixed (6.111), (6.112). This proves that the true YM vacuum is below the perturbative vacuum and that there is a nonzero energy gap gap > 0 between perturbative and non-perturbative vacua. The true vacuum is characterised by the nonzero value of the chromomagnetic field strength tensor (1.14), (6.99) and the energy density gap gap = | vac | (1.16) [13,32]. the forth and fifth sections the analyses will be extended to the Yang-Mills theory and the proof of the formation of the chromomagnetic gluon condensation will be reexamined. We shall present the absence of the imaginary part in the YM effective Lagrangian in chromomagnetic field and the stability of the chromomagnetic gluon condensate by Niels Bohr theory group and by Kurt Flory.

Heisenberg-Euler Effective Lagrangian in Massless Limit
The effective action Γ and the effective Lagrangian L in gauge field theories can be represented as a sum of the one-particle irreducible loop diagrams W (n) , n = 1, 2... : (2.25) where S is the Yang-Mills or Maxwell action: We shall analyse the behaviour of the effective actions in both theories and shall consider first the Quantum Electrodynamic in massless limit ( see, in particular, [76]). The Heisenberg-Euler Lagrangian in QED [1,2,3,4,5,6,75] is a sum of the one-loop diagrams with a vacuum electronpositron pair running in the loop: (2.27) and can be expressed through the functional determinant of the Dirac operator [5]: where Π µ = −i∂ µ − eA µ . The general expression for the one-loop effective Lagrangian L (1) has therefore the following form: In the case of the constant electromagnetic field strength tensor F µν the matrix element of the operator (x|U (s)|x) can be calculated exactly and has the following form [5]: The Lagrangian will take the form By the rotation of the integration contour in complex plane s as s → −is one can get The traces in (2.30) can be evaluated by using the eigenvalues of the field strength tensor matrix The characteristic equation is (2.32) and has the solutions The Lagrangian (2.30) will take the following form [4,5]: and with real eigenvalues In pure magnetic field configurations G = 0, F > 0 one can get f 1 = (2F) 1/2 , f 2 = 0 and in pure We shall consider the QED in the massless limit and impose the following renormalisation condition on the effective Lagrangian introduced in [13, 14] § : where µ 2 is the renormalisation scale parameter. This condition defines the renormalisation of the effective Lagrangian in a covariant gauge L r = L un − Z F. In the case of pure magnetic field (2.35) the Lagrangian L (1) has the following form: , § This renormalisation scheme is alternative to the standard MS and other schemes, see, in particular, [32]. and diverges at the boundaries of the proper time integration region. With the use of the renormalisation condition (2.36) one can handle both divergences [13,14]. This leads to the following renormalisation of the Heisenberg-Euler Lagrangian in the massless limit: One can get convinced that this expression is well defined in both limits, in the ultraviolet s → 0 and in the infrared s → ∞ regions. One can calculate this integral exactly. The integrals appearing in this expression can be expressed in terms of the Riemann zeta function and its extension (see the Appendix for details). The Lagrangian (2.37) will take the following form: and in the limit k → −1 we shall get where we used the identity [67] lim k→−1 (k + 1)Γ(k)ζ(k) = lim k→−1 Thus, in terms of Lorentz and gauge invariant F = 1 4 F 2 µν , the exact expression of the one-loop Lagrangian in massless QED is: where 2F = H 2 − E 2 > 0, G = E H = 0 and the effective Lagrangian will take the following form: As it follows from this expression, the QED vacuum responds to the background magnetic field by inducing a vacuum current of the electron-positron pairs, which attenuates the magnetic field imposed on the vacuum. The magnetic induction B of the QED vacuum is [11]: 41) and the QED vacuum responds to the background magnetic field as a diamagnet with the magnetic permeability of the following form: The diamagnetism of the QED vacuum means that it repels the magnetic fields by forming induced magnetic field in the direction opposite to that of the applied magnetic field. This phenomenon is similar to the Landau diamagnetism of free electron gas when the counteracting field is formed when the electron trajectories are curved due to the Lorentz force. This also can be seen from the vacuum energy expression ( see Fig.1 ): In the case of pure electric field the one-loop Lagrangian has the following form: and has singularities at s = s n = πn/b = µ 2 πn/eE. The integration path is considered to lie above the real axis, therefore we shall obtain a large positive imaginary contribution to L (1) ¶ : The real part of the Lagrangian in the electric field is The formulas (2.40), (2.47) and (2.46) prove that the effective Lagrangian is the analytical function of the variable F and has the general form (2.40). The corresponding energy density takes the following form: and its behaviour is similar to the one shown on Fig.1. .
In the next section we shall consider the renormalisation group invariant derivation of the all-loop effective Lagrangian (2.25) and the generalised expressions for the magnetic induction (2.41) and permeability (2.42) as well as the electromagnetic energy-momentum tensor and its trace. ¶ The universal character of the electric instability of the vacuum was discussed in the recent article [77].
Let us derive the exact expression of the effective Lagrangian by using the renormalisation group equation [13,14]. The effective action Γ is renormalisation group invariant quantity: because the vertex functions and gauge fields transforms as follows: The renormalisation group equation takes the form where β(g) is the Callan-Symanzik beta function, the γ(g) is the anomalous dimension. When where in the covariant background gauge β = −gγ [11]. By introducing a dimensionless quantity and (2.36) plays the role of the boundary condition: From equations (3.51) and (3.53) it follows that The solution of the renormalisation group equation (3.51) in terms of effective coupling constant g(g, t), with the boundary conditionḡ(g, 0) = g, has the following form [13,14]: The behaviour of the effective Lagrangian at large fields is similar to the behaviour of the gauge theory at large momentum. It follows that M(g, t) is completely determined for all t in terms of its first derivative (3.54) at t = 0. To define the effective Lagrangian L one should perform additional integration, which we shall do in the next section.
The above results allow to obtain renormalisation group expressions for the physical quantities considered above in one-loop approximation. Indeed, with these expressions in hand we can calculate different observables of physical interest, that will include the effective energy momentum tensor, vacuum energy density, the magnetic permeability, the effective coupling constants and their behaviour as functions of the external fields.

Massless QED
By using the one loop expression (2.40) derived above one can calculate the derivative and the Callan-Symanzik beta function (3.54) takes the following form: (4.57) The effective coupling constant (3.55) in the one-loop approximation is and tends to infinity at the magnetic field In order to estimate the value of the critical field one can consider the mass parameter µ to be of the order of the electron mass m. Then one can get The perturbation expansion breaks down at the "Moscow zero" shown on Fig.1.
As far as the derivative (3.55) of the effective Lagrangian (3.50) has transparent expression in terms of the effective coupling constant (3.55) one can obtain the effective Lagrangian by integration over F: By using the relation (3.50) to express the differential g 2 dF = µ 4 e 2τ dτ through dτ one can represent the Lagrangian in the form: In massless QED the magnetic induction (2.41) will take the following form: Therefore the vacuum permeability (2.42) can be expressed through the effective coupling constant e 2 (t): . (4.63) The effective Lagrangian approach allows to calculate the quantum-mechanical corrections to the energy momentum tensor by using the formula derived by Schwinger in [5]: In our case, when G = 0, we shall find all-loop expression for T µν by using (4.61): And for the vacuum energy density we shall get: where the trace of the energy momentum tensor T µµ is not equal to zero and characterises the breaking of conformal symmetry in massless QED: where G = 0, t = 1 2 ln(2e 2 F/µ 4 ). It is also useful to obtain the derivative of T expressed in terms of the effective coupling constant: and by using (3.55) we shall get The integration of (4.68) over F provides the alternative forms of (4.67): . (4.69) The last two formulas (4.67) and (4.68) provide the all-loop expressions for the conformal anomaly in QED in the massless limit. If one considers the approximation in whichē(t) is field independent e(t) ≡ e then (4.67), (4.68) will reduce to the expression T = 2 β(g) g F given in literature [52,53,54]. For the one-loop energy momentum tensor (4.64) we shall get where we used the expressions (4.56). The energy density and the trace of the energy momentum tensor are and they represent the one-loop approximation of (4.66) and (4.67). In the next section we shall consider the behaviour of the effective Lagrangian in Yang-Mills theory and QCD.

Effective Lagrangian of Yang-Mills theory
The loop expansion of the effective action in Yang-Mills theory has the following form: and the one-loop effective Lagrangian has the form [7,8,9,10,11,12] 75) and for the effective Lagrangian the following expression: where U (s) = e −iH(α)s , U 0 (s) = e −iH F P s . The Green function in the background field has the following form: As is was proven in [11,12], the L ef f is α independent functional on the solutions of the YM classical equations. On the covariantly constant gauge field solution [9, 10, 12, 11] the matrix elements can be calculated and have the following form [12,11]: where the corresponding matrices are By substituting the matrix elements and calculating the traces one can get [11,12]: where The first integral here coincides, up to the coefficient 2, with the expression of the one-loop Lagrangian in the scalar electrodynamics. The doubling of this expression is associated with the additional degrees of freedom due to the vector bosons isospin. The second term is due to the spin contribution −2gG µν in the operator H µν . We introduced the mass parameter µ 2 in order to control the infrared singularities and to make the integrals convergent at infinity [11]. Still, this is not enough to make integrals convergent at infinity. By using the real eigenvalues one can observe that the second term in the square bracket will take the form sinh(gf 2 s) sin(gf 1 s) and the integral diverges exponentially in the infrared region at infinity. We shall choose the integration counter in the complex plane s so as to guarantee the convergence of the last integral. For that one should rotate the integration counter in the third integral by the substitution s → −is. The same rotation of the counter can be performed in the first integral as far it is convergent in any way. Thus we shall get [12,11] The integrals are still diverging in the ultraviolet region at the s = 0. In order to renormalise the Lagrangian we have to identify the ultraviolet divergences in the above integrals. These are (gf 1 s) (gf 2 s) sinh(gf 1 s) sin(gf 2 s) Subtracting these terms, which are quadratic in the field strength tensor, we shall get the renormalised effective Lagrangian [11]: Now the integrals are convergent in both regions, in the infrared and in the ultraviolet. First let us consider a pure chromomagnetic field: The Lagrangian (5.87) will take the form At strong chromomagnetic fields g 2 H 2 a ≥ µ 4 the asymptotic behaviour of the real part is [11] where the first term represents the diamagnetism, which counteracts to the external field caused by the quantum current induced by the charged vector bosons in the vacuum and the second term represents the paramagnetism, an effect associated with the polarisation of the gluon spins, which, as one can see, dominates the asymptotic behaviour [11]. The imaginary part of the effective Lagrangian (5.88) in background chromomagnetic field in our regularisation scheme vanishes [11]: where f 1 = H a 2 = H. A similar conclusion was derived in [46] by using alternative regularisation. The significance of the absence/presence of the imaginary part in the effective Lagrangian connected with the fact that it defines the quantum-mechanical stability of a given field configuration. The above conclusion on the absence of the imaginary part is not conclusive due to the tachyonic eigenmode in the operator (5.74) [30,31], and we shall discuss the stability of the background field configurations in the subsequent eighth section. Here we shall refer to the articles of Leutwyler [42,43] and Flory [44,45], where they come to the same conclusion that there is no imaginary part in the effective Lagrangian in chromomagnetic field. Leutwyler was considering the self-dual chromomagnetic background field configurations and demonstrated that there is no imaginary part in the effective Lagrangian [42,43] and that the real part of the effective Lagrangian has the form identical to (6.95). Flory included the quartic self-interaction of eigenmodes and also came to the same conclusion [48]. The physical reason behind this universality lies in the fact that even when the background field depends on space-time coordinates the part of the effective LagrangianL(F, G) which depends only on the field strength tensor but not of its covariant derivatives has a universal form [11,13,14]. In other words, as far as the wavelength of the fluctuating fields is very long, the effective action is not sensitive to the fine structure of the fluctuating fields.
Let us now consider pure chromoelectric fields G = 0, F < 0 and f 2 1 = 0, f 2 2 = −2F: The Lagrangian has singularities on the real axis at s n = nπ/eE, and the integration path is considered to lie above the real axis [12,11]: This is the probability per unit time and per unit volume that gluons are created by the chromoelectric field. Due to the masslessness of QCD gluons the above formula does not contain the Sauter-Schwinger exponentially small tunnelling factor exp (−π m 2 c 3 e E ) and even a weak chromoelectric field will break down creating a cloud of soft gluons from the vacuum neutralising the imposed colour electric field. While the exact results for the imaginary part of the effective action (5.91) depend on the details of the background field, it was argued in [77,78] that the threshold singularity is universal. The physical reason for this universality lies in the fact that the onset of pair production is dominated by the long-range fluctuations of the particles created from the vacuum and becomes insensitive to the details of the field profile [77,78].
The infrared long wavelength gluons created from the vacuum are strongly interacting. In the one-loop approximation the interaction between the produced gluons is not considered. In [28,29] the authors considered a semiclassical corrections to the production rates due to the interactions between created pairs and suggested a mechanism of colour neutralisation.
By considering a cylindrical non-homogeneous chromoelectic field configurations in [27] Flory demonstrated a possible formation of a chromoelectic flux tube of a finite radius between quarkantiquark pairs which are embedded into the chromomagnetic gluon condensate [13] (see the next section). The alternative mechanisms of formation of chromoelectric and chromomagnetic flux tubes were considered in [26] and in [25].
Let us now apply the renormalisation condition (2.36) to the Yang-Mills effective Lagrangian (5.87) and (5.88). This leads to the following expression for the renormalised effective Lagrangian in chromomagmentic fields [13,11]: and the effective Lagrangian in SU(N) gauge theory will take the following form [13]: The QCD vacuum responds to the background magnetic field as paramagnet with the magnetic permeability of the following form [11]: The paramagnetism of the QCD vacuum means that it amplifies the applied chromomagnetic field by generating induced chromomagnetic field in the direction of the applied field. This phenomenon is similar to the Pauli paramagnetism, an effect associated with the polarisation of the electron spins. In QCD the polarisation of the vector boson spins is responsible for this amplification of the background field (5.89). This also can be seen from the vacuum energy density ( see Fig.2 ) with its new minimum outside of the perturbative vacuum < F >= 0, at the renormalisation group invariant field strength [13] < 2g 2 F > vac = µ 4 exp (− understood as average over these field configurations (see also [63,64,51]).
For the energy momentum tensor (4.64) we shall get The trace of the energy momentum tensor is not equal to zero and characterises the breaking of conformal symmetry in QCD: The vacuum energy density is given in (6.98): T 00 = (F) with its minimum at (6.99) [13].
Substituting this value into the expression for the energy momentum tensor (6.100) we shall get the expression which is proportional to the metric tensor g µν : and is therefore a relativistically invariant characterisation of the vacuum with its negative energy density vac =< T 00 > vac and the pressure P vac = − vac . This is an important result because the vacuum state should be Lorentz invariant and its stress tensor T µν should be the same in The Lorentz invariant average < > over the covariantly constant field Gµν orientations can be performed as in [58], [13]. In [58] the invariant measure was taken in the form all frames [15,16]. As a result, its vacuum average value can only be of the cosmological type < 0|T µν |0 >= vac g µν , and indeed it is.
Let us consider the behaviour of the effective Lagrangian from the renormalisation group point of view and compare it with the behaviour of the effective coupling constant. The equations derived above are universally true for the non-Abelian field as well. Thus when G = E a H a = 0 we have The vacuum magnetic permeability introduced in (2.42) will take the following form [11]: , The Callan-Symanzik beta function can be calculated by using (6.95): and the effective coupling constant as a function of the field has the form where we introduced the Casimir operator C 2 (G) = N for the gauge group G = SU (N ).
Let us consider the value of the field strength tensor F 0 at which the vacuum energy density (6.98) vanishes (F 0 ) = 0, as it is shown on Fig.2: The effective coupling constant (6.106) at this field strength has the valuē It follows that the effective coupling constant at the intersection point F 0 is small: The energy density curve (F) (6.98) intersects the horizontal zero energy line at the nonzero angle θ (see Fig.2): tan θ = 11g 2 N 96π 2 > 0. (6.110) This means that the true vacuum state is below the perturbative vacuum, and we proved that i) the true vacuum is below the perturbative vacuum, ii) that there is a nonzero chromomagnetic field in the vacuum and iii) that there is a nonzero energy gap between perturbative and nonperturbative vacuum.
Now the question is, how far into the infrared region one can continue the energy density curve by using the perturbative result? Let us consider the fields which are approaching the infrared pole. This can be done, in particular, by using the following parametrisation: (6.111) where the parameter n is less than one, and we have F n → < F > vac when n tends to unity from below. At these fields values the effective coupling constant (6.106) tends to zero: if the product N (1 − n) → ∞ is large and the t'Hooft coupling constant g 2 N = λ is fixed and small. It follows then that the effective coupling constant can be made small to justify the use of the perturbative result and the energy density curve can be continuously extended infinitesimally close to the value of the vacuum field < F > vac , as it is shown on Fig. 2. In an alternative way one can study the ratio of quantum mechanical correction L (1) to the classical Lagrangian At the intersection point F 0 this ratio is 11g 2 N 96π 2 1, that is, one should have which is fulfilled for large N and small coupling constant g 2 or in terms of t'Hooft coupling N . It is interesting to know if the energy density curve is a continuous function of the field strength F in the region [0, < F > vac ] which is outside of the validity of the perturbative calculations and that the energy density curve is a convex function.
A nonperturbative functional method developed by Zwanziger in [47] is answering to these questions affirmatively. It seems that further development of his approach can shed even more light to the behaviour of the effective Lagrangian in the nonpertubative region. We already obtained the first derivative of the energy density curve (6.103) and can calculate its second derivative as well: The sign of the second derivative depends on the sign of the ratio of the beta function to the cube of the effective coupling constant. In QCD, in the perturbative regime this ratio is negative and the second derivative (6.115) is positive: Thus the energy density curve is convex (see Fig.2). In QED the overall sign is negative and the energy density curve is concave (see Fig.1). Any non-perturbative information about the ratiō β(ḡ)/ḡ 3 can be translated into the information about property of the energy density curve. As far as the beta functionβ(g) has no zeros, is a negative analytical function of the coupling constant and ∞ g dḡ β(g) < ∞ (6.117) the minimum of the energy density curve is defined by the extremum, where the derivative (6.103) vanishes. By using the expression (6.103) one can derive the value of the chromomagnetic condensate [13]: To all orders in the perturbation theory the derivative of the energy momentum tensor trace can be obtained by using the renormalisation group invariant result (4.68): Integration over F gives the trace If one considers the approximation in which the effective coupling constant (3.55) is field independentḡ(t) ≡ g then this formula after integration over F will reduce to the one given in literature [52,53,54]: Otherwise the field dependence of the energy momentum trace is defined through the beta functions and effective coupling constant and has more complicated dependence on field strength tensor F.
As is follows from (6.102), in the ground state the following relation between energy density vac and pressure P vac takes place: vac = −P vac . It is a relativistically invariant characterisation of the vacuum [15,16], and it represents a field-theoretical contribution into the effective cosmological The chromomagnetic condensate (6.99) is of order Λ 4 QCD , and the vacuum energy density is negative and is about vac ≈ −b 10 −8 GeV 4 . The value of the cosmological constant measured in the observation of the high-z Type Ia supernovae [19,20,21,22] and by the Plank Collaboration [23,24] Λ = c 4 Λ obser /8πG ≈ 10 −47 GeV 4 is about 39 decimal places smaller and positive. It is important to mention that the energy gap depends on a gauge group and a matter content, the beta function (the b parameter in one-loop approximation), as well as of the temperature of the universe [49]. At high temperatures the curve of the effective potential moves upward, the value of the chromomagnetic gluon condensate tends to zero, as well as the Λ ef f , and the scaling invariance get restored. The phase transition is of the second-order [49].
In the recent article [68] there was suggested a possible cancelation mechanism between chromomagnetic and its "mirror chromoelectric" condensates. In this proposal, which involves adding to the SM particles a mirror world (dark matter) [69,70,71,72,73,74], the entire SM is replicated in a mirror world. The new Z 2 symmetry interchanges SM with the mirror SM, ensuring identical particles and interactions. It is conjectured that the quantum vacua of the "Mirror SM" contribute to the cosmological constant on the same footing as the SM, since mirror particles are expected to gravitate in the same way as the usual ones, and that the mirror chromoelectric gluon condensate contributes to the energy density of the universe with a positive sign and thus may, in principle, eliminates the negative QCD vacuum effect by yielding a cosmological constant small.
In electroweak theory the Higgs vacuum field also generates a negative contribution to the effective cosmological constant [17,18]. In both field theories the value of the Λ ef f is of many orders of magnitude larger than the observed value.
A careful inspection of the charged vector bosons spectrum in a chromomagnetic field by Nielsen, Olesen [30] and Skalozub [31] demonstrated that due to the unstable mode k 2 0 = k 2 || − gf 1 (k 2 || ≤ gf 1 ) there is an imaginary part in the effective Lagrangian: In the subsequent publications [32,33,34,38,36,39,40,47,48,49,51,94] the theoretical groups at the Niels Bohr Institute, New York University, SLAC and Bari University came to the conclusion that due to the quartic self-interaction term in the YM action there is a hidden Higgs mechanism, which stabilises the system so that the effective Lagrangian remains a real function in the background chromomagnetic field [33,38,40,47,48]. This is the reflection of the fact that calculations are performed in the approximation in which only the quadratic term Lagrangian where the zero field configuration is unstable and a field configuration at the bottom of the potential provides a stable field configuration due to the quartic term. In the pure YM theory the higgs-like action for the unstable mode was derived in [26,33,38,40]. It was proposed to search the stable solutions of the classical Yang-Mills field equations in a fixed background chromo-magnetic field which plays the role of an external order parameter. Without the presence of the order mass parameter, the gH in the given case, the conformal invariance of the pure classical Yang-Mills equations prevents the existence of localised solutions [25]. This program was successfully realised with the discovery of the field configurations which are varying in space due to the development of the unstable mode, the colour magnetic flux tubes and the spaghetti magnetic tubes forming the domain-like field configurations [38]. The configurations are supported by the external chromomagnetic field gH. The difficulty here is to calculate the quantum-mechanical fluctuations around these classical field configurations and to see if they remain localised when the external field is switched off. The important conclusion of the investigation was that it pointed out to the fact that the stability of the chromomagnetic field configurations is a natural consequence of the quartic self interaction of the Yang-Mills field.
This result became the initial point for the investigation initiated by Curt Flory in his article devoted to the resolution of the higgs-like mode problem [48]. His breakthrough idea was to integrate exactly the functional integral over the higgs-like mode from the start in order to get the quantum-mechanical contribution to the effective Lagrangian of that mode instead of searching the corresponding classical field configurations. Presenting the amplitude of the higgs-like mode and of the corresponding action in terms of dimensionless variables one can get * * : is the dimensionless amplitude of the higgs-like mode. The action of the higgslike mode will take the following form: What is essential in this representation is that the dependence on the chromomagnetic field does not show up in the Lagrangian (8.124) and appears only in front of the higgs-like field amplitude (gH) 1/2 in (8.123). The factor coming from the integration of the action (8.124) over the field Φ k 2 is background field independent and does not show up in the renormalised effective Lagrangian.
Thus the contribution of the higgs-like mode to the effective Lagrangian is only through the integration measure and its degeneracy: where C is the gH independent value of the functional integral over Φ k 2 (x 0 , x 3 ). This contribution is a real function of chromomagnetic field [48,49]. After taking into account the contributions from all other modes the effective Lagrangian takes the form which identically coincides with (6.94). This confirms the expression (6.94) being without imaginary part (5.90).
One can consider the above approach of calculating the effective action as an alternative to a standard loop expansion in the following sense: The expansion is organised by rearranging the perturbative expansion (5.72) in a background field A so that the quartic self-interactions of eigenmodes are included into the propagator of the gauge field G(x, y; A) and the loop expansion is performed in terms of the remaining cubic and cross-mode quartic vertices of the YM action. * * In the SU(2) case Wµ = 1 √ 2 (A 1 µ + A 2 µ ), Aµ = A 3 µ and W = W1 = −iW2, as it is defined in [26,33,38,40].
The short overview of the publications devoted to the chromomagnetic gluon condensation and QCD vacuum are given below. The confinement problem from the point of view of the QCD vacuum and chromomagnetic gluon condensate was considered in the articles of Mandelshatam [85,86], Nambu [87], Adler and Piran [88] and Nielsen and Olesen [37]. The thermodynamics of the Yang-Mills gas by Linde [18]. The publication on generation of galactic magnetic field due to the condensation of vector field was considered in [81]. The induced gravity was considered by Adler [56] and the effective cosmological constant in [68]. The magnetostatics was considered in [89]. The phenomenology of hadrons and the properties of the QCD vacuum by Shuryak [96]. The mechanism of dynamical supersymmetry breaking and string compactification to four dimension due to the properties of the non-Abelian effective action was suggested by Veneziano and Taylor [95]. The dynamical mass generation in QCD and glueballs spectrum by Cornwall [91,92]. The string-like solution of pure YM equations stabilised by the presence of the condensate by Faddeev and Niemi [93].
The discovery of the chromomagnetic gluon condensation in [13] initiated series of publications by the ITEP group where they used the gluon condensate to improve their perturbative sum rule equations [65,66] † † . Modern determination of the gluon condensate numerical value from hadronic τ decay data and from the charmonium sum rules can be found in the review article of Ioffe [99]. The best values of condensates, extracted from QCD sum rules from experimental data, are given in Table 1 in [99]. These data do not allow to exclude the zero value for the gluon condensate [97,98,99]. The separation of perturbative and nonperturbative contributions has some arbitrariness [66], as it was pointed out by Ioffe in [99].
Here we reexamined the proof of the existence of the YM condensate [13]. It is of the chromomagnetic type and it has a numerical value Λ 4 QCD which is of the order of few hundred M eV 4 ) = Λ 4 QCD > 0. (9.125) or in terms of the strong coupling constant (9.126) † † In 1977 the author gave a theoretical seminar on the chromomagnetic gluon condensation [13] in ITEP. At end of the seminar one of the participants, Victor Novikov, on our way back to the metro station by tram, remarked to the author that the theoretical prediction of the chromomagnetic condensate presented at the seminar [13] can be crucial in improving the naive sum rule equations published earlier in [65] by introducing the chromomagnetic condensate in the form of power corrections. A year later, the proposal was realised in [66].
The integrals appearing in the effective Lagrangian (2.37) have the following form: where k can be considered as a dimensional regularisation parameter and the integrals should be calculated in the limit k → −1 [11]. As far as b = a, the second integral does not completely coincide with the one appearing in the effective Lagrangian (2.37) and we have to consider its extension. In order to calculate the integral we shall take b = a − and consider the limit → 0: By the definition the Riemann zeta function ζ(k, q) is [67] ζ(k, q) = The Lagrangian should be calculated in the limit k → −1.