{'fraction_non_alphanumeric': 0.033043542326677965, 'fraction_numerical': 0.023183280532429577, 'mean_word_length': 4.456183840278708, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 16, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Gaia-ESO Survey: Spectroscopic-asteroseismic analysis of K2 stars in Gaia-ESO

The extensive stellar spectroscopic datasets that are available for studies in Galactic Archeaology thanks to, for example, the Gaia-ESO Survey, now benefit from having a significant number of targets that overlap with asteroseismology projects such as Kepler, K2 and CoRoT. Combining the measurements from spectroscopy and asteroseismology allows us to attain greater accuracy with regard to the stellar parameters needed to characterise the stellar populations of the Milky Way. The aim of this Gaia-ESO Survey special project is to produce a catalogue of self-consistent stellar parameters by combining measurements from high-resolution spectroscopy and precision asteroseismology. We carried out an iterative analysis of 90 K2@Gaia-ESO red giants. The spectroscopic values of Teff were used as input in the seismic analysis to obtain log(g) values. The seismic estimates of log(g) were then used to re-determine the spectroscopic values of Teff and [Fe/H]. Only one iteration was required to obtain parameters that are in good agreement for both methods and thus, to obtain the final stellar parameters. A detailed analysis of outliers was carried out to ensure a robust determination of the parameters. The results were then combined with Gaia DR2 data to compare the seismic log(g) with a parallax-based log(g) and to investigate instances of variations in the velocity and possible binaries within the dataset. This analysis produced a high-quality catalogue of stellar parameters for 90 red giant stars observed by both K2 and Gaia-ESO that were determined through iterations between spectroscopy and asteroseismology. We compared the seismic gravities with those based on Gaia parallaxes to find an offset which is similar to other studies that have used asteroseismology. Our catalogue also includes spectroscopic chemical abundances and radial velocities, as well as indicators for possible binary detections.


Introduction
The characterisation of the Milky Way stellar populations in studies of Galactic Archaeology has been greatly advanced with the recently acquired wealth of high-quality spectroscopic data that are now available for hundreds of thousands of stars in our Galaxy. The recent Data Release 2 for the European Space Agency (ESA) space mission, Gaia, has led to the publication of highly accurate astrometry for over a billion stars in the Milky Way (Gaia Collaboration et al. 2018).
In anticipation of this astrometric mecca, a suite of stellar spectroscopic surveys of high-resolution, such as Gaia-ESO (Gilmore et al. 2012), APOGEE (Majewski et al. 2017) and GALAH (De Silva et al. 2015) were created, initiating a new era of large databases containing the spectra and scientific measurements for hundreds of thousands of stars. The next wave of surveys, such as WEAVE (Dalton et al. 2012), 4MOST (de Jong et al. 2019), MOONS (Cirasuolo & MOONS Consortium 2016) and MSE (Szeto et al. 2018) will expand the coverage of these databases into the millions.
However, the robustness of the stellar parameters determined for these large spectroscopic datasets depends, in particular, on the accuracy of the parameters of just a few small samples of reference stars (see discussion in e.g. Jofré et al. 2019). A stellar reference set commonly used for the validation and verification of automated stellar parameterisation pipelines is the Gaia FGK Benchmark Stars Heiter et al. 2015). These are very bright stars but they number only 36 (Jofré et al. 2018, for the latest list) and so, they sparsely sample the FGK stellar parameter space. Reference sets are also drawn from stellar spectral libraries, in particular the ELODIE library (Prugniel & Soubiran 2001), compilations of high quality (but inhomogeneous) literature values such as PASTEL (Soubiran et al. 2016), and stars that are well-known members of open or globular clusters.
These relatively small reference sets are being used to define the parameter scale zero-point upon which the large scale analyses are then based (e.g., Pancino et al. 2017;Kunder et al. 2017). They are therefore crucial for determining the accurate absolute values of the reported stellar parameters (effective temperature T eff , surface gravity log g, metallicity [Fe/H], microturbulence ξ) for these large-scale surveys. They are also essential for effective comparison of the survey datasets with stellar and galactic evolution models as this is key to making a straightforward comparison and combination of data from multiple surveys analysed by Article number, page 1 of 33 arXiv:2007.10189v1 [astro-ph.SR] 20 Jul 2020 A&A proofs: manuscript no. k2ges different pipelines. In the future, parameters for these reference sets derived within Gaia-ESO and other surveys can be compared to results drawn from more sophisticated models of stellar atmospheres that may include advances in non-LTE, 3D and dynamical atmospheres for example.
However, the relatively small number of reference stars stands as a problem given the required stellar parameter space is not well-sampled for comprehensive analysis by the automated pipelines. Greater coverage is needed all the way from cool premain sequence stars to hot OB stars. New stellar reference sets of independently determined stellar parameters are required to keep up with the demands of upcoming large scale surveys. The work presented here seeks to define a new sample of reference stars generated by combining the stellar parameters of T eff and [Fe/H] from spectroscopy with log g from asteroseismology.
Today, the availability of datasets with a wide sky coverage has resulted in many overlapping targets between research fields and presents, thus, an opportunity to make simultaneous use of the collective strengths of multiple types of analyses. This is the case here, among the thousands of targets with high-resolution and high-quality stellar spectra, where there is also asteroseismic information for some thanks to the dedicated monitoring of their oscillations by Kepler (Borucki et al. 2010), CoRoT (Baglin et al. 2006) and, recently, by K2 (Howell et al. 2014). These asteroseismic measurements of the interior of stars combined with the spectroscopic 'exterior' measurements of the same stars and the direct measurement of the parallax of each star, are part of the unfolding revolution in Galactic Archaeology (Miglio et al. 2017).
To take advantage of and further develop this multi-analysis approach, the Gaia-ESO Survey observed several hundreds of stars in the K2 Campaign 3 (C3) field, located towards the South Galactic pole.
The Gaia-ESO Survey is an ESO Large Public Spectroscopic survey designed to target over 100,000 in the key stellar populations of the Milky Way (Gilmore et al. 2012). It has made use of a large range of analysis methods and thus developed key homogenisation procedures to bring all the results together in a robust single star catalogue, the final data release of which currently underway. Gaia-ESO observed stars using the medium-resolution spectrograph GIRAFFE (R∼20,000) and the high-resolution spectrograph UVES (R∼47,000) on the VLT.
K2 is a re-purposing of the Kepler satellite for which the science goals are focussed on the detection of the variations in the light curves of stars in 19 fields along the ecliptic to look for transiting exo-planets. At the time this special project began, 90 giant stars with oscillations detected by K2 were identified within the sample of Gaia-ESO stars observed in C3 with medium-and high-resolution spectroscopy. This sample is referred to hereafter as K2@Gaia-ESO.
In this paper, which is part of the series of the K2 Galactic Caps Project (see also Rendle et al. 2019, hereafter Paper I) and is, in particular, a collaboration with Gaia-ESO, we describe the process of obtaining accurate atmospheric parameters that are consistent with the results from asteroseismology. We use our results to explore age dependencies with abundance ratios and examine potential binary stars in our sample. This sample will further provide a good opportunity for the study of any possible offsets between standard spectroscopic parameters and future sets of parameters derived from more sophisticated models of stellar atmospheres (e.g. non-LTE, 3D, dynamical atmospheres).
In Sect. 2 we describe in more detail the data we use for this work. In Sect. 3 we describe the process we used to determine atmospheric parameters iteratively between spectroscopy Fig. 1: a) HR Diagram of iDR5 FGK stars (grey points), K2 C3 stars analysed in Gaia-ESO iDR5 (black stars) and PARSEC Solar track. b) Kiel diagram of the final set of K2@Gaia-ESO stars for which we obtained spectroscopic parameters as analysed in this paper, with PARSEC Solar track. and asteroseismology. In Sect. 4 we present our chemical abundance results. The comparison to log g based on Gaia parallaxes is presented in Sect. 5. Our findings regarding binary stars are described in Sect. 6 and the final discussion and conclusions are presented in Sect. 7.

The Gaia-ESO Survey sample of K2 stars
Targets within the K2 C3 field were prepared as part of the Gaia-ESO observing programme. This resulted in 496 observations that were available for analysis. The initial set of targets were observed in May and June 2016 and, thus, added into the internal Data Release (iDR) 5 data analysis cycle of Gaia-ESO that had begun at the start of May 2016. The rest of the fields were observed in October 2016 and are part of iDR6.
In total there were 231 targets observed using UVES 580 (blue and red arms) and 265 targets observed using the HR10 and HR21 setups of GIRAFFE. Of these, 182 UVES targets and 133 GIRAFFE targets had been included as part of the iDR5 analysis. them how many would ultimately have asteroseismic detections. At the start of this spectroscopic analysis, 90 of the 496 targets were identified as having K2 asteroseismic detections; 28 of these were observed with GIRAFFE, 62 were observed with UVES. It is possible that more of the full sample will have asteroseismic detections as the K2 analysis advances, but we leave this consideration to future works. Figure. 1a shows the HR diagram of the Gaia-ESO iDR5 FGK stars (S/N>30) and the 224 K2 targets (out of 496) observed by Gaia-ESO that were analysed in iDR5 with this signalto-noise ratio (S/N) cut. A PARSEC stellar track of Solar metallicity and age is also shown (Bressan et al. 2012). Figure. 1b shows the final stellar parameters of the 90 K2@Gaia-ESO sample analysed in this work with the [Fe/H] colourmap. Details of the seismic analysis for these targets can be found in Paper I.

