Simulations of Dust in Interacting Galaxies I: Dust Attenuation

A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation.