We provide a Mathematica package that evaluates the QCD analytic couplings (in the complex domain) A ν (Q 2 ), which are analytic analogs of the powers a(Q 2 ) ν of the underlying perturbative QCD (pQCD) coupling a(Q 2 ) ≡ α s (Q 2 )/π, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2δanQCD), and Massive Perturbation Theory (MPT). The analytic (holomorphic) running couplings A ν (Q 2 ), in contrast to the corresponding pQCD expressions a(Q 2 ) ν , reflect correctly the analytic properties of the spacelike observables 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 | ∼ 1 GeV 2 .
Program Summary
Title of program: anQCDThe main program ( anQCD.m) and supplementary modules ( Li nu.m and s0r.m), and the zipped file containing all three files ( anQCD Mathematica.zip), available from the web page:No. of bytes in distributed program including test data etc.: 63 kB (main module anQCD.m), 2 kB (supplementary module Li nu.m), 18 kB (supplementary module s0r.m);Nature of the physical problem: Evaluation of the values for analytic couplings A ν (Q 2 ; N f ) in analytic QCD [the analytic analog of the power (α s (Q 2 ; N f )/π) ν ] based on the dispersion relation; A ν represents a physical (holomorphic) function in the plane of complex squared momenta −q 2 ≡ Q 2 . In anQCD.m we collect the formulas for three different analytic models depending on the energy scale, Q 2 , number of flavors N f , the QCD scale Λ N f , and the (nonpower) index ν. The considered models are: Analytic Perturbation theory (APT), Two-delta analytic QCD (2δanQCD) and Massive Perturbation Theory (MPT).
Method of solution:anQCD uses Mathematica functions to perform numerical integration of spectral function for each analytic model, in order to obtain the corresponding analytic images A ν (Q 2 ) via dispersion relation.Restrictions on the complexity of the problem: It could be that for an unphysical choice of the input parameters the results are meaningless. Typical running time: For all operations the running time does not exceed a few seconds.