BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
J. A. M. Vermaseren appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We report results for the lowest even-$N$ moments of the flavor-nonsinglet structure functions $F_2$ and $F_L$ in QCD at the fourth order in the perturbative expansion in the strong coupling constant $α_s$. Our results are presented in numerical form and we compare them with the leading and subleading terms of the threshold expansion for large values of $N$, which corresponds to the limit $x \to 1$.
We have computed the four lowest even-N moments of all four splitting functions for the evolution of flavour-singlet parton densities of hadrons at the fourth order in the strong coupling constant alpha_s. The perturbative expansion of these moments, and hence of the splitting functions for momentum fractions x >~ 0.1, is found to be well behaved with relative alpha_s-coefficients of order one and sub-percent effects on the scale derivatives of the quark and gluon distributions at alpha_s ~< 0.2. More intricate computations, including other approaches such as the operator-product expansion, are required to cover the full x-range relevant to LHC analyses. Our results are presented analytically for a general gauge group for detailed checks and validations of such future calculations.
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 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 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 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 discuss a number of FORM features that are essential in the automatic processing of very large numbers of diagrams as used in the Forcer program for 4-loop massless propagator diagrams. Most of these features are new.
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 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 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.
By using the GRACE-Loop system, we calculate the full $\mathcal{O}(α)$ electroweak radiative corrections to the process $e^+e^- \rightarrow e^+e^- γ$, 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.
The full $\mathcal{O}(α)$ electroweak radiative corrections to $t \bar{t} γ$ and $e^-e^+γ$ 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 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}.
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}(α)$ electroweak radiative corrections to the process $e^+e^- \rightarrow t \bar{t} γ$ 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} γ$ 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 $α$ scheme and the $G_μ$ scheme ($δ_{W}^{G_μ}$). We obtain numerical values for $δ_{W}^{G_μ}$ 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.
Optimizing the cost of evaluating a polynomial is a classic problem in computer science. For polynomials in one variable, Horner's method provides a scheme for producing a computationally efficient form. For multivariate polynomials it is possible to generalize Horner's method, but this leaves freedom in the order of the variables. Traditionally, greedy schemes like most-occurring variable first are used. This simple textbook algorithm has given remarkably efficient results. Finding better algorithms has proved difficult. In trying to improve upon the greedy scheme we have implemented Monte Carlo tree search, a recent search method from the field of artificial intelligence. This results in better Horner schemes and reduces the cost of evaluating polynomials, sometimes by factors up to two.
Over the past few years considerable progress has been made on the resummation of double-logarithmically enhanced threshold (large-x) and high-energy (small-x) higher-order contributions to the splitting functions for parton and fragmentation distributions and to the coefficient functions for inclusive deep-inelastic scattering and semi-inclusive e^+e^- annihilation. We present an overview of the methods which allow, in many cases, to derive the coefficients of the highest three logarithms at all orders in the strong coupling from next-to-next-to-leading order results in massless perturbative QCD. Some representative analytical and numerical results are shown, and the present limitations of these resummations are discussed.
We confirm a conjecture about the construction of basis elements for the multiple zeta values (MZVs) at weight 27 and weight 28. Both show as expected one element that is twofold extended. This is done with some lengthy computer algebra calculations using TFORM to determine explicit bases for the MZVs at these weights.
I give an overview of FORM development based on a few pilot projects, explaining how they have influenced the FORM capabilities. Next I explain what is happnening right now in the field of Open Sourcing and the FORM Forum.
Some of the new features of the symbolic manipulation system FORM are discussed. Then some recent results running its multithreaded version TFORM are shown. Finally the plans for the future are presented.
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 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.
2nd workshop on the implications of HERA for LHC physics. Working groups: Parton Density Functions Multi-jet final states and energy flows Heavy quarks (charm and beauty) Diffraction Cosmic Rays Monte Carlos and Tools
We provide a data mine of proven results for multiple zeta values (MZVs) of the form $ζ(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\{(ε_1^{n_1} >...ε_1 ^{n_k})/ (n_1^{s_1} ... n_k^{s_k}) \}$ with signs $ε_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 present TFORM, the version of the symbolic manipulation system FORM that can make simultaneous use of several processors in a shared memory architecture. The implementation uses Posix threads, also called pthreads, and is therefore easily portable between various operating systems. Most existing FORM programs will be able to take advantage of the increased processing power, without the need for modifications. In some cases some minor additions may be needed. For a computer with two processors a typical improvement factor in the running time is 1.7 when compared to the traditional version of FORM. In the case of computers with 4 processors a typical improvement factor in the execution time is slightly above 3.
We extend the threshold resummation exponents G^N in Mellin-N space to the fourth logarithmic (N^3LL) order collecting the terms alpha_s^2 (alpha_s ln N)^n to all orders in the strong coupling constant as. Comparing the results to our previous three-loop calculations for deep-inelastic scattering (DIS), we derive the universal coefficients B_q and B_g governing the final-state jet functions to order alpha_s^3, extending the previous quark and gluon results by one and two orders. A curious relation is found at second order between these quantities, the splitting functions and the large-angle soft emissions in Drell-Yan type processes. We study the numerical effect of the N^3LL corrections using both the fully exponentiated form and the expansion of the coefficient function in towers of logarithms.
We study the photon-quark-quark and Higgs-gluon-gluon form factors for on-shell massless quarks and gluons in perturbative QCD. Previous third-order results for the quark case are extended by calculating the fermion-loop contributions up to the finite terms in dimensional regularization. For the gluon case the complete set of infrared poles at three loops is derived. Using the exponentiation of the form factor, the latter results are employed to extract a function entering the infrared factorization of general third-order amplitudes. We evaluate the infrared finite absolute ratio of the time-like and space-like gluon form factors up to the fourth order in the strong coupling constant. The result supports previous indications that the perturbative expansion of the Higgs boson production rate at the LHC is under control.
We present the results of an analytic next--to--next--to leading order QCD calculation of the non--singlet anomalous dimension $γ_{\rm NS}^+(N)$ and the coefficient functions $C_{2,L}(N)$ associated to the deeply inelastic structure functions $F_2(x,Q^2)$ and $F_L(x,Q^2)$ for the Mellin moment N=16. Comparisons are made with results in the literature.
We compute the complete third-order contributions to the coefficient functions for the longitudinal structure function F_L, thus completing the next-to-next-to-leading order (NNLO) description of unpolarized electromagnetic deep-inelastic scattering in massless perturbative QCD. Our exact results agree with determinations of low even-integer Mellin moments and of the leading small-x terms in the flavour-singlet sector. In this letter we present compact and accurate parametrizations of the results and illustrate the numerical impact of the NNLO corrections.
We compute the next-to-next-to-leading order (NNLO) contributions to the splitting functions governing the evolution of the unpolarized flavour-singlet parton densities in perturbative QCD. The exact expressions are presented in both Mellin-N and Bjorken-x space. We also provide accurate parametrizations for practical applications. Our results agree with all partial results available in the literature. As in the non-singlet case, the correct leading logarithmic predictions for small momentum fractions x do not provide good estimates of the respective complete splitting functions. We investigate the size of the corrections and the stability of the NNLO evolution under variation of the renormalization scale. The perturbative expansion appears to converge rapidly at x >~ 10^-3. Relatively large third-order corrections are found at smaller values of x.
We compute the fermionic (n_f) contributions to the flavour non-singlet structure functions in unpolarized electromagnetic deep-inelastic scattering at third order of massless perturbative QCD. Complete results are presented for the corresponding nf-parts of the three-loop anomalous dimension and the three-loop coefficient functions for the structure functions F_2 and F_L. Our results agree with all partial and approximate results available in the literature. The present calculation also facilitates a complete determination of the threshold-resummation parameters B_2 and D_2^DIS of which only the sum was known so far, thus completing the information required for the next-to-next-to-leading logarithmic resummation. We find that D_2^DIS vanishes in the MSbar scheme.
The analytic calculation of the three-loop QCD $β$-function and anomalous dimensions within the minimal subtraction scheme in an arbitrary covariant guage is presented. The result for the $β$-function coincides with the previous calculation \cite{tvz}.