Far-infrared and dust properties of present-day galaxies in the EAGLE simulations

TL;DR

This study uses the EAGLE simulations and radiative transfer modeling to analyze dust properties of present-day galaxies, comparing results with observations and exploring model parameter effects.

Contribution

It introduces a detailed method for deriving dust properties in EAGLE galaxies and compares simulated dust scaling relations with observed data, highlighting model limitations.

Findings

EAGLE galaxies show similar dust-to-stellar mass relations as observed.

Discrepancies in submm color ratios suggest insufficient dust heating in simulations.

Adjusting dust parameters helps constrain dust-related properties for galaxy attenuation modeling.

Abstract

The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a diffuse dust distribution from the gas particles and we re-sample the star-forming gas particles and the youngest star particles into star-forming regions that are assigned dedicated emission templates. We select a set of redshift-zero EAGLE galaxies that matches the K-band luminosity distribution of the galaxies in the Herschel Reference Survey (HRS), a volume-limited sample of about 300 normal galaxies in the Local Universe. We find overall agreement of the EAGLE dust scaling relations with those observed in the HRS, such as the dust-to-stellar mass ratio versus stellar mass and versus NUV-r colour relations. A discrepancy in the f_250/f_350 versus f_350/f_500 submm colour-colour relation implies that part of the simulated dust is insufficiently heated, likely because of limitations in our sub-grid model for star-forming regions. We also investigate the effect of adjusting the metal-to-dust ratio and the covering factor of the photodissociation regions surrounding the star-forming cores. We are able to constrain the important dust-related parameters in our method, informing the calculation of dust attenuation for EAGLE galaxies in the UV and optical domain.