A scenario for the formation of the bi-modality in the chemical space [$\alpha$/Fe] vs [Fe/H] of the Milky Way was recently proposed in which $\alpha$-enhanced stars are produced early and quickly in clumps. Besides accelerating the enrichment of the medium with $\alpha$-elements, these clumps scatter the old stars, converting in-plane to vertical motion, forming a geometric thick disc. In this paper, by means of a detailed analysis of the data from smooth particle hydrodynamical simulations, we investigate the geometric properties (in particular of the chemical thick disc) produced in this scenario. For mono-age populations we show that the surface radial density profiles of high-[$\alpha$/Fe] stars are well described by single exponentials, while that of low-[$\alpha$/Fe] stars require broken exponentials. This break is sharp for young populations and broadens for older ones. The position of the break does not depend significantly on age. The vertical density profiles of mono-age populations are characterized by single exponentials, which flare significantly for low-[$\alpha$/Fe] stars but only weakly (or not at all) for high-[$\alpha$/Fe] stars. For low-[$\alpha$/Fe] stars, the flaring level decreases with age, while for high-[$\alpha$/Fe] stars it weakly increases with age (although with large uncertainties). All these properties are in agreement with observational results recently reported for the Milky Way, making this a plausible scenario for the formation of the Galactic thick disc.

We introduce the study of box/peanut (B/P) bulges in the action space of the initial axisymmetric system. We explore where populations with different actions end up once a bar forms and a B/P bulge develops. We find that the density bimodality due to the B/P bulge (the X-shape) is better traced by populations with low radial, JR,0, or vertical, Jz,0, actions, or high azimuthal action, J{\phi},0. Generally populations separated by JR,0 have a greater variation in bar strength and vertical heating than those separated by Jz,0. While the bar substantially weakens the initial vertical gradient of Jz,0, it also drives a strikingly monotonic vertical profile of JR,0. We then use these results to guide us in assigning metallicity to star particles in a pure N-body model. Because stellar metallicity in unbarred galaxies depends on age as well as radial and vertical positions, the initial actions are particularly well suited for assigning metallicities. We argue that assigning metallicities based on single actions, or on positions, results in metallicity distributions inconsistent with those observed in real galaxies. We therefore use all three actions to assign metallicity to an N-body model by comparing with the actions of a star-forming, unbarred simulation. The resulting metallicity distribution is pinched on the vertical axis, has a realistic vertical gradient and has a stronger X-shape in metal-rich populations, as found in real galaxies.

The direction of the spin vectors of disk galaxies change over time. We present the tilting rate of a sample of galaxies in the NIHAO suite of cosmological hydrodynamical simulations. Galaxies in our sample have been selected to be isolated and to have well determined spins. We compare the tilting rates to the predicted observing limit of Gaia, finding that our entire sample lies above the limit, in agreement with previous work. To test the role of dark matter and of gas we compare the weighted Pearson's correlation coefficients between the tilting rates and various properties. We find no correlation between the dark halo's tilting rate, shape, or misalignment with respect to the disc, and the tilting rate of the stellar disc. Therefore, we argue that, in the presence of gas, the dark halo plays a negligible role in the tilting of the stellar disc. On the other hand, we find a strong correlation between the tilting rate of the stellar disc and the misalignment of the cold gas warp. Adding the stellar mass fraction improves the correlation, while none of the dark matter's properties together with the cold gas misalignment improves the correlation to any significant extent. This implies that the gas cooling onto the disc is the principal driver of disc tilting.