Evolution of galaxy attenuation curves driven by evolving dust mass and grain size distributions

Abstract: We investigate the impacts of the evolution of dust mass and grain size distribution within a Milky Way-like (MW-like) galaxy simulation on global attenuation curves, focusing on the optical-UV slope and the 2175 $AA$ bump. We discuss the contributions of star-dust geometry, scattering, and dust properties. Post-processing dust radiative transfer was performed using SKIRT based on the MW-like galaxy simulation. The simulation was carried out with GADGET4-OSAKA, which models the evolution of grain size distributions. For lower inclination angles (closer to face-on), the attenuation curve flattens over time up to t=1 Gyr, then becomes progressively steeper. This steeper slope arises from the interplay between scattering and the dust disk becoming more extended over time (changes in star-dust geometry). At higher inclination, scattering is suppressed, and the attenuation curves slightly steepen over time due to small-grain formation and the bias of observed UV light toward older stars. The bump strengthens on a timescale of ~250 Myr due to the formation of small carbonaceous grains. The bump strength is affected not only by the abundance of small grains but also by star-dust geometry. At higher $A_V$ or higher inclination, the bump weakens. These results may help interpret flatter attenuation curves and weaker bumps in high-redshift galaxies. Variations in star-dust geometry alter the amount of scattered photons escaping the galaxy, driving the anti-correlation between the slope and $A_V$. Scatter in this relation arises from differences in dust optical depth along and perpendicular to the line of sight, reflecting inclination and star-dust geometry. Additional contributions come from variations in grain size distribution and the fraction of obscured young stars.