The Sun Also Rises: the High-Temperature Expansion of the Thermal Sunset

We give a prescription for calculating the high-temperature expansion of the thermal sunset integral to arbitrary order. Results up to $\mathcal{O}(T^0)$ are given for both bosonic and fermionic thermal sunsets in dimensional regularisation, and for all odd powers of $T$ up to order $\epsilon^0$. The methods used generalize to non-zero external momentum. We verify the results with sundry analytical and numerical cross-checks. Intermediate steps involve integrals over three Bessel functions.