Search

We construct a QCD coupling ${\mathcal{A}}(Q^2)$ in the Effective Charge (ECH) scheme of the canonical part $d(Q^2)$ of the (inelastic) polarised Bjorken Sum Rule (BSR) ${\overline \Gamma}_1^{{\rm p-n}}(Q^2)$. In the perturbative domain, the coupling ${\mathcal{A}}(Q^2)$ practically coincides with the perturbative coupling $a(Q^2)$ [$\equiv \alpha_s(Q^2)/\pi$] in the four-loop ECH renormalisation scheme. In the deep infrared (IR) regime, ${\mathcal{A}}(Q^2)$ behaves as suggested by the Holographic Light-Front QCD up to the second derivative. Furthermore, in contrast to its perturbative counterpart $a(Q^2)$, the coupling ${\mathcal{A}}(Q^2)$ is holomorphic in the entire complex $Q^2$-plane with the exception of the negative semiaxis, reflecting the holomorphic properties of the BSR observable $d(Q^2)$ [or: ${\overline \Gamma}_1^{{\rm p-n}}(Q^2)$] as dictated by the general principles of the Quantum Field Theory. It turns out that the obtained coupling, used as ECH, reproduces quite well the experimental data for ${\overline \Gamma}_1^{{\rm p-n}}(Q^2)$ in the entire $N_f=3$ regime $0 < Q^2 \lesssim 5 \ {\rm GeV}^2$., Comment: v2, 15 pages, 7 figures; additional Figs.4 and 5; version as it appears in JHEP; relations Eqs.(35)-(39) of our arXiv 2312.13134v2, and the corresponding explanations, are summarized here as Appendix B for convenience

Approximate knowledge of the renormalon structure of the Bjorken polarised sum rule (BSR) ${\overline \Gamma}_1^{{\rm p-n}}(Q^2)$ leads to the corresponding BSR characteristic function that allows us to evaluate the leading-twist part of BSR. In our previous work \cite{pPLB}, this evaluation (resummation) was performed using perturbative QCD (pQCD) coupling $a(Q^2) \equiv \alpha_s(Q^2)/\pi$ in specific renormalisation schemes. In the present paper, we continue this work, by using instead holomorphic couplings [$a(Q^2) \mapsto {\mathcal A}(Q^2)$] that have no Landau singularities and thus require, in contrast to the pQCD case, no regularisation of the resummation formula. The $D=2$ and $D=4$ terms are included in the Operator Product Expansion (OPE) of inelastic BSR, and fits are performed to the available experimental data in a specific interval $(Q^2_{\rm min},Q^2_{\rm max})$ where $ Q^2_{\rm max}=4.74 \ {\rm GeV}^2$. We needed relatively high $Q^2_{\rm min} \approx 1.7 \ {\rm GeV}^2$ in the pQCD case since the pQCD coupling $a(Q^2)$ has Landau singularities at $Q^2 \lesssim 1 \ {\rm GeV}^2$. Now, when holomorphic (AQCD) couplings ${\mathcal A}(Q^2)$ are used, no such problems occur: for the $3 \delta$AQCD and $2 \delta$AQCD variants the preferred values are $Q^2_{\rm min} \approx 0.6 \ {\rm GeV}^2$. The preferred values of $\alpha_s$ in general cannot be unambiguously extracted, due to large uncertainties of the experimental BSR data. At a fixed value of $\alpha_s^{{\overline {\rm MS}}}(M_Z^2)$, the values of the $D=2$ and $D=4$ residue parameters are determined in all cases, with the corresponding uncertainties., Comment: 21 pages, 8 figures; v4: improved presentation, improved expression for D=2 OPE terms used and thus the results slightly changed, references added; to appear in Nucl.Phys.B

