The origin of strong $α$-element bimodalities in FIRE simulations of Milky Way-mass galaxies
One of the Milky Way's characteristic features is a strongly bimodal distribution of $α$-process elements, such as Mg, at fixed [Fe/H] in stellar abundances. We examine patterns in [Mg/Fe] versus [Fe/H] in FIRE-2 simulations of Milky Way-mass galaxies. Out of 16 galaxies, 4 are capable of producing a strongly bimodal distribution. In all four galaxies, the high-$α$ population corresponds to an older, radially-compact, thick disk, and the low-$α$ population corresponds to a younger, radially-extended, thin disk, similar to the MW.The transition from high- to low-$α$ took $0.3-1.2\Gyr$ and began $5.5-6.5\Gyr$ ago. [Mg/Fe] decreased at relatively fixed [Fe/H], both in the galaxy overall and at fixed radii: Fe enrichment nearly balanced gas accretion (and therefore dilution), but Mg enrichment was weaker. Importantly, this transition occurred during a period of relatively low gas fraction ($5-15\%$), immediately after a rapid decline in star formation (halving within a few hundred Myr), which caused an increase in Fe-producing white-dwarf supernovae relative to Mg-producing core-collapse supernovae. Only one case coincided with a major merger coalescence. We find similar trends in measuring stars by their current radius and by their birth radius, therefore, radial redistribution did not play a dominant role in the formation of a bimodality or its spatial dependence today. Overall, in FIRE-2, strong $α$-element bimodalities are relatively uncommon ($\sim25\%$), often not associated with a major merger, and arise primarily from a rapid decline in star formation during relatively low gas fraction.