REBELS-IFU: Dust attenuation curves of 12 massive galaxies at $z\simeq7$

Abstract: We present measurements of the dust attenuation curves of 12 massive ($9~<~\log$($M_{\star}/{M}{\odot})$ $<~10$) Lyman-break galaxies at $z=6.5-7.7$ derived from James Webb Space Telescope (JWST) NIRSpec integral field unit (IFU) spectroscopy. The galaxies are drawn from the Atacama Large Millimeter/submillimeter Array (ALMA) Reionization Era Bright Emission Line Survey (REBELS) large program. The dust attenuation curves were obtained by fitting spectral energy distribution (SED) models with a flexible dust law to the full galaxy spectra over observed wavelengths $0.6-5.3$ $\mu$m. These attenuation curves show a range of recovered slopes ($-0.39\leq\delta\leq0.08$) that are on average slightly flatter than seen in local sources of the same stellar masses, with none exhibiting very steep slopes. Three galaxies exhibit evidence for a 2175 \r{A} dust bump ($>4\sigma$) and we find SED fitting excluding the bump can overestimate derived stellar masses by up to $0.4$ dex. Correcting for the dust attenuation with our best-fit attenuation curves we recover a range of intrinsic UV-slopes ($-2.5\leq\beta_0\leq-2.2$). The galaxies show moderate reddening ($A_V~=~0.1-0.6$ mag) and the $A_V$ to stellar mass relation is consistent with local sources. The attenuation law slope is found to correlate with $A_V$, while we see no strong correlation with stellar mass, ${M{\rm UV}}$, or gas-phase metallicity. Overall, our results show little evolution in dust properties in the REBELS sources compared to the local Universe. Comparing our recovered trends to empirical models suggests that the most important factor driving the variation in the attenuation curves in our sample is the dust-star geometry, not the properties of the dust grains themselves.