The Diversity and Evolution of Dust Attenuation Curves from Redshift z ~ 1 to 9 (2509.01795v1)

Abstract: The UV-optical dust attenuation curve is key to interpreting the intrinsic properties of galaxies and provides insights into the nature of dust grains and their geometry relative to stars. In this work, we constrain the UV-optical slope of the stellar attenuation curve using a spectroscopic-redshift sample of ~3300 galaxies at z~1-9, to characterize the diversity and redshift evolution of stellar attenuation curves and to gain insight into dust production and evolution at high redshifts. The sample is constructed from three JWST/NIRCam grism surveys in GOODS and A2744 fields, with a wealth of JWST/NIRCam and HST photometry. With constraints from spectroscopic redshifts and emission line fluxes, we use the Prospector SED fitting code with a flexible dust model. We find that the attenuation curve slope varies strongly with Av at all redshifts, becoming flatter at higher attenuation. We find no strong correlation between attenuation curve slope and size or axis ratio, and the trends with stellar mass and star-formation rate are largely driven by their correlation with Av. We find strong evidence that at fixed Av, the curve becomes flatter with increasing redshift. On average, the attenuation curves derived here are shallower than those at z~0 and than the SMC curve. The highest redshift galaxies at z=7-9 (124 galaxies, a significantly larger sample than in previous studies) show slopes even flatter than the Calzetti curve, implying reduced UV obscuration and lower IR luminosities than expected from an SMC dust curve, by as large as an order of magnitude. Hydrodynamical simulations that couple dust growth to gas chemical enrichment successfully reproduce the different loci of high- and low-redshift galaxies in the slope-Av diagram, suggesting that dust in high-redshift galaxies is increasingly dominated by large grains produced in supernova ejecta with limited ISM processing at early times.