• We present the complete set of vertex, wave function and charge renormalisation constants in QCD in a general simple gauge group and with the complete dependence on the covariant gauge parameter $\xi$ in the minimal subtraction scheme of conventional dimensional regularisation. Our results confirm all already known results, which were obtained in the Feynman gauge, and allow the extraction of other useful gauges such as the Landau gauge. We use these results to extract the Landau gauge five-loop anomalous dimensions of the composite operator $A^2$ as well as the Landau gauge scheme independent gluon, ghost and fermion propagators at five loops.
  • We have computed the self-energies and a set of three-particle vertex functions for massless QCD at the four-loop level in the MSbar renormalization scheme. The vertex functions are evaluated at points where one of the momenta vanishes. Analytical results are obtained for a generic gauge group and with the full gauge dependence, which was made possible by extensive use of the Forcer program for massless four-loop propagator integrals. The bare results in dimensional regularization are provided in terms of master integrals and rational coefficients; the latter are exact in any space-time dimension. Our results can be used for further precision investigations of the perturbative behaviour of the theory in schemes other than MSbar. As an example, we derive the five-loop beta function in a relatively common alternative, the minimal momentum subtraction (MiniMOM) scheme.
  • We explain the construction of Forcer, a FORM program for the reduction of four-loop massless propagator-type integrals to master integrals. The resulting program performs parametric IBP reductions similar to the three-loop Mincer program. We show how one can solve many systems of IBP identities parametrically in a computer-assisted manner. Next, we discuss the structure of the Forcer program, which involves recognizing reduction actions for each topology, applying symmetries, and transitioning between topologies after edges have been removed. This part is entirely precomputed and automatically generated. We give examples of recent applications of Forcer, and study the performance of the program. Finally we demonstrate how to use the Forcer package and sketch how to prepare physical diagrams for evaluation by Forcer.
  • We have computed the five-loop corrections to the scale dependence of the renormalized coupling constant for Quantum Chromodynamics (QCD), its generalization to non-Abelian gauge theories with a simple compact Lie group, and for Quantum Electrodynamics (QED). Our analytical result, obtained using the background field method, infrared rearrangement via a new diagram-by-diagram implementation of the R* operation and the Forcer program for massless four-loop propagators, confirms the QCD and QED results obtained by only one group before. The numerical size of the five-loop corrections is briefly discussed in the standard MSbar scheme for QCD with n_f flavours and for pure SU(N) Yang-Mills theory. Their effect in QCD is much smaller than the four-loop contributions, even at rather low scales.
  • We have computed the fourth-order nf^2 contributions to all three non-singlet quark-quark splitting functions and their four nf^3 flavour-singlet counterparts for the evolution of the parton distributions of hadrons in perturbative QCD with nf effectively massless quark flavours. The analytic form of these functions is presented in both Mellin N-space and momentum-fraction x-space; the large-x and small-x limits are discussed. Our results agree with all available predictions derived from lower-order information. The large-x limit of the quark-quark cases provides the complete nf^2 part of the four-loop cusp anomalous dimension which agrees with two recent partial computations.
  • We present a new FORM program for analytically evaluating four-loop massless propagator-type Feynman integrals in an efficient way. Our program Forcer implements parametric reductions of the aforementioned class of Feynman integrals into a set of master integrals and can be considered as a four-loop extension of the three-loop Mincer program. Since the program structure at the four-loop level is highly complicated and the equations easily become lengthy, most of the code was generated in an automatic way or with computer-assisted derivations. We have checked correctness of the program by recomputing already-known quantities in the literature.
  • We have calculated the coefficient functions for the structure functions F_2, F_L and F_3 in nu-nubar charged-current deep-inelastic scattering (DIS) at the third order in the strong coupling alpha_s, thus completing the description of unpolarized inclusive W^(+-) exchange DIS to this order of massless perturbative QCD. In this brief note, our new results are presented in terms of compact approximate expressions that are sufficiently accurate for phenomenological analyses. For the benefit of such analyses we also collect, in a unified notation, the corresponding lower-order contributions and the flavour non-singlet coefficient functions for nu+nubar charged-current DIS. The behaviour of all six third-order coefficient functions at small Bjorken-x is briefly discussed.
  • We present version two of the Latex graphical style file Axodraw. It has a number of new drawing primitives and many extra options, and it can now work with \program{pdflatex} to directly produce output in PDF file format (but with the aid of an auxiliary program).
  • We present results on the fourth-order splitting functions and coefficient functions obtained using Forcer, a four-loop generalization of the Mincer program for the parametric reduction of self-energy integrals. We have computed the respective lowest three even-N and odd-N moments for the non-singlet splitting functions and the non-singlet coefficient functions in electromagnetic and nu+nu(bar) charged-current deep-inelastic scattering, and the N=2 and N=4 results for the corresponding flavour-singlet quantities. Enough moments have been obtained for an LLL-based determination of the analytic N-dependence of the nf^3 and nf^2 parts, respectively, of the singlet and non-singlet splitting functions. The large-N limit of the latter provides the complete nf^2 contributions to the four-loop cusp anomalous dimension. Our results also provide additional evidence of a non-vanishing contribution of quartic group invariants to the cusp anomalous dimension.
  • We present Forcer, a new FORM program for the calculation of four-loop massless propagators. The basic framework is similar to that of the Mincer program for three-loop massless propagators: the program reduces Feynman integrals to a set of master integrals in a parametric way. To overcome an ineludible complexity of the program structure at the four-loop level, most of the code was automatically generated or made with computer-assisted derivations. Correctness of the program has been checked with the recomputation of some quantities in the literature.
  • We discuss the prescription for the Dirac matrix gamma_5 in dimensional regularization used in most second- and third-order QCD calculations of collider cross sections. We provide an alternative implementation of this approach that avoids the use of an explicit form of gamma_5 and of its (anti-) commutation relations in the most important case of no more than one gamma_5 in each fermion trace. This treatment is checked by computing the third-order corrections to the structure functions F_2 and g_1 in charged-current deep-inelastic scattering with axial-vector couplings to the W-bosons. We derive the so far unknown third-order helicity-difference splitting function Delta P_ns^(2)s that contributes to the next-to-next-to-leading order (NNLO) evolution of the polarized valence quark distribution of the nucleon. This function is negligible at momentum fractions x >~ 0.3 but relevant at x << 1.
  • We present the next-to-next-to-leading order (NNLO) contributions to the main splitting functions for the evolution of longitudinally polarized parton densities of hadrons in perturbative QCD. The quark-quark and gluon-quark splitting functions have been obtained by extending our previous all Mellin-N calculations to the structure function g_1 in electromagnetic deep-inelastic scattering (DIS). Their quark-gluon and gluon-gluon counterparts have been derived using third-order fixed-N calculations of structure functions in graviton-exchange DIS, relations to the unpolarized case and mathematical tools for systems of Diophantine equations. The NNLO corrections to the splitting functions are small outside the region of small momentum fractions x where they exhibit a large double-logarithmic enhancement, yet the corrections to the evolution of the parton densities can be unproblematic down to at least x about 10^{-4}.
  • By using the GRACE-Loop system, we calculate the full $\mathcal{O}(\alpha)$ electroweak radiative corrections to the process $e^+e^- \rightarrow e^+e^- \gamma$, which is important for future investigations at the International Linear Collider (ILC). With the GRACE-Loop system, the calculations are checked numerically by three consistency tests: ultraviolet finiteness, infrared finiteness, and gauge-parameter independence. The results show good numerical stability when quadruple precision is used. In the phenomenological results, we find that the electroweak corrections to the total cross section range from $\sim -4\%$ to $\sim -21\%$ when $\sqrt{s}$ varies from $250$ GeV to $1$ TeV. The corrections also significantly affect the differential cross sections, which are a function of the invariant masses and angles and the final-particle energies. Such corrections will play an important role for the high-precision program at the ILC.
  • We have calculated the complete matrix of three-loop helicity-difference (`polarized') splitting functions Delta P_ik^(2), i,k = q,g, in massless perturbative QCD. In this note we briefly discuss some properties of the polarized splitting functions and our non-standard determination of the hitherto missing lower-row quantities Delta P_gq^(2) and Delta P_gg^(2). The resulting next-to-next-to-leading order (NNLO) corrections to the evolution of polarized parton distributions are illustrated and found to be small even at rather large values of the strong coupling constant alpha_s.
  • The full $\mathcal{O}(\alpha)$ electroweak radiative corrections to $t \bar{t} \gamma$ and $e^-e^+\gamma$ productions at the International Linear Collider (ILC) are presented in this paper. The computation is performed with the help of GRACE-Loop system. In the physical results, we discuss on the cross section, electroweak corrections, and the top quark forward-backward asymmetry ($A_{FB}$) which are the function of the center-of-mass energy.
  • We describe the implementation of output code optimization in the open source computer algebra system FORM. This implementation is based on recently discovered techniques of Monte Carlo tree search to find efficient multivariate Horner schemes, in combination with other optimization algorithms, such as common subexpression elimination. For systems for which no specific knowledge is provided it performs significantly better than other methods we could compare with. Because the method has a number of free parameters, we also show some methods by which to tune them to different types of problems.
  • We present the full $\mathcal{O}(\alpha)$ electroweak radiative corrections to the process $e^+e^- \rightarrow t \bar{t} \gamma$ at the International Linear Collider (ILC). The computation is performed with the help of the GRACE-Loop system. We present the total cross-section and the top quark forward-backward asymmetry ($A_{FB}$) as a function of the center-of-mass energy and compare them with the process $e^+e^- \rightarrow t \bar{t}$. We find that the value of $A_{FB}$ in $t \bar{t} \gamma$ production is larger than $A_{FB}$ in $t\bar{t}$ production. It is an important result for the measurement of the top quark forward-backward asymmetry at the ILC. Applying a structure function method, we also subtract the QED correction to gain the genuine weak correction in both the $\alpha$ scheme and the $G_{\mu}$ scheme ($\delta_{W}^{G_{\mu}}$). We obtain numerical values for $\delta_{W}^{G_{\mu}}$ which are changing from 2% to -24% when we vary the center-of-mass energy from 360 GeV to 1 TeV.
  • We present version 4.0 of the symbolic manipulation system FORM. The most important new features are manipulation of rational polynomials and the factorization of expressions. Many other new functions and commands are also added; some of them are very general, while others are designed for building specific high level packages, such as one for Groebner bases. New is also the checkpoint facility, that allows for periodic backups during long calculations. Lastly, FORM 4.0 has become available as open source under the GNU General Public License version 3.
  • We discuss the large-x behaviour of the splitting functions P_qg and P_gq and of flavour-singlet coefficient functions, such as the gluon contributions C_2,g and C_L,g to the structure functions F_2,L, in massless perturbative QCD. These quantities are suppressed by one or two powers of 1-x with respect to the 1/(1-x) terms which are the subject of the well-known threshold exponentiation. We show that the double-logarithmic contributions to P_qg, P_gq and C_L at order alpha_s^4 can be predicted from known third-order results and present, as a first step towards a full all-order generalization, the leading-logarithmic large-x behaviour of P_qg, P_gq and C_2,g at all orders in alpha_s.
  • The symbolic manipulation program FORM is specialized to handle very large algebraic expressions. Some specific features of its internal structure make FORM very well suited for parallelization. We have now two parallel versions of FORM, one is based on POSIX threads and is optimal for modern multicore computers while another one uses MPI and can be used to parallelize FORM on clusters and Massive Parallel Processing systems. Most existing FORM programs will be able to take advantage of the parallel execution without the need for modifications.
  • We have studied the physical evolution kernels for nine non-singlet observables in deep-inelastic scattering (DIS), semi-inclusive e^+e^-annihilation and the Drell-Yan (DY) process, and for the flavour-singlet case of the photon- and heavy-top Higgs-exchange structure functions (F_2, F_phi) in DIS. All known contributions to these kernels show an only single-logarithmic large-x enhancement at all powers of 1-x. Conjecturing that this behaviour persists to (all) higher orders, we have predicted the highest three (DY: two) double logarithms of the higher-order non-singlet coefficient functions and of the four-loop singlet splitting functions. The coefficient-function predictions canbe written as exponentiations of 1/N-suppressed contributions in Mellin-N space which, however, are less predictive than the well-known exponentiation of the ln^k N terms.
  • We provide a data mine of proven results for multiple zeta values (MZVs) of the form $\zeta(s_1,s_2,...,s_k)=\sum_{n_1>n_2>...>n_k>0}^\infty \{1/(n_1^{s_1} >... n_k^{s_k})\}$ with weight $w=\sum_{i=1}^k s_i$ and depth $k$ and for Euler sums of the form $\sum_{n_1>n_2>...>n_k>0}^\infty t\{(\epsilon_1^{n_1} >...\epsilon_1 ^{n_k})/ (n_1^{s_1} ... n_k^{s_k}) \}$ with signs $\epsilon_i=\pm1$. Notably, we achieve explicit proven reductions of all MZVs with weights $w\le22$, and all Euler sums with weights $w\le12$, to bases whose dimensions, bigraded by weight and depth, have sizes in precise agreement with the Broadhurst--Kreimer and Broadhurst conjectures. Moreover, we lend further support to these conjectures by studying even greater weights ($w\le30$), using modular arithmetic. To obtain these results we derive a new type of relation for Euler sums, the Generalized Doubling Relations. We elucidate the "pushdown" mechanism, whereby the ornate enumeration of primitive MZVs, by weight and depth, is reconciled with the far simpler enumeration of primitive Euler sums. There is some evidence that this pushdown mechanism finds its origin in doubling relations. We hope that our data mine, obtained by exploiting the unique power of the computer algebra language {\sc form}, will enable the study of many more such consequences of the double-shuffle algebra of MZVs, and their Euler cousins, which are already the subject of keen interest, to practitioners of quantum field theory, and to mathematicians alike.
  • We notice that the existing literature about the reaction e+e- --> hh suffers from a mistake in the relative sign between the t-channel and u-channel graphs. Correcting this lowers the crosssections by about an order of magnitude.
  • We report on the first calculation of the structure function g_1 in polarised deep-inelastic scattering to the third order in massless perturbative QCD. The calculation follows the dispersive approach already used for the corresponding unpolarised cases of F_2,L, but additionally involves higher tensor integrals and the Dirac matrix gamma_5 in D unequal 4 dimensions. Our results confirm all known two-loop expressions including the coefficient functions of Zijlstra and van Neerven not independently verified before. At three loops we extract the helicity-difference next-to-next-to-leading order (NNLO) quark-quark and gluon-quark splitting functions Delta P_qq and Delta P_qg. The results exhibit interesting features concerning sum rules and the momentum-fraction limits x to 1 and x to 0.
  • The FORM project (0806.4080)

    June 25, 2008 hep-ph
    The necessity of the FORM project is discussed. Then the evolutionary needs in particle physics are considered, looking at the trends over the years. A guess is made at what will be needed in the (near) future. The whole is concluded with some critical remarks concerning the publication of results and programs.