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, $\rm J_{R,0}$, or vertical, $\rm J_{z,0}$, actions, or high azimuthal action, $\rm J_{\phi ,0}$. Generally populations separated by $\rm J_{R,0}$ have a greater variation in bar strength and vertical heating than those separated by $\rm J_{z,0}$. While the bar substantially weakens the initial vertical gradient of $\rm J_{z,0}$, it also drives a strikingly monotonic vertical profile of $\rm J_{R,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.

Gaia DR2 has revealed breathing motions in the Milky Way, with stars on both sides of the Galactic mid-plane moving coherently towards or away from it. The generating mechanism of these breathing motions is thought to be spiral density waves. Here we test this hypothesis. Using a self-consistent, high-resolution simulation with star formation, and which hosts prominent spirals, we first study the signatures of breathing motions excited by spirals. In the model, the breathing motions induced by the spiral structure have an increasing amplitude with distance from the mid-plane, pointing to an internal cause for them. We then show that, at fixed height, the breathing motion amplitude decreases with age. Next, we investigate the signature of the breathing motions in the Gaia DR2 dataset. We demonstrate that, at the location with a consistently large breathing motion, the corresponding amplitude increases monotonically with distance from the mid-plane, in agreement with the model. Furthermore, we show that at the same location, the breathing motion amplitude decreases with age, again similar to what we find in the model. This strengthens the case that the observed breathing motions are driven by spiral density waves.

Context. In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey, the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. Aims. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. Methods. We made use of the information recorded and archived in the observing blocks; during the observing runs; in a number of relevant documents; in the spectra and master catalogue of spectra; in the parameters delivered by the analysis nodes and the working groups; in the final catalogue; and in the science papers. Based on these sources, we critically analyse and discuss the output and products of the Survey, including science highlights. We also determined the average metallicities of the open clusters observed as science targets and of a sample of clusters whose spectra were retrieved from the ESO archive.Results. The Gaia-ESO Survey has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 32 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. Conclusions. The final catalogue will be released through the ESO archive in the first half of 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come.

A scenario for the formation of the bi-modality in the chemical space [α/Fe] versus [Fe/H] of the Milky Way was recently proposed in which α-enhanced stars are produced early and quickly in clumps. Besides accelerating the enrichment of the medium with α-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-[α/Fe] stars are well described by single exponentials, while that of low-[α/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-[α/Fe] stars but only weakly (or not at all) for high-[α/Fe] stars. For low-[α/Fe] stars, the flaring level decreases with age, while for high-[α/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.

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.