Preliminary spectroscopic parameters
It was important to initiate the iterations between the spectroscopic and seismic analyses from the best starting point possible in T eff . As not all of the K2@Gaia-ESO targets were observed in time to be included in iDR5, it was necessary to compile the rest of the initial stellar parameter set from a variety of other sources including: a photometric T eff ; parameters associated with the synthetic template used in the radial velocity determination in the reduction pipelines; and parameters derived using an available Gaia-ESO node analysis. These values are provided for just the 90 K2@Gaia-ESO stars with asteroseismic detections in Table 1 along with the Gaia-ESO CNAME, EPIC identifier and the instrument with which the spectrum was observed for Gaia-ESO.
For both samples the infrared flux method (IRFM) calibrations of Ramírez & Meléndez (2005) were used to estimate a photometric T eff using the APASS V magnitude and the 2MASS K s magnitude. The [Fe/H] from the iDR5 recommended parameters were used as input to the IRFM calibration equations where possible. Otherwise the [Fe/H] from the GIRAFFE radial velocity determination (Gilmore & et al. 2020) were used for GI-RAFFE, and the [Fe/H] from the Gaia-ESO Nice Node iDR5 analysis (see Smiljanic et al. 2014;Worley et al. 2016, for description) were used for UVES.

Outliers
The sample was then investigated for outliers and discrepancies within this range of parameters, which are discussed below.

Signal-to-Noise
The signal-to-noise (S/N) is a good indicator of the quality of the observed spectra. The distribution of the S/N for the UVES and GIRAFFE samples for the 90 K2 stars are shown in Fig. 2. The majority of the spectra have S/N above 50. Those spectra with the lowest S/N may potentially suffer from insufficient signal causing deviations in the derived stellar parameters. This is considered in Sect. 3.4 in light of the stellar parameters determined by the two analysis teams.

IRFM not applicable
There were three cases for which the magnitudes of the stars did not lie within the range of acceptable values to which the IRFM calibration relations can be applied. They are CNAMEs: 22072768-1440392, 22092416-0610474, and 22105015-1119135. The preliminary parameters were however complete for each and were thus used without an assessment against an IRFM T eff . Otherwise the median difference between the IRFM and preliminary T eff is 20 ± 75 K showing good agreement in general.

No iDR5 [Fe/H]
For the target with CNAME 22034179-0815421 and EPIC ID K2_206298620, no recommended [Fe/H] was reported for iDR5, although T eff and log g were provided. The radial velocity determination provided an associated [Fe/H] of −2.55, however the T eff associated with the radial velocity determination was greater than the iDR5 T eff by ∼300 K and the log g was lower by ∼0.3. The IRFM T eff was in better agreement with the iDR5 T eff . An inspection of the spectrum using iSpec (Blanco-Cuaresma et al. 2014) was carried out, based on a comparison with synthetic spectra generated at both the iDR5 parameters and at the radial velocity determination parameters. In both cases, the [Fe/H] was re-derived, obtaining −2.08 and −1.80, respectively. Based on this and the agreement with the IRFM T eff , the initial parameter estimate for this star was supplemented by [Fe/H]=−2.08 derived from the iDR5 parameters using iSpec.
Based on the preliminary set being complete in T eff and in reasonable agreement with the available photometric T eff , the preliminary T eff from spectroscopy were used to derive the preliminary log g from seismology using the seismic analysis described in Paper I carried out by the Birmingham team (hereafter referred to as BHAM).

Initial spectroscopic parameters
The set of preliminary T eff and [Fe/H] compiled in Table 1 are comprehensive but were unavoidably inconsistent in their source because not all of them were processed previously by the Gaia-ESO Survey. Thus, for the 90 stars found to have asteroseismic Article number, page 3 of 33 A&A proofs: manuscript no. k2ges detections an initial analysis using just iSpec was carried out solely on these stars. For this analysis, the surface gravity was fixed to the preliminary Seismic log g that was based on the preliminary T eff . This now homogeneous set from iSpec was used to explore any further inconsistencies between the iSpec spectroscopic and preliminary seismic results and fill in missing values as in Sect. 2.2.1. The iSpec spectroscopic T eff was used to determine an associated seismic log g. These comprise the initial set of parameters used as the starting point for the following iterative process between the two Gaia-ESO spectroscopic analyses and the seismic analysis. These initial parameters are listed in Table 1.

Iterative determination of parameters
The goal of this process was to iterate between the spectroscopic effective temperature (T eff ,Spec ) and the seismic surface gravity (log g Seis ) to converge on a final set of independently-confirmed stellar parameters.
The seismic log g was determined considering the parameters determined from seismology, namely the frequency of maximum power (ν max ) from the p-mode pulsation analysis and the spectroscopic parameter T eff . We follow the scaling relation of as in Morel & Miglio (2012). Details of how this analysis works can be further found in Morel et al. (2014) for a sample of CoRoT targets and in Pinsonneault et al. (2018) for stars observed with Kepler and APOGEE. As discussed there, the seismic analysis of the p-modes is model-independent and so, the main source of uncertainty stems from the input temperature. Therefore, as the T eff ,Spec determination improves, the log g Seis determination also improves. We point to the discussion of Morel et al. (2014) that a change of 100 K in T eff only affects log g by about 0.005. Therefore, significant improvement in log g by, say, a change of the order of 0.1, requires a change in T eff that is much larger than typical uncertainties of T eff . Nonetheless, by fixing log g to a value that is primarily only affected by ν max allows us to set a spectroscopic T eff scale that is consistent with the seismic log(g) and leads to improvement in other stellar quantities such as chemical composition, masses and ultimately ages (See also discussion in Paper I).
Two Gaia-ESO analysis nodes, EPINARBO and Lumba (see Smiljanic et al. 2014, for further details of these and other nodes), carried out the spectroscopic analysis of the GIRAFFE and UVES spectra of the K2 stars in the iterative procedure. The analysis methods are based on equivalent widths for EPINARBO and spectrum synthesis for Lumba (see Sects. 3.2 and 3.3). Following the requirements of the Gaia-ESO Survey, both methods use MARCS stellar atmosphere models (Gustafsson et al. 2008) and the Gaia-ESO linelist ). These two nodes were selected because they represent two of the most widely used methods for parameter determination in stellar spectroscopy (equivalent widths and syntheses, see Jofré et al. 2019). Figure 3 illustrates the iterative process between spectroscopic and seismic parameter determinations that was followed in this analysis.
In summary, the initial spectroscopic parameters (T eff,iniSpec , [Fe/H] iniSpec ) were used to determine the initial seismic log g iniSeis . This set of T eff,iniSpec , [Fe/H] iniSpec and log g iniSeis were given to both spectroscopic nodes. EPINARBO and Lumba were asked to fix the log g of their analysis to log g iniSeis and could otherwise use T eff,iniSpec and [Fe/H]  These three values: T eff,finSpec , [Fe/H] finSpec and log g finSeis comprise the final stellar parameters of the K2@Gaia-ESO sample. In the final phase EPINARBO and Lumba were then asked to derive chemical abundances for each star based on these stellar parameters.
This was the defined procedure and goal of the K2@Gaia-ESO special project. However the homogenisation and combination of the results into a single final set of parameters per star required a detailed investigation of individual results. This was to ensure that each result was well-understood in an informed manner, which allows for reproducibility. We note that a careful homogenisation of the different node results has been a crucial focus of the Gaia-ESO Survey (Hourihane & et al. 2020;Worley & et al. 2020).
Article number, page 4 of 33 C. C. Worley et al.: The Gaia-ESO Survey: Spectroscopic-asteroseismic analysis of K2 stars in Gaia-ESO

Homogenisation strategy
There were six sets of T eff , six sets of [Fe/H] and four sets of log g produced in the iterative process. Within these sets are the high-resolution (UVES: 62 targets) and medium-resolution (GIRAFFE: 28 targets) subsamples. The high-and mediumresolution results were homogenised separately, as the lower resolution and smaller wavelength range of the GIRAFFE observations required more detailed quality assessment. The individual node results and the analysis undertaken to homogenise them are explained in the subsequent sections.

