A method for getting a finite $\alpha$ in the IR region from an all-order beta function

The analytical method of QCD running coupling constant is extended to a model with an all-order beta function which is inspired by the famous Novikov-Shifman-Vai\-n\-s\-htein-Zakharov beta function of N=1 supersymmetric gau\-g\-e theories. In the approach presented here, the running coupling is determined by a transcendental equation with non-elementary integral of the running scale $\mu$. In our approach $\alpha_{an}(0)$, which reads 0.30642, does not rely on any dimensional parameters. This is in accordance with results in the literature on the analytical method of QCD running coupling constant. The new"analytically im\-p\-roved"running coupling constant is also compatible with the property of asymptotic freedom.

The analytical method of QCD running coupling constant is extended to a model with an all-order beta function which is inspired by the famous Novikov-Shifman-Vainshtein-Zakharov beta function of N=1 supersymmetric gauge theories. In the approach presented here, the running coupling is determined by a transcendental equation with non-elementary integral of the running scale µ. In our approach αan(0), which reads 0.30642, does not rely on any dimensional parameters. This is in accordance with results in the literature on the analytical method of QCD running coupling constant. The new "analytically improved" running coupling constant is also compatible with the property of asymptotic freedom. The renormalization scale dependence or evolution of running coupling constant is described by the beta function of a theory concerned. The explicit expression for beta function can be obtained by means of perturbation methods, and the validity of its use depends on whether its resultant coupling constant is excessively large or not. Solving a given beta function, we can obtain an expression for running coupling constant, which may blow up or become enough large somewhere, indicating that the perturbation breaks down and the invalidity of using this given beta function. This problem leads us to expect that, as has been done in Ref. [1], running coupling constant must cease to increase somewhere (i.e. freeze somewhere).
The freezing of running coupling constant (FORCC) is closely related to the explicit form of beta function. Using beta function to study the FORCC is meaningful and direct; the freezing of running coupling constant might appear as a result of beta function's vanishing-FORCC is related to infrared or ultraviolet zero of beta function. In Ref. [2] infrared zero of beta function was studied at three-loop order. In Ref. [3][4][5][6][7][8][9] this was carried out to four-loop order. Recently, in Ref. [10] this has been carried out to five-loop order.
Also, there are theoretical reasons in favor of FORCC [2,22]. Even though the behaviour of running coupling of QCD in IR region is nonperturbative in usually sense, in Ref. [16] it was shown that perturbation theory provided some theoretical evidence for FORCC.
Many methods have been used to investigate FORCC. For example, in Ref. [23] Schwinger-Dyson equation was used to investigate FORCC. Another well known way to study FORCC is to study coupling constant directly-we can obtain the knowledge of running coupling in IR region from the knowledge of it in UV region [1]. In Ref. [1], through the use of "analytization procedure", an "analytically-improved" expression for running coupling was obtained, which is free of ghost-pole and has a universal limit value. The "analytization procedure" elaborated in Ref. [24][25][26] and used in Ref. [1] includes three steps: (I) Finding an explicit expression for α s (µ 2 ) in the Euclidean region; (II)Performing analytical continuation into the Minkowski region.
In some sense these three steps are perturbative, since α is obtained by perturbation methods usually. An allorders beta function was first proposed in Ref. [27]. A proof of this is given in Ref. [28]. In our work, this beta function is used to extract α an (µ) by a similar way as the way used in Ref. [1].
The remainder of this paper is organised as follows. In section 2, we start with an all-orders beta function and concentrate on eliminating the ghost pole. In section 3, after having removed singularity, we get an equation from which we define a new "analytically-improved" running coupling which is free of ghost pole and does respects asymptotic freedom. In section 4, the issue of scheme dependence is discussed. As has been said before, an all-order beta function was first proposed in Ref. [27] and proved in Ref. [28], which takes the form For pure Yang-Mills theory this can be simplified to . ( In SU(N) gauge theory, C 2 (G) is just N. Carrying out the obvious transformation α 2π → α converts Eq.
(2) to a new compact beta function in the form where C 1 = 11C 2 (G)/3 and C 2 = 17C 2 (G)/11. The general solution of this Equation is of the form where Λ is just a integral constant. For later convenience, inverting Eq. (4) and making a replacement For simplicity, we denote by y 0 (µ 2 E ) the right-hand side of Eq. (5) and by y 0 (α(µ)) the left-hand side of Eq. (5).
The analyticity property of y 0 (α(µ)) is almost the same as that of running coupling α(µ) with the exception that the coupling constant is non-positive or may be a pole of y 0 (α(µ)). Therefore, in this work we shall concentrate on studying the analyticity properties of y 0 (α(µ)) and y 0 (µ 2 E ). For simplicity, we set α(Λ) = 1 and λ = C 2 , then .
It may be shown that the denominator of y 0 (α(µ)) can be zero only if the running coupling α is a real number.
For the continuation of this work, we make some reasonable assumptions here; self-consistency of some of them can be verified by later check.
Let us make the first assumption: for the running coupling to be real-valued, what is necessary is that µ 2 must be real-valued. Self-consistence of this assumption can be checked later.
The second assumption is the absence of Landau-pole. Therefore, according to asymptotic freedom, which will be discussed later, there is no point on the real axis at which α s is negative. When µ 2 E goes to−∞, according to asymptotic freedom, α(µ) tends to 0, which leads to the denominator of y 0 (α(µ)) having a non-zero value. Therefore, the distance between points at which the denominator of y 0 (α(µ)) vanishes and the origin must reach its maximum (if exist) at a point on the negative real axiswe set it to be −µ 2 0 . As is well known, if a function is analytical in a connected region surrounded by a contour, we can use Cauchy theorem which in this paper is used to subtract singularity terms. The third assumption is that the integral along infinite contour is summed to zero. Note that y 0 (α(µ)) is analytical except µ 2 E is a pole of it on the real axis. Therefore, in this paper, the integral is taken over the contour shown in Fig. 1. Thus we immediately arrive at We note that Im ln Thus Eq. (8), can be further simplified into

III. EVALUATION OF COUPLING CONSTANT
It's almost always the case that the values of y 0 (µ 2 E ) are related to Λ. Two cases Λ = 200M eV and Λ = 400M eV are shown in Fig. 2. However, there is one obvious exception that y an (0) is obviously independent of the value of Λ. Therefore, the value for α an (0) is independent of Λ.
Since the first part of the integrand in Eq. (10) is a monotone decreasing function of energy scale µ, it's obvious that, as can be seen form Fig. 2, y 0 (µ 2 E ) is a monotone decreasing function of energy scale.
By straightforward calculation, we obtain an explicit expression for Eq. (10) in the form which decreases monotonically from 2/C 1 to zero as µ varies from zero to infinity, and is free of ghost-pole compared with the original expression 2/(C 1 ln µ 2 /Λ 2 ). Acquisition of α an (µ 2 ), in our approach, is determined by equation y an (α an (µ)) = α an (µ) [1 + α an (µ)(λ ln α an (µ) − 1)] which is an "analytically-improved" equation corresponding to y 0 (α(µ))=y 0 (µ 2 E ). As can be seen from Fig. 3 the new "analyticallyimproved" running coupling α an as a function of µ is two-valued: each value of the y an (µ) corresponds to two different values of α an . The difference between this two choices of value for running coupling, as can be seen from Fig. 3, is obviously large and varies monotonously with energy scale. Therefore, the determination of running coupling is a little more complicated.
It may be shown that, of the three regions represented by A, B and C respectively in Fig. 3, the middle one (represented by B) corresponds to un-physical region and therefore can be discarded, since y an (µ), with an explicit expression shown in Eq. (11), is always positive. Since y an (µ) is a monotone decreasing function of energy scale, of the remainder region A and C, the left region A in Fig.  3 corresponds to region obeying asymptotic freedom (as energy scale increases, y an (µ) decreases which leads to the decrement of α an (µ)), and the right region C in Fig.  3 corresponds to region violating asymptotic freedom (as energy scale increases coupling constant increases, what is worse, in UV region the coupling constant blows up). Taking our assumption of asymptotic freedom, physical intuition and mathematical continuity into consideration region C also can be discarded. The distinction between regions A and C also may be traced back to the explicit form of all-orders beta function in Eq. (3). As energy scale changes from IR region to UV region, evolution of coupling constant is essentially determined by the initial value of running coupling in IR region. If this value is smaller than 1 C2 , asymptotic freedom is obeyed, otherwise is violated dramatically. Now for a fixed µ, taking all the discussion given above into consideration, we can fix the corresponding α an (µ). Even though, it's still impossible to acquire a expression for α an (µ) in terms of µ explicitly, numerical value for α an can be obtained for any µ. Letting µ vanish, for SU(3) gauge theory, we get a numerical value 0.0487688 in the region A of Fig. 3 (region obeying asymptotic freedom) and a numerical value 3.8413 in the region C of Figure.3 (region violating asymptotic freedom dramatically). Therefore, we choose 0.0487688 as our result. Since a transformation α 2π → α has been made early, the running coupling constant in IR region freezes at 0.30642.

IV. DISCUSSION OF SCHEME DEPENDENCE
As is well known, for a given theory to be renormalizable, we impose the requirement that the divergent part of any one-particle irreducible diagram can be cancelled by a local counterterm generated from the bare Lagrangian of this given theory. Therefore, once the divergences of a given theory have been subtracted out by a certain renormalization scheme, one is still free to perform further finite renormalization, thus resulting in the may-be scheme dependence of quantities in field theory. Here in this section, our discussion is restricted within the class of mass-independent renormalization scheme.
As far as we are concerned here, for beta function, usually only the first two none-vanishing coefficients of the expansion for beta function are scheme independent. There is one obvious may-be exception to this statement just given above. If, for a given theory, an all-orders beta function can be obtained in a certain renormalization scheme and written in the form and if this form can be maintained in anther scheme, it's obvious that the all-orders beta functions in this two schemes are equal to each other. The expression for beta function expanded in terms of running coupling is usually of the form β(α) = dα d ln µ = β 0 α n0 + β 1 α n1 + β 2 α n2 + . . . , (15) which can be rewritten as These new coefficients of Eq. (16) can be obtained by matching with Eq. (15) order by order. Along the same procedures carried out above, we may obtain an equation between the new "analytically-improved" coupling α an and energy scales in the form where Λ is just an integral constant and C also a constant. It's obvious that, if α an is in the perturbation region (i.e. small enough), of all the terms in the denominator of the expression in the left-hand side of Eq. (17), the term with a lower index i contribute more (the first two none-vanishing term with scheme independent coefficients contribute most).
To conclude this section, we may make the following remarks concerning the scheme dependence of the new "analytically-improved" freezing coupling. First, in our approach the existence of freezing coupling is independent of renormalization scheme we choose, though the explicit value is scheme dependent. Second, if the running coupling is in perturbation region, the freezing coupling is effectively (within the accuracy required) determined by a finite number of none-vanishing coefficients of beta function.

V. CONCLUSION
The appearance of ghost pole in the expression for running coupling of QCD or QED, leads to the hypothesis that the running coupling must freeze somewhere. In this article, starting with an all-orders beta function, we have obtained an 'analytically-improved" coupling by some analyticity procedure to remove singularity. This new running coupling is free of ghost pole and does respect asymptotic freedom.
Many works have been done to investigate FORCC theoretically or phenomenologically. It is clear that much more work has to be done in order to get a profound understanding of this amazing property which is essential to physics. Now we can set − π 2 ≤ θ ≤ π 2 . From this condition we have 2 π ≤ λr ≤ 1, (A5) λ ln r + 1 r ≥ 1. (A6) Note that under the condition π 2 ≤ λr ≤ 1 there is no solution to λ ln r + 1 r > 1-this two conditions contradict. Thus this statement is proved.
Secondly, we check the assumption-when µ 2 is not real, coupling is not real.
From Eq.10, we can extract the imaginary part of y an (µ 2 ), which can be written as It's obvious that Im(µ 2 E ) = 0 implies Imy(µ 2 E ) = 0, which leads to the appearance of an imaginary part in y an (α an ). This can happen only when α is not real number. Thus we finish our check.