Observational abundances of CEMP stars with patterns in between those produced by the rapid and slow nucleosynthesis processes (CEMP-rs stars) are currently invoked as evidence of synthesis via the intermediate process in the early AGB evolutionary phase of metal-poor low mass stars. Nevertheless, discriminating between r+s- and i-process hypotheses requires high-precision abundances obtained through advanced spectral modelling techniques. Theoretical models of the i-process have become more robust, incorporating refined stellar modelling and nuclear reaction physics, providing ranges of probable elemental abundances and isotopic ratios predictions to be confronted with observational determinations. We performed a new analysis of a high resolution and high S/N UVES spectrum of TYC 6044-714-1. We derived accurate effective temperature and highly precise atmospheric parameters, element abundances, and isotopic ratios using state-of-the-art 1D non-LTE and 3D non-LTE spectral line modelling. Using the latest AGB nucleosynthesis models, we assessed the possibility of the i-process to act aside the s-process. We find that TYC 6044-714-1 was likely born as a normal in-situ halo star about 13 Gyr ago, pre-enriched by the r-process through a standard Galactic chemical-evolution pathway. The s+r model provides the best overall reproduction of the observed heavy-element abundance pattern and Ba isotopic ratios, yielding excellent agreement across all three s-process peaks. While i+s+r models with increasing overshooting efficiency improve the fit for specific elements, they do not consistently reproduce the full abundance pattern. The i+s+r models require extreme and physically implausible conditions, and predict s-process Ba fractions inconsistent with those inferred from isotopic ratios of the 4934 Å resonance line. We conclude that the pure s+r scenario is the most plausible explanation.

Context. The Gaia-ESO survey (GES) is a large public spectroscopic survey that acquired spectra for more than 100 000 stars across all major components of the Milky Way. In addition to atmospheric parameters and stellar abundances that have been derived in previous papers of this series, the GES spectra allow us to detect spectroscopic binaries with one (SB1), two (SB2), or more (SBn >= 3) components. Aims. The present paper discusses the statistics of GES SBn >= 2 after analysing 160 727 GIRAFFE HR10 and HR21 spectra, amounting to 37 565 unique Milky Way field targets. Methods. Cross-correlation functions (CCFs) have been re-computed thanks to a dozen spectral masks probing a range of effective temperatures (3900 K < T-eff < 8000 K), surface gravities (1.0 < log g < 4.7), and metallicities (-2.6 < [Fe/H] < 0.3). By optimising the mask choice for a given spectrum, the newly computed, so-called NACRE (NArrow CRoss-correlation Experiment) CCFs are narrower and allow more stellar components to be unblended than standard masks. The DOE (Detection Of Extrema) extremum-finding code then selects the individual components and provides their radial velocities. Results. From the sample of HR10 and HR21 spectra corresponding to 37 565 objects, the present study leads to the detection of 322 SB2, ten SB3 (three of them being tentative), and two tentative SB4. In particular, compared to our previous study, the NACRE CCFs allowed us to multiply the number of SB2 candidates by approximate to 1.5. The colour-magnitude diagram reveals, as expected, the shifted location of the SB2 main sequence. A comparison between the SB identified in Gaia DR3 and the ones detected in the present work was performed and the complementarity of the two censuses is discussed. An application to the mass-ratio determination is presented, and the mass-ratio distribution of the GES SB2 is discussed. When accounting for the SB2 detection rate, an SB2 frequency of approximate to 1.4 % is derived within the present stellar sample of mainly FGK-type stars. Conclusions. As primary outliers identified within the GES data, SBn spectra produce a wealth of information and useful constraints for the binary population synthesis studies.

Chemical clocks based on [s-process elements/alpha-elements] ratios are widely used to estimate ages of Galactic stellar populations. However, the [s/alpha] vs. age relations are not universal, varying with metallicity, location in the Galactic disc, and specific s-process elements. Current Galactic chemical evolution models struggle to reproduce the observed [s/alpha] increase at young ages. We provide chemical evolution models for the Milky Way disc to identify the conditions required to reproduce the observed [s/H], [s/Fe], and [s/alpha] vs. age relations. We adopt a multi-zone chemical evolution model including state-of-the-art nucleosynthesis prescriptions for neutron-capture elements (AGB stars, rotating massive stars, neutron star mergers, magneto-driven supernovae). We explore variations in gas infall, AGB yield dependencies on progenitor stars, and rotational velocity distributions for massive stars. Results are compared with open cluster data from the Gaia-ESO survey. A three-infall scenario for disc formation captures the rise of [s/alpha] with age in the outer regions but fails in the inner ones, especially for second s-process peak elements. Ba production in the last 3 Gyr of chemical evolution would need to increase by half to match observations. S-process contributions from low-mass AGB stars improve predictions but require increases not supported by nucleosynthesis calculations, even with potential i-process contribution. Variations in the metallicity dependence of AGB yields show inconsistent effects across elements. Distributions of massive star rotational velocities fail to improve results due to balanced effects on elements. We confirm that there is no single relationship [s/alpha] vs. age, but that it varies along the MW disc. Current prescriptions for neutron-capture element yields cannot fully capture the complexity of evolution, particularly in the inner disc.