EPINARBO analysis
The EPINARBO analysis measures the equivalent widths (EW) with the DOOp code (Cantat-Gaudin et al. 2014), which automatically measures equivalent widths with DAOSPEC (Stetson & Pancino 2008). It then derives the stellar parameters and abundances with FAMA (Magrini et al. 2014), which calls spectrum synthesis code MOOG (Sneden et al. 2012). The initial microturbulence parameter (ξ) was computed with the Gaia-ESO relation for stars with different T eff and log g (Smiljanic et al. 2014).
For the K2@Gaia-ESO analysis the surface gravity was fixed to the provided seismic value and EPINARBO iterated to converge on the equilibrium T eff , [Fe/H] and ξ. The stars for which the analysis found a lack of sufficient Fe i and Fe ii lines were flagged EPINARBO and, in particular, it was noted that the blended lines in the medium-resolution GIRAFFE spectra were not ideal for EW methods. Figure 4 shows the progression of the EPINARBO results iteration. Figure 4a- The median and median absolute difference (MAD) of the difference between iteration sets for each parameter and each iteration are also shown. There is little variation from Spec1 to Spec2 for both UVES and GIRAFFE. GIRAFFE shows more scatter in the results while the UVES results seem more stable. The plots of seismic log g are included for completeness showing that any large variation in T eff between iterations does not result in much variation in the log g. This agrees with the findings of Morel et al. (2014). See Sect. 3.4 for more discussion on this.

Lumba analysis
The Lumba analysis (Gavel et al. 2020) performs spectrum synthesis using Spectroscopy Made Easy (SME: Valenti & Piskunov 1996;Piskunov & Valenti 2017). For the K2@Gaia-ESO analysis the surface gravity was fixed to the provided seismic value and Lumba iterated to converge on the equilibrium T eff , [Fe/H] and ξ (Smiljanic et al. 2014). Figure 5 shows the iterative process for the Lumba UVES and GIRAFFE analyses as for Fig. 4. There are some distinct outliers for each instrument set, however the Lumba results are generally very stable between iterations as expected (Morel et al. 2014).

Comparison of node parameters
Within the Gaia-ESO Survey, all of the nodes perform a fully spectroscopic analysis of the Gaia-ESO spectra. While the determination of unconstrained spectroscopic T eff , log g and [Fe/H] was not the goal of this study, it is interesting to compare the unconstrained spectroscopic parameters to those determined by iteration between spectroscopy and astereoseismology. For the purposes of this paper, the unconstrained spectroscopic T eff , log g and [Fe/H] determined by each node are referred to as Spec0. Figure 6 directly compares the EPINARBO and Lumba parameters derived for Spec0 (top row) and for the final iteration, Spec2 (bottom row), for GIRAFFE (blue) and UVES (red) respectively. Based on Fig. 5 and Fig. 4, there was little movement between iterations for each node. Therefore inspecting the final parameters (Spec2) from each node was deemed sufficient. The median and MAD for each parameter between the two nodes are also given. The greater spread of the difference in the GIRAFFE parameters compared to the UVES parameters is clearly seen. The UVES results are in good agreement between the nodes for both Spec0 and Spec2.
As stated in Sect. 3.2, EPINARBO reported that the EW method found inconsistencies in Fe i and Fe ii abundances for the GIRAFFE spectra due to the blending of spectral features at that resolution and, indeed, a total of nine stars are not included at all for GIRAFFE Spec0 as the EW method did not converge on a result. The particular stars with converged results but that were flagged by EPINARBO are highlighted in all panels in Fig. 6 as data points with a yellow circle as those rejected by EPINARBO in Spec0, and data points with central white dots (GIRAFFE), or white triangle (UVES) as those rejected by EPINARBO in Spec2. The use of the EPINARBO flagged stars for assessing the quality of the results is explained in detail in Sect. 3.4.2.
For Spec0 even those that are not rejected have a wide spread, while for Spec2 the stars that are not rejected are tightly distributed and in good agreement with the Lumba results. Certainly this comparison shows that for the GIRAFFE spectra (lower resolution and smaller wavelength range that UVES) fixing the log g using astereoseismology has allowed the spectroscopic T eff to be better constrained and, thus, there is better agreement between the methods and more of the sample is available for abundance analysis. The high-resolution and greater wavelength range of the UVES sample produces good agreement between the nodes results for both Spec0 and Spec2.
The differences in node parameters in Fig. 6 are shown against S/N. There is no obvious indication that as the S/N decreases, the differences between the node parameters increase. Indeed, for the UVES sample, the spread in the differences is fairly consistent and minimal across the S/N range for both Spec0 and Spec2. For the GIRAFFE sample there is a large scatter generally for this smaller sample of 28 stars. As there was no obvious trend with S/N, the differences between the node parameters were used directly in the assessment of the quality of the results.
Considering the Spec2 results in particular, which are the set from which the final stellar parameters will be determined, there is, overall, an offset in [Fe/H] between the nodes for both the UVES and GIRAFFE analyses (∆ [Fe/H] −0.15). There is an offset in T eff (∆ T eff −80 K). This is discussed further in Sect. 3.4.7.
The log g in all cases are those from the seismic analysis, based on the respective spectroscopic T eff . There are some high discrepancies found between the nodes (∆ T eff > 500 K, ∆ [Fe/H] > 0.5) for certain stars, most particularly in the GIRAFFE analysis reflected in the high value of the MAD (MAD T eff > 308 K, MAD [Fe/H]    largest disagreement in T eff equates to a very small shift in seismic log g (∆ T eff ≈ −2200 K corresponds to ∆ log g ≈ −0.09). Figure 7 shows the difference between the two sets of node results for the final seismic log g against the difference in final spectroscopic T eff upon which the final seismic log g values are based. In cyan we show the typical uncertainties for T eff and log g when determined spectroscopically. Differences in spectroscopic T eff greater than 500 K equate to less than 0.1 difference in seismic log g which agrees with the discussion in Sect.3. While asteroseismology pinpoints the log g, it does so from a large potential range in T eff . Therefore complementary methods are needed, such as spectroscopy, to accurately converge on all stellar parameters.
The goal at this point was to combine the EPINARBO and Lumba Spec2 results to produce a final spectroscopic T eff (and [Fe/H]) from which a final seismic log g could be calculated. However, prior to this, it was important to understand the differences between the node analyses, particularly with regard to stars for which there was large disagreement, as we did not want to blindly assume that a mean of the parameters from the two nodes was sufficient as a best final value.
To assess the goodness of the results, ancillary information was compared to the differences in the parameters to ensure that the results that were selected were of the best quality possible.
As stated above, there was no obvious trend with S/N so it was not considered in the following quality assessment.
The IRFM and node final parameters are given in Table 2 for the UVES and GIRAFFE samples. Figure 8 shows the uncertainty distributions for the spectroscopic parameters T eff and [Fe/H] for both EPINARBO and Lumba for the GIRAFFE and UVES analyses. There is a systematic offset between the node distributions and EPINARBO reports a tighter distribution (with strict upper limit for T eff ) than Lumba for both the UVES and GIRAFFE samples.

Uncertainties on parameters
Each node has reported the uncertainties as best suits their analysis pipelines and, thus, they are internally consistent. However the uncertainties are not calculated in the same way between the nodes. Due to this difference, they cannot be used to compare the node results (i.e. as a weight) to determine the best value for individual stars.
The difficulty of making comparisons between node results based on the respective uncertainty distributions, when the uncertainties are not defined in a standard way, has been a challenge throughout the lifetime of the Gaia-ESO Survey and for studies in Galactic Archeaology in general. Further discussions on uncertainty analysis can be found in Jofré et al. (2019).
Article number, page 7 of 33 A&A proofs: manuscript no. k2ges Fig. 7: Node difference in final seismic log g against the difference in final T eff for the high (red) and medium (blue) resolution data. Typical spectroscopic uncertainties are shown in cyan.

Quality assessment
EPINARBO provided a list of stars for which the sample of Fe i and Fe ii lines was inadequate for a robust analysis. These are indicated with a '*' in the CNAME column of Table 2. There are eight of these stars within the GIRAFFE analysis and 1 star within the UVES analysis.
The uncertainty distribution of the Lumba results shows high tails inferring less confidence in the parameter determination for those stars. Objects with T eff uncertainty > 294 K (244 K), which is the median T eff uncertainty + 3σ, and [Fe/H] uncertainty > 0.28 (0.28), which is the median [Fe/H] uncertainty + 3σ, were assumed to indicate lower confidence by Lumba in the parameters for the GIRAFFE (UVES) samples.

Microturbulence (ξ)
The ξ was provided by the nodes for both the GIRAFFE and UVES analyses. The ξ relation derived based on iDR1 was used often as a starting point for iterations or as a derived value depending on the respective node procedures. The difference between node ξ values were in some cases very high, particularly for the GIRAFFE analysis. Figure 9 compares the node ξ values for both the a) UVES and b) GIRAFFE analyses directly with a colourmap of ∆T eff,LM−EP . For UVES there was one star with a large difference in ξ, but otherwise there was no strong evidence of ∆ξ correlating with ∆T eff .
For GIRAFFE the extreme differences in ξ typically were accompanied by other extremes in parameters which are described below. A difference of ξ> 1 km/s was used as a threshold when considered alongside the other indicators of goodness of fit.
This comparison of the ξ values between the two nodes and the two resolutions for the K2@Gaia-ESO sample revealed an error in the microturbulence relation made available for iDR1.
Thanks to this work, the relation was re-derived based on iDR5 values for the iDR6 analysis (see Worley & et al. 2020).

