Inspired by recent lattice measurements, we determine the short-distance ( a ≪ r ≪ 1 / π T ) as well as large-frequency (1 /a ≫ ω ≫ π T ) asymptotics of scalar (trace anomaly) and pseudoscalar (topological charge density) correlators at 2-loop order in hot Yang-Mills theory. The results are expressed in the form of an Operator Product Expansion. We confirm and refine the determination of a number of Wilson coefficients; however some discrepancies with recent literature are detected as well, and employing the correct values might help, on the qualitative level, to understand some of the features observed in the lattice measurements. On the other hand, the Wilson coefficients show slow convergence and it appears uncertain whether this approach can lead to quantitative comparisons with lattice data. Nevertheless, as we outline, our general results might serve as theoretical starting points for a number of perhaps phenomenologically more successful lines of investigation.