The ηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c$$\end{document} decays into light hadrons using the principle of maximum conformality

In the paper, we analyze the ηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c$$\end{document} decays into light hadrons at the next-to-leading order QCD corrections by applying the principle of maximum conformality (PMC). The relativistic correction at the O(αsv2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{{O}}(\alpha _s v^2)$$\end{document}-order level has been included in the discussion, which gives about 10%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10\%$$\end{document} contribution to the ratio R. The PMC, which satisfies the renormalization group invariance, is designed to obtain a scale-fixed and scheme-independent prediction at any fixed order. To avoid the confusion of treating nf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_{f}$$\end{document}-terms, we transform the usual MS¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm{MS}}$$\end{document} pQCD series into the one under the minimal momentum space subtraction scheme. To compare with the prediction under conventional scale setting, RConv,mMOM-r=4.12-0.28+0.30×103\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_\mathrm{{Conv,mMOM}-r}= \left( 4.12^{+0.30}_{-0.28}\right) \times 10^3$$\end{document}, after applying the PMC, we obtain RPMC,mMOM-r=6.09-0.55+0.62×103\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_\mathrm{PMC,mMOM-r}=\left( 6.09^{+0.62}_{-0.55}\right) \times 10^3$$\end{document}, where the errors are squared averages of the ones caused by mc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_c$$\end{document} and ΛmMOM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _\mathrm{mMOM}$$\end{document}. The PMC prediction agrees with the recent PDG value within errors, i.e. Rexp=6.3±0.5×103\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^\mathrm{exp}=\left( 6.3\pm 0.5\right) \times 10^3$$\end{document}. Thus we think the mismatching of the prediction under conventional scale-setting with the data is due to improper choice of scale, which however can be solved by using the PMC.