Identifying the nucleosynthesis processes behind heavy-element enrichment in stellar atmospheres is challenging. It typically relies on comparing observed abundance-to-iron ratios with theoretical predictions relative to the Sun, but this method is prone to uncertainty due to limitations of classical 1D hydrostatic models. One promising but still underexplored approach is to measure the isotopic composition of stellar atmospheres by focusing on elements that have both slow (s)-process and rapid (r)-process contributions. While the study of total elemental abundances offers a simplified view, isotopic ratios are directly linked to the underlying nucleosynthesis processes. Our aim is to provide a reliable method for quantifying the contributions of the s- and r-processes to barium in stellar atmospheres. This is achieved by determining barium isotopic ratios using 1D atmospheric models in combination with a carefully calibrated microturbulence, based on the comparison between subordinate and resonance Ba lines. In this initial study, we use member stars of the globular cluster NGC 6752 to calibrate the microturbulence (v_mic) value for both subordinate and resonance barium lines across different stellar evolutionary stages. This allows us to provide a reliable estimate of v_mic that can be applied to accurately determine barium abundances and isotopic ratios in stars ranging from the main sequence to the upper red giant branch. The v_mic scale adapted for barium subordinate lines is consistent with that derived from 3D model atmospheres, and thus the T_eff-log g dependent relations of the later can be used safely. The v_mic for the resonance line at λ4934 Angstrom – for the determination of the isotopic ratio – is higher, and depends on the equivalent width (EW). We provide calibrated relations between v_mic and EW for measuring isotopic ratios.

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further.

We present the first large-scale 3D kinematic study of ~2000 spectroscopically-confirmed young stars (<20 Myr) in 18 star clusters and OB associations (hereafter groups) from the combination of Gaia astrometry and Gaia-ESO Survey spectroscopy. We measure 3D velocity dispersions for all groups, which range from 0.61 to 7.4 km/s (1D velocity dispersions of 0.35 to 4.3 km/s). We find the majority of groups have anisotropic velocity dispersions, suggesting they are not dynamically relaxed. From the 3D velocity dispersions, measured radii and estimates of total mass we estimate the virial state and find that all systems are super-virial when only the stellar mass is considered, but that some systems are sub-virial when the mass of the molecular cloud is taken into account. We observe an approximately linear correlation between the 3D velocity dispersion and the group mass, which would imply that the virial state of groups scales as the square root of the group mass. However, we do not observe a strong correlation between virial state and group mass. In agreement with their virial state we find that nearly all of the groups studied are in the process of expanding and that the expansion is anisotropic, implying that groups were not spherical prior to expansion. One group, Rho Oph, is found to be contracting and in a sub-virial state (when the mass of the surrounding molecular cloud is considered). This work provides a glimpse of the potential of the combination of Gaia and data from the next generation of spectroscopic surveys.

Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way.

The Gaia -ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted 10 5 stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ( R ~ 47 000) and GIRAFFE medium ( R ~ 20 000) resolution spectra. The analysis of the Gaia -ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a Bayesian Inference method developed specifically for the Gaia -ESO Survey by WG11. The method was also used for the chemical abundance homogenisation within WG11, however, the WG10 chemical abundance data set was too sparsely populated so basic corrections for each node analysis were employed for the homogenisation instead. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6 231 with typical uncertainties for T eff , log g and [Fe/H] of 32 K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6 188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76 675 with typical uncertainties for T eff , log g and [Fe/H] of 64 K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64177 of the stars.