Normalised χ 2
It was necessary to independently assess the goodness of fit of the node solutions for the GIRAFFE analyses due to large differences between the reported node parameters. To this end synthetic spectra for both HR10 and HR21 wavelength ranges were generated for each star using Turbospectrum (Plez 2012) and the MARCS stellar atmosphere models (Gustafsson et al. 2008). Interpolating between the models (Masseron 2006), the spectra were generated at the parameters derived by each node and then a normalised χ 2 was calculated between the synthetic spectrum and the observed spectrum. The observed spectrum was normalised to the synthetic spectrum in each case. In this way the fit was optimised to the solution of each node. The goal was to look for obvious discrepancies in χ 2 to discard extreme outliers.
Two measures were used for assessing the goodness indication of the χ 2 : a) absolute difference between the EPINARBO χ 2 and the Lumba χ 2 (|∆χ 2 |); b) absolute values of each node χ 2 .
The set of parameters that represent the best-fit by χ 2 BF between EPINARBO and Lumba was used as the comparison set to avoid assuming that either EPINARBO or Lumba parameters were the best. Figure 10a compares |∆χ 2 | with χ 2 BF also showing log |∆χ 2 | as a colourmap. This log |∆χ 2 | colourmap is used in Fig. 10b-d to explore the trend with other goodness of fit indicators. When considered in combination with other indicators of goodness of fit, typically all |∆χ 2 | > 1.0 indicated poor agreement between the nodes and so in these cases the node parameter set with the highest χ 2 value was rejected. This threshold is shown as a dashed line in Fig. 10a, and the rejected points have a white dot at the centre. Figure 10b compares the difference in [Fe/H] between the nodes with the difference in T eff between the nodes. Typically, when one is large the other is large, and the log |∆χ 2 | is also large, as expected. Figures 10c and d explore whether there is any trend of log |∆χ 2 | and ∆T eff ,LM−EP with T eff ,LM , and any trend of log |∆χ 2 | and ∆[Fe/H] LM−EP with [Fe/H] LM , in case a particular part of the parameter space was particularly susceptible. The Lumba T eff and [Fe/H] were used as reference. There were no obvious trends in either case.

Photometric T eff
The GIRAFFE sample showed, in particular, large disagreements between the two sets of node results. A comparison to the photometric T eff provided another useful indicator as the photometric and spectroscopic T eff should be in relatively good agreement. Interstellar extinction is not expected to play a significant role in the photometric colours as these stars are towards the Galactic pole. When considered in combination with the other indicators a difference of T eff > 250 K between the photometric and spectroscopic T eff was empirically defined as being too great.

Final GIRAFFE spectroscopic parameters
The final spectroscopic stellar parameters (T eff , [Fe/H], ξ) of the 28 stars in the GIRAFFE sample were defined by first inspecting the node results using the criteria as described above to reject any poor parameter determinations (see Table 2 for the photometric and node parameters, and the outlier assessment, as well as Ta ble 3 for the final spectroscopic parameters. In these cases the remaining node parameter set was used as indicated in column 4 of Table 3. There were 12 cases where only one of the node parameter sets was used. For eight of these cases EPINARBO had reported large discrepancies between the Fe i and Fe ii abundances, or a lack of suitable Fe ii lines, for these stars. The χ 2 test confirms these as poor fits to the data, with the Lumba node providing significantly better fitting results.
The remaining four cases were considered as follows: 1. CNAME 22082566-1532383: The large discrepancies between the EPINARBO T eff and the IRFM T eff (-288 K), and also the Lumba T eff (-360 K), considered with the large discrepancy between the EPINARBO ξ and the Lumba ξ (-1.29 km/s) indicated that there were issues with the EPINARBO result and so the Lumba result was taken as the final parameter set. 2. CNAME 22114679-1126477: Similarly the large discrepancies between the EPINARBO T eff and the IRFM T eff (-455 K), and also the Lumba T eff (-497 K), considered with the large discrepancy between the EPINARBO [Fe/H] and the Lumba [Fe/H] (-0.55) indicated that there were issues with the EPINARBO result and so the Lumba result was taken as the final parameter set. 3. CNAME 22032304-0754111: The large discrepancy between the EPINARBO χ 2 and the Lumba χ 2 (-1.81) indicated an issue with the fit to the data, in this case EPINARBO had the worst fit. This was borne out by inspecting the individual HR10 and HR21 χ 2 values, for which the HR21 χ 2 values were comparable between the nodes, but the EPINARBO χ 2 for HR10 was significantly worse. The difference in T eff reflected this. The Lumba result was taken as the final parameter set. 4. CNAME 22154067-0627110: The discrepancy above the threshold of 1.0 between the EPINARBO χ 2 and the Lumba χ 2 (-1.06) indicated an issue with the fit to the data, also in this case EPINARBO had the worst fit. Inspection of the individual HR10 and HR21 χ 2 values also showed that the HR21 χ 2 values were comparable between the nodes, but the EPINARBO χ 2 for HR10 was significantly worse. The difference in T eff reflected this. The Lumba result was therefore taken as the final parameter set.

Final UVES spectroscopic parameters
The final spectroscopic stellar parameters (T eff , [Fe/H], ξ) of the 62 stars in the UVES sample were defined also by first inspecting the node results using the criteria as described above to reject any poor parameter determinations. Overall there was much better agreement between Lumba and EPINARBO for the UVES sample than for the GIRAFFE sample and so it was not necessary to derive the independent χ 2 as an additional measure of goodness of fit. The difference in T eff and [Fe/H] between nodes is seen in Fig. 6. While there is a mean difference that is comparable to the scatter, we are not in a position to favour one node over the other since we do not have reference values for stars of this sample. For Gaia-ESO, both nodes have been independently cal-Article number, page 9 of 33 A&A proofs: manuscript no. k2ges ibrated with the common set of reference parameters of the Gaia FGK benchmark stars. Therefore, it is not a straightforward task to decide which of the nodes in this particular case, where the pipelines were adapted to iterate with seismology, might be less accurate. In this light, the average of the node results were taken as the final spectroscopic parameters. See Table 2 for the photometric and node parameters and the outlier assessment as well as Table 3 for the final spectroscopic parameters.
There are two remaining cases where the differences are significantly larger and are considered as follows: 1. CNAME 22000793-1203412: The large discrepancies between the Lumba T eff and the IRFM T eff (-263 K) considered with the large discrepancy between the EPINARBO ξ and the Lumba ξ (-0.90 km/s) and the very large uncertainty on the Lumba T eff (512 K) and [Fe/H] (0.50) indicated that there were issues with the Lumba result and so the EPINARBO result was taken as the final parameter set.
2. CNAME 22032202-0829154: EPINARBO flagged this star as the only UVES star that they found to have a large discrepancy between the abundances derived from the Fe i and Fe ii lines. The EPINARBO and Lumba T eff disagree by ∆T eff,LM−EP = −679 K. Both do not agree well with the IRFM T eff with differences of -296 K and 383 K respectively. The initial spectroscopic T eff derived from iSpec was 4643 K, which also does not agree with the IRFM T eff nor EPINARBO. There is, however, no substantial difference between the seismic log g derived for each node (0.03). As EPINARBO flags an issue between the Fe i and Fe ii, the Lumba parameters were taken as the final set. In the subsequent chemical analysis the final Fe i and Fe ii are in good agreement (Fe i-Fe ii=0.1), confirming this choice of parameters.

Fig
. 10: Assessment of rejection criteria based on absolute difference of χ 2 between EPINARBO and Lumba. In each case the colourmap is log |∆χ 2 |: a) |∆χ 2 | in a log scale against the best fit (BF) χ 2 . Rejection limit at 1.0 is shown as a dashed line; The data points with a white centre are those rejected when |∆χ 2 | > 1.0

Final K2@Gaia-ESO stellar parameters
The homogenisation process carried out above resulted in the final catalogue of stellar parameters for the 90 K2@GaiaESO red giant stars. These are shown in Fig. 1 with the [Fe/H] as a colourmap. The final stellar parameters are list in Table 3 and are provided in the final K2@Gaia-ESO catalogue of parameters and abundances which is available online. The columns of the full catalogue are listed in Table B.1.
Article number, page 10 of 33 C. C. Worley et al.: The Gaia-ESO Survey: Spectroscopic-asteroseismic analysis of K2 stars in Gaia-ESO

