  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.