The Li abundance observed in pre-main sequence and main sequence late-type stars is strongly age-dependent, but also shows a complex pattern depending on several parameters, such as rotation, chromospheric activity and metallicity. The best way to calibrate these effects, with the aim of studying Li as an age indicator for FGK stars, is to calibrate coeval groups of stars, such as open clusters (OCs) and associations. We present a considerable target sample of 42 OCs and associations, ranging from 1 Myr to 5 Gyr, observed within the Gaia-ESO survey (GES), and using the latest data provided by GES iDR6 and the most recent release of Gaia that was then available, EDR3. As part of this study, we update and improve the membership analysis for all 20 OCs presented in our previous article. We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES and based on numerous criteria: from radial velocity distributions, to the astrometry and photometry provided by Gaia, to gravity indicators, [Fe/H] metallicity, and Li content. We obtain updated lists of cluster members for the whole target sample, as well as a selection of Li-rich giant contaminants obtained as an additional result of the membership process. Each selection of cluster candidates was thoroughly contrasted with numerous existing membership studies using data from Gaia to ensure the most robust results. These final cluster selections will be used in the third and last paper of this series, which reports the results of a comparative study characterising the observable Li dispersion in each cluster and analysing its dependence on several parameters, allowing us to calibrate a Li-age relation and obtain a series of empirical Li envelopes for key ages in our sample.

Context. The light element (anti-)correlations shown by globular clusters (GCs) are the main spectroscopic signature of multiple stellar populations. These internal abundance variations provide us with fundamental constraints on the formation mechanism of stellar clusters. Aims. Using Gaia-ESO, the largest and most homogeneous survey of open clusters (OCs), we intend to check whether these stellar aggregates display the same patterns. Based on previous studies of many GCs, several young and massive clusters in the Magellanic Clouds, as well as a few OCs, we do not expect to find any anti-correlation, given the low mass of Milky Way OCs. Methods. We used the results based on UVES spectra of stars in Gaia-ESO to derive the distribution of Na and O abundances and see whether they show an unexplained dispersion or whether they are anti-correlated. By selecting only high-probability members with high-precision stellar parameters, we ended up with more than 700 stars in 74 OCs. We examined the O-Na distribution in 28 OCs with at least 4 stars available as well as the Na distribution in 24 OCs, with at least 10 stars available. Results. We find that the distribution of Na abundances is compatible with a single-value population, within the errors. The few apparent exceptions can be explained by differences in the evolutionary phase (main sequence and giant post first dredge-up episode) or by difficulties in analysing low gravity giants. We did not find any indication of an Na-O anti-correlation in any of the clusters for which O has been derived. Conclusions. Based on the very small spread we find, OCs maintain the status of single stellar populations. However, a definitive answer requires studying more elements and larger samples covering different evolutionary phases. This will be possible with the next generation of large surveys.

Context. Recent evidence from spectroscopic surveys points towards the presence of a metal-poor, young stellar population in the low- alpha, chemically thin disk. In this context, the investigation of the spatial distribution and time evolution of precise, unbiased abundances is fundamental to disentangle the scenarios of formation and evolution of the Galaxy. Aims. We study the evolution of abundance gradients in the Milky Way by taking advantage of a large sample of open star clusters, which are among the best tracers for this purpose. In particular, we used data from the last release of the Gaia-ESO survey. Methods. We performed a careful selection of open cluster member stars, excluding those members that may be affected by biases in spectral analysis. We compared the cleaned open cluster sample with detailed chemical evolution models for the Milky Way, using well-tested stellar yields and prescription for radial migration. We tested different scenarios of Galaxy evolution to explain the data, namely, the two-infall and the three-infall frameworks, which suggest the chemical thin disk is formed by one or two subsequent gas accretion episodes, respectively. Results. With the performed selection in cluster member stars, we still find a metallicity decrease between intermediate-age (1 < Age/Gyr < 3) and young (Age < 1 Gyr) open clusters. This decrease cannot be explained in the context of the two-infall scenario, even by accounting for the effect of migration and yield prescriptions. The three-infall framework, with its late gas accretion in the last 3 Gyr, is able to explain the low metallic content in young clusters. However, we have invoked a milder metal dilution for this gas infall episode relative to previous findings. Conclusions. To explain the observed low metallic content in young clusters, we propose that a late gas accretion episode triggering a metal dilution would have taken place, extending the framework of the three-infall model for the first time to the entire Galactic disk.

Context. Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims. This study aims to investigate correlations between the Mean Interdistance (), Mean Closest Interdistance () and Median Weighted Central Interdistance () with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods. We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as age indicators for star clusters. Results. Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions. By conservatively restraining the distance to a maximum of similar to 200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages.