K2@Gaia-ESO chemical abundances
The final phase of the analysis was the measurement of chemical abundances for those elements typically measured for Gaia-ESO by each node. Based on the final set of stellar parameters (T eff , log g, [Fe/H], and ξ) each node determined their typical chemical abundances for the GIRAFFE and UVES spectra. The nodes provided abundances per star and also per spectral line. The homogenisation was undertaken using the line-by-line abundances in order for the results to be combined in a more informed way such that the effect of each spectral feature is considered separately. Even with the detailed work requiredto produce a refined, high-quality line list , some of the spectral features used here are less well-modelled than others, which can be due to less robust line information or a strong influence coming from surrounding features or less robust normalisation of the continuum about that feature. The larger pool of values for each spectral feature across the two nodes compared to a single global value from each node allows for the outlying spectral features to be identified and discarded if needed, and features that are in good agreement between the nodes to be more highly weighted.
This was carried out only for those elements for which at least one of the nodes measured more that 50% of the sample, that is, more than 14 stars for GIRAFFE and more than 31 stars for UVES. The threshold of 50% was chosen to meet the requirement for a reference set to be as complete as possible in the reported measurements across all the stars in the set. The node results for any element with less than these thresholds were discarded. The remaining elements and number of K2 stars that were measured by each node and the median number of spectral lines each node measured per element are listed in Table 4.
The compilation of spectral lines used for Gaia-ESO is described in Heiter et al. (2019)

Line-by-line cleaning
Lines were accepted or rejected following the rules outlined in Smiljanic et al. (2014) but with appropriate modifications as this project has results from only two nodes, rather than upward of three nodes for the full Gaia-ESO WG11 homogenisation. This meant Rules 1, 2 and 3 of the full Gaia-ESO WG11 homogenisation process were not applicable here. Also we did not apply the weights based on the parameters used in Smiljanic et al. (2014) so Rule 7 was also not applicable. The remaining rules (4, 5, 6, 8) were applied as follows: -Rule 4: When information of the EWs was available, only lines with 5≤EW(mÅ)≤120 were used. Exceptions were sodium (5≤EW (mÅ)≤140) and barium (5≤EW(mÅ)≤250). This could only be applied for the EW node EPINARBO and no cases outside the limits were found for these stars.
-Rule 5: If, for a given species at a given star, abundances from 20 or more different spectral lines were available, we removed the ones that are flagged as blended in the Gaia-ESO line list (U). This was applied as specified.
-Rule 6: If [...] the total number of spectral lines with abundances (for a given species of a given star) is more than 20, a 2σ clipping from the mean value was applied. (The total number of lines is counted across all nodes, therefore if eight nodes provide abundances for 5 lines each, it counts as 40 lines for the clipping). This was adapted to the total lines from the two nodes.
-Rule 8: The median value of multiple lines is adopted as the recommended abundance. This was adapted to be the median of all lines from both nodes.
Article number, page 11 of 33 A&A proofs: manuscript no. k2ges The reported spectral lines from each node were 'cleaned' so to include only those lines flagged as Y or U in the two Gaia-ESO line list flags gf_flag and synflag, except in the case of Rule 4 above. The definition of the Gaia-ESO line list flags are taken from Heiter et al. (2019) and are as follows: 1. gf_flag -Y: Data come from a trusted source (mainly laboratory measurements with excellent accuracy). -U: Quality of data is not decided (lower accuracy laboratory) -N: Data are expected to have low accuracy data and advanced theoretical calculations. 2. synthflag -Y: Line is unblended or only blended with line of same species in both the Sun and Arcturus. -U: Line may be inappropriate in at least one of the stars.
-N: Line is strongly blended with line(s) of different species in both stars.
The reported abundances are thus the median of the remaining line-by-line abundances once the above rules and cleaning were applied. The associated uncertainty was then the MAD converted to a standard deviation by the scale factor of 1.48 (σ=MAD×1.48).
In the cases of Cr ii and Zr ii for UVES, and Ca i, Cr i and Y ii for GIRAFFE, only one node provided results for a sufficient number of stars and these were based on one spectral line only. This resulted in an associated uncertainty of zero by the above prescription. In these cases the uncertainty reported by the node on the spectral line was taken as the uncertainty.
Also, in some cases, the associated uncertainty for a particular star for an element was calculated to be zero while the remaining stars for that element had a non-zero associated uncertainty by the above prescription. In those cases, the median non-zero uncertainty and the MAD converted to sigma on the non-zero uncertainty was calculated from the rest of the sample. A lower threshold was set as the median non-zero uncertainty minus twice the sigma. All associated uncertainties lower than this were set instead to this threshold.

Intra-node corrections
There were seven elements (Mg i, Al i, Si i, Ca i, Ti ii, Ni i, Y ii,) for which combining the results proved difficult due to systematic offsets between the GIRAFFE and UVES datasets of each node. The seven sets of per instrument and per node results are shown in Fig. A.1.
As there were no stars in common between the GIRAFFE and UVES datasets the way to assess an offset was to look at the overlapping region in [Fe/H] [-0.8 to -0.4] between the two datasets per node. Two assumptions were made to assess the magnitude and direction of the corrections: 1) Results from a single node should be continuous in [Fe/H] between the UVES and GIRAFFE results; 2) the Node abundance values around solar metallicity should be similar to the solar abundance.
These allowed the assessment of an offset (difference of the medians across the [Fe/H] range) and the direction of the correction. For instance, for Mg i the EPINARBO GIRAFFE and UVES medians disagree by 0.3 in the overlapping [Fe/H] region. However, we expect solar Mg i at solar [Fe/H], hence the correction was applied to the EPINARBO UVES results. This assessment was carried out resulting in the offsets in Table 5 which were then applied at the per-line level. In most cases there is an offset between species on the order of the scatter on the difference. For Zr, for which the median offset is 0.31±0.11, the largest uncertainties correspond to the stars with the greatest disagreement. Typically for all elements the individual stars with largest difference have a large uncertainty on one of the species measurements. There are no obvious trends with any of the parameters. Figure 12 compares Fe i with Fe ii, then the difference of these against [Fe/H], T eff and log g. These are the abundances derived for each species of Fe using the final set of stellar parameters. They are distinct from, but based on, the stellar parameter [Fe/H] which is the global metallicity determined in the parameter round. There is excellent agreement between Fe i and [Fe/H] for both the UVES (-0.03±0.02) and GIRAFFE (0.03±0.05) samples. The agreement for Fe ii is less good, showing a larger offset with [Fe/H] for UVES (-0.11±0.04), as well as a large offset and much larger scatter for GIRAFFE (0.13±0.14).
In particular, for GIRAFFE, 10 of the 28 stars have a Fe i -Fe ii > 0.2, which is the median combined error on the abundance differences of the species for the GIRAFFE spectra. Two of these have a difference greater than 0.5.
At the parameter determination stage, EPINARBO flagged stars for which the Fe i and Fe ii abundances were inconsistent. These are indicated with white dots for GIRAFFE and as a yellow triangle for the single UVES spectrum in Fig. 12. While some of these are still outliers, some have been reclaimed as now being consistent between Fe i and Fe ii.
Of the two with Fe i and Fe ii abundance differences greater than 0.5, only CNAME 22063424-1530038 was flagged by EPINARBO in the parameter round as having inconsistent Fe i and Fe ii abundances. There was also a difference in Spec2 T eff between the nodes of 504 K. Consequently only the Lumba parameters were used for the final parameters (see Table 3). This star has an outlying rotational velocity which is discussed further in Sect. 6.2. CNAME 22071427-1431390, also with Fe i -Fe ii > 0.5, was not flagged as an outlier during the parameter stage although there is a difference in the Spec2 [Fe/H] values between Lumba and EPINARBO of -0.24. It has, however, been flagged as a potential binary (see Sect. 6). Two other spectra with Fe i -Fe ii > 0.5 were also flagged as potential binaries: CNAMEs 22003290-0808595, 22031541-0753433.
Article number, page 12 of 33 To further understand the possible source of inconsistencies for the remaining seven GIRAFFE stars, for GIRAFFE, the Fe ii value is based on a maximum of two spectral lines and only comes from the EPINARBO analysis, compared to at least nine lines for the Fe i value from both EPINARBO and Lumba (see Table 4). Additionally the Fe ii lines are also typically weaker and thus more difficult to measure in medium-compared to highresolution spectra, resulting in less certain values as indicated by the larger uncertainty bars.
Article number, page 13 of 33 A&A proofs: manuscript no. k2ges For UVES, 3 of the 62 stars have a Fe i -Fe ii > 0.16, which is the median combined error on the abundance differences of the species for the UVES spectra, although none of these have a difference greater than 0.3. None of these three (CNAMEs 22065112-1504580, 22195215-1234594, 22021848-1139147) were identified as an outlier in the parameter homogenisation. CNAME 22065112-1504580 has the largest difference between Fe i and Fe ii of -0.28 and it is the only one of these three to be flagged as a potential binary (see Sect. 6).
We note that the ionisation equilibrium (consistent abundance values between neutral and ionised species) may not have been achieved in these cases because the log g value was constrained by non-spectroscopic data, which may reflect shortcomings of the atmospheric models used in the spectroscopic analysis, such as the LTE and 1D prescriptions.

