Abstract: We provide three Fortran programs which evaluate the QCD analytic (holomorphic) couplings A(v)(Q(2)) for complex or real squared momenta Q(2). These couplings are holomorphic analogs of the powers a(Q(2))(v) of the underlying perturbative QCD (pQCD) coupling a(Q(2)) equivalent to alpha(s)(Q(2))/pi, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2 delta anQCD), and Massive Perturbation Theory (MPT). The index v can be noninteger. The provided programs do basically the same job as the Mathematica package anQCD.m published by us previously (Ayala and Cvetic, 2015), but are now written in Fortran. Program summary Program title: AanQCDext Catalogue identifier: AEYKv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEYICv1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12105 No. of bytes in distributed program, including test data, etc.: 98822 Distribution format: tar.gz Programming language: Fortran. Computer: Any work-station or PC where Fortran 95/200312008 (gfortran) is running. Operating system: Operating system Linux (Ubuntu and Scientific Linux), Windows (in all cases using gfortran). Classification: 11.1, 11.5. Nature of problem: Calculation of values of the running analytic couplings A(v)(Q(2); N-f) for general complex squared momenta Q(2) equivalent to -q(2), in three analytic QCD models, where A(v)(Q(2); N-f) is the analytic (holomorphic) analog of the power (alpha(s)(Q(2); N-f)/pi)(v). Here, A(v)(Q(2); N-f) is a holomorphic function in the Q(2) complex plane, with the exception of the negative semiaxis (-infinity, -M-thr(2)], reflecting the analyticity properties of the spacelike renormalization invariant quantities D(Q(2)) in QCD. In contrast, the perturbative QCD power (alpha(s)(Q(2); N-f)/pi)(v) has singularities even outside the negative semiaxis (Landau ghosts). The three considered models are: Analytic Perturbation theory (APT); Two-delta analytic QCD (2 delta anQCD); Massive Perturbation Theory (MPT). We refer to Ref. [1] for more details and literature. Solution method: The Fortran programs for FAPT and 2 delta anQCD models contain routines and functions needed to perform two-dimensional numerical integrations involving the spectral function, in order to evaluate A(v)(Q(2)) couplings. In MPT model, one-dimensional numerical integration involving A(1)(Q(2)) is sufficient to evaluate any A(v)(Q(2)) coupling. Restrictions: For unphysical choices of the input parameters the results are meaningless. When Q(2) is close to the cut region of the couplings (Q(2) real negative), the calculations can take more time and can have less precision. Running time: For evaluation of a set of about 10 related couplings, the times vary in the range t similar to 10(1)-10(2) s. MPT requires less time, t similar to 1-10(1) s. References: [1] C. Ayala and G. Cvetic, anQCD: a Mathematica package for calculations in general analytic QCD models, Comput. Phys. Commun. 190 (2015) 182.

Abstract: We present the first comprehensive analysis of the unitarity thresholds and anomalous thresholds of scattering amplitudes at two loops and beyond based on the loop- tree duality, and show how non-causal unphysical thresholds are locally cancelled in an efficient way when the forest of all the dual on-shell cuts is considered as one. We also prove that soft and collinear singularities at two loops and beyond are restricted to a compact region of the loop three-momenta, which is a necessary condition for implementing a local cancellation of loop infrared singularities with the ones appearing in real emission; without relying on a subtraction formalism.

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.

Abstract: We present a novel benchmark application of a quantum algorithm to Feynman loop integrals. The two on-shell states of a Feynman propagator are identified with the two states of a qubit and a quantum algorithm is used to unfold the causal singular configurations of multiloop Feynman diagrams. To identify such configurations, we exploit Grover's algorithm for querying multiple solutions over unstructured datasets, which presents a quadratic speed-up over classical algorithms when the number of solutions is much smaller than the number of possible configurations. A suitable modification is introduced to deal with topologies in which the number of causal states to be identified is nearly half of the total number of states. The output of the quantum algorithm in IBM Quantum and QUTE Testbed simulators is used to bootstrap the causal representation in the loop-tree duality of representative multiloop topologies. The algorithm may also find application and interest in graph theory to solve problems involving directed acyclic graphs.

Abstract: A measurement of the Z/gamma* transverse momentum (p(T)(Z)) distribution in proton-proton collisions at sqrt(s) = 7 TeV is presented using Z/gamma* -> e(+)e(-) and Z/gamma* -> mu(+)mu(-) decays collected with the ATLAS detector in data sets with integrated luminosities of 35 pb(-1) and 40 pb(-1), respectively. The normalized differential cross sections are measured separately for electron and muon decay channels as well as for their combination up to pi of 350 GeV for invariant dilepton masses 66 GeV < m(ll) < 116 GeV. The measurement is compared to predictions of perturbative QCD and various event generators. The prediction of resummed QCD combined with fixed order perturbative QCD is found to be in good agreement with the data.