
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 fiveloop 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 selfenergies and a set of threeparticle vertex
functions for massless QCD at the fourloop 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 fourloop propagator integrals. The bare results in
dimensional regularization are provided in terms of master integrals and
rational coefficients; the latter are exact in any spacetime 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 fiveloop 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
fourloop massless propagatortype integrals to master integrals. The resulting
program performs parametric IBP reductions similar to the threeloop Mincer
program. We show how one can solve many systems of IBP identities
parametrically in a computerassisted 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 fiveloop corrections to the scale dependence of the
renormalized coupling constant for Quantum Chromodynamics (QCD), its
generalization to nonAbelian 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
diagrambydiagram implementation of the R* operation and the Forcer program
for massless fourloop propagators, confirms the QCD and QED results obtained
by only one group before. The numerical size of the fiveloop corrections is
briefly discussed in the standard MSbar scheme for QCD with n_f flavours and
for pure SU(N) YangMills theory. Their effect in QCD is much smaller than the
fourloop contributions, even at rather low scales.

We have computed the fourthorder nf^2 contributions to all three nonsinglet
quarkquark splitting functions and their four nf^3 flavoursinglet
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 Nspace and momentumfraction
xspace; the largex and smallx limits are discussed. Our results agree with
all available predictions derived from lowerorder information. The largex
limit of the quarkquark cases provides the complete nf^2 part of the fourloop
cusp anomalous dimension which agrees with two recent partial computations.

We present a new FORM program for analytically evaluating fourloop massless
propagatortype 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 fourloop
extension of the threeloop Mincer program. Since the program structure at the
fourloop level is highly complicated and the equations easily become lengthy,
most of the code was generated in an automatic way or with computerassisted
derivations. We have checked correctness of the program by recomputing
alreadyknown quantities in the literature.

We have calculated the coefficient functions for the structure functions F_2,
F_L and F_3 in nunubar chargedcurrent deepinelastic 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 lowerorder contributions and the flavour
nonsinglet coefficient functions for nu+nubar chargedcurrent DIS. The
behaviour of all six thirdorder coefficient functions at small Bjorkenx 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 fourthorder splitting functions and coefficient
functions obtained using Forcer, a fourloop generalization of the Mincer
program for the parametric reduction of selfenergy integrals. We have computed
the respective lowest three evenN and oddN moments for the nonsinglet
splitting functions and the nonsinglet coefficient functions in
electromagnetic and nu+nu(bar) chargedcurrent deepinelastic scattering, and
the N=2 and N=4 results for the corresponding flavoursinglet quantities.
Enough moments have been obtained for an LLLbased determination of the
analytic Ndependence of the nf^3 and nf^2 parts, respectively, of the singlet
and nonsinglet splitting functions. The largeN limit of the latter provides
the complete nf^2 contributions to the fourloop cusp anomalous dimension. Our
results also provide additional evidence of a nonvanishing contribution of
quartic group invariants to the cusp anomalous dimension.

We present Forcer, a new FORM program for the calculation of fourloop
massless propagators. The basic framework is similar to that of the Mincer
program for threeloop 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 fourloop level, most of
the code was automatically generated or made with computerassisted
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 thirdorder 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 thirdorder
corrections to the structure functions F_2 and g_1 in chargedcurrent
deepinelastic scattering with axialvector couplings to the Wbosons. We
derive the so far unknown thirdorder helicitydifference splitting function
Delta P_ns^(2)s that contributes to the nexttonexttoleading 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 nexttonexttoleading order (NNLO) contributions to the main
splitting functions for the evolution of longitudinally polarized parton
densities of hadrons in perturbative QCD. The quarkquark and gluonquark
splitting functions have been obtained by extending our previous all MellinN
calculations to the structure function g_1 in electromagnetic deepinelastic
scattering (DIS). Their quarkgluon and gluongluon counterparts have been
derived using thirdorder fixedN calculations of structure functions in
gravitonexchange 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 doublelogarithmic 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 GRACELoop 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 GRACELoop system, the calculations are checked
numerically by three consistency tests: ultraviolet finiteness, infrared
finiteness, and gaugeparameter 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 finalparticle
energies. Such corrections will play an important role for the highprecision
program at the ILC.

We have calculated the complete matrix of threeloop helicitydifference
(`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 nonstandard determination of the
hitherto missing lowerrow quantities Delta P_gq^(2) and Delta P_gg^(2). The
resulting nexttonexttoleading 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 GRACELoop system. In the physical results, we discuss on the cross
section, electroweak corrections, and the top quark forwardbackward asymmetry
($A_{FB}$) which are the function of the centerofmass 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
GRACELoop system. We present the total crosssection and the top quark
forwardbackward asymmetry ($A_{FB}$) as a function of the centerofmass
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 forwardbackward 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 centerofmass 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 largex behaviour of the splitting functions P_qg and P_gq and
of flavoursinglet 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 1x with respect to the
1/(1x) terms which are the subject of the wellknown threshold exponentiation.
We show that the doublelogarithmic contributions to P_qg, P_gq and C_L at
order alpha_s^4 can be predicted from known thirdorder results and present, as
a first step towards a full allorder generalization, the leadinglogarithmic
largex 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 nonsinglet
observables in deepinelastic scattering (DIS), semiinclusive
e^+e^annihilation and the DrellYan (DY) process, and for the flavoursinglet
case of the photon and heavytop Higgsexchange structure functions (F_2,
F_phi) in DIS. All known contributions to these kernels show an only
singlelogarithmic largex enhancement at all powers of 1x. Conjecturing that
this behaviour persists to (all) higher orders, we have predicted the highest
three (DY: two) double logarithms of the higherorder nonsinglet coefficient
functions and of the fourloop singlet splitting functions. The
coefficientfunction predictions canbe written as exponentiations of
1/Nsuppressed contributions in MellinN space which, however, are less
predictive than the wellknown 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 BroadhurstKreimer 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 doubleshuffle 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 tchannel and uchannel 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
deepinelastic 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 twoloop expressions including the coefficient functions of Zijlstra and
van Neerven not independently verified before. At three loops we extract the
helicitydifference nexttonexttoleading order (NNLO) quarkquark and
gluonquark splitting functions Delta P_qq and Delta P_qg. The results exhibit
interesting features concerning sum rules and the momentumfraction limits x to
1 and x to 0.

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.