A first empirical derivation of the average dust attenuation law at 2<z<7

Dust attenuation strongly affects the observed spectral energy distributions of galaxies, introducing significant uncertainties in the derivation of key physical properties such as star formation rates, stellar masses, and metallicities. While attenuation curves have been extensively studied in the local Universe and at intermediate redshift, direct spectroscopic constraints at earlier cosmic epochs have remained limited prior to JWST. We aim to derive the average dust attenuation law of star-forming galaxies over the redshift range 2<z<7. We combine NIRSpec spectroscopy from the JADES survey with deep multi-wavelength photometry from the ASTRODEEP catalogs. Using a mass-selected sample (log(M_\star/M_\odot) > 9) of 120 galaxies with reliable Balmer decrement (Ha/Hb), we construct stacked spectral energy distributions in bins of Balmer optical depth and derive the selective attenuation curve following the empirical methodology introduced by Calzetti et al. (2000). The wavelength coverage is further extended toward the near-infrared using MIRI photometry. The resulting attenuation curve spans the rest-frame range 0.16-1.14mu and is well described by a smooth function. We derive a normalization factor R_V=3.98, finding that the average attenuation law is consistent with the local starburst relation in both slope and normalization. Compared to several determinations at intermediate redshift, however, our curve appears systematically flatter in the ultraviolet. We find no significant evidence for a 2175A UV bump in the average attenuation curve. Our results provide the first empirical determination of the average dust attenuation law for star-forming galaxies at 2<z<7 based on JWST spectroscopy. Despite the diversity of attenuation properties observed in individual systems, the ensemble-average behavior remains consistent with the local starburst relation.