2022ApJ...930..142D

Some recent literature has claimed there to be an evolution in galaxies' dust temperatures toward warmer (or colder) spectral energy distributions (SEDs) between low and high redshift. These conclusions are driven by both theoretical models and empirical measurement. Such claims sometimes contradict one another and are prone to biases in samples or SED fitting techniques. What has made direct comparisons difficult is that there is no uniform approach to fitting galaxies' infrared/millimeter SEDs. Here we aim to standardize the measurement of galaxies' dust temperatures with a python-based SED fitting procedure, MCIRSED. We draw on reference data sets observed by Infrared Astronomical Satellite, Herschel, and SCUBA-2 to test for redshift evolution out to z ∼ 2. We anchor our work to the L_IR-λ_peak plane, where there is an empirically observed anticorrelation between IR luminosity and rest-frame peak wavelength (an observational proxy for luminosity-weighted dust temperature) such that \left\langle {\lambda }_{\mathrm{peak}}\right\rangle ={\lambda }_{{\rm{t}}}{({L}_{\mathrm{IR}}/{L}_{{\rm{t}}})}^{\eta } where η = -0.09 ± 0.01, L_t = 10¹²L_☉, and λ_t = 92 ± 2 µm. We find no evidence for redshift evolution of galaxies' temperatures, or λ_peak, at fixed L_IR from 0 < z < 2 with >99.99% confidence. Our finding does not preclude evolution in dust temperatures at fixed stellar mass, which is expected from a nonevolving L_IR-λ_peak relation and a strongly evolving SFR-M_* relation. The breadth of dust temperatures ({\sigma }_{\mathrm{log}{\lambda }_{\mathrm{peak}}}) at a given L_IR is likely driven by variation in galaxies' dust geometries and sizes, and it does not evolve. Testing for L_IR-λ_peak evolution toward higher redshift (z ∼ 5-6) requires better sampling of galaxies' dust SEDs near their peaks (observed ∼200-600 µm) with <=1 mJy sensitivity. This poses a significant challenge to current instrumentation.