We use the known renormalon structure of Bjorken polarised sum rule (BSR) ${\overline \Gamma}_1^{p-n}(Q^2)$ to evaluate the leading-twist part of that quantity. In addition, we include $D=2$ and $D=4$ Operator Product Expansion (OPE) terms and fit this expression to available experimental data for inelastic BSR. Since we use perturbative QCD (pQCD) coupling, which fails at low squared spacelike momenta $Q^2 \lesssim 1 \ {\rm GeV}^2$ due to Landau singularities, the fit is performed for $Q^2 \geq Q^2_{\rm min}$ where $Q^2_{\rm min} \approx (1.7 \pm 0.3) \ {\rm GeV}^2$. Due to large BSR experimental uncertainties, the extracted value of the pQCD coupling has very large uncertainties, especially when $Q^2_{\rm min}$ is varied. However, when we fix the pQCD coupling to the known world average values, the $D=2$ and $D=4$ residue parameters can be determined within large but reasonable uncertainties., Comment: 10 pages, 3 figures; v2: minor changes and references added; v3: improved text, version as published in Physics Letters B

We obtain an improved determination of the normalization of the leading infrared renormalon of the Bjorken sum rule: $Z_B (n_f=3)= -0.407\pm 0.119 $. Estimates of higher order terms of the perturbative series are given. We compute the Bjorken sum rule with hyperasymptotic precision by including the leading terminant, associated with the first infrared renormalon. We fit the experimental data to the operator product expansion theoretical prediction with $\hat f_{3}^{\rm PV}$ as the free parameter. We obtain a good agreement with the experiment with $\hat f_{3}^{\rm PV}\times 10^3=32^{+187}_{-196}\;{\rm GeV}^{2}$ for $Q^2 \geq 1$ GeV$^2$., Comment: 25 pages, 8 Figures; v2: small numerical modification due to finite charm mass corrections. Physical results unchanged. References added. It matches the PRD version

We perform numerical analysis of double-pinched Borel-Laplace QCD sum rules for the strangeless semihadronic $\tau$-decay data. The $D=0$ contribution to the theoretical contour integral in the sum rules is evaluated by the (truncated) Fixed Order perturbation theory method (FO) and by the Principal Value (PV) of the Borel integration. We use for the full Adler function the Operator Product Expansion (OPE) with the terms $\sim \langle O_D \rangle$ of dimension $D=2 n$ where $2 \leq n \leq 5$ for the (V+A)-channel, and $2 \leq n \leq 7$ for the V-channel data. In our previous works [1,2], only the (V+A)-channel data was analysed. In this work, the analysis of a new set of V-channel data is performed as well. Further, a renormalon-motivated construction of the $D=0$ part of the Adler function is improved in the $u=3$ infrared renormalon sector, by involving the recently known information on the two principal noninteger values $k^{(j)}=\gamma^{(1)}(O_6^{(j)})/\beta_0$ of the effective leading-order anomalous dimensions. Additionally, the OPE of the Adler function has now the D=6 contribution with the principal anomalous dimension ($\sim \alpha_s^{k^{(1)}}$), and terms of higher dimension (with zero anomalous dimension). Cross-checks of the obtained extracted values of $\alpha_s$ and of the condensates were performed by reproduction of the (central) experimental values of several double-pinched momenta $a^{(2,n)}$. The averaged final extracted values of the (MSbar) coupling are: $\alpha_s(m_{\tau}^2) = 0.3169^{+0.0070}_{-0.0096}$, corresponding to $\alpha_s(M_Z^2)=0.1183^{+0.0009}_{-0.0012}$., Comment: 29 pages, 10 figures; v3: minor changes in the text, version as published in J.Phys.G