Context. The third data release of Gaia, has provided stellar parameters, metallicity [M/H], [alpha/Fe], individual abundances, broadening parameter from its Radial Velocity Spectrograph (RVS) spectra for about 5.6 million objects thanks to the GSP-Spec module, implemented in the Gaia pipeline. The catalogue also publishes the radial velocity of 33 million sources. In recent years, many spectroscopic surveys with ground-based telescopes have been undertaken, including the public survey Gaia-ESO, designed to be complementary to Gaia, in particular towards faint stars. Aims. We took advantage of the intersections between Gaia RVS and Gaia-ESO to compare their stellar parameters, abundances and radial and rotational velocities. We aimed at verifying the overall agreement between the two datasets, considering the various calibrations and the quality-control flag system suggested for the Gaia GSP-Spec parameters. Methods. For the targets in common between Gaia RVS and Gaia-ESO, we performed several statistical checks on the distributions of their stellar parameters, abundances and velocities of targets in common. For the Gaia surface gravity and metallicity we considered both the uncalibrated and calibrated values. Results. Overall, there is a good agreement between the results of the two surveys. We find an excellent agreement between the Gaia and Gaia-ESO radial velocities given the uncertainties affecting each dataset. Less than 25 out of the approximate to 2100 Gaia-ESO spectroscopic binaries are flagged as non-single stars by Gaia. For the effective temperature and in the bright regime (G <= 11), we found a very good agreement, with an absolute residual difference of about 5 K (+/- 90 K) for the giant stars and of about 17 K (+/- 135 K) for the dwarf stars; in the faint regime (G >= 11), we found a worse agreement, with an absolute residual difference of about 107 K (+/- 145 K) for the giant stars and of about 103 K (+/- 258 K) for the dwarf stars. For the surface gravity, the comparison indicates that the calibrated gravity should be preferred to the uncalibrated one. For the metallicity, we observe in both the uncalibrated and calibrated cases a slight trend whereby Gaia overestimates it at low metallicity; for [M/H] and [alpha/Fe], a marginally better agreement is found using the calibrated Gaia results; finally for the individual abundances (Mg, Si, Ca, Ti, S, Cr, Ni, Ce) our comparison suggests to avoid results with flags indicating low quality (XUncer = 2 or higher). These remarks are in line with the ones formulated by GSP-Spec. We confirm that the Gaia vbroad parameter is loosely correlated with the Gaia-ESO v sin i for slow rotators. Finally, we note that the quality (accuracy, precision) of the GSP-Spec parameters degrades quickly for objects fainter than G approximate to 11 or GRVS approximate to 10. Conclusions. We find that the somewhat imprecise GSP-Spec abundances due to its medium-resolution spectroscopy over a short wavelength window and the faint G regime of the sample under study can be counterbalanced by working with averaged quantities. We extended our comparison to star clusters using averaged abundances, using not only the stars in common, but also the members of clusters in common between the two samples, still finding a very good agreement. Encouraged by this result, we studied some properties of the open-cluster population, using both Gaia-ESO and Gaia clusters: our combined sample traces very well the radial metallicity and [Fe/H] gradients, the age-metallicity relations in different radial regions, and allows us to place the clusters in the thin disc.

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate

In this paper, we describe the wide-field spectroscopic survey telescope (WST) project. WST is a 12-metre wide-field spectroscopic survey telescope with simultaneous operation of a large field-of-view (3 sq. degree), high-multiplex (20,000) multi-object spectrograph (MOS), with both a low and high-resolution modes, and a giant 3 ' 3 arcmin2 integral field spectrograph (IFS). In scientific capability, these specifications place WST far ahead of existing and planned facilities. In only 5 years of operation, the MOS would target 250 million galaxies and 25 million stars at low spectral resolution, plus 2 million stars at high resolution. Without need for preimaged targets, the IFS would deliver 4 billion spectra offering many serendipitous discoveries. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work in synergy with future ground and space-based facilities. We show how it can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; the origin of stars and planets; and time domain and multimessenger astrophysics. WST's uniquely rich dataset may yield unforeseen discoveries in many of these areas. The telescope and instruments are designed as an integrated system and will mostly use existing technology, with the aim to minimise the carbon footprint and environmental impact. We will propose WST as the next European Southern Observatory (ESO) project after completion of the 39-metre ELT.

