Measuring the total infrared light from galaxy clusters at z=0.5-1.6: connecting stellar populations to dusty star formation

Massive galaxy clusters undergo strong evolution from z similar to 1.6 to z similar to 0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present `total light' stacking techniques spanning 3.4-500 mu m aimed at revealing the total cluster emission, including low-mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our techniques to 232 well-studied log M-200/M-circle dot similar to 13.8 clusters in multiple redshift bins, recovering extended cluster emission at all wavelengths. We measure the averaged IR radial profiles and spectral energy distributions (SEDs), quantifying the total stellar and dust content. The near-IR profiles are well described by anNFWmodel with a high (c 7) concentration. Dust emission is similarly concentrated, albeit suppressed at r less than or similar to 0.3Mpc. The measured SEDs lack warm dust, consistent with the colder SEDs of low-mass galaxies. We derive total stellar masses consistent with the theoretical M-halo -M-* relation and specific star formation rates that evolve strongly with redshift, echoing that of log M*/M-circle dot greater than or similar to 10 cluster galaxies. Separating out the massive population reveals the majority of cluster far-IR emission (similar to 70-80 per cent) is provided by the low-mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low-mass cluster galaxies.