Double-pinched Borel-Laplace sum rules are applied to ALEPH $\tau$-decay data. For the leading-twist ($D=0$) Adler function a renormalon-motivated extension is used, and the 5-loop coefficient is taken to be $d_4=275 \pm 63$. Two $D=6$ terms appear in the truncated OPE ($D \leq 6$) to enable cancellation of the corresponding renormalon ambiguities. Two variants of the fixed order perturbation theory, and the inverse Borel transform, are applied to the evaluation of the $D=0$ contribution. Truncation index $N_t$ is fixed by the requirement of local insensitivity of the momenta $a^{(2,0)}$ and $a^{(2,1)}$ under variation of $N_t$. The averaged value of the coupling obtained is $\alpha_s(m_{\tau}^2)=0.3235^{+0.0138}_{-0.0126}$ [$\alpha_s(M_Z^2)=0.1191 \pm 0.0016$]. The theoretical uncertainties are significantly larger than the experimental ones., Comment: 4 pages, 1 figure, based on a presentation at: "alphas-2022: Workshop on precision measurements of the QCD coupling constant," January 31 - February 4, 2022, ECT* Trento, Italy; written for the Snowmass-2022 White Paper (The strong coupling constant: State of the art and the decade ahead)

We use improved truncated Operator Product Expansion (OPE) for the Adler function, involving two types of terms with dimension $D=6$, in the double-pinched Borel-Laplace Sum Rules and Finite Energy Sum Rules for the V+A channel strangeless semihadronic $\tau$ decays. The generation of the higher order perturbative QCD terms of the $D=0$ part of the Adler function is carried out using a renormalon-motivated ansatz incorporating the leading UV renormalon and the first two leading IR renormalons. The trunacted $D=0$ part of the Sum Rules is evaluated by two variants of the fixed-order perturbation theory (FO), by Principal Value of the Borel resummation (PV), and by contour-improved perturbation theory (CI). For the experimental V+A channel spectral function we use the ALEPH $\tau$-decay data. We point out that the truncated FO and PV evaluation methods account correctly for the renormalon structure of the Sum Rules, while this is not the case for the truncated CI evaluation. We extract the value of the ${\overline {\rm MS}}$ coupling $\alpha_s(m_{\tau}^2) = 0.3235^{+0.0138}_{-0.0126}$ [$\alpha_s(M_Z^2)=0.1191 \pm 0.0016$] for the average of the two FO methods and the PV method, which we consider as our main result. If we included in the average also CI extraction, the value would be $\alpha_s(m_{\tau}^2) = 0.3299^{+0.0232}_{-0.0225}$ [$\alpha_s(M_Z^2)=0.1199^{+0.0026}_{-0.0028}$]. This work is an extension and improvement of our previous work [Eur.Phys.J.C81 (2021) 10, 930] where we used for the truncated OPE a more naive (and widely used) form and where the extracted values for $\alpha_s(M_Z^2)$ were somewhat lower., Comment: v3: 20 pages, 6 figures; corrected list of references

We present a determination of the perturbative QCD (pQCD) coupling using the V+A channel ALEPH $\tau$-decay data. The determination involves the double-pinched Borel-Laplace Sum Rules and Finite Energy Sum Rules. The theoretical basis is the Operator Product Expansion (OPE) of the V+A channel Adler function in which the higher order terms of the leading-twist part originate from a model based on the known structure of the leading renormalons of this quantity. The applied evaluation methods are contour-improved perturbation theory (CIPT), fixed-order perturbation theory (FOPT), and Principal Value of the Borel resummation (PV). All the methods involve truncations in the order of the coupling. In contrast to the truncated CIPT method, the truncated FOPT and PV methods account correctly for the suppression of various renormalon contributions of the Adler function in the mentioned sum rules. The extracted value of the ${\overline {\rm MS}}$ coupling is $\alpha_s(m_{\tau}^2) = 0.3116 \pm 0.0073$ [$\alpha_s(M_Z^2)=0.1176 \pm 0.0010$] for the average of the FOPT and PV methods, which we regard as our main result. On the other hand, if we include in the average also the CIPT method, the resulting values are significantly higher, $\alpha_s(m_{\tau}^2) = 0.3194 \pm 0.0167$ [$\alpha_s(M_Z^2)=0.1186 \pm 0.0021$]., Comment: v3: 36 pages, 11 figures. Matched published version