K2@Gaia-ESO final chemical abundances
The final chemical abundances for the K2@Gaia-ESO stars with the above corrections applied are shown in full in Fig. C.1. The full catalogue table is available online but the columns are provided in Table B.1. The UVES (high-resolution) and GIRAFFE (medium-resolution) measurements are shown separately. In general the two resolutions track each other for each element distribution, with the improvements due to the corrections applied above. Overall, there appears to be greater scatter in the GIRAFFE datasets, which is not unexpected, particularly for those elements where the measured spectral lines are weak or more blended at lower resolutions (see figures in e.g. Jofré et al. 2019).
The final stage of this work is focussed on the derivation of ages and masses within the greater K2 Galactic Caps Project (see Paper I). Specifically for this paper, the ages of the K2@Gaia-ESO sample were derived, and are shown as a colour map on the chemical distribution of key elements in Fig. 13.
This set of elements is arranged in five key nucleosynthetic channels: Light odd-Z, α, Fe-Peak, s-process and r-process. Each of these channels reflect different nucleosynthetic origins in stellar evolution. Elements created in the same channel show similar behaviour in chemical distribution morphology. These are interpreted in consideration of the recent release of the GALAH DR2 abundances presented in (Buder et al. 2019) which used a pipeline which was developed based on the experience with the Lumba pipelines. We note that this is a qualitative comparison only. On the one hand, there are no stars in common between this sample and GALAH DR2. In fact, Gaia-ESO and GALAH only have a small overlap in observed targets (the Gaia benchmark stars, some targets in M67 and the CoRoT fields). On the other hand, our methodology differs from the standard Gaia-ESO parameters so direct comparison between surveys is not the aim here.
Considering the two Light odd-Z elements, Na i shows the expected trend of depleted Na at low metallicity increasing to ∼0.15 at super-solar metallicity. Its distribution is similar to Ni i and Cu i as expected from previous studies (Buder et al. 2019). On the other hand, Al i shows enhanced Al at low metallicity, although no non-LTE corrections have been applied contrary to Buder et al. (2019). Al behaves similarly to Mg and the other α elements, also noted in Buder et al. (2019).
The α elements presented here (Mg i, Ca i, Ti i 1 ) all bear the same morphology of showing enhancement at low metallicity decreasing to solar at solar metallicity, which is the typical chemical distribution of α elements (Buder et al. 2019).
The Fe-peak elements (V i, Mn i, Co i, Ni i, Cu i) show at least four types of morphology. V i is unique in that it stays generally constant with a large scatter over the presented metallicity range, which is consistent with other studies. Mn i also is consistent with other studies showing a clear trend of increasing Mn i with metallicity. Co i, Ni i and Cu i can be grouped with similar if not the same morphology of slightly increasing at the metal-poor end, then slightly decreasing then slightly increasing again at solar. This is also not dissimilar to Na i as noted above. The overall increasing trend with metallicity follows that presented in other studies (Buder et al. 2019) although the behaviour around solar begs further investigation.
The s-process is separated into the light (Y ii and Zr ii) and heavy (Ba ii, La ii, Nd ii) s-process peaks. Y ii shows a large scatter with there being enhancements and depletions for both old and young stars. This was similarly found in Buder et al. (2019). There are fewer measurements for Zr ii, particularly at the metalpoor end, but a similarly large spread is indicated.
For the heavy s-process, there are only few measurements at the metal-poor end but these mostly indicate enhanced levels for metal-poor stars for Ba ii, La ii and Na ii. Despite showing a larger scatter for Ba ii, all are fairly constant around solar.
Finally, the single r-process element, Eu ii, shows enhancements at low metallicity, decreasing to a constant solar values at solar and super-solar metallicities, not unlike the α distributions. Although the sample is small compared to other samples where spectroscopy and seismology is available (e.g. Pinsonneault et al. 2018), even for this sample of 90 stars, the morphology of the different nucleosynthetic channels can be seen.
The key outliers in this sample are the two most metal-poor stars which appear to be quite young (∼2 Gyr). This is discussed further in Sect. 6.

K2@Gaia-ESO@Gaia: parallaxes
The iterative determination of the seismic log g and spectroscopic T eff and [Fe/H] was completed prior to the release of Gaia The log g based on the Gaia parallax was calculated using the following additional inputs: the final K2@Gaia-ESO T eff and [Fe/H]; the reddening, E(B-V), assigned to each target as part of the Gaia-ESO dataset was taken from Schlegel et al. (1998); and the 2MASS K band photometry (Skrutskie et al. 2006).
Of the 90 stars, only 58 are compared here due to the Gaia parallax either not being available or not good enough (parallax uncertainty was greater than 15%).
An iterative procedure using isochrones to determine the optimal stellar mass and surface gravity based on the Gaia parallaxes was carried out (log g iso ). Further details on the calculation are provided in Appendix D. For comparison, the surface gravity based on the Gaia parallax assuming only solar mass was also calculated (log g M ). Figure 14a compares log g M with the final K2@Gaia-ESO log g and Fig. 14b compares log g iso with the final K2@Gaia-ESO log g. The mean offset for each are ∆log g M −K2 -0.20±0.14 and ∆log g iso−K2 -0.08±0.08 respectively. Overall log g iso are in better agreement with the K2@Gaia-ESO log g than log g M . The typical reported log g uncertainties of K2 are 0.02 and for the Gaia analysis are 0.13. Spectroscopic log g typical uncertainties are 0.10 to 0.25 for the Gaia-ESO high-resolution UVES spectra (Smiljanic et al. 2014) and 0.15 to 0.40 for the Gaia-ESO medium-resolution GIRAFFE spectra (Worley & et al. 2020).
It is known that the parallaxes provided by Gaia DR2 are affected by systematics (Lindegren et al. 2018;Zinn et al. 2019). A comparison of Gaia parallaxes and K2 seismic parallaxes for stars in the C3 field in Khan et al. (2019) found an offset ranging between -45 to -55 µas. The authors found that the correction in the parallax between Gaia and seismic values seems to depend, to varying degrees, on the position on the sky, the magnitude, and, potentially, the colour. Also, specifically for the K2 C3 and C6 fields, the correction seems to be smaller than in the Kepler field. Following the recommendation in Khan et al. (2019), no ad hoc correction was attempted here but the possible systematics are noted as follows.
Using the equations in Appendix D, the possible range in offsets in parallax corresponds to a range in offsets in log g Gaia-K2 of -0.05 to -0.07. This is in the same direction and of similar magnitude to that measured here, showing the offset between parallaxes from Gaia DR2 and asteroseismology found in this study is in agreement with the previous studies (Zinn et al. 2019;Khan et al. 2019).
The main outlier is CNAME 22235003-1422417 with ∆log g Gaia-K2 = −0.30 ± 0.16. This star was observed with UVES and had not been flagged as particularly peculiar otherwise. The seismic log g values for each node were in good agreement which is not unexpected as shown above. The difference between the node T eff values is LM-EP=-130 K and [Fe/H] values is LM-EP=-0.39. These are within the acceptable limits imposed in this study for classifying outliers. Three other stars have |∆log g Gaia-K2 | > 0.20.

