Origin of neutron capture elements with the Gaia-ESO survey: the evolution of s- and r-process elements across the Milky Way

We study the abundance patterns and the radial gradients of s-process elements (Y, Zr, Ba, La and Ce), r-process elements (Eu) and mixed-process elements (Mo, Nd and Pr) in the Galactic thin disc by means of a detailed two-infall chemical evolution model for the Milky Way with state-of-the-art nucleosynthesis prescriptions. We consider r-process nucleosynthesis from merging neutron stars (MNS), magneto-rotational supernovae (MR-SNe) and s-process synthesis from low- and intermediate- mass stars (LIMS) and rotating massive stars. The predictions of our model are compared with data from the sixth data release of the Gaia-ESO survey, from which we consider 62 open clusters with age > 0.1 Gyr and 1300 Milky Way disc field stars. We conclude that: i) the [Eu/Fe] vs. [Fe/H] is reproduced by both a prompt and a delayed source, but the quick source completely dominates the Eu production; ii) rotation in massive stars contribute substantially to the s-process elements of the first peak, but MNS and MR-SNe are necessary in order to reproduce the observations; iii) due to the adopted yields, our model overpredicts Pr and underpredicts Nd, while the [Mo/Fe] vs. [Fe/H] is nicely reproduced. For the radial gradients, we conclude that: i) our predicted slope of the [Fe/H] gradient is in agreement with the one observed in open clusters by Gaia-ESO and other high-resolution spectroscopic surveys. ii) The predicted slope of the [Eu/H] radial gradient is steeper than the observed one, independently on how quick the production of Eu is. We discuss the possible causes of this discrepancy in terms of both different Galaxy formation scenarios and stellar radial migration effects. iii) For all the elements belonging to the second s-process peak (Ba, La, Ce) as well as for Pr, we predict a plateau at low Galactocentric distances, which is probably due to the enhanced enrichment from LIMS in the inner regions.