mMOM . The PMC prediction agrees with the recent PDG value within errors, i.e. R exp = (6.3 ± 0.5) × 10 3 . Thus we think the mismatching of the prediction under conventional scale-setting with the data is due to improper choice of scale, which however can be solved by using the PMC.
The heavy quark mass provides a natural hard scale for the heavy quarkonium decays into light hadrons or photons. Calculations of their decay rates are considered as one of the earliest applications of pQCD. The charmonium has become a popular field since the discovery of J/ψ resonance at SLAC and Brookhaven in 1974. There are lots of successful experimental studies about charmonium, including the precise measurements of spectrum, lifetimes and branch ratios, cf. a comprehensive review given in the PDG [1]. At the same time, many theoretical efforts have been tried for an appropriate description of charmonium. As an important breakthrough, a systematic pQCD analysis of the heavy quarkonium inclusive annihilation and production has been given within the nonrelativistic QCD theory (NRQCD) in 1995 [2]. a e-mail: wuxg@cqu.edu.cn According to the NRQCD framework, the quarkonium decay rate can be factored into a sum of products of the shortdistance coefficients and the long-distance matrix elements (LDMEs). The short-distance coefficients are perturbatively calculable in a power series of α s . The LDMEs can be estimated by means of the velocity power counting rule, i.e. the LDMEs can be classified in terms of the relative velocity between the constituent quarks of the heavy quarkonium. Especially, the color-singlet ones can be directly related to the wavefunction (derivative of the wavefunction) at the origin, which then can be calculated via proper potential models.
The decay rates of the pseudoscalar quarkonium into light hadrons and photons have been calculated at the next-toleading order (NLO) level [3,4]. The relativistic corrections at the O(α s v 2 )-order have been given in Refs. [5,6]. Within the NRQCD factorization framework, the decay rate of the η c into light hadrons or photons can be expressed as and where F 1 , G 1 , F γ γ and G γ γ are short distance coefficients. The symbol · · · stands for the contributions from highdimensional LDMEs which are at least at the level of O(v 4 ). m c is the c-quark pole mass. 1 v 2 is the squared heavy quark or antiquark velocity in the meson rest frame. For the case of η c , it can be calculated by To suppress the uncertainty from the LDMEs, one usually calculates the ratio where an arbitrary renormalization scale, and β 0 = 11 − 2 3 n f (n f being the active flavor number) is the leading β-term of the renormalization group function. It is noted that the factorization scale dependence is missing at this level, which is the case even at the NNLO level [8], we are thus free of the factorization scale-setting problem.
It is conventional to take the renormalization scale as the typical momentum flow of the process or the one to eliminate the large logs of the pQCD series, we call this conventional scale-setting approach. As will be shown later, such a simple treatment on scale introduces large scale uncertainty and makes the lower-order prediction unreliable. At present, the η c decays into light hadrons or photons have been calculated up to NNLO level, which however still shows a poor pQCD convergence [8][9][10]. Thus by simply pursuing higherand-higher order terms may not be the solution for those high-energy processes. In fact, even if we obtain a small scale uncertainty for global quantities such as the total crosssection or the decay rate at a certain fixed order, it is due to cancelations among different orders; the scale uncertainty for each order is still uncertain and could be very large. Two such examples for Higgs boson decay and the hadronic production of Higgs boson can be found in Refs. [11,12]. When one applies conventional scale-setting, the renormalization scheme-and initial renormalization scale-dependence are introduced at any fixed order. Thus, a proper scale-setting approach is important for the fixed-order predictions.
We should point out that those predictions are different from the value derived from the new experimental measurements, which shows R exp = (6.3 ± 0.5)×10 3 [1]. As will be shown later, the BLM prediction given in Ref. [13] is questionable. Thus it is interesting to show whether an improved pQCD analysis could be done and could explain the new R exp , as is the purpose of this paper. Especially, it is important to show whether the mismatching of the data and the pQCD prediction is caused by improper choice of scale or by some other reasons.
A novel scale-setting approach, the Principle of Maximum Conformality (PMC) [18][19][20][21], has been developed in recent years. The PMC satisfies renormalization group invariance [22] and it reduces in N C → 0 Abelian limit [23] to the standard Gell-Mann-Low method [24]. A more convergent pQCD series without factorial renormalon divergence can be obtained. The PMC scales are physical in the sense that they reflect the virtuality of the gluon propagators at a given order, as well as setting the effective number (n f ) of active flavors. The resulting resummed pQCD expression thus determines the relevant "physical" scales for any physical observable, thereby providing a solution to the renormalization scalesetting problem. Because all the scheme-dependent {β i }terms in pQCD series have been resummed into the running couplings with the help of renormalization group equation, the PMC predictions are renormalization scheme independent at every order. Such scheme independence can be demonstrated by using commensurate scale relations [25] among different observables. A number of PMC applications have been summarized in the review [26][27][28]. The PMC provides the underlying principal for the BLM, and in the following, we shall adopt the PMC to set the renormalization scale.
Up to NLO level, the expression of R can be rewritten as where the MS-coefficients r i, j can be read from Eq. (4), in which r i,0 are conformal ones. Following the standard PMC procedures, we get where ln Q 2 1 /μ 2 = −r 2,1 /r 1,0 . Here, we have set the unknown PMC scale Q 2 = Q 1 such that to ensure the scheme-independence of R under any renormalization schemes via proper commensurate scale relations [25], whose exact value can be determined by the NNLO term which is not available at present.
If directly using the MS-scheme expression (4), we shall obtain a small PMC scale Q 1 = 0.86 or 0.78 GeV for the prediction with or without relativistic correction. It is already close to the low-energy region, this explains why a large R BLM is obtained in Ref. [13]. [At the NLO level, the BLM prediction is the same as the PMC prediction if all n f -terms are pertained to α s -running.] For this case, a reliable prediction can only be obtained by using certain low-energy α s -model, which however will introduce extra model dependence for the prediction.
Following the idea of PMC, only those {β i }-terms that are pertained to the renormalization of running coupling should be absorbed into the running coupling. For the processes involving three-gluon or four-gluon vertex, the scalesetting problem is more involved [29]. The MOM scheme is a physical scheme which is based on the renormalization of the triple-gluon vertex at some symmetric off-shell momentum. The MOM scheme carries information about the vertex at a specific momentum configuration. This external momentum configuration is non-exceptional and there are no infrared issues, thus avoiding the confusion of distinguishing {β i }-terms. Thus to avoid the ambiguity of applying the PMC on R, similar to the case of QCD BFKL Pomeron [30][31][32], we shall first transform the results from the MS-scheme to the momentum space subtraction scheme (MOM-scheme) [33,34] and then apply the PMC. Another reasons for choosing the MOM scheme lie in that a better pQCD convergence can be obtained by using the MOM scheme than using the MS-scheme, and a more reasonable PMC scale in perturbative region can be achieved.
For the purpose, we adopt the perturbative relation between the MS-scheme running coupling and the mMOMscheme one as [35] a MS (μ) = a mMOM (μ) 1 − 4D 1 a mMOM + · · · , where for the Landau gauge, D 1 = d 1,0 + d 1,1 n f , d 1,0 = 169 144 N C , and d 1,1 = − 5 18 . We then obtain where R mMOM After applying the PMC, we obtain a new PMC scale Q 1 = exp(−3d 1,1 )Q 1 , which equals to 1.99 or 1.80 GeV for the prediction with or without relativistic correction. Such a larger PMC scale indicates a reliable pQCD prediction can be achieved by using the mMOM scheme.
we obtain m c = 1.49 ± 0.03 GeV. To be consistent, we adopt the two-loop α s -running, whose behavior is fixed by using the reference point α s (m Z ) = 0.1181 ± 0.0011 [1]. And we adopt v 2 η c = 0.430 GeV 2 /m 2 c [40,41]. Numerical results of the QCD asymptotic scales MS and mMOM under Landau gauge are listed in Table 1, where the errors are dominantly caused by the uncertainty α s (m Z ) = ±0.0011. The asymptotic scales for different schemes satisfy the relation [11,35], mMOM / MS = exp(2D 1 /β 0 ).
As a cross-check, by using the same input parameters, we obtain the same MS-scheme prediction on R under conventional scale-setting as that of Ref. [13]. Due to the reasons listed above, we shall adopt the mMOM-scheme to do our following discussions.
We present the PMC prediction on R at the NLO level versus the initial choice of μ in Fig. 1, which is under the mMOM scheme and both the results before and after applying the PMC are presented. Under conventional scale setting, R shows a strong scale dependence which decreases with the increment of μ. More explicitly, by varying μ from m c to 4m c , the ratio R will change from ∼ 9 × 10 3 to ∼ 3 × 10 3 .  μ = 2m c . The symbol "−r " stands for relativistic corrections. κ = R NLO /R LO , which shows the relative importance of the NLO-term and the LO-term of R After applying the PMC, the PMC scale Q 1 is the same for any choice of μ, leading to scale independent prediction. The relativistic correction brings an extra ∼ 2% contribution to the conventional prediction and ∼ 14% contribution to the PMC prediction. Thus the relativistic correction is important, especially for the PMC predictions. Figure 1 shows that if choosing μ = Q 1 , the values of R under conventional scale setting shall be equal to the PMC ones. After applying the PMC, due to the elimination of divergent renormalon terms as n!β n 0 α n s , the pQCD series shall be more convergent. We present the LO and NLO terms of R before and after applying the PMC in Table 2. We define a parameter κ = R NLO /R LO to show the relative importance of the NLO-term and the LO-term. Table 2 confirms that a better pQCD convergence can be achieved by applying the PMC. A larger κ and a larger scale uncertainty for each term indicate that one cannot get the exact value for each term by using a guessed scale suggested by conventional scale-setting.
Analyzing the pQCD series in detail, we observe that the scale errors for conventional scale-setting are rather large for each term, and a possible net small scale error for a pQCD approximant is due to correlations/cancelations among different orders. On the other hand, due to the fact that the running of α s at each order has its own {β i }-series governed by the renormalization group equation, the β-pattern for the pQCD series at each order is a superposition of all the {β i }terms which govern the evolution of the lower-order α s contributions at this particular order. Thus, inversely, the PMC scale at each order is determined by the known β-pattern, and the individual terms of R at each order shall be well determined.
We present the theoretical uncertainties for the conventional and the PMC scale settings in Fig. 2, in which the errors are squared averages of the ones from the choices of the cquark pole mass m c and the asymptotic scale mMOM . As a comparison, the experimental prediction of Ref. [1] is also presented. Under conventional scale-setting, Fig. 2 shows that the errors caused by m c and mMOM is smaller than the case of PMC scale-setting, which is however diluted by the quite large scale uncertainty. For example, the value of R with R Conv,mMOM−r = 4.12 +0.12+0.28 −0.11−0.26 × 10 3 , R PMC,mMOM = 7.09 +0.32+0.75 −0.29−0.66 × 10 3 , R PMC,mMOM−r = 6.09 +0.21+0.58 −0.19−0.52 × 10 3 , where the first error is for m c ∈ [1.46, 1.52] GeV and the second one is caused by taking mMOM to be the values listed in Table 1. Figure 2 shows that the conventional prediction of R with or without relativistic correction is about 3.6σ deviation from the data. This discrepancy becomes even larger by including the NNLO term [8], thus the authors there even doubt the validity of NRQCD theory for this particular observable. 2 However, Fig. 2 shows that after applying the PMC, the pQCD prediction and the data are consistent with each other within reasonable errors even at the NLO level. The condition of the branching ratio 1/R is similar. This indicates that the large discrepancy between the data and the pQCD predic-tion is caused by improper choice of renormalization scale, and a simple guessed scale may lead to false prediction or false conclusion. Thus a proper setting of the renormalization scale is important for lower-order predictions.
As a final remark, one may also calculate R by using the determined predictions on the decay widths (1) and (2) separately. If taking all input parameters as the central value of our present choices, we obtain R Conv.,mMOM−r ∼ 2.57 × 10 −3 and R PMC,mMOM−r ∼ 2.64 × 10 −3 , both of which are quite different from our above predictions (11,13). Thus there are large differences for those two treatments on R, which starts at α 4 s -order level. Such large differences can be explained by the weaker pQCD convergence as can be seen from Table 2, which shall be suppressed by including more-and-more loop terms. We prefer the usually adopted way of using Eq. (4) to calculate R, in which the uncertainty from the LDME is suppressed and there is no factorization scale dependence up to the NNLO level.
As a summary, in this paper, we have studied the ratio of the η c (1S) decay rate into hadrons over its decay rate into photons by applying the PMC. The PMC provides a systematic way to set the optimal renormalization scale for high energy process, whose prediction is free of initial renormalization scale dependence at any fixed order. A more convergent pQCD series can be achieved and the residual scale dependence due to unknown high-order terms are highly suppressed. Figure 2 shows that the large discrepancy between the data and the pQCD prediction by using a guessed scale suggested by conventional scale-setting can be cured by applying the PMC. The PMC, with its solid physical and theoretical background, greatly improves the precision of standard model tests, and it can be applied to a wide variety of perturbatively calculable processes.