K2@Gaia-ESO@Gaia: binary stars
Both the spectra of Gaia-ESO and the astrometric measurements of Gaia offer important information about a star's movement. These can be combined to study potential binary stars in the sample. This information is given in Table 6.
We compared the radial velocities (RV) of Gaia-ESO with those from Gaia Katz et al. 2019), from the Gaia Radial Velocity Spectrometer (RVS: Cropper et al. 2018). We generally find a very good agreement, with a zero point offset of −0.44 km/s for UVES and of −0.95 km/s for GI-RAFFE, based on which we infer that Gaia-ESO provides systematically lower velocities. The offset is determined from the mean of the RV differences between Gaia and Gaia-ESO. We note the relatively large offset for the GIRAFFE case is based on very few measurements. Indeed, the faint magnitudes of the GIRAFFE sample make the overlap with Gaia RVS very small.
Individual comparisons are shown in Fig. 15, in which we plot the difference of the reported velocities from Gaia-ESO and Gaia as a function of the Gaia uncertainty for all stars that are contained in both datasets. While the majority of the stars have uncertainties in Gaia that are below 1 km/s, there are a few stars that are more uncertain. There is, however, no relation between Gaia RV uncertainty and apparent magnitude. The values used for the figure can be found in Table 6.
In addition, by taking a conservative cut of (1) an uncertainty reported by Gaia larger than 1.8 km/s or (2) a difference in radial velocities between Gaia and Gaia-ESO larger than 1.8 km/s, we can select potential binaries. In Fig. 15, lines marking an uncertainty or difference of 1.8 km/s are shown vertically and horizontally. It is clear that some stars do not agree that well, falling outside the box marked by the lines within 1.8 km/s. These measurements might be attributed to binary stars. We note that the reported Gaia RV corresponds to the median of each of various RV measurements at different epochs, and its uncertainty may be attributed as the variability in RV for each target. The number of transits varies a lot between each Gaia target, and is indicated in Table 6. A large uncertainty might be related to a large variability, which can be an indication of binary systems (see e.g. recent paper on comparisons of RV of Gaia and RAVE by Birko et al. 2019).
We selected our binaries as those with differences in RV of at least 1.8 km/s to account for the intrinsic variation in RV due to jitter, time-variable winds, spot visibility and so forth, which can have a notable effect in evolved giant stars (see e.g. Carney et al. 2003, for a reference about jitter on metal-poor giants), and for differences of zero point between Gaia-ESO and Gaia, which in our case is for both datasets below 1 km/s. This gives a total of 14 potential binaries in our sample, corresponding to 15%. This value agrees with the estimate of 17% for metal-poor giants found by Carney et al. (2003) but the comparison should be taken as a reference only since both samples are selected entirely different. Indeed, the fact that we have been able to provide parameters for the stars implies that no peculiarity in the spectra has been found. It is important to remark that this comparison is also not fully representative as Gaia-DR2 delivered RVs only for stars that are not double-lined spectroscopic binaries. Gaia-ESO, on the other hand, determines RVs for every star, regardless of their binary status. Indeed, a special working group within Gaia-ESO (WG14) analyses double-line spectroscopic binaries and peculiar stars, flagging every target prior to the spectral analysis (Merle et al. 2017(Merle et al. , 2020. None of the stars in this sample have been flagged as a double-line spectroscopic binary or as a peculiar star by WG14. Fig. 16: Age-metallicity relation of the stars in our sample. Grey indicates the stars that are binaries (B) and coloured symbols indicate stars that are constant (C) according to our classification (see text). Error bars in age represent the mean between upper and lower age uncertainty estimate (see Paper I).
It is worth commenting here that we aimed to identify and investigate the potential binaries that this sample might contain. Yet, variable RV is not the only indication that the star may be a binary. If stars have relatively small RV variation and are binaries at the same time, the companion is likely to be less massive, and, som it exhibits larger RV excursions. As a result, the lines from the companion fall in the wings of the lines of the primary star so that the spectrum appears to be that of a fast rotator, even though in reality it is a set of blended lines. This is seen also in RAVE where the spectrum of a cool dwarf appeared to be a fast rotator, but was ultimately found to be a binary. A similar effect might exist for giants.
In Fig. 16, we plot the ages and metallicities of the stars, coloured according to the respective instrument and indicating in grey the stars that are binaries according to our classification. In general, the expected age-metallicity relation is found, in which metal-poor stars are old and metal-rich stars can have a wide range of ages. We note that the oldest stars in our sample show ages older than the Universe. As discussed in Paper I, this reflects that the prior used for the age determination was set to a maximum of 20 Gyr. When uncertainties of the results are taken into account, the oldest stars in our sample have ages within the age of the Universe. Here we plot the absolute values with the intention to study their distribution but full details can be found in Paper I.
We see that two metal-poor stars show notably young ages. For one of them we have further RV information from Gaia and we note that this star is classified as a possible binary. This star, along with other binary candidates that stand out, are discussed further below.

EPIC 206101493: potentially an evolved blue straggler
As in several other recent works, some metal-poor (and alphaenhanced) stars show unexpectedly young ages. The case of EPIC 206101493 is one example found in this K2@Gaia-ESO Article number, page 17 of 33 A&A proofs: manuscript no. k2ges sample. Its CNAME is 22000793-1203412 and has been observed with UVES.
Such apparently young stars have been reported by Chiappini et al. (2015), and were called 'young alpha-rich stars' (YAR). The observation that their high masses might imply they are young (Martig et al. 2015) and that their chemical composition is more consistent with that of thick disk stars (Matsuno et al. 2018) challenges our belief that the thick disk formed a long time ago and very rapidly. It is thus very important to first test if they are (or were) binaries, as this determines whether we can or can not apply stellar models of isolated objects to constrain the age of YAR stars. This is discussed in Jofré et al. (2016) with the help of RV monitoring of the YAR stars listed by Martig et al. (2015), and then studied in Izzard et al. (2018) with the help of stellar synthesis models which include the interaction of binary stars. It is thus interesting to see that in our K2@Gaia-ESO sample, EPIC 206101493, as another identified YAR star, is flagged as a binary due to its RV inconsistencies between Gaia-ESO and Gaia.
EPIC 206101493 has a UVES spectrum but since it is one of the faintest stars in our sample, its spectrum has a very low S/N (∼ 29), this does not allow us to perform a very detailed chemical analysis to investigate its properties further, like Matsuno et al. (2018) did for a sample of these stars. In any case, at a first glance, its parameters are well-determined and its abundances are consistent with a thick disk star. It presents no evident enhancements of s-process or Li. This means that EPIC 206101493 is as typical a thick disk star as the other YAR stars studied in Matsuno et al. (2018), and its high mass agrees with what has been postulated by Jofré et al. (2016), namely that it is likely a product of mass transfer.

EPIC 206322094 and EPIC 205977363: induced rotation
In Fig. 17 we show the mass and the projected rotational velocity (v sin i) for our stars, following the same symbols as in previous figures. The projected rotational velocity is determined by Gaia-ESO alongside the radial velocity through the cross-correlation with templates. This is similar to the procedure reported by Carney et al. (2008) for their study of line broadening and rotation of metal-poor giants.
We can see that the giants in our sample have typical v sin i of 2-4 km/s, regardless of their binary status or stellar mass. There are however two stars that stand out, with v sin i > 6 km/s (EPIC 206322094 v sin i = 7.85 km/s; EPIC 205977363 v sin i = 6.25 km/s). We note there are other stars in our sample with even higher v sin i but they are excluded from this discussion since they do not have RVs from Gaia DR2. While this value is higher than the rest of the stars in the sample, an abundance analysis is still possible and could be done safely.
Giants with these values of rotation have been classified as 'anomalous rotators' by Tayar et al. (2015) in the APOKASC sample, and are attributed to be the cause of a recent interaction.
As extensively discussed by Carney et al. (2008), the higher broadening of these stars is not likely to be due to macroturbulence. There is a trend of macroturbulence with stellar parameters which is not followed for the outlier cases. This suggests that something else might be causing the broadening. Carney et al. (2003) offer other interesting explanations for outstanding line broadening which are worthy of note.
EPIC 206322094 (CNAME 22031541-0753433) is marked as a binary, and its line broadening might be reflecting an induced rotation due to tidal interactions. Gaia reports an uncertainty in RV of 2.5 km/s. If the star is tidally interacting with a companion, then its orbit must be circular, which can be tested with individual RV measurements. Unfortunately, this information is still not available from Gaia, but a follow-up campaign on this star is ongoing and will be reported on in a separate paper.
The other star with outlying rotational velocity is EPIC 205977363 (CNAME 22063424-1530038), however, this star does not show indication of having a binary companion from its radial velocity behavior. The star could be the result of a merger (Tayar et al. 2015). Another interesting possibility for its increase in v rot is the capture of a giant planet of few Jupiter masses due to which the star grew in size and became a red giant. It has been shown in the literature that this process can spin up giants a few km/s depending on the planet mass and the stellar radius (Privitera et al. 2016). Siess & Livio (1999) have studied this scenario extensively, discussing possible effects that can be observed. One effect of planet capture is mass loss, which can be detected by the presence of an emission feature or asymmetry in the Hα line. Alternatively, recent discussion in the literature connects Li-rich giants to planet engulfment, where the Li abundance now unexpectedly present in the star was initially material from the planet (see e.g. discussions in Casey et al. 2016, using Gaia ESO data, and references therein). Although it is not entirely clear that engulfment might explain all Li-rich stars (see e.g. Aguilera-Gómez et al. 2016). Since this star is observed with the GIRAFFE setup, we do not have information about the Li abundance or the Hα profile unless a spectrum covering this wavelength domain is undertaken in the future.

EPIC 206070270 and its unseen companion
In our sample we have one star (EPIC 206070270 -CNAME 22182719-1252466) with a notably large uncertainty in RV as determined by Gaia (above 8 km/s). The radial velocities measured by Gaia-ESO and Gaia differ by 25 km/s, which is consistent with the large RV uncertainty. While this is probably due to it being a binary system, without single RV measurements it Article number, page 18 of 33 C. C. Worley et al.: The Gaia-ESO Survey: Spectroscopic-asteroseismic analysis of K2 stars in Gaia-ESO is difficult to know for sure. However it is possible that this star is similar to the binaries identified in APOGEE by Price-Whelan et al. (2018). In their catalogue, they use the RV variation (and single epoch RV visit) of APOGEE giants to find the properties of possible binaries. The spectrum is further used to determine the stellar parameters and hence the mass of the giant. The stellar masses they determined cover a range that encompasses the mass of EPIC 206070270 as determined by seismology. Furthermore, they show examples of giants with high variations in RV, in some cases reaching the 25 km/s we find here. This value is also a typical value for giants found by the study of mutiplicity of giants in APOGEE by Badenes et al. (2018). If we consider that this difference is the maximum amplitude of the RV variation, and that the giant has a mass of 1.2 M (see Paper I), we can use the information in Table 5 of Price-Whelan et al. (2018) to estimate that the companion of EPIC 206070270 must be a star with approximately half a solar mass.
The fact that we can constrain the mass and age of the primary with asteroseismology makes this system a potentially interesting target for follow-up studies of RV to constrain the orbits and masses of the secondary. This campaign is taking place already with the HERMES telescope in La Palma and the results will be presented elsewhere once enough RV measurements and a sufficient time span has been achieved to robustly determine the parameters of the binary system.

Discussion and conclusion
This work presents a detailed study to determine the stellar parameters of 90 red giant stars by combining results from spectroscopy and asteroseismology. The spectroscopy comes from data taken as part of the Gaia-ESO Survey in a project designed to expand the available set of calibration samples. The asteroseismology comes from data taken by K2 in the C3 field. The parameters could therefore be determined with higher accuracy due to the iterative procedure of deriving T eff from the spectra to be used as input for deriving log g from oscillation frequencies which is then used to re-derive T eff . A substantial change in T eff results in just a small change in the seismic log g and so both methods are required to ensure the required accuracy in parameters.
Understanding the origin of the Galactic disks is based on a fundamental question and the answer is sought in the exploitation of stellar surveys. However, there are underlying challenges due to the inaccuracy of the stellar parameters and abundances derived for metal-poor and ancient stars. Having as many reference stars as possible in this regime is key to finding the answer to this fundamental question.
The C3 field was chosen as it points towards the South Galactic Pole, thus favouring the inclusion of thick disk and metalpoor stars. These stars are under-represented in other seismic fields, notably the Kepler field (see Rendle et al. (2019) for the space distribution of the stars of K2C3 and Kepler).
In addition, a study of chemical abundance evolution and binary detection was performed given the additional information pertaining to the masses and ages from K2, as well as the information of Gaia DR2.
The spectroscopic analysis was based on the combination of two distinct methods, both of which are part of the analysis procedure of the Gaia-ESO data. One method was based on spectrum synthesis and the other one was based on equivalent widths. Together, they allowed us to assess systematic uncertainties and reliability of results, thus providing a catalogue of spectral parameters and abundances that is accurate and not biased towards a given methodology. However, based on the intra-node abundance assessment and the corrections that were needed for seven of the elements, some possible systematic offsets in individual element abundances may remain.
The combination of seismology and spectroscopy allows us to analyse trends of abundances as a function of age, finding, as expected, that more metal-poor and α-capture enhanced stars are older than stars with solar abundances. The addition of the Gaia DR2 data also allows us to identify possible binary stars: the ones that are suspiciously more massive for their metal-content, those rotating unusually faster, or those having differences in radial velocities that are larger than expected.
Catalogues of parameters and abundances for which both seismology and spectroscopy are available are few. Pioneering works in this direction are the APOKASC sample (Pinsonneault et al. 2014(Pinsonneault et al. , 2018 and the CoRoGEE sample (Anders et al. 2017). Hawkins et al. (2016) published a catalogue with abundances for this sample considering the infrared APOGEE spectra based on temperatures determined photometrically and surface gravities from scaling relations of Kepler.
Another group at the forefront of this work is led by Valentini et al. (2017), who published a study similar to the one presented here (i.e. iterating between spectroscopy and seismology to reach to consistent parameters). They analysed 87 giants observed by RAVE. The spectra are of intermediate resolution (R 7000) and cover a short wavelength range around the Ca ii triplet, equivalent to Gaia RVS. In Valentini et al. (2017) it is possible to see how the overall stellar parameters improve when seismology is available. The group further performed a followup, detailed, high-resolution study of four metal-poor stars of their sample, which is presented in Valentini et al. (2019). The full RAVE catalogue will soon become public (Steinmetz et al. 2020). Targets with seismology have proven to be fundamental for the calibration of surface gravities of RAVE (Kunder et al. 2017). Similarly, Gaia-ESO is in the process of using CoRoT data as part of the calibration plan (Pancino et al. 2017), with the intention to publish the CoRoT@Gaia-ESO catalogue (Masseron & et al. 2020). In Valentini et al. (2016), the UVES and GIRAFFE stars were analysed following the same procedure as for the RAVE study.
Recently, Sharma et al. (2019) presented a study of the age and metallicity gradient of a large sample of K2 stars observed in GALAH. The study includes the same fields we analyse here, but GALAH spectra are currently not public to the entire community. The parameters and abundances are part of GALAH DR2 (Buder et al. 2018), although they were not derived using the seismic parameters as a prior.
In addition, Nissen et al. (2017) published a catalogue of parameters and abundances derived self consistently with seismology using high-resolution optical spectra of HARPS-N for a sample of solar twins in Kepler. Finally, Morel et al. (2014) presents the study of high resolution optical spectra of 20 red giants that were targeted by CoRoT.
Our work is a novel contribution as the largest highresolution optical catalogue with parameters and abundances derived homogeneously. We have taken advantage of forefront spectroscopic and astereoseismic methodologies, along with the sample targets of the outer regions of the disk to capture soughtafter old and metal-poor stars. Our catalogue of parameters and abundances is available online, including the results at a line-byline and method-by-method level to ensure reproducibility (see Tables B.1 and B.2 for list of columns). The spectra are public via the ESO archives. In addition, the chosen field of K2 was selected so that our catalogue would include more thick disk stars Article number, page 19 of 33 A&A proofs: manuscript no. k2ges than other samples for which spectroscopy and astereoseismology are available.
While the interior of stars gives us clues as to their evolutionary phase (based on the derived mass and age), the exterior of stars tell us about their formation site (chemistry) and stellar motion (radial velocities). Individually, these latter two types of data can provide the means to calibrate the analysis of large samples of stars, however the combination of both provides an even more powerful basis upon which to make a more effective and robust analysis that is key to constraining models of galactic formation and evolution.
Bolometric corrections were interpolated from the tables of Houdashelt et al. (2000). These corrections are tabulated for the K magnitude in the CIT/CTIO system (Elias et al. 1982). Thus, first the Ks magnitude of 2MASS were converted to K CIT/CTIO using the relations provided by Carpenter (2001).
The bolometric corrections (BC) are mostly independent of log g. First the two tables of metallicity values closest to the metallicity of the star (usually [Fe/H] = 0.00 and [Fe/H] = -0.50) are found. These tables are used to linearly interpolate the corrections to the T eff of the star. Then a second linear interpolation is carried out to the metallicity of the star.
The E(B-V) is converted to A Ks using the coefficients of Mc-Call (2004). In this case, A Ks /E(B-V) = 0.350.
The absolute magnitude (Mabs) is then: where +0.024 is the conversion to K CIT/CTIO and the parallax (π) has to be in arcsec (hence the 0.001) This implicitly assumes that distance = 1/parallax which is only then calculated if the relative uncertainty of the parallax is better than 15%.
Then the bolometric magnitude (Mbol) and luminosity of the star relative to the Sun log(L * /L ) are: The stellar masses were then estimated using the UniDAM code (see http://www2.mps.mpg.de/homes/mints/unidam.html Mints & Hekker 2017). It interpolates masses from PARSEC isochrones (Bressan et al. 2012) using a Bayesian scheme. A first estimate of mass was computed with this code using the spectroscopic parameters (T eff , log g, [Fe/H]), 2MASS magnitudes, parallaxes and extinction.
With this first mass estimate, we compute a first estimate of log g using Equation D. The new log g values are then given again to UniDAM, to recompute the mass estimates. These iterations were repeated (a total of four iterations was sufficient) until the estimated mass was consistent with the estimated log g.
UniDAM gives an estimate of the error of the mass, and this was also used in the computation of the total error budget in log g. Uncertainties are from 10000 Monte Carlo simulations, assuming Gaussian uncertainties in π, K 2MASS , stellar mass (with error coming from UniDAM), and ±0.02 mag for the bolometric correction. No uncertainty was assumed in E(B-V), nor for T eff , or log g .     of the 90 stars analysed in the K2 Gaia-ESO special project. Listed are: the Gaia-ESO CNAME; EPIC identifier; the spectroscopic instrument used; the IRFM T eff ; the preliminary set of spectroscopic parameters (Prelim. Spec.) sourced from the iDR5 recommended parameters where available (noted as *) or otherwise parameters from the Nice Node analysis or the GIRAFFE reduction pipeline radial velocity determination; the preliminary seismic log g; the initial spectroscopic T eff and [Fe/H] from iSpec; the initial seismic log g.