NNLO phase space master integrals for two-to-one inclusive cross sections in dimensional regularization

We evaluate all phase space master integrals which are required for the total cross section of generic 2 -> 1 processes at NNLO as a series expansion in the dimensional regulator epsilon. Away from the limit of threshold production, our expansion includes one order higher than what has been available in the literature. At threshold, we provide expressions which are valid to all orders in terms of Gamma functions and hypergeometric functions. These results are a necessary ingredient for the renormalization and mass factorization of singularities in 2 -> 1 inclusive cross sections at NNNLO in QCD.