Spectroscopic surveys of the Milky Way’s stars have revealed spatial, chemical, and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at z ∼ 1.5, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[α/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[α/Fe] stars. This event leads to a bimodality in [α/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since z ∼ 1 is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner disc to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scale heights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H] profiles above the mid-plane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, α-enhancements, and ages of stars; e.g. the most metal-poor stars in the low-[α/Fe] sequence are found to have a scale height comparable to old high-[α/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [α/Fe]–[Fe/H] cannot be captured.

Using the VINTERGATAN cosmological zoom simulation, we explore the contributions of the in situ and accreted material, and the effect of galaxy interactions and mergers in the assembly of a Milky Way-like galaxy. We find that the initial growth phase of galaxy evolution, dominated by repeated major mergers, provides the necessary physical conditions for the assembly of a thick, kinematically hot disc populated by high-[α/Fe] stars, formed both in situ and in accreted satellite galaxies. We find that the diversity of evolutionary tracks followed by the simulated galaxy and its progenitors leads to very little overlap of the in situ and accreted populations for any given chemical composition. At a given age, the spread in [α/Fe] abundance ratio results from the diversity of physical conditions in VINTERGATAN and its satellites, with an enhancement in [α/Fe] found in stars formed during starburst episodes. Later, the cessation of the merger activity promotes the in situ formation of stars in the low-[α/Fe] regime, in a radially extended, thin and overall kinematically colder disc, thus establishing chemically bimodal thin and thick discs, in line with observations. We draw links between notable features in the [Fe/H]-[α/Fe] plane with their physical causes, and propose a comprehensive formation scenario explaining self-consistently, in the cosmological context, the main observed properties of the Milky Way.

Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 10 −5 cm −3 to 10 −3 cm −3 . This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.

Using the cosmological zoom simulation VINTERGATAN, we present a new scenario for the onset of star formation at the metal-poor end of the low-[α/Fe] sequence in a Milky Way-like galaxy. In this scenario, the galaxy is fuelled by two distinct gas flows. One is enriched by outflows from massive galaxies, but not the other. While the former feeds the inner galactic region, the latter fuels an outer gas disc, inclined with respect to the main galactic plane, and with a significantly poorer chemical content. The first passage of the last major merger galaxy triggers tidal compression in the outer disc, which increases the gas density and eventually leads to star formation, at a metallicity 0.75 dex lower than the inner galaxy. This forms the first stars of the low-[α/Fe] sequence. These in situ stars have halo-like kinematics, similar to what is observed in the Milky Way, due to the inclination of the outer disc that eventually aligns with the inner one via gravitational torques. We show that this tilting disc scenario is likely to be common in Milky Way-like galaxies. This process implies that the low-[α/Fe] sequence is populated in situ, simultaneously from two formation channels, in the inner and the outer galaxy, with distinct metallicities. This contrasts with purely sequential scenarios for the assembly of the Milky Way disc and could be tested observationally.

We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-thefly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e. 𝜌 < 10 5 cm −3 ), including photoheating, photoelectric heating, H 2 heating/cooling, cosmic-ray heating, H/He cooling, metal-line cooling, CO cooling, and dust cooling (recombination and gas-grain collisions). We validate the model by comparing 1D equilibrium simulations across six dex in metallicity to existing 1D ISM models in the literature. We apply PRISM to high-resolution (4.5 pc) isolated dwarf galaxy simulations that include state-of-the-art models for star formation and stellar feedback to take an inventory of which cooling and heating processes dominate each different gas phase of a galaxy and to understand the importance of non-equilibrium effects. We show that most of the ISM gas is either close to thermal equilibrium or exhibits a slight cooling instability, while from a chemical perspective, the non-equilibrium electron fraction is often more than three times higher or lower than the equilibrium value, which impacts cooling, heating, and observable emission lines. Electron enhancements are attributed to recombination lags while deficits are shown to be due to rapid cosmic-ray heating. The PRISM model and its coupling to RAMSES-RTZ is applicable to a wide variety of astrophysical scenarios, from cosmological simulations to isolated giant molecular clouds, and is particularly useful for understanding how changes to ISM physics impact observable quantities such as metallic emission lines.