We give the hyperasymptotic expansion of the plaquette with a precision that includes the terminant associated to the leading renormalon. Subleading effects are also considered. The perturbative series is regulated using the principal value prescription for its Borel integral. We use this analysis to give a determination of the gluon condensate in SU(3) pure gluodynamics that is independent of the scale and renormalization scheme used for the coupling constant: $\langle G^2 \rangle_{\rm PV} (n_f=0)=3.15(18)\; r_0^{-4}$., Comment: 20 pages, 9 figure. To meet journal version

We give the hyperasymptotic expansion of the energy of a static quark-antiquark pair with a precision that includes the effects of the subleading renormalon. The terminants associated to the first and second renormalon are incorporated in the analysis when necessary. In particular, we determine the normalization of the leading renormalon of the force and, consequently, of the subleading renormalon of the static potential. We obtain $Z_3^F(n_f=3)=2Z_3^V(n_f=3)=0.37(17)$. The precision we reach in strict perturbation theory is next-to-next-to-next-to-leading logarithmic resummed order both for the static potential and for the force. We find that the resummation of large logarithms and the inclusion of the leading terminants associated to the renormalons are compulsory to get accurate determinations of $\Lambda_{\overline{\rm MS}}$ when fitting to short-distance lattice data of the static energy. We obtain $\Lambda_{\overline{\rm MS}}^{(n_f=3)}=338(12)$ MeV and $\alpha(M_z)=0.1181(9)$. We have also found strong consistency checks that the ultrasoft correction to the static energy can be computed at weak coupling in the energy range we have studied., Comment: 37 pages, 14 figures. References added. Minor modification. To meet journal version

We evaluate QCD effects in the neutrinoless double beta ($0\nu\beta\beta$) decay, originating from new physics short-range mechanism in the form of five dimension-9 operators. For this, we employ the one-loop and two-loop renormalization group equations (RGEs) for the corresponding Wilson coefficients, performing the RGE-evolution from the new physics scales (estimated as $\Lambda ~ \sim 10^2$ GeV) to the typical spacelike $0\nu\beta\beta$-scale $Q \sim 0.1$ GeV. Since the latter scale is clearly nonperturbative, we apply various infrared-safe (IR-safe) variants of QCD where the running coupling has no Landau singularities at low spacelike $Q$. We point out that the correct treatment of the IR-safe analogs of the (noninteger) powers of the couplings is important. It turns out that in most cases of the considered operators the resulting QCD effects can be significant in this process, i.e., can be stronger than the effects of the present uncertainties in the nuclear matrix elements., Comment: 34 pages, 11 figures. v3: discussions added, principally in the new Sec. VI; results unchanged; version published in Phys. Rev. D; App.B is based partly on App.B of arXiv:1812.01580

These proceedings review recent work on hyperasymptotic constructions to the operator product expansion. Quantities we consider are the static potential and the pole mass., Comment: 9 pages, 7 figures, Proceedings for the 22th International Conference in Quantum Chromodynamics (QCD 19), 2-5 July 2019, Montpellier, France

We construct hyperasymptotic expansions for the heavy quark pole mass regulated using the principal value (PV) prescription. We apply such hyperasymptotic expansions to the $B/D$ meson masses, and $\bar \Lambda $ computed in the lattice. The issue of the uncertainty of the (top) pole mass is critically reexamined. The present theoretical uncertainty in the relation between ${\bar m}_t$, the $\bar{\rm MS}$ top mass, and $m_{t, \rm PV}$, the top pole mass regulated using the PV prescription, is numerically assessed to be $\delta m_{t,\rm PV}= 28\;{\rm MeV}$ for ${\bar m}_t =163$ GeV., Comment: 52 pages, 27 figures; v2: 56 pages, extra clarifications, one extra appendix, added references. Journal version

Given an observable and its operator product expansion (OPE), we present expressions that carefully disentangle truncated sums of the perturbative series in powers of $\alpha$ from the non-perturbative (NP) corrections. This splitting is done with NP power accuracy. Analytic control of the splitting is achieved and the organization of the different terms is done along an super/hyper-asymptotic expansion. As a test we apply the methods to the static potential in the large $\beta_0$ approximation. We see the superasymptotic and hyperasymptotic structure of the observable in full glory., Comment: 60 pages, 11 figures, 2 Tables. Minor changes to meet journal version. Conclusions unchanged

