The impact of binary stars on the dust and metal evolution of galaxies

• Published in: Monthly notices of the Royal Astronomical Society 2024-01, Vol.527 (3), p.6292-6311
• Main Authors: Yates, Robert M, Hendriks, David, Vijayan, Aswin P, Izzard, Robert G, Thomas, Peter A, Das, Payel
• Format: Article
• Language: English
• Subjects: Binary stars; Cosmic dust; Dust; Galactic evolution; Galaxies; Massive stars; Metallicity; Novae; Simulation; Star & galaxy formation; Star formation; Stars & galaxies; Stellar evolution
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1093/mnras/stad3419

ABSTRACT We present detailed implementations of (i) binary stellar evolution (BSE; using binary_c) and (ii) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, l-galaxies. This new version of l-galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts (z). We find that binaries have a negligible impact on the stellar masses, gas masses, and star formation rates of galaxies if the total mass ejected by massive stars is unchanged. This is because massive stars determine the strength of supernova (SN) feedback, which in turn regulates galaxy growth. Binary effects, such as common envelope ejection and novae, affect carbon and nitrogen enrichment in galaxies; however, heavier alpha elements are more affected by the choice of SN and wind yields. Unlike many other simulations, the new l-galaxies reproduces observed dust-to-metal (DTM) and dust-to-gas (DTG) ratios at z ∼ 0–4. This is mainly due to shorter dust accretion time-scales in dust-rich environments. However, dust masses are under-predicted at z ≳ 4, highlighting the need for enhanced dust production at early times in simulations, possibly accompanied by increased star formation. On sub-galactic scales, there is very good agreement between l-galaxies and observed dust and metal radial profiles at z = 0. A drop in DTM ratio is also found in diffuse low-metallicity regions, contradicting the assumption of a universal value. We hope that this work serves as a useful template for BSE implementations in other cosmological simulations in future.