Unveiling the trends between dust attenuation and galaxy properties at z ∼ 2 − 12 with the James Webb Space Telescope

A large variety of dust attenuation and/or extinction curves has been observed in high-redshift galaxies. Some studies investigated their correlations with fundamental galaxy properties, which yielded mixed results. These variations are likely driven by underlying factors such as the intrinsic dust properties, the total dust content, and the spatial distribution of dust relative to stars. Aims. We investigate the correlations between the shape of dust attenuation curves, defined by the UV-optical slope (S) and the UV bump parameter (B), and fundamental galaxy properties. Our goal is to identify the key physical mechanisms that shape the dust attenuation curves through cosmic time in the broader context of galaxy formation and evolution. Methods. We extended the analysis of 173 dusty high-redshift (z ∼ 2 − 11.5) galaxies, whose dust attenuation curves were inferred by fitting James Webb Space Telescope (JWST) data with a modified version of the spectral energy distribution (SED) fitting code BAGPIPES. We investigate the trends between the dust attenuation parameters and different galaxy properties as inferred from the SED fitting: V-band attenuation (A V ), star formation rate (SFR), stellar mass (M *), specific SFR (sSFR = SFR/M∗), mass-weighted stellar age (⟨a⟩*m), ionization parameter (log U), and metallicity (Z). For a subset of sources, we additionally explored the trends with oxygen abundance (12 + log(O/H)), which we derived using the direct T e -based method. Results. We report moderate correlations between S and A V , and B and A V . Galaxies characterized by lower (higher) A V exhibit steeper (flatter) slopes and stronger (weaker) UV bumps. These results agree with radiative transfer (RT) predictions that account for the total dust content and the relative spatial distribution of dust with respect to stars. Additionally, we find that S flattens with decreasing ⟨a⟩*m and increasing sSFR. These two trends can be explained if the strong radiation fields associated with young stars (low ⟨a⟩*m) and/or bursty galaxies (high sSFR) preferentially destroy small dust grains, which would shift the size distribution toward larger grains. Finally, the positive correlation between B and 12 + log(O/H) that emerged from our analysis might be driven by variations in the intrinsic dust properties with the gas metallicity. Conclusions. The shape of the dust attenuation curves primarily correlates with four key galaxy properties: (1) With the redshift, which traces variations in the intrinsic and/or reprocessed dust properties, (2) with A V , which reflects RT effects, (3) with the mass-weighted stellar age or sSFR, which might be driven by the radiation field strength, and (4) with the oxygen abundance, which might be linked to intrinsic dust properties. The overlap between some of these mechanisms makes it difficult to isolate their contributions, however. Further progress requires a combination of observations of a larger galaxy sample, point-like sources, spatially resolved galaxy studies, and theoretical models incorporating dust evolution in cosmological simulations.