2021ApJ...910...89B

The underlying distribution of galaxies' dust spectral energy distributions (SEDs) (i.e., their spectra reradiated by dust from rest-frame ∼3 µm to 3 mm) remains relatively unconstrained owing to a dearth of far-IR/(sub)millimeter data for large samples of galaxies. It has been claimed in the literature that a galaxy's dust temperature-observed as the wavelength where the dust SED peaks (λ_peak)-is traced most closely by its specific star formation rate (sSFR) or parameterized "distance" to the SFR-M_* relation (the galaxy "main sequence"). We present 0.''24 resolved 870 µm ALMA dust continuum observations of seven z = 1.4-4.6 dusty star-forming galaxies chosen to have a large range of well-constrained luminosity-weighted dust temperatures. We also draw on similar-resolution dust continuum maps from a sample of ALESS submillimeter galaxies from Hodge et al (2016). We constrain the physical scales over which the dust radiates and compare those measurements to characteristics of the integrated SED. We confirm significant correlations of λ_peak with both L_IR (or SFR) and Σ_IR (∝SFR surface density). We investigate the correlation between log₁₀(λ_peak) and log₁₀(Σ_IR) and find the relation to hold as would be expected from the Stefan-Boltzmann law, or the effective size of an equivalent blackbody. The correlations of λ_peak with sSFR and distance from the SFR-M_* relation are less significant than those for Σ_IR or L_IR; therefore, we conclude that the more fundamental tracer of galaxies' luminosity-weighted integrated dust temperatures are indeed their star formation surface densities in line with local universe results, which relate closely to the underlying geometry of dust in the interstellar medium.