Experimental data obtained for the polarized Bjorken sum rule (BSR) $\Gamma_1^{p-n}(Q^2)$ are fitted by using predictions derived within a truncated operator product expansion (OPE) approach to QCD. Four QCD versions are considered: perturbative QCD (pQCD) in the MSbar scheme, Analytic Perturbation Theory (APT), and 2$\delta$ and 3$\delta$ analytic QCD versions. In contrast to pQCD, these QCD variants do not have Landau singularities at low positive $Q^2$, which facilitates the fitting procedure significantly. The fitting procedure is applied first to the experimental data of the inelastic part of BSR, and the known elastic contributions are added after the fitting. In general, when 2$\delta$ and 3$\delta$ QCD coupling is used the fitted curves give the best results, within the $Q^2$-range of the fit as well as in extended $Q^2$-intervals. When the fitting procedure is applied to the total BSR, i.e., to the sum of the experimental data and the elastic contribution, the quality of the results deteriorates significantly., Comment: 29 pages, 8 figures; to appear in EPJ C

We consider a phenomenologycal parametrization of the QCD running coupling which arises from the dispersion relation respecting the holomorphic properties of the physical QCD observables in the complex momentum plane. The parameters are fixed by the following requirements: 1) at enough high energies, it reproduces the underlying perturbative coupling, 2) at intermediate energy momenta, it reproduces the experimental semihadronic tau decay ratio, and 3) in the deep IR regime, it satisfies the qualitative properties coming from recent lattice results. Finally, we apply this new coupling to low-energy available experimental data. In particular, to Borel sum rules for {\tau}-decay, extracting the values of the dimension 4 and 6 condensates, to the V-channel Adler function, and to polarized Bjorken Sum Rule., Comment: 6 pages, 1 table, 5 figures. Contribution to the proceedings of the 20th International Conference in Quantum Chromodynamics (QCD 17), 3 July - 7 July 2017, Montpellier, France

We investigate the Bjorken polarized sum rule (BSR) in three approaches to QCD with analytic (holomorphic) coupling: Analytic Perturbation Theory (APT), Two-delta analytic QCD (2$\delta$anQCD), and Three-delta lattice-motivated analytic QCD in the three-loop and four-loop MOM scheme (3l3$\delta$anQCD, 4l3$\delta$anQCD). These couplings do not have unphysical (Landau) singularities, and have finite values when the transferred momentum goes to zero, which allows us to explore the infrared regime. With the exception of APT, these theories at high momenta practically coincide with the underlying perturbative QCD (pQCD) in the same scheme. We apply them in order to verify the Bjorken sum rule within the range of energies available in the data collected by the experimental JLAB collaboration, i.e., $0.05 GeV^2

The product of the gluon dressing function and the square of the ghost dressing function in the Landau gauge can be regarded to represent, apart from the inverse power corrections 1/Q^{2n}, a nonperturbative generalization A(Q^2) of the perturbative QCD running coupling a(Q^2)=alpha_s(Q^2)/pi. Recent large volume lattice calculations for these dressing functions indicate that the coupling defined in such a way goes to zero as A(Q^2)~Q^2 when the squared momenta Q^2 go to zero (Q^2<<1 GeV^2). In this work we construct such a QCD coupling A(Q^2) which fulfills also various other physically motivated conditions. At high momenta it becomes the underlying perturbative coupling a(Q^2) to a very high precision. And at intermediate low squared momenta Q^2~1 GeV^2 it gives results consistent with the data of the semihadronic tau lepton decays as measured by OPAL and ALEPH. The coupling is constructed in a dispersive way, resulting as a byproduct in the holomorphic behavior of A(Q^2) in the complex Q^2-plane which reflects the holomorphic behavior of the spacelike QCD observables. Application of the Borel sum rules to tau-decay V+A spectral functions allows us to obtain values for the gluon (dimension-4) condensate and the dimension-6 condensate, which reproduce the measured OPAL and ALEPH data to a significantly better precision than the perturbative MSbar coupling approach., Comment: 43 pages, 15 figures, v4: improved presentation in the text, results unchanged; new references added ([7-8,28-29,40,44-45,47,52-53,127-128]); to appear in J.Phys.G

