Abstract: We present a detailed study of the dynamics associated with the ghost sector of quenched QCD in the Landau gauge, where the relevant dynamical equations are supplemented with key inputs originating from large-volume lattice simulations. In particular, we solve the coupled system of Schwinger-Dyson equations that governs the evolution of the ghost dressing function and the ghost-gluon vertex, using as input for the gluon propagator lattice data that have been cured from volume and discretization artifacts. In addition, we explore the soft gluon limit of the same system, employing recent lattice data for the three-gluon vertex that enters in one of the diagrams defining the Schwinger-Dyson equation of the ghost-gluon vertex. The results obtained from the numerical treatment of these equations are in excellent agreement with lattice data for the ghost dressing function, once the latter have undergone the appropriate scale-setting and artifact elimination refinements. Moreover, the coincidence observed between the ghost-gluon vertex in general kinematics and in the soft gluon limit reveals an outstanding consistency of physical concepts and computational schemes.

Abstract: We present a novel method for computing the nonperturbative kinetic term of the gluon propagator from an ordinary differential equation, whose derivation hinges on the central hypothesis that the regular part of the three-gluon vertex and the aforementioned kinetic term are related by a partial Slavnov-Taylor identity. The main ingredients entering in the solution are projection of the three-gluon vertex and a particular derivative of the ghost-gluon kernel, whose approximate form is derived from a Schwinger-Dyson equation. Crucially, the requirement of a pole-free answer determines the initial condition, whose value is calculated from an integral containing the same ingredients as the solution itself. This feature fixes uniquely, at least in principle, the form of the kinetic term, once the ingredients have been accurately evaluated. In practice, however, due to substantial uncertainties in the computation of the necessary inputs, certain crucial components need be adjusted by hand, in order to obtain self-consistent results. Furthermore, if the gluon propagator has been independently accessed from the lattice, the solution for the kinetic term facilitates the extraction of the momentum-dependent effective gluon mass. The practical implementation of this method is carried out in detail, and the required approximations and theoretical assumptions are duly highlighted.

Abstract: The numerical evaluation of multi-loop scattering amplitudes in the Feynman representation usually requires to deal with both physical (causal) and unphysical (non-causal) singularities. The loop-tree duality (LTD) offers a powerful framework to easily characterise and distinguish these two types of singularities, and then simplify analytically the underling expressions. In this paper, we work explicitly on the dual representation of multi-loop Feynman integrals generated from three parent topologies, which we refer to as Maximal, Next-to-Maximal and Next-to-Next-to-Maximal loop topologies. In particular, we aim at expressing these dual contributions, independently of the number of loops and internal configurations, in terms of causal propagators only. Thus, providing very compact and causal integrand representations to all orders. In order to do so, we reconstruct their analytic expressions from numerical evaluation over finite fields. This procedure implicitly cancels out all unphysical singularities. We also interpret the result in terms of entangled causal thresholds. In view of the simple structure of the dual expressions, we integrate them numerically up to four loops in integer space-time dimensions, taking advantage of their smooth behaviour at integrand level.