(abridged) The third data release of Gaia, has provided stellar parameters, metallicity [M/H], [α/Fe], individual abundances, broadening parameter from its RVS spectra for about 5.6 million objects thanks to the GSP-Spec module. The catalogue publishes the radial velocity of 33 million sources. We took advantage of the intersections between Gaia RVS and Gaia-ESO to compare their stellar parameters, abundances and radial and rotational velocities. We aimed at verifying the overall agreement between the two datasets, considering the various calibrations and the quality-control flag system suggested for the Gaia GSP-Spec parameters. For the targets in common between Gaia RVS and Gaia-ESO, we performed several statistical checks on the distributions of their stellar parameters, abundances and velocities of targets in common. For the Gaia surface gravity and metallicity we considered both the uncalibrated and calibrated values. We find an excellent agreement between the Gaia and Gaia-ESO radial velocities given the uncertainties affecting each dataset. Less than 25 of  2100 Gaia-ESO spectroscopic binaries are flagged as non-single stars by Gaia. The temperature scales are in good agreement. The calibrated GSP-Spec gravity should be preferred. We note that the quality (accuracy, precision) of the GSP-Spec parameters degrades quickly for objects fainter than G 11. We find that the somewhat imprecise GSP-Spec abundances due to its medium-resolution spectroscopy over a short wavelength window and the faint G regime of the sample under study can be counterbalanced by working with averaged quantities. We studied some properties of the open-cluster population: our combined sample traces very well the radial [Fe/H] and [Ca/Fe] gradients, the age-metallicity relations in different radial regions, and it places the clusters in the thin disc.

The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G similar to 21 mag. Due to the nature of the Gaia onboard-selection algorithms, these are mostly point-source-like objects. Using data provided by the satellite, we have identified quasar and galaxy candidates via supervised machine learning methods, and estimate their redshifts using the low resolution BP/RP spectra. We further characterise the surface brightness profiles of host galaxies of quasars and of galaxies from pre-defined input lists. Here we give an overview of the processing of extragalactic objects, describe the data products in Gaia DR3, and analyse their properties. Two integrated tables contain the main results for a high completeness, but low purity (50 70%), set of 6.6 million candidate quasars and 4.8 million candidate galaxies. We provide queries that select purer sub-samples of these containing 1.9 million probable quasars and 2.9 million probable galaxies (both similar to 95% purity). We also use high quality BP /RP spectra of 43 thousand high probability quasars over the redshift range 0.05 4.36 to construct a composite quasar spectrum spanning restframe wavelengths from 72 1000 nm.

We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some $800\,000$ astrometric, spectroscopic and eclipsing binaries. More than $150\,000$ Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged)

Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples of the stars of interest. We validate our results by using the Gaia catalogue itself and by comparison with external data. We have produced six homogeneous samples of stars with high quality astrophysical parameters across the HR diagram for the community to exploit. We first focus on three samples that span a large parameter space: young massive disk stars (~3M), FGKM spectral type stars (~3M), and UCDs (~20K). We provide these sources along with additional information (either a flag or complementary parameters) as tables that are made available in the Gaia archive. We furthermore identify 15740 bone fide carbon stars, 5863 solar-analogues, and provide the first homogeneous set of stellar parameters of the Spectro Photometric Standard Stars. We use a subset of the OBA sample to illustrate its usefulness to analyse the Milky Way rotation curve. We then use the properties of the FGKM stars to analyse known exoplanet systems. We also analyse the ages of some unseen UCD-companions to the FGKM stars. We additionally predict the colours of the Sun in various passbands (Gaia, 2MASS, WISE) using the solar-analogue sample.

A renewed interest about the origin of \emph{r}-process elements has been stimulated by the multi-messenger observation of the gravitational event GW170817, with the detection of both gravitational waves and electromagnetic waves corresponding to the merger of two neutron stars. Such phenomenon has been proposed as one of the main sources of the \emph{r}-process. However, the origin of the \emph{r}-process elements at different metallicities is still under debate. We aim at investigating the origin of the \emph{r}-process elements in the Galactic thin disc population. From the sixth internal data release of the \emph{Gaia}-ESO we have collected a large sample of Milky Way thin- and thick-disc stars for which abundances of Eu, O, and Mg are available. The sample consists of members of 62 open clusters, located at a Galactocentric radius from $\sim 5$ kpc to $\sim 20$ kpc in the disc, in the metallicity range $[-0.5, 0.4]$ and covering an age interval from 0.1 to 7 Gy, and about 1300 Milky Way disc field stars in the metallicity range $[-1.5, 0.5]$. We compare the observations with the results of a chemical evolution model, in which we varied the nucleosynthesis sources for the three considered elements. Our main result is that Eu in the thin disc is predominantly produced by sources with short lifetimes, such as magneto-rotationally driven SNe. There is no strong evidence for additional sources at delayed times. Our findings do not imply that there cannot be a contribution from mergers of neutron stars in other environments, as in the halo or in dwarf spheroidal galaxies, but such a contribution is not needed to explain Eu abundances at thin disc metallicities.