We compute the right top quark coupling in the aligned two-Higgs-doublet model. In the Standard Model the real part of this coupling is dominated by QCD-gluon-exchange diagram, but the imaginary part, instead, is purely electroweak at one loop. Within this model we show that values for the imaginary part of the coupling up to one order of magnitude larger than the electroweak prediction can be obtained. For the real part of the electroweak contribution we find that it can be up to three orders of magnitude larger than the standard model one. We also present detailed results of the one loop analytical computation.

Newer lattice results indicate that, in the Landau gauge at low spacelike momenta, the gluon propagator and the ghost dressing function are finite nonzero. This leads to a definition of the QCD running coupling, in a specific scheme, that goes to zero at low spacelike momenta. We construct a running coupling which fulfills these conditions, and at the same time reproduces to a high precision the perturbative behavior at high momenta. The coupling is constructed in such a way that it reflects qualitatively correctly the holomorphic (analytic) behavior of spacelike observables in the complex plane of the squared momenta, as dictated by the general principles of Quantum Field Theories. Further, we require the coupling to reproduce correctly the nonstrange semihadronic decay rate of tau lepton which is the best measured low-momentum QCD observable with small higher-twist effects. Subsequent application of the Borel sum rules to the V+A spectral functions of tau lepton decays, as measured by OPAL Collaboration, determines the values of the gluon condensate and of the V+A 6-dimensional condensate, and reproduces the data to a significantly higher precision than the usual MSbar running coupling., Comment: v4: 23 pages, 9 figures; improved discussion; new Sec.V and new references added; version to appear in J.Phys.G

We update our perturbative determination of MSbar bottom quark mass mb(mb), by including the recently obtained four-loop coefficient in the relation between the pole and MSbar mass. First the renormalon subtracted (RS or RS') mass is determined from the known mass of the Upsilon(1S) meson, where we use the renormalon residue Nm obtained from the asymptotic behavior of the coefficient of the 3-loop static singlet potential. MSbar mass is then obtained using the 4-loop renormalon-free relation between the RS (RS') and MSbar mass. We argue that the effects of the charm quark mass are accounted for by effectively using Nf=3 in the mass relations. The extracted value is mb(mb) = 4222(40) MeV, where the uncertainty is dominated by the renormalization scale dependence., Comment: 6 pages; presented at ACAT 2016, 18-22 January 2016, UTFSM, Valparaiso, Chile

We implement our past investigations in the quark-antiquark interaction through a non-perturbative running coupling defined in terms of a gluon mass function, similar to that used in some Schwinger-Dyson approaches. This coupling leads to a quark-antiquark potential, which satisfies not only asymptotic freedom but also describes linear confinement correctly. From this potential, we calculate the bottomonium and charmonium spectra below the first open flavor meson-meson thresholds and show that for a small range of values of the free parameter determining the gluon mass function an excellent agreement with data is attained., Comment: 13 pages, 4 figures; v3: version accepted for publication in Journal of Physics G

We provide three Fortran programs which evaluate the QCD analytic (holomorphic) couplings $\mathcal{A}_{\nu}(Q^2)$ for complex or real squared momenta $Q^2$. These couplings are holomorphic analogs of the powers $a(Q^2)^{\nu}$ of the underlying perturbative QCD (pQCD) coupling $a(Q^2) \equiv \alpha_s(Q^2)/\pi$, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2$\delta$anQCD), and Massive Perturbation Theory (MPT). The index $\nu$ can be noninteger. The provided programs do basically the same job as the Mathematica package anQCD.m in Mathematica published by us previously, Ref.[1], but are now written in Fortran., Comment: 7 pages; the programs can be downloaded from the web page gcvetic.usm.cl

