Related papers: Difference Equations in Massive Higher Order Calcu…
A system of inhomogeneous second-order difference equations with linear parts given by noncommutative matrix coefficients are considered. Closed form of its solution is derived by means of newly defined delayed matrix sine/cosine using the…
In this paper we prove a weighted sum formula for multiple harmonic sums modulo primes, thereby proving a weighted sum formula for finite multiple zeta values. Our proof utilizes difference equations for the generating series of multiple…
We propose a novel method to compute multi-loop master integrals by constructing and numerically solving a system of ordinary differential equations, with almost trivial boundary conditions. Thus it can be systematically applied to problems…
We investigate conditions for logarithmic complete monotonicity of a quotient of two products of gamma functions, where the argument of each gamma function has different scaling factor. We give necessary and sufficient conditions in terms…
We present the two-mass QCD contributions to the polarized pure singlet operator matrix element at three loop order in $x$-space. These terms are relevant for calculating the polarized structure function $g_1(x,Q^2)$ at $O(\alpha_s^3)$ as…
We perform the analytic calculation of the Mellin moments of the structure functions F_2 and F_L in perturbative QCD up to second order corrections and in leading twist approximation. We calculate the 2-loop contributions to the anomalous…
We present a complete description of the calculation of the spin-dependent next-to-leading order splitting functions. The calculation is performed in the light-cone gauge. We give results for different prescriptions for the Dirac matrix…
We calculate the Laurent series expansion up to ${\cal O}(\epsilon^2)$ for all scalar one-loop one-, two-, three- and four-point integrals that are needed in the calculation of hadronic heavy flavour production. The Laurent series up to…
Following the Mellin and inverse Mellin transform techniques presented in our paper arXiv:1606.02150 (NT), we have established close forms of Laurent series expansions of products of bi- and trigamma functions /psi(z)*/psi(-z) and…
The usefulness of recursive equations to compute scattering matrix elements for arbitrary processes is discussed. Explicit results at tree and one-loop order, obtained by the HELAC/PHEGAS package that is based on the Dyson-Schwinger…
The infrared divergences of QCD scattering amplitudes can be derived from an anomalous dimension \Gamma, which is a matrix in color space and depends on the momenta and masses of the external partons. It has recently been shown that in…
Determining the number of factors in high-dimensional factor modeling is essential but challenging, especially when the data are heavy-tailed. In this paper, we introduce a new estimator based on the spectral properties of Spearman sample…
The iteration procedure of supersymmetric transformations for the two-dimensional Schroedinger operator is implemented by means of the matrix form of factorization in terms of matrix 2x2 supercharges. Two different types of iterations are…
We develop summation by parts (SBP) approach for generating high-order finite-difference schemes on the interval and propose new sets of schemes up to the 12th order. The coefficients of the schemes are governed by values of grid spacing…
We extend several well-known tools from the theory of second-order divergence-form elliptic equations to the case of higher-order equations. These tools are the Caccioppoli inequality, Meyers's reverse Holder inequality for gradients, and…
Due to its highly oscillating solution, the Helmholtz equation is numerically challenging to solve. To obtain a reasonable solution, a mesh size that is much smaller than the reciprocal of the wavenumber is typically required (known as the…
A summary of recent results obtained for higher-order corrections to precision observables in the Standard Model and the Minimal Supersymmetric Standard Model is given. In the Standard Model electroweak two-loop results valid for arbitrary…
We derive the form of the variance-covariance matrix for any affine equivariant matrix-valued statistics when sampling from complex elliptical distributions. We then use this result to derive the variance-covariance matrix of the sample…
Critical exponents are computed for a variety of twist-2 composite operators, which occur in polarized and unpolarized deep inelastic scattering, at leading order in the 1/N_f expansion. The resulting d-dimensional expressions, which depend…
We compute the two-loop heavy quark jet-function in the heavy quark limit. This is one of the key ingredients in next-to-next-to-leading order (NNLO) and next-to-next-to-leading-log order (NNLL) computations of the invariant mass…