We present a program for the reduction of large systems of integrals to master integrals. The algorithm was first proposed by Laporta; in this paper, we implement it in MAPLE. We also develop two new features which keep the size of intermediate expressions relatively small throughout the calculation. The program requires modest input information from the user and can be used for generic calculations in perturbation theory.
We discuss the prospects of performing high-order perturbative calculations in systems characterized by a vanishing temperature but finite density. In particular, we show that the determination of generic Feynman integrals containing fermionic chemical potentials can be reduced to the evaluation of three-dimensional phase space integrals over vacuum on-shell amplitudes - a result reminiscent of a previously proposed naive real-time formalism for vacuum diagrams. Applications of these rules are discussed in the context of the thermodynamics of cold and dense QCD, where it is argued that they facilitate an extension of the Equation of State of cold quark matter to higher perturbative orders.
With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% confidence intervals in the Bayesian approach using the new public code MiHO.
We consider two approaches to estimate and characterise the theoretical uncertainties stemming from the missing higher orders in perturbative calculations in Quantum Chromodynamics: the traditional one based on renormalisation and factorisation scale variation, and the Bayesian framework proposed by Cacciari and Houdeau. We estimate uncertainties with these two methods for a comprehensive set of more than thirty different observables computed in perturbative Quantum Chromodynamics, and we discuss their performance in properly estimating the size of the higher order terms that are known. We find that scale variation with the conventional choice of varying scales within a factor of two of a central scale gives uncertainty intervals that tend to be somewhat too small to be interpretable as 68% confidence-level-heuristic ones. We propose a modified version of the Bayesian approach of Cacciari and Houdeau which performs well for non-hadronic observables and, after an appropriate choice of the relevant expansion parameter for the perturbative series, for hadronic ones too.
We describe in detail the implementation of a systematic perturbative approach to observables in the QCD gradient-flow formalism. This includes a collection of all relevant Feynman rules of the five-dimensional field theory and the composite operators considered in this paper. Tools from standard perturbative calculations are used to obtain Greens functions at finite flow time $t$ at higher orders in perturbation theory. The three-loop results for the quark condensate at finite $t$ and the conversion factor for the ringed quark fields to the $overline{mbox{MS}}$ scheme are presented as applications. We also re-evaluate an earlier result for the three-loop gluon condensate, improving on its accuracy.
The recently developed algorithm FIRE performs the reduction of Feynman integrals to master integrals. It is based on a number of strategies, such as applying the Laporta algorithm, the s-bases algorithm, region-bases and integrating explicitly over loop momenta when possible. Currently it is being used in complicated three-loop calculations.