We apply (Fractional) Analytic Perturbation Theory (FAPT) to the QCD analysis of the nonsinglet nucleon structure function $F_2(x,Q^2)$ in deep inelastic scattering up to the next leading order and compare the results with ones obtained within the standard perturbation QCD. Based on a popular parameterization of the corresponding parton distribution we perform the analysis within the Jacobi Polynomial formalism and under the control of the numerical inverse Mellin transform. To reveal the main features of the FAPT two-loop approach, we consider a wide range of momentum transfer from high $Q^2\sim 100 {\rm GeV}^2$ to low $Q^2\sim 0.3 {\rm GeV}^2$ where the approach still works., Comment: 14 pages, 4 figures, 1 table; v3: Improve clarity and some details, references are also added. Version accepted for publication in PRD

We outline here the motivation for the existence of analytic QCD models, i.e., QCD frameworks in which the running coupling $A(Q^2)$ has no Landau singularities. The analytic (holomorphic) coupling $A(Q^2)$ is the analog of the underlying pQCD coupling $a(Q^2) \equiv \alpha_s(Q^2)/\pi$, and any such $A(Q^2)$ defines an analytic QCD model. We present the general construction procedure for the couplings $A_{\nu}(Q^2)$ which are analytic analogs of the powers $a(Q^2)^{\nu}$. Three analytic QCD models are presented. Applications of our program (in Mathematica) for calculation of $A_{\nu}(Q^2)$ in such models are presented. Programs in both Mathematica and Fortran can be downloaded from the web page: gcvetic.usm.cl., Comment: 7 pages, 2 figures; uses JPConf style; based on the presentation given by G.C. at the 16th International workshop on Advanced Computing and Analysis Techniques in physics research (ACAT 2014), Prague, Czech Republic, September 1-5, 2014; to appear in the proceedings by the IOP Conference Series publishing

We provide a Mathematica package that evaluates the QCD analytic couplings (in the Euclidean domain) $\mathcal{A}_{\nu}(Q^2)$, which are analytic analogs of the powers $a(Q^2)^{\nu}$ of the underlying perturbative QCD (pQCD) coupling $a(Q^2) \equiv \alpha_s(Q^2)/\pi$, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2$\delta$anQCD), and Massive Perturbation Theory (MPT). The analytic (holomorphic) running couplings $\mathcal{A}_{\nu}(Q^2)$, in contrast to the corresponding pQCD expressions $a(Q^2)^{\nu}$, reflect correctly the analytic properties of the spacelike observables ${\cal D}(Q^2)$ in the complex $Q^2$ plane as dictated by the general principles of quantum field theory. They are thus more suited for evaluations of such physical quantities, especially at low momenta $|Q^2| \sim 1 \ {\rm GeV}^2$., Comment: v3: 35 pp., 7 figs.; improved programs which work in Mathematica 9 and 10; improved presentation of the text (Sec.2 shorter, Sec.4 enlarged); to appear in Computer Physics Communications

We obtain an improved determination of the normalization constant of the first infrared renormalon of the pole mass (and the singlet static potential). For $N_f=3$ it reads $N_m=0.563(26)$. Charm quark effects in the bottom quark mass determination are carefully investigated. Finally, we determine the bottom quark mass using the NNNLO perturbative expression for the $\Upsilon(1S)$ mass. We work in the renormalon subtracted scheme, which allows us to control the divergence of the perturbation series due to pole mass renormalon. Our result for the ${\overline {\rm MS}}$ mass reads ${\overline m}_{b}({\overline m}_{b})=4201(43)$ MeV., Comment: 37 pages, 11 figures, minor changes in presentation, version to appear in JHEP

We consider the gluon propagator in the Landau gauge at low spacelike momenta and with the dressing function $Z(Q^2)$ at the two-loop order. We incorporate the nonperturbative effects by making the (noninteger) powers of the QCD coupling in the dressing function $Z(Q^2)$ analytic (holomorphic) via the Fractional Analytic Perturbation Theory (FAPT) model, and simultaneously introducing the gluon dynamical mass in the propagator as motivated by the previous analyses of the Dyson-Schwinger equations. The obtained propagator has behavior compatible with the unquenched lattice data ($N_f=2+1$) at low spacelike momenta $0.4 \ {\rm GeV} < Q \lesssim 10$ GeV. We conclude that the removal of the unphysical Landau singularities of the powers of the coupling via the (F)APT prescription, in conjunction with the introduction of the dynamical mass $M \approx 0.62$ GeV of the gluon, leads to an acceptable behavior of the propagator in the infrared regime., Comment: 9 pages, 3 figures

We extract MSbar quark masses mqbar (q=b,c) from the evaluation of the masses of quarkonia Upsilon(1S) and J/\psi(1S), performed in two analytic QCD models, and in perturbative QCD in two renormalization schemes. In analytic QCD the running coupling has no unphysical singularities in the low-momentum regime. We apply the analytic model of Shirkov et al.[Analytic Perturbation Theory (APT)], extended by Bakulev et al.[Fractional Analytic Perturbation Theory (FAPT)], and the two-delta analytic model (2delanQCD). The latter, in contrast to (F)APT, at higher energies basically coincides with the perturbative QCD (in the same scheme). We use the renormalon-free mass mqbar as input. The separation of the soft and ultrasoft parts of the binding energy E is performed by the requirement of the cancellation of the leading infrared renormalon. The analysis in the 2delanQCD model indicates that the low-momentum ultrasoft regime is important for the extraction of the masses mqbar, especially mbbar. The 2delanQCD model gives us clues on how to estimate the influence of the ultrasoft sector on mqbar in general. These effects lead to relatively large values mbbar(mbbar) =approx (4.35 +- 0.08) GeV in the 2delanQCD model, which, however, are compatible with recent lattice calculations. In perturbative QCD in MSbar scheme these effects are even stronger and give larger uncertainties in mbbar. The (F)APT model gives small ultrasoft effects and the extracted values of mbbar agree with those in most of the literature: mbbar(mbbar) =approx 4.2 GeV. The extracted values of mcbar(mcbar) in all four models are about 1.26-1.27 GeV and agree well with those in the literature., Comment: 31 pages, 4 figures; in comparison with v1: improved presentation; (F)APT calculation performed in the beta_2=beta_3=...=0 scheme; comparison with the results of other methods in the literature included (cf.Table IV); part of the text moved into the new Appendix A; the basic conclusions unchanged; extended acknowledgments; Refs.~[12,38-39,41,44,80-83,89-99,110-11] are new

In contrast to perturbative QCD, the analytic QCD models have running coupling whose analytic properties correctly mirror those of spacelike observables. The discontinuity (spectral) function of such running coupling is expected to agree with the perturbative case at large timelike momenta; however, at low timelike momenta it is not known. In the latter regime, we parametrize the unknown behavior of the spectral function as a sum of (two) delta functions; while the onset of the perturbative behavior of the spectral function is set to be 1.0-1.5 GeV. This is in close analogy with the "minimal hadronic ansatz" used in the literature for modeling spectral functions of correlators. For the running coupling itself, we impose the condition that it basically merges with the perturbative coupling at high spacelike momenta. In addition, we require that the well-measured nonstrange semihadronic (V+A) tau decay ratio value be reproduced by the model. We thus obtain a QCD framework which is basically indistinguishable from perturbative QCD at high momenta (Q > 1 GeV), and at low momenta it respects the basic analyticity properties of spacelike observables as dictated by the general principles of the local quantum field theories., Comment: 15 pages, 6 figures; in v2 Sec.IV is extended after Eq.(48) and refs.[51-52] added; v2 published in Phys.Rev.D85,114043(2012)