# TOPoS VI. The metal-weak tail of the metallicity distribution functions of the Milky Way and of the Gaia-Sausage-Enceladus structure

P Bonifacio, L Monaco, S Salvadori, E Caffau, M Spite, L Sbordone, F Spite, H.-G Ludwig, P Di Matteo, M Haywood, et al.

► To cite this version:

P Bonifacio, L Monaco, S Salvadori, E Caffau, M Spite, et al., TOPoS VI. The metal-weak tail of the metallicity distribution functions of the Milky Way and of the Gaia-Sausage-Enceladus structure. Astronomy and Astrophysics - A&A, 2021, 10.1051/0004-6361/202140816. obspm-03226191v3

HAL Id: obspm-03226191

https://hal-obspm.ccsd.cnrs.fr/obspm-03226191v3

Submitted on 15 Jul 2021

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Public Domain
# TOPoS VI. The metal-weak tail of the metallicity distribution functions of the Milky Way and the Gaia-Sausage-Enceladus structure

P. Bonifacio<fnref n="1" />, L. Monaco², S. Salvadori³,⁴, E. Caffau<fnref n="1" />, M. Spite<fnref n="1" />, L. Sbordone⁵, F. Spite<fnref n="1" />, H.-G. Ludwig⁶, P. Di Matteo<fnref n="1" />, M. Haywood<fnref n="1" />, P. François<fnref n="1" />,⁷, A.J. Koch-Hansen⁸, N. Christlieb⁶, and S. Zaggia⁹

¹ GEPI, Observatoire de Paris, Université PSL, CNRS, Place Jules Janssen, F-92195 Meudon, France

² Departamento de Ciencias Fisicas, Universidad Andres Bello, Fernandez Concha 700, Las Condes, Santiago, Chile

³ Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, Via G. Sansone 1, I-50019 Sesto Fiorentino, Italy

⁴ INAF/Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, I-50125 Firenze, Italy

⁵ European Southern Observatory, Casilla 19001, Santiago, Chile

⁶ Zentrum für Astronomie der Universität Heidelberg, Landessternwarte, Königstuhl 12, 69117 Heidelberg, Germany

⁷ UPJV, Université de Picardie Jules Verne, 33 rue St Leu, 80080 Amiens, France

⁸ Zentrum für Astronomie der Universität Heidelberg, Astronomisches Rechen-Institut, Mönchhofstr. 12, 69120 Heidelberg, Germany

⁹ INAF—Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, Padova, I-35122, Italy

Received; accepted

## ABSTRACT

**Context.** The goal of the Turn-Off Primordial Stars survey (TOPoS) project is to find and analyse turn-off (TO) stars of extremely low metallicity. To select the targets for spectroscopic follow-up at high spectral resolution, we relied on low-resolution spectra from the Sloan Digital Sky Survey (SDSS).

**Aims.** In this paper, we use the metallicity estimates we obtained from our analysis of the SDSS spectra to construct the metallicity distribution function (MDF) of the Milky Way, with special emphasis on its metal-weak tail. The goal is to provide the underlying distribution out of which the TOPoS sample was extracted.

**Methods.** We made use of SDSS photometry, Gaia photometry, and distance estimates derived from the Gaia parallaxes to derive a metallicity estimate for a large sample of over 24 million TO stars. This sample was used to derive the metallicity bias of the sample for which SDSS spectra are available.

**Results.** We determined that the spectroscopic sample is strongly biased in favour of metal-poor stars, as intended. A comparison with the unbiased photometric sample allows us to correct for the selection bias. We selected a sub-sample of stars with reliable parallaxes for which we combined the SDSS radial velocities with Gaia proper motions and parallaxes to compute actions and orbital parameters in the Galactic potential. This allowed us to characterise the stars dynamically, and in particular to select a sub-sample that belongs to the Gaia-Sausage-Enceladus (GSE) accretion event. We are thus also able to provide the MDF of GSE.

**Conclusions.** The metal-weak tail derived in our study is very similar to that derived in the H3 survey and in the Hamburg/ESO Survey. This allows us to average the three MDFs and provide an error bar for each metallicity bin. Inasmuch as the GSE structure is representative of the progenitor galaxy that collided with the Milky Way, that galaxy appears to be strongly deficient in metal-poor stars compared to the Milky Way, suggesting that the metal-weak tail of the latter has been largely formed by accretion of low-mass galaxies rather than massive galaxies, such as the GSE progenitor.

**Key words.** Stars: Population II - Stars: abundances - Galaxy: abundances - Galaxy: halo

## 1. Introduction

Considerable information on a galaxy's formation and evolution is imprinted in its metallicity distribution function (MDF). To our knowledge, the first attempt to determine the MDF of the Galaxy was that of van den Bergh (1962). Since Wallerstein (1962) robustly established the correlation between metallicity and ultraviolet excess, $\delta(U-B)$, van den Bergh (1962) proceeded to study the distribution of the ultraviolet excess of a sample of

142 stars, finding that metal-poor stars are very rare, and that no star in his sample had a metallicity below -0.85¹. It was quickly appreciated that the ultraviolet excess, as a metallicity indicator, saturates around a metallicity of -2.0, thus calling for spectroscopic metallicities, or photometric indices based on intermediate or narrow-band photometry. Simple models of Galactic chemical evolution and the long life time of low-mass stars led to the expectation that it should be possible to identify stars below -2.0 and -3.0 in dedicated observational programmes,

<footnote>* Spectroscopic and photometric metallicities derived and discussed in this paper as well as orbital actions computed and discussed this paper are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/</footnote>

<footnote n="1">¹ One should be aware that in the literature 'metallicity' may refer to slightly different quantities; often it is [Fe/H], especially in the older papers, but not always. In Sect. 2, we provide the definition of metallicity used in this paper.</footnote>
and hence considerable effort was devoted specifically to the discovery of metal-poor stars.

Extensive objective-prism surveys (Bond 1970; Bidelman & MacConnell 1973; Bond 1980) failed to detect stars below –3.0 and led Bond (1981) to conclude that 'very few stars with [Fe/H] < –3 may exist. Clearly one must go to quite faint objects before encountering such Population III stars'<fnref n="2" />, in spite of the fact that Bessell (1977) had already claimed [Ca/H] = –3.5 for the giant CD –38° 245, based on a low-resolution spectrum.

Hartwick (1976) made use of spectroscopic metallicities determined for a few individual stars in several globular clusters and showed that their MDF extends at least down to a metallicity of –2.2. Cayrel (1986) showed that if star formation started on a gas cloud of $10^6 M_\odot$, the O-type stars that formed would rapidly evolve and enrich the cloud to a level of $Z \approx 0.2 Z_\odot$, thus explaining the paucity of metal-poor stars. This scenario is still considered valid today, modulo the introduction of a dark matter halo, which allows us to relax the hypothesis on the baryonic mass. The question remains, however, as to whether low-mass stars were also formed at the same time as O-stars.

Beers et al. (1985) demonstrated the efficiency of coupling a wide-field objective-prism survey with low-resolution (R ~ 2 000) spectroscopic follow-up. From the low-resolution spectra, they were able to extract metallicities using a spectral index centred on the Ca II H&K lines, which resulted in a catalogue of 134 metal-poor stars. This sample extended down to metallicity –4.0 (star CS 22876-32). The MDF was not explicitly shown in that paper, but it was concluded that the slope of the metal-poor tail of the MDF is approximately constant. Beers (1987) presented the MDF of that sample of 134 stars. The success of Beers et al. (1985) in finding extremely metal-poor stars on one hand lies in the fainter magnitude limit of their survey (B ~ 15, compared to B ~ 11.5 of Bond 1980), and on the other hand in the choice of shortening the spectra on the plate by using an interference filter that selected only the UV-blue wavelengths, reducing the fraction of overlapping spectra, and also minimising the diffuse sky background, hence allowing for a fainter magnitude limit.

Carney et al. (1987) demonstrated how metallicities could be reliably derived from high-resolution (R~ 30 000), low signal-to-noise spectra, resulting in a catalogue of 818 stars that provided an MDF extending to –3.0 (Laird et al. 1988). Presumably that was the first time at which the derived MDF clearly displayed a bi-modality, although Laird et al. (1988) did not comment on it.

Ryan & Norris (1991a) used low-resolution spectroscopy and defined a spectral index, centred on Ca II H&K lines, slightly different from that defined by Beers et al. (1985), and they obtained metallicities for a sample of 372 kinematically selected halo stars. The MDF was presented and discussed in Ryan & Norris (1991b). The MDF extends down to ~ –4.0 and is not bi-modal, this can however be attributed to the kinematical selection that excluded thick disc stars.

In the last twenty years, there have been many attempts to define the MDF of the Galactic halo, a non-exhaustive list includes Schörck et al. (2009), Yong et al. (2013b), Allende Prieto et al. (2014), Youakim et al. (2020), Naidu et al. (2020), and Carollo & Chiba (2020), and we compare our results with some of theirs. Our group has long been exploiting the spectra of the Sloan Digital Sky Survey (hereafter SDSS, York et al. 2000; Adelman-McCarthy et al. 2008; Ahn et al. 2012; Alam et al. 2015) to select extremely metal-poor stars (Caffau et al. 2011; Bonifacio et al. 2011; Caffau et al. 2012; Bonifacio et al. 2012; Bonifacio et al. 2013a,b; Bonifacio et al. 2015; Caffau et al. 2016; Bonifacio et al. 2018; François et al. 2018a,b, 2020). The TOPoS project (Caffau et al. 2013b) was designed to exploit the high success rate of the selection of extremely metal-poor stars to study their detailed abundances. The high-resolution follow-up observations were carried out with XShooter and UVES at the ESO VLT telescope (Bonifacio et al. 2015; Caffau et al. 2016; Bonifacio et al. 2018; François et al. 2018b). The selection technique is described in Ludwig et al. (2008) and Caffau et al. (2011, 2013b). It made use of spectra of turn off (TO) stars from SDSS DR6 (Adelman-McCarthy et al. 2008), DR9 (Ahn et al. 2012), and DR12 (Alam et al. 2015), selected as described in Appendix A. A version of the raw MDF from SDSS DR6 has been shown in Ludwig et al. (2008), Figure 3. Some interesting features were already apparent, including the bi-modality of the MDF, which is similar to that of Laird et al. (1988) identified at the time as 'halo' and 'thick disc'. However, it was clear that the metal-weak tail of the raw MDF was unrealistic, providing a flat MDF at low metallicity. Visual inspection of the spectra allowed us to identify many white dwarfs with a Ca II K line absorption (likely due to circumstellar material or material accreted from a planet orbiting the white dwarf). However, we were unable at that time to robustly discriminate against these objects by an automated analysis of the spectra, even when also using photometric data. In addition, we were well aware that the selection function of the SDSS spectra we were using was difficult to define. For these reasons, we have previously not published an analysis of the MDF derived from our analysis of SDSS spectra.

The situation changed dramatically with the Gaia (Gaia Collaboration et al. 2016) Data Release 2 (Gaia Collaboration et al. 2018; Arenou et al. 2018). The Gaia photometry and astrometry allowed us to remove the white dwarfs contaminating our sample, and to clean the sample from stars with inaccurate SDSS photometry. Distances based on Gaia parallaxes allowed us to also derive a photometric metallicity estimate from SDSS photometry. This was used to quantify the metallicity bias present in our SDSS spectroscopic sample. With these improvements, we now present our estimate of the MDF. In Sect. 6, we also provide a detailed comparison with some of the previously published MDFs.

## 2. Selection and analysis of SDSS spectra

We selected from SDSS data release 12 (Alam et al. 2015) stars with TO colours, obtaining a set of 312 009 spectra of 266 442 unique stars (see Appendix A). The average S/N of the spectra at the effective wavelength of the g band is 24 ± 17. Both the SDSS and BOSS spectrographs have a red and a blue arm. The resolving power, $\lambda/\Delta\lambda$, is 1850 in the blue and 2200 in the red for the SDSS spectrograph, and 1560 in the blue and 2300 in the red for the BOSS spectrograph. We used the reddening provided in the SDSS catalogue, which is derived from the maps of Schlegel et al. (1998). At the time the data were retrieved from the SDSS archive, as detailed in Appendix A, Gaia distances were not yet available, nor were 3D maps based on these distances. Our sample spans a distance range of roughly 1 to 5 kpc. Available 3D maps (e.g. Lallement et al. 2019) are limited to 3 kpc, and in many directions at high galactic latitude, only to less than 1 kpc. Usage of those maps would therefore introduce a distance- and coordinate-dependent bias to the sample. All the aforementioned spectra were analysed with MyGIsFOS (Sbordone et al. 2014), using a fixed $T_{\text{eff}}$ derived from $(g-z)_0$ with a third-order polynomial.

<footnote n="2">² Bond (1981) refers to any star with [Fe/H] < –3 as Population III.</footnote>
nomial fit to the theoretical colour for metallicity −4.0, and an
α-enhancement of +0.4.<fnref n="3" />

The surface gravity was kept fixed at log g = 4.0, the micro-
turbulence was set to 1.5 km s⁻¹. The first guess for each spec-
trum was a metallicity of −1.5. The α-element enhancement was
kept fixed at [α/Fe] = +0.4.

MyGIsFOS used a special grid of synthetic spectra. The metal-poor part was computed explicitly using the OSMARCS (Gustafsson et al. 2008) TOPoS grid (see Caffau et al. 2013b, for details). This was complemented with MARCS interpolated models at higher metallicity. The synthetic spectra were computed with version 12.1.1 of turbospectrum (Alvarez & Plez 1998; Plez 2012).

The synthetic spectra contained two wavelength inter-
vals: 380–547 nm and 845–869 nm. They were broadened by
162 km s⁻¹ and 136 km s⁻¹ for the red and blue SDSS spectra,
and 192 km s⁻¹ and 132 km s⁻¹ for the red and blue BOSS spec-
tra, respectively. This corresponds to the spectral resolution of
the spectra.

A special line list and continuum files were used, in which
many features, including the Ca II K and Ca II IR triplet were la-
belled as 26.00, thus treating them as if they were iron. We de-
liberately exclude the features due to C-bearing molecules from
our estimate (G-band and Swan band), although the C abun-
dance is estimated from them. This was done because the reso-
lution and S/N of the SDSS spectra are not sufficient for a multi-
element analysis, especially at low metallicity; in many spectra,
the only measurable feature is the Ca II K line. So what MyGIs-
FOS provides in this case is an overall ‘metallicity’ that is essen-
tially a mean of Mg, Ca, Ti, Mn, and Fe. The choice of [α/Fe]
only impacts the lower metallicities, where the metallicity esti-
mate is dominated by the Ca II line. At higher metallicities, the
Fe feature holds more weight, and our estimate is not biased. We
verified this by looking at our estimates based on the SDSS spec-
tra for a sample of seven stars from Caffau et al. (2018) with low
[α/Fe]. At low metallicity, there is a large scatter, but at metallic-
ity > −1.0 there is a tight correlation with no offset. Throughout
this paper, we refer to this quantity as metallicity, and it should
not be confused with Z, the mass fraction of metals. Our metal-
licity estimate is closer to [Fe/H] than to Z, and the difference
can be very large in the case of CEMP stars, which, as discussed
below, we cannot detect in a complete and reliable manner.

It should be kept in mind that our procedure makes a spe-
cific assumption of [α/Fe] as a function of [Fe/H]. We assumed
[α/Fe] = +0.4 for [Fe/H] ≤ −0.5, and [α/Fe] = 0.0 otherwise.
Stars that do not follow this simplistic law will be assigned a
wrong metallicity. We stress that our analysis was specifically
designed to detect low-metallicity stars for follow-up at higher
resolution, and in this respect it has been extremely successful.
In Sect. 3, we show that it nevertheless provides a useful metal-
licity estimate at higher metallicities as well.

The combination of the low resolution and low S/N of the SDSS spectra implies that C can be measured reliably only when its spectral features are very strong. We did not make any attempt to remove CEMP stars from our sample, since this would result in a strong bias. Many CEMP stars would be missed simply because the S/N of the spectrum is too low. The drawback is that for some (but not all) CEMP stars our metallicity estimate is far too high, due to the presence of strong lines of C-bearing molecules, as we show in Sect. 3

<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Metallicities derived from SDSS spectra (LR) compared to those of independent analysis of the same stars (ext). The black dots are the TOPoS sample (see text); the data are provided in Table B.1. The red dots are the PASTEL sample (see text); the data is provided in Table B.2. The green dots are data from Caffau et al. (2018), and open symbols are stars with uncertain abundance according to Caffau et al. (2018). The blue dots are data from Caffau et al. (2020). The stars circled in black are removed from the average (see text).</figcaption></figure>

**3. Precision and accuracy of the MyGIsFOS**
**metallicities**

**3.1. Precision**

We took advantage of the fact that 15 496 stars of our sample
have multiple spectra (i.e. up to 33) to estimate the precision
of our metallicities. The root-mean-square error of the derived
metallicity for stars with multiple spectra is as large as 1.3 dex;
however, its mean value is 0.10 dex, with a standard deviation of
0.11 dex. We hence conclude that the internal precision of our
metallicities is 0.1 dex.

3.2. Accuracy

We now consider the accuracy of our metallicities. The com-
parison is based on three external datasets. The TOPoS sample
provides an analysis of the same stars based on higher resolution
spectra obtained by our group. We have a total of 70 stars (we
excluded stars with upper limits in the high-resolution analysis,
as well as stars for which MyGIsFOS did not converge), and they
are summarised in Table B.1. There are 88 stars extracted from
the PASTEL catalogue (Soubiran et al. 2016) in common with
our sample. We excluded the stars already present in the TOPoS
sample. The data of this sample is summarised in Table B.2. The
low-resolution sample features the stars observed with FORS at
VLT and GMOS at Gemini by Caffau et al. (2018), with resolv-
ing power *R* ~ 5000 and high S/N (typically above 100, 30 stars
in common), and by Caffau et al. (2020), with FORS at VLT, a
resolving power of *R* ~ 2000, and a high S/N (again above 100,
three stars in common). In the Caffau et al. (2018) sample, one
star is analysed independently using two spectra. The results are
very similar, and in the discussion we keep both points.

We stress that the comparison with the low-resolution sample
is relevant, since on one hand the spectra have a very high S/N,
and on the other hand they have been analysed with a full multi-

<footnote n="3">³ T<sub>eff</sub> = 7.2165 × 10³ − 3.2366 × 10³(g − z)<sub>0</sub> + 1.1578 × 10³(g − z)<sub>0</sub>² −
1.8493 × 10²(g − z)<sub>0</sub>³</footnote>
element approach; hence, in that analysis Fe and the α elements
are measured independently of one another.

In Fig. 1, we show the literature [Fe/H] as a function of
the metallicity derived from the SDSS spectra by MyGIsFOS.
It is obvious that our abundances are tightly correlated with
the external abundances. In some cases, our abundances are
very high compared to the literature value. We highlight eight
highly deviating stars in the plot. Three of them are well-
known CEMP stars: SDSS J029+0238 (Bonifacio et al. 2015;
Caffau et al. 2016, A(C)=8.50), SDSS J1349-0229 (Aoki et al.
2013, A(C)=8.20), and SDSS 1349-0229 (Caffau et al. 2018,
A(C)=6.18). Three other stars (i.e. SDSS J014828+150221,
SDSS J082511+163459, and SDSS J220121+010055) have a
strong G-band that is clearly visible in the SDSS spectrum and
is consistent with a CEMP classification. We suspect that the re-
maining two stars are also carbon-enhanced, although the qual-
ity of the available spectra has not allowed us to confirm that.
The mean difference (i.e. the external sample minus ours) of
the whole comparison sample (192 stars) is −0.21 dex with a
standard deviation of 0.55 dex. If we exclude the eight strongly
deviating stars, reducing the sample to 184 stars, the mean is
−0.12 dex with a standard deviation of 0.35 dex.

We hence adopted 0.35 dex as the accuracy of our spectro-
scopic metallicities. This is probably an upper limit to the actual
accuracy, because it does not take into account the errors of the
external metallicities. There is also an indication that our metal-
licity estimates are 0.1–0.2 dex too high.

4. Cross-match with Gaia and down selection

As soon as the metallicities from our processing of the SDSS
spectra became available, we began to look at the MDF result-
ing from this sample. It was obvious that there was a clear ex-
cess of metal-poor stars, which we suspect is due to the presence
of cool white dwarfs, some of which show Ca II K absorption,
as suggested by visual inspection of a few spectra. We there-
fore cross-matched our catalogue with Gaia to use parallaxes to
remove white dwarfs from the sample. To construct the Gaia
colour-magnitude diagram, we need to correct the colours for
the effect of reddening. To be on the same scale as SDSS, we
used the extinction in the g band, A<sub>g</sub>, provided in the SDSS
catalogue. The reddening on the Gaia colours was estimated
in three steps: A<sub>V</sub> = A<sub>g</sub>/1.161 ; E(B − V) = A<sub>V</sub>/3.1; correct
E(B − V) as in Bonifacio et al. (2000), A<sub>VC</sub> = 3.1 × E(B − V)<sub>C</sub>;
(G<sub>BP</sub> − G<sub>RP</sub>)<sub>0</sub> = 1.289445 E(B − V)<sub>C</sub> ; G<sub>0</sub> = G − 0.85926 A<sub>VC</sub>. The
extinction coefficients are adopted from the PARSEC web site<sup>4</sup>
(Bressan et al. 2012) and correspond to a G2V star, assuming the
extinction laws by Cardelli et al. (1989) plus O’Donnell (1994),
and a total-to-selective extinction ratio of 3.1. The cross-match
with Gaia left us 290 498 spectra for 245 838 unique stars. For a
large part of the sample, the parallaxes were very poor, leading
to implausibly high luminosities for such stars.

The Gaia photometry, with a typical intrinsic precision of
a few mmag, gives us the opportunity to clean the sample of
stars with inaccurate SDSS photometry. We started from the
(g − z)₀ versus (G<sub>BP</sub> − G<sub>RP</sub>)₀ plane, shown in Fig. 2. We fitted a
straight line<fnref n="5" /> and selected only the stars for which the observed
(G<sub>BP</sub> − G<sub>RP</sub>)₀ was within ±3σ of the fitted value. We assume
that the stars that fall out of this strip have either systematic er-
rors in their g − z colours or are highly reddened. Either way, the

<figure><img src="image_2.png" /><figcaption>Fig. 2. (g − z)₀ versus (G<sub>BP</sub> − G<sub>RP</sub>)₀. The full dataset is shown as grey dots; the green dots are the selected stars.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 3. CMD with the selection using Gaia photometry and parallaxes. Two sets of Parsec isochrones (Bressan et al. 2012) of ages 1–14 Gyr and metallicities −0.5 (red) and −1.5 (blue) are shown. The absolute magnitudes were computed assuming the distances of Bailer-Jones et al. (2018).</figcaption></figure>

T<sub>eff</sub> used in the MyGIsFOS analysis is wrong; therefore, we dis-
carded these stars. To exclude all the very-low-luminosity stars,
which we assume to be mostly cool white dwarfs, we also cut
on absolute magnitude, excluding stars with G<sub>0,abs</sub> > 7.0. Since
many of the stars in our sample have poor parallaxes, we de-
cided to use the distance estimates of Bailer-Jones et al. (2018)
to perform this cut. These geometric distances are the result of
a Bayesian estimate, based on the prior of an exponentially de-
creasing space density. The exponential length scale varies as a
function of Galactic coordinates (l, b). For each star, an estimated
distance is provided, r<sub>est</sub>, as well as a ‘highest distance’ r<sub>hi</sub> and
‘lowest distance’ r<sub>lo</sub>. Over the whole catalogue (1.33×10<sup>9</sup> stars),
the 50th percentile of the difference r<sub>hi</sub> − r<sub>lo</sub> is about 3.6 kpc.
This selection results in 171 645 spectra for which MyGIsFOS
has converged. This resulted in 139 943 unique stars.

The colour-magnitude diagram (CMD) of the selection is
shown in Fig. 3. The cleaned sample thus appears to be compat-
ible with a sample of old TO stars, as is our intention. In Fig. 4,
we show the histogram of the semi-dispersion of the geometric
distance estimates from Bailer-Jones et al. (2018).

<footnote n="4">⁴ http://stev.oapd.inaf.it/cgi-bin/cmd</footnote>

<footnote n="5">⁵ (G<sub>BP</sub> − G<sub>RP</sub>)<sub>0</sub> = 0.43742296099663 + 0.63260418176651(g − z)<sub>0</sub>;
rms = 0.06694397</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 4.** Histogram of the distance errors in pc, from the geometric distances of Bailer-Jones et al. (2018), estimated as $0.5 \times (D_{max} - D_{min})$</figcaption></figure>

<table><caption>**Table 1.** Target types and SDSS programs for which the spectra have been observed</caption><thead><tr><th>Type</th><th># of spectra</th></tr></thead><tbody><tr><td></td><td>6425</td></tr><tr><td>GALAXY</td><td>736</td></tr><tr><td>NA</td><td>6221</td></tr><tr><td>NONLEGACY</td><td>69917</td></tr><tr><td>QA</td><td>1308</td></tr><tr><td>QSO</td><td>5202</td></tr><tr><td>REDDEN_STD</td><td>10106</td></tr><tr><td>ROSAT_D</td><td>601</td></tr><tr><td>SERENDIPIT</td><td>31</td></tr><tr><td>SPECTROPHO</td><td>47930</td></tr><tr><td>STANDARD</td><td>4</td></tr><tr><td>STAR</td><td>20973</td></tr><tr><td>STAR_BHB</td><td>3032</td></tr><tr><td>STAR_WHTE</td><td>159</td></tr></tbody><thead><tr><th>Programme</th><th># of spectra</th></tr></thead><tbody><tr><td>boss</td><td>45250</td></tr><tr><td>sdss</td><td>42543</td></tr><tr><td>segue1</td><td>58298</td></tr><tr><td>segue2</td><td>26554</td></tr></tbody></table>

## 5. Bias in the observed spectroscopic metallicity distribution function

Our sample consists of 139 493 unique stars for which we derived metallicities from at least one SDSS spectrum. For stars with multiple spectra, each spectrum was independently analysed, and then the derived metallicities were averaged (see Sect. 3.1). Our spectroscopic sample contains stars that have been observed in four different SDSS programs, comprising 39 different surveys. It seems arduous to extract the bias in our sample starting from the selection criteria of the different surveys convoluted with our own selection criteria. In Table 1, we assembled the programmes and the types of our sample of spectra, with the numbers for each. We note that 6 425 spectra have a blank type, for 6221 targets the type is 'not available' (NA), and for 20 973 it is STAR. The largest fraction of our targets are NONLEGACY. We were unable to find any reference to this type of target in the SDSS documentation, and so we may guess that these are essentially 'fill fibres' to complete a plate, with no particular selection criterion. The sample is in fact dominated by these targets and spectro-photometric standards. We therefore employed a different route to estimate the bias in our sample.

<figure><img src="image_5.png" /><figcaption>**Fig. 5.** *p* index as a function of metallicity for stars with $5875\ K \le T_{eff} \le 6125\ K$ compared with theoretical predictions for $log\ g = 4.0$ (red) and $log\ g = 2.0$ (blue), derived from the Castelli & Kurucz (2003) grid of ATLAS 9 model atmospheres. For the observed stars, we show the spectroscopic metallicities.</figcaption></figure>

### 5.1. Photometric metallicities

We began by determining metallicities using a photometric index that is metallicity sensitive. The obvious choice would be *u* − *g*; however, we preferred to use a reddening-free index that can be formed combining *u* − *g* and *g* − *z*:

$$p = (u - g) - 0.5885914(g - z). \qquad (1)$$

The choice of the colours to form the *p* index is dictated by the need to have a metallicity-sensitive colour, the fact that *u* − *g* is the most metallicity-sensitive colour in *ugriz* photometry, and the need to factor out the temperature sensitivity of the same index. The choice of *g* − *z* has been done to be consistent with the colour chosen to determine effective temperatures. The coefficient in front of *g* − *z* is nothing other than the ratio of the two reddening coefficients as provided by Schlegel et al. (1998) in order to make *p* reddening free. In fact, *p* is nothing other than the *ugz* equivalent of the *Q* index introduced by Johnson & Morgan (1953) for *UBV* photometry.

The *p* index varies with effective temperature, surface gravity, and metallicity, as illustrated in Fig. 5, where we show the spectroscopic metallicities on the x-axis. It is clear that a spread in *p* at any given metallicity remains, even by fixing the effective temperature of the sample. This is because *p* is also sensitive to log *g*, and the set of a few stars that are around metallicity −1.5 and *p* ≈ 0.9 also contains horizontal-branch stars. In order to obtain a metallicity, we adopted a very simple approach: we derived $T_{eff}$ from the (*g* − *z*)₀ colour, and once the $T_{eff}$ is known, we derived the surface gravity using the parallax and the Stefan-Boltzmann equation, assuming a mass of 0.8M⊙. Once $T_{eff}$ and log *g* were known, we only needed to interpolate in a grid of theoretical values to obtain a metallicity from the observed *p* index.

We computed the theoretical *p* index for the Castelli & Kurucz (2003) grid of ATLAS 9 model atmospheres. The procedure needs to be iterated, since to determine the surface gravity one needs to know the bolometric magnitude. The latter can be derived from the Gaia $G_0$ magnitude and a theoretical bolometric correction that we computed for all the Castelli & Kurucz (2003) grid. To do this, one needs to know the metallicity. In this sense, the procedure is iterative: one starts from a guess of the metallicity to determine the bolometric correction and then computes the metallicity. At that point, with the new metallicity, one derives a new bolometric correction and iterates the procedure.
<figure><img src="image_6.png" /><figcaption>**Fig. 6.** Spectroscopic metallicity as a function of the metallicity derived from the *p* index. The blue line is a linear fit we made to the sub-set of stars along the highest density ridge.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 7.** Spectroscopic (red) and photometric (blue) samples in distance from the Galactic centre and height above the Galactic plane.</figcaption></figure>

In the vast majority of cases, a single iteration was sufficient. Bolometric corrections for the *G* magnitude are small for TO stars anyway. In order to accelerate the procedure, we did not interpolate in the grid of theoretical *p* values for the current *T*<sub>eff</sub> and log *g*, but we simply used the grid values that are nearest. In order to be able to use the method, we need parallaxes for all the stars. For the majority of the stars in our spectroscopic and photometric samples, the Gaia parallaxes are too uncertain to be used. Therefore, to derive the metallicities we used the distance estimates provided by Bailer-Jones et al. (2018), which were derived from Gaia parallaxes using a Bayesian estimate and a prior of exponentially decreasing stellar density. The procedure to determine the metallicities from the *p* index is explained in detail in Appendix C.

In Fig. 6, we show the comparison of photometric and spectroscopic metallicities. Aside from the large scatter, which we attribute mainly to the large uncertainty in the photometric metallicity estimate, the two variables are tightly correlated. We interactively selected a sub-set of 118 158 stars that lie along the highest density ridge and fitted a straight line (shown in blue) to derive the spectroscopic metallicity from the photometric metallicity.

<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Distribution of photometric metallicities for the photometric and spectroscopic samples. Each sample has been normalised so that the area under the histogram is equal to one.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Bias function (black line) and smoothed bias function (red line)</figcaption></figure>

By construction, these ‘adjusted’ photometric metallicities (i.e. not those derived directly from *p* using the procedure in Appendix C, but those derived from the latter and the straight line fit shown in Fig. 6) have a distribution that is similar to that of the spectroscopic metallicities.

## 5.2. Comparison with a larger photometric sample

In order to estimate the bias in our spectroscopic sample, we make the bold assumption that a sample selected from the SDSS DR12 photometry, using the same colour and absolute magnitude cuts that have defined our spectroscopic sample, is bias-free. Having selected the sample, we applied to it our photometric calibration. The distances were adopted from the Bailer-Jones et al. (2018) catalogue, and we discarded about one million stars whose colours fall outside our synthetic grid. The final sample, which we refer to in the following as ‘photometric’, consists of 24 037 008 unique stars.
<figure><img src="image_10.png" /><figcaption>**Fig. 10.** Raw metallicity distribution function (black line) and bias-corrected (red line).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Normalised metallicity distribution of our full spectroscopic sample (black) compared to the subset with good parallaxes (violet). The normalisation is such that the area under each histogram is equal to one.</figcaption></figure>

In Fig. 7, we show the two samples in the diagram of height above the Galactic plane and distance from the Galactic centre. We adopted a distance of 8.43 kpc for the Galactic centre (Reid et al. 2014) and a distance of the Sun above the Galactic plane of 27 pc (Chen et al. 2001). The photometric sample reaches higher values of $|z|$ and larger distances from the Sun. This is a result of the fact that the photometric sample has many more faint stars than the spectroscopic sample. We computed the bias using the full sample (see Sect. 5.3) and also using a sub-set of the photometric sample that covers the range in $|z|$ and $R_{GC}$ occupied by the spectroscopic sample more closely. The difference in the derived bias is negligible.

We can now compare the two metallicity distributions. In Fig. 8, the distributions of the photometric metallicities of the two samples are compared. It is clear that the spectroscopic sample is boosted in metal-poor stars. An obvious reason for this is the large number of spectrophotometric standards that are part of the sample. These stars are selected to have colours consistent

<figure><img src="image_12.png" /><figcaption>**Fig. 12.** Normalised MDFs for our sample with good parallaxes (black line), the sub-sample selected according to the criteria of H3 (red line) and the H3 (Naidu et al. 2020, green).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 13.** Metal-weak tail of the raw metallicity distribution function (black line), normalised at -1.95, compared with the MDFs of the Hamburg/ESO survey (Schörck et al. 2009, red), the Pristine survey (Youakim et al. 2020, blue), and the H3 survey (Naidu et al. 2020, green).</figcaption></figure>

with metal-poor F dwarfs to allow a reliable flux calibration of the spectra on the same plate. We suspect that many of the NONLEGACY targets were selected according to similar criteria. These two contributions alone could explain the boost of metal-poor stars in the spectroscopic sample.

### 5.3. Correcting for the bias

A straightforward way to correct our sample for the selection bias is to take the ratio of the two histograms in Fig. 8, (photometric sample)/(spectroscopic sample), to derive a 'bias function'. This function is shown in Fig. 9 and relies on photometric metallicities in both samples. In the same plot, we also show a smoothed version of the bias function (median filter with a 0.25 dex kernel) that avoids the oscillations seen at the extremes
<figure><img src="image_14.png" /><figcaption>**Fig. 14.** Metal-weak tail of the corrected metallicity distribution function (black line), normalised at [Fe/H] = -1.95, compared with the corrected Hamburg/ESO (Schörck et al. 2009, red) and the corrected Pristine (Youakim et al. 2020, blue) MDFs.</figcaption></figure>

of the metallicity range that are due to the small number of stars contained in each bin in either of the two histogram functions. We can then produce a bias-corrected version of the histogram of our spectroscopic metallicities, simply by multiplying the observed histogram with this bias function. This rests on the assumption that the bias in spectroscopic metallicities is the same as that in photometric metallicities. Or to put it another way, if we had an SDSS spectrum for each of the stars in the photometric sample and derived a spectroscopic metallicity for it, the bias derived dividing the histogram of these spectroscopic metallicities by that of the spectroscopic metallicities of the current spectroscopic sample would be the same as that shown in Fig. 9.

The result of this exercise is shown in Fig. 10, where the observed and bias-corrected MDFs are shown for our spectroscopic sample. Two features in the bias-corrected MDF are obvious: i) the bi-modality disappears and the metal-poor peak is now detectable only as a change in the slope of the MDF; and ii) the slope of the MDF from a metallicity ~ -1.6 is much more gentle than in the observed MDF.

# 6. The MDF

Previous investigations of the MDF used different techniques. To keep the comparison to our results manageable, we selected just three MDFs: Schörck et al. (2009, Hamburg/ESO Survey, hereafter HES), Naidu et al. (2020, H3 Survey), and Youakim et al. (2020, Pristine Survey). The first two are based on spectra, while the third uses a photometric metallicity estimate. Our MDF is provided in Table 2.

The HES MDF is derived from a sample of metal-weak candidates selected from the objective prism spectra (see Christlieb et al. 2008), then observed at low resolution ($R \sim 2000$) with several telescopes and spectrographs, to provide a sample of 1638 unique stars.

The H3 sample presented in Naidu et al. (2020) is selected with the following conditions: i) $15 < r < 18$; $\bar{\omega} - 2 \times \sigma_{\bar{\omega}} < 0.5$; $|b| > 40^{\circ}$. Naidu et al. (2020) claim that this sample is unbiased with respect to metallicity. The stars selected by H3 are observed with the multi-fibre Hectoechelle spectrograph (Szent-

gyorgyi et al. 2011) on the 6m MMT telescope (Conroy et al. 2019). The spectra cover the range 513–530 nm at a resolving power $R \sim 23 000$. Using the MINESweeper code (Cargile et al. 2020), an analysis of the spectra combined with photometry provides radial velocities, spectrophotometric distances, and abundances ([Fe/H] and [α/Fe]).

The Pristine sample is aimed at selecting the tip of the main-sequence turn-off of the old populations of the Galaxy using SDSS photometry. A magnitude cut was applied to select stars with $19 < g_0 < 20$ and colour $0.15 < (g-i) < 0.4$, combined with the exclusion of a region in the $(u-g)_0$ versus $(g-i)_0$ plane that is dominated by young stars (see Fig. 3 of Youakim et al. 2020). This sample amounts to about 80 000 stars. A bias correction of the sample is proposed to correct for the selection effects and for the photometric errors.

We did not compare to the MDF of Allende Prieto et al. (2014), because its metal-weak tail is essentially identical to ours, as is expected, since there is a substantial overlap between the two sets of spectra, although the spectra were analysed in different way. No bias correction was attempted in Allende Prieto et al. (2014). With respect to our sample, the Allende Prieto et al. (2014) sample lacks the metal-rich peak, and this is just a matter of selection criteria.

## 6.1. Comparison with H3

We carried out a more detailed comparison with H3, since it is the one that is more easily directly comparable to our results. Both are based on spectroscopy, and an MDF is available up to super-solar metallicities.

Here, and in Section 8, we need to define sub-samples that make use of parallaxes. For this reason, we need to define a sub-set with ‘good’ parallaxes. We define this sample as stars for which $\bar{\omega} > 3 \times \Delta\bar{\omega}$, which leaves us with a sample of 57 386 unique stars. In terms of metallicity distribution, this sub-sample favours the metal-rich component, so that now the most prominent peak is the metal-rich one, as shown in Fig. 11. On the other hand, the metal-weak tail of the two distributions is pretty much the same. In the following, we refer to this sub-sample as the good parallax sample. To have more meaningful comparison between the H3 MDF and our sample, we selected a sub-sample of our good parallax sample, which fulfils the same selection criteria as Naidu et al. (2020). This is composed of 28 165 stars, which is 49 % of our good parallax sample. We refer to this sub-sample as the H3-like sample. In Fig. 12, we compare the normalised MDFs of our good parallax and of our H3-like samples to that of Naidu et al. (2020).

From this comparison four facts are obvious: our sample is boosted in metal-poor stars with respect to an unbiased sample, confirming what was deduced from the comparison with our photometric sample. A H3-like selection, although it does not make any explicit reference to metallicity, is likely to introduce a bias in metallicity. This can be inferred by the comparison of our good parallax and H3-like samples. It is, however, obvious that this bias is minor, especially in the low-metallicity tail, and we ignore it in the following. The H3 MDF has a much more gentle slope of the metal-weak tail than our sample. The metal-rich peak in our full sample seems very similar to that in the H3 sample.
## 6.2. The metal-weak tail of the MDF

While the metal-rich portion of the MDF is the result of the complex chemical evolution of the Galaxy, we may assume that its metal-weak tail is much simpler, since many of the stars are the descendants of only a few generations of stars beyond the first one. For more than a decade, theoretical models have shown that the properties of the first stars affect the metal-weak tail of the MDF (e.g. Hernandez & Ferrara 2001; Prantzos 2003; Oey 2003; Karlsson 2006; Tumlinson 2006; Salvadori et al. 2007). Furthermore, this should hold true even if this metal-weak tail contains stars that have been formed in Milky Way satellites that later merged to form the present-day Milky Way (e.g. Salvadori et al. 2015).

To study the metal-weak tail of the MDF, Schörck et al. (2009) normalised their distributions at metallicity −1.95. This normalised distribution provides, for any lower metallicity, the fraction of stars in that metallicity bin, with respect to those in the bin centred at −1.95. We use the same approach, which also allows us to directly compare our results with theirs.

In Fig. 13, we compare the normalised metal-weak tail of the MDF, as observed, for *i*) our whole sample; *ii*) the Hamburg/ESO sample (Schörck et al. 2009); *iii*) the Pristine sample (Youakim et al. 2020); and *iv*) the H3 sample (Naidu et al. 2020). The H3 Survey adopts a target selection that is based only on apparent magnitudes, parallaxes, and Galactic latitude, thus they claim that their sample is unbiased with respect to metallicity. It is surprising that the metal-weak tail of their MDF is essentially identical to that of the Hamburg-ESO Survey down to metallicity −3.0, since the Hamburg/ESO Survey has been considered to be biased in favour of metal-poor stars (Schörck et al. 2009). The abrupt drop of the H3 MDF at −3.0 is not easy to understand, but it may be due to the fact that the size of their sample is too small to adequately sample the populations of lower metallicity.

The selection function for the four samples is clearly different. Schörck et al. (2009) provide three possible bias corrections for their observed MDF, while Youakim et al. (2020) provide two. In Fig. 14, we show the comparison with all the corrected MDFs. What is obvious is that while in the raw distributions ours is always *below* the Hamburg/ESO, once bias corrections are taken into account ours is *above* the Hamburg-ESO. As noticed in Sect. 5.3, this is due to the fact that in our corrected MDF, the slope is much more gentle than in our raw MDF; as a consequence, once we normalise at −1.95, the corrected MDF lies above the raw MDF. The reverse is true if we look at the non-normalised MDF, as shown in Fig. 10.

Above metallicity −3.0 the H3, the *uncorrected* HES, and our *corrected* MDF are essentially identical. If we make the assumption that our bias-corrected MDF is ground truth, then the metallicity bias of the H3 and HES MDFs must be small, and in fact negligible. We thus conclude that it is legitimate to use the average of these three MDFs for comparison with theoretical models. The assumption we are making is equivalent to stating that the bias corrections provided by Schörck et al. (2009) and Youakim et al. (2020) are not adequate to describe the actual bias in their stellar sample. The bias corrections proposed by Schörck et al. (2009) are based on models, while those of Youakim et al. (2020) are based on the statistical properties of the sample. One advantage of using the average is that in each metallicity bin the standard deviation can be used as an error estimate on the MDF. A further advantage is that the averaging process makes the statistical errors smaller. The H3 MDF is essentially undefined below −3.0; therefore, we suggest only using the average of our MDF and that of HES in this metallicity range. In Fig. 15,

<table><caption>**Table 2.** Raw and corrected metallicity distribution functions for the 139 493 unique stars of our sample.</caption><thead><tr><th>[M/H]</th><th>raw</th><th>corrected</th><th>[M/H]</th><th>raw</th><th>corrected</th></tr></thead><tbody><tr><td>-4.43</td><td>0</td><td>0</td><td>-2.12</td><td>1771</td><td>586</td></tr><tr><td>-4.38</td><td>0</td><td>0</td><td>-2.08</td><td>1990</td><td>640</td></tr><tr><td>-4.32</td><td>0</td><td>0</td><td>-2.03</td><td>2368</td><td>703</td></tr><tr><td>-4.28</td><td>0</td><td>0</td><td>-1.98</td><td>2656</td><td>787</td></tr><tr><td>-4.22</td><td>1</td><td>2</td><td>-1.92</td><td>2955</td><td>873</td></tr><tr><td>-4.18</td><td>0</td><td>0</td><td>-1.88</td><td>3196</td><td>944</td></tr><tr><td>-4.12</td><td>1</td><td>2</td><td>-1.83</td><td>3403</td><td>1005</td></tr><tr><td>-4.07</td><td>1</td><td>2</td><td>-1.77</td><td>3666</td><td>1083</td></tr><tr><td>-4.03</td><td>1</td><td>2</td><td>-1.73</td><td>4069</td><td>1202</td></tr><tr><td>-3.97</td><td>0</td><td>0</td><td>-1.67</td><td>4311</td><td>1273</td></tr><tr><td>-3.92</td><td>1</td><td>1</td><td>-1.62</td><td>4520</td><td>1383</td></tr><tr><td>-3.88</td><td>1</td><td>1</td><td>-1.58</td><td>4703</td><td>1733</td></tr><tr><td>-3.83</td><td>1</td><td>1</td><td>-1.52</td><td>4696</td><td>1902</td></tr><tr><td>-3.78</td><td>2</td><td>3</td><td>-1.48</td><td>4592</td><td>1896</td></tr><tr><td>-3.72</td><td>3</td><td>4</td><td>-1.42</td><td>4441</td><td>1957</td></tr><tr><td>-3.67</td><td>1</td><td>1</td><td>-1.38</td><td>4063</td><td>2069</td></tr><tr><td>-3.62</td><td>2</td><td>2</td><td>-1.33</td><td>3849</td><td>2198</td></tr><tr><td>-3.58</td><td>9</td><td>9</td><td>-1.27</td><td>3653</td><td>2529</td></tr><tr><td>-3.53</td><td>3</td><td>3</td><td>-1.23</td><td>3289</td><td>2772</td></tr><tr><td>-3.47</td><td>4</td><td>4</td><td>-1.17</td><td>3191</td><td>3393</td></tr><tr><td>-3.42</td><td>7</td><td>7</td><td>-1.12</td><td>3189</td><td>3698</td></tr><tr><td>-3.38</td><td>10</td><td>9</td><td>-1.08</td><td>3310</td><td>4283</td></tr><tr><td>-3.33</td><td>8</td><td>8</td><td>-1.02</td><td>3334</td><td>4678</td></tr><tr><td>-3.28</td><td>17</td><td>16</td><td>-0.98</td><td>3487</td><td>5230</td></tr><tr><td>-3.22</td><td>19</td><td>18</td><td>-0.93</td><td>3760</td><td>5874</td></tr><tr><td>-3.17</td><td>26</td><td>23</td><td>-0.88</td><td>3862</td><td>6761</td></tr><tr><td>-3.12</td><td>29</td><td>26</td><td>-0.82</td><td>3856</td><td>7215</td></tr><tr><td>-3.08</td><td>32</td><td>29</td><td>-0.77</td><td>3849</td><td>7354</td></tr><tr><td>-3.03</td><td>49</td><td>44</td><td>-0.73</td><td>3748</td><td>7161</td></tr><tr><td>-2.97</td><td>70</td><td>59</td><td>-0.68</td><td>3357</td><td>6414</td></tr><tr><td>-2.92</td><td>74</td><td>60</td><td>-0.62</td><td>3329</td><td>6361</td></tr><tr><td>-2.88</td><td>107</td><td>83</td><td>-0.57</td><td>3025</td><td>5780</td></tr><tr><td>-2.83</td><td>152</td><td>107</td><td>-0.52</td><td>2691</td><td>5142</td></tr><tr><td>-2.78</td><td>145</td><td>100</td><td>-0.47</td><td>2495</td><td>4728</td></tr><tr><td>-2.72</td><td>208</td><td>142</td><td>-0.43</td><td>2160</td><td>4006</td></tr><tr><td>-2.67</td><td>244</td><td>158</td><td>-0.38</td><td>1949</td><td>3516</td></tr><tr><td>-2.62</td><td>303</td><td>189</td><td>-0.32</td><td>1534</td><td>2768</td></tr><tr><td>-2.58</td><td>401</td><td>241</td><td>-0.28</td><td>1302</td><td>2349</td></tr><tr><td>-2.53</td><td>440</td><td>262</td><td>-0.22</td><td>1012</td><td>1826</td></tr><tr><td>-2.47</td><td>564</td><td>295</td><td>-0.17</td><td>728</td><td>1313</td></tr><tr><td>-2.42</td><td>656</td><td>296</td><td>-0.12</td><td>559</td><td>992</td></tr><tr><td>-2.38</td><td>822</td><td>364</td><td>-0.08</td><td>451</td><td>814</td></tr><tr><td>-2.33</td><td>988</td><td>395</td><td>-0.03</td><td>353</td><td>662</td></tr><tr><td>-2.28</td><td>1152</td><td>438</td><td>+0.03</td><td>287</td><td>551</td></tr><tr><td>-2.22</td><td>1300</td><td>489</td><td>+0.08</td><td>244</td><td>468</td></tr><tr><td>-2.17</td><td>1557</td><td>539</td><td>+0.12</td><td>185</td><td>355</td></tr></tbody></table>

we show the two average values with the standard deviation in each bin. The values are provided in Table 3.

## 7. 3D bias and metallicity maps

We can now take advantage of the precise measurements of metallicity and position for our large stellar sample to study the spatial distribution of stars with different metallicities as has been suggested by cosmological models (e.g. Salvadori et al. 2010). To achieve this goal, we need to derive the 3D bias function, meaning that we also need to study how the bias function depends on the position of the stars. To this aim, we adopt the
<figure><img src="image_15.png" /><figcaption>**Fig. 15.** Average metal-weak tail of the MDF, normalised at metallicity −1.95: black present paper (corrected) and Hamburg/ESO (uncorrected) red; the H3 is also included in the average.</figcaption></figure>

cylindrical coordinates $|z|$ and $\rho = \sqrt{(R_{GC}^2 - z^2)}$ and divide the space in annuli of radial width $dR = 0.5$ kpc and height $dz = 0.3$ kpc. In each annulus, the photometric MDFs from the spectroscopic and photometric samples have been derived and used to compute the local bias function by adopting metallicity bins of 0.25 dex. In the following, we use Fe to denote the metallicity, as defined in Sect. 2. From a global point of view, this means that we have constructed a 3D bias function, Bias = Bias(Fe, $\rho$, |z|), that we can use to compute the correct local MDF and to construct stellar metallicity maps.

In Fig. 16, we show a simplified view of the 3D bias function obtained by dividing the stellar populations into different metallicity ranges: extremely metal-poor stars (−4 ≤ [Fe/H] < −3, upper left panel), very metal-poor stars (−3 ≤ [Fe/H] < −2, upper right), metal-poor stars (−2 ≤ [Fe/H] < −1, lower left), and metal-rich stars (−1 ≤ [Fe/H] < 0, lower right). Colours in the maps illustrate the value of the bias obtained in each annulus, which is identified by a single pixel in the ($\rho$, |z|) space. A bias equal to zero means that there are no stars from the spectroscopic sample within the pixel, which is indeed the case for extremely metal-poor stars (upper left) located in the outermost regions (|z| > 3 kpc, $\rho$ > 10 kpc).

Bearing in mind that the colour-bar scale is different in each panel, we can notice several features in Fig. 16. First, the contrast between the minimum and maximum values of the bias is higher for very and extremely metal-poor stars ([Fe/H] < −2, upper panels), where it reaches values of Bias > 5 in several small and isolated regions. Our analysis shows that these extreme bias values are a result of the small statistic in the considered pixel. Conversely, the bias is much more uniform for metal-poor stars (lower left panel), only spanning an interval of two in spite of the more extended regions covered by the map. The contrast between minimum and maximum bias increases again for metal-rich stars, as seen in the lower right panel of Fig. 16, where we can notice three zones with Bias ≥ 2. The high contrast between the bias in these isolate small regions, and the rest of the space is at the origin of the oscillations seen in the (1D) bias function analysed in Sect. 5.3 (Fig. 9).

<table><caption>**Table 3.** Average MDF, normalised at metallicity −1.95. Aver2 is the average of the present MDF (corrected) and the Hamburg/ESO (uncorrected), $\sigma_2$ is the standard deviation, Aver3 is the average of the previous two and the H3 MDF, and $\sigma_3$ is the standard deviation.</caption><thead><tr><th>[M/H]</th><th>Aver2</th><th>σ<sub>2</sub></th><th>Aver3</th><th>σ<sub>3</sub></th></tr></thead><tbody><tr><td>−4.30</td><td>0.000</td><td></td><td></td><td></td></tr><tr><td>−4.25</td><td>0.000</td><td>0.001</td><td></td><td></td></tr><tr><td>−4.20</td><td>0.000</td><td>0.001</td><td></td><td></td></tr><tr><td>−4.15</td><td>0.000</td><td>0.001</td><td></td><td></td></tr><tr><td>−4.10</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−4.05</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−4.00</td><td>0.000</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.95</td><td>0.000</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.90</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.85</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.80</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.75</td><td>0.002</td><td>0.002</td><td></td><td></td></tr><tr><td>−3.70</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.65</td><td>0.001</td><td>0.001</td><td></td><td></td></tr><tr><td>−3.60</td><td>0.004</td><td>0.004</td><td></td><td></td></tr><tr><td>−3.55</td><td>0.004</td><td>0.005</td><td></td><td></td></tr><tr><td>−3.50</td><td>0.002</td><td>0.005</td><td></td><td></td></tr><tr><td>−3.45</td><td>0.003</td><td>0.007</td><td></td><td></td></tr><tr><td>−3.40</td><td>0.005</td><td>0.009</td><td></td><td></td></tr><tr><td>−3.35</td><td>0.005</td><td>0.011</td><td></td><td></td></tr><tr><td>−3.30</td><td>0.007</td><td>0.015</td><td></td><td></td></tr><tr><td>−3.25</td><td>0.010</td><td>0.019</td><td></td><td></td></tr><tr><td>−3.20</td><td>0.012</td><td>0.028</td><td></td><td></td></tr><tr><td>−3.15</td><td>0.015</td><td>0.037</td><td></td><td></td></tr><tr><td>−3.10</td><td>0.017</td><td>0.044</td><td></td><td></td></tr><tr><td>−3.05</td><td>0.022</td><td>0.058</td><td></td><td></td></tr><tr><td>−3.00</td><td>0.031</td><td>0.060</td><td></td><td></td></tr><tr><td>−2.95</td><td>0.108</td><td>0.036</td><td>0.108</td><td>0.036</td></tr><tr><td>−2.90</td><td>0.120</td><td>0.035</td><td>0.123</td><td>0.034</td></tr><tr><td>−2.85</td><td>0.142</td><td>0.031</td><td>0.147</td><td>0.029</td></tr><tr><td>−2.80</td><td>0.162</td><td>0.045</td><td>0.170</td><td>0.040</td></tr><tr><td>−2.75</td><td>0.187</td><td>0.055</td><td>0.199</td><td>0.046</td></tr><tr><td>−2.70</td><td>0.206</td><td>0.059</td><td>0.223</td><td>0.040</td></tr><tr><td>−2.65</td><td>0.217</td><td>0.049</td><td>0.233</td><td>0.029</td></tr><tr><td>−2.60</td><td>0.274</td><td>0.046</td><td>0.289</td><td>0.029</td></tr><tr><td>−2.55</td><td>0.333</td><td>0.030</td><td>0.334</td><td>0.030</td></tr><tr><td>−2.50</td><td>0.340</td><td>0.022</td><td>0.347</td><td>0.013</td></tr><tr><td>−2.45</td><td>0.342</td><td>0.068</td><td>0.365</td><td>0.041</td></tr><tr><td>−2.40</td><td>0.415</td><td>0.069</td><td>0.437</td><td>0.042</td></tr><tr><td>−2.35</td><td>0.509</td><td>0.054</td><td>0.514</td><td>0.052</td></tr><tr><td>−2.30</td><td>0.537</td><td>0.053</td><td>0.551</td><td>0.042</td></tr><tr><td>−2.25</td><td>0.577</td><td>0.048</td><td>0.592</td><td>0.031</td></tr><tr><td>−2.20</td><td>0.673</td><td>0.054</td><td>0.675</td><td>0.053</td></tr><tr><td>−2.15</td><td>0.761</td><td>0.083</td><td>0.761</td><td>0.083</td></tr><tr><td>−2.10</td><td>0.799</td><td>0.060</td><td>0.800</td><td>0.060</td></tr><tr><td>−2.05</td><td>0.835</td><td>0.057</td><td>0.852</td><td>0.039</td></tr><tr><td>−2.00</td><td>0.906</td><td>0.025</td><td>0.914</td><td>0.016</td></tr><tr><td>−1.95</td><td>1.000</td><td>0.000</td><td>1.000</td><td>0.000</td></tr></tbody></table>

A more intuitive view of the bias function showing its dependency on the single variables $\rho$ and $|z|$ is presented in Fig. 17 (upper panels), while we refer the reader to Fig. 9 (Sec. 5.3) for the metallicity dependence. The 1D bias functions shown in Fig. 17 were obtained starting from the spectroscopic and photometric samples from which we derived the radial and height stellar distributions displayed in the same figure (lower panels). The upper left panel of Fig. 17 shows how the bias function varies with
<figure><img src="image_16.png" /><figcaption>**Fig. 16.** Bias maps in the cylindrical coordinate plane $(\rho, |z|)$. Each panel shows the results for different stellar populations (see labels).</figcaption></figure>

the radius $\rho$. We notice that Bias($\rho$) is maximum at the lowest radii, and it rapidly declines towards larger $\rho$, reaching its minimum at $\rho \approx 8.5$ kpc, and then it rises again. This behaviour can be easily understood: the position of the minimum corresponds to the position of the Sun, $\rho = 8.43$ kpc (see Sect. 5.2), and hence it identifies the radius where the number of stars in the spectroscopic sample is larger, $N^* \approx 10^4$, as can be seen in the lower left panel of Fig. 17. In the same panel, we see that both at smaller and larger radii the number of stars in the spectroscopic sample decreases, reaching the value $N^* \le 400$ at the extremes. Conversely, the number of stars in the photometric sample remains roughly at the constant value of $N^* > 10^6$ for 6.5 kpc $< \rho \le 10$ kpc and slowly declines at larger $\rho$, following the same trend of the spectroscopic sample. This behaviour drives the shape of the bias function, which we remind the reader is obtained by dividing the photometric and the spectroscopic metallicity distributions (see Sect. 5.2 for more details).

The variation of the bias function with $|z|$ is shown in the right upper panel of Fig. 17. We see that Bias($|z|$) is maximum on the Galactic plane and then it rapidly decreases at higher $|z|$, reaching values close to zero for $|z| \approx 4$ kpc. We can also notice a little bump at $|z| \approx 5$ kpc. Once again, these trends can be understood by looking at the effective number of stars from the spectroscopic and photometric samples in each $|z|$ bin (lower

right panel of Fig. 17). We see that the number of stars in the spectroscopic sample gradually increases with $|z|$, reaching the maximum $N^* \approx 10^4$, at $|z| \approx 1$ kpc, and then rapidly declining at higher $z$. Conversely, the number of stars in the photometric sample is maximum close to the Galactic plane, $|z| \le 0.5$ kpc, and then it starts to decline. This explains the trend of the bias function with $|z|$.

We note that at $|z| > 4$ kpc the spectroscopic sample is very poorly represented, $N^* < 100$, which means that the bump observed in Bias($z$) at $|z| \approx 5$ kpc is driven by a very small number of stars, $N^* < 10$. To understand the origin of this bump we should inspect the complete 3D bias maps shown in Fig. 16 again. With the exception of metal-poor stars (left lower panel), all stellar populations show zero bias at $|z| \ge 4$ kpc for all $\rho$. This means that there are no stars in the spectroscopic sample in these regions and at these metallicities. At $|z| \approx 5$ kpc ($\rho \approx 7.5$ kpc) instead, very metal-poor stars show a small isolated ‘island’ with Bias $\approx 1.5$. It is likely this small region, counting < 10 stars, which originates the small bump, observed Bias($z$) at $|z| = 5$ (right upper panel).

In Fig. 18, we show the spatial distribution of extremely metal-poor, very metal-poor, metal-poor, and metal-rich stars corrected for the bias. In other words, these maps were obtained by multiplying the raw spatial distributions of the spec-
<figure><img src="image_17.png" /><figcaption>**Fig. 17.** Variation of the bias function with the cylindrical radius, $\rho$ (left upper panel), and the height above the Galactic plane, $|z|$ (right upper panel). The lower panels show the distribution of stars from the photometric (blue) and spectroscopic (red) samples as functions of $\rho$ (left) and $|z|$ (right).</figcaption></figure>

troscopic sample for the local 3D bias function. In Fig. 18, we can first notice that the spatial distribution of extremely metal-poor stars (EMPs, left upper panel) is completely different with respect to those of the more metal-rich stellar populations (i.e. [Fe/H] > −3). EMPs are clustered in a smaller region of the ($\rho, |z|$) plane, and their spatial distribution does not appear uniform like the others do. This is most likely a consequence of the low number of statistics: in the three main regions filled by EMPs, we can indeed see that the number of stars per pixel is < 15. Conversely, more metal-rich stars cover almost the same regions of the ($\rho, |z|$), and they are all much better represented in each pixel, where $N^* > 30$ most of the time. However, we can notice that the maxima of the stellar distributions (red regions) reside in three different locations, which implies that there are spatial gradients in the metallicity distribution function. A more in-depth and quantitative analysis of these results, which includes a comparison with cosmological models, is beyond the scope of the present paper and will be presented in our future work (Salvadori et al., in prep).

# 8. Dynamical properties of the SDSS spectroscopic sample

In order to discuss the dynamical properties of the sample, we only make use of the good parallax sample. In terms of Galactic location, the two samples occupy approximately the same space in the $R_{GC}$-$z$ plane, as shown in Fig. 7. We corrected the Gaia parallaxes for the global zero point of −0.029 mas (Lindegren et al. 2018). With these parallaxes, the Gaia proper motions, and the SDSS radial velocities, we used the galpy code<fnref n="6" /> (Bovy 2015) to evaluate orbital parameters and actions under the Staeckel approximation (Binney 2012; Mackereth & Bovy 2018) and using the MWPotential2014 potential (Bovy 2015).

## 8.1. The Gaia-Sausage-Enceladus structure

From the studies of Belokurov et al. (2018), Haywood et al. (2018), Helmi et al. (2018) and Di Matteo et al. (2019), we know that the Galactic halo<fnref n="7" /> is dominated by stars that were accreted from a galaxy that we refer to as Gaia-Sausage-Enceladus (GSE). The precise dynamical definition of this object is still not

<footnote n="6">⁶ http://github.com/jobovy/galpy</footnote>

<footnote n="7">⁷ We loosely define the Galactic halo as the stars with a total space velocity higher than 180km s⁻¹ (Di Matteo et al. 2019).</footnote>
<figure><img src="image_18.png" /><figcaption>**Fig. 18.** Distribution of extremely metal-poor (upper left), very meta-poor (upper right), metal-poor (lower left), and metal-rich (lower right) stars in the cylindrical coordinate plane $(\rho, |z|)$. The distributions have been corrected using the 3D bias function.</figcaption></figure>

known. We refrain from identifying, in dynamical spaces, multiple accretion events (Naidu et al. 2020) and small sub-structures, since we are aware that, from the theoretical point of view, a single massive accretion event can give rise to several dynamical sub-structures (Jean-Baptiste et al. 2017). We therefore concentrate on isolating GSE, although we are aware that it is practically impossible to identify a 'pure' sample of GSE stars, since one of the outcomes of the accretion is that stars that come from both the accreted satellite galaxy and the accreting galaxy will end up in the same region of any dynamical space. Chemical information may be helpful to disentangle the accreted from the accreting, but our information is limited to an average metallicity and is not sufficient for this purpose.

Based on dynamical data and photometric metallicities, Feuillet et al. (2020) suggested that the best way to select a pure GSE sample is in the space of angular momentum $L_z$ and radial action ($J_r$, or its square root). In Fig. 19, we show such a diagram for our good parallax sample. The points in blue correspond to the GSE selection of Feuillet et al. (2020), namely $-500 \text{ kpc km s}^{-1} < L_z < 500 \text{ kpc km s}^{-1}$ and $30 (\text{kpc km s}^{-1})^{1/2} < \sqrt{J_r} < 50 (\text{kpc km s}^{-1})^{1/2}$.

The existence of two sequences of $[\alpha/\text{Fe}]$ in the solar neighbourhood as a function of metallicity was initially discovered by

<figure><img src="image_19.png" /><figcaption>**Fig. 19.** Diagram of angular momenta and square root of the radial action. Our good parallax sample is shown via red dots, the blue rectangle is the GSE selection proposed by Feuillet et al. (2020), the yellow dots are the high-$\alpha$ stars of Nissen & Schuster (2010), while the black dots are their low-$\alpha$ stars.</figcaption></figure>

Fuhrmann (1998, and subsequent papers) and later confirmed by Nissen & Schuster (2010), who introduced the terms low-$\alpha$ and
<figure><img src="image_20.png" /><figcaption>**Fig. 20.** Observed metallicity histogram of GSE, for the whole sample (black), and the bias-corrected histogram (red). The histogram of Feuillet et al. (2020) is shown in blue, while that of Naidu et al. (2020) is shown in green.</figcaption></figure>

high-$\alpha$ stars. Haywood et al. (2018) demonstrated convincingly that the low-$\alpha$ stars should be associated with GSE, while high-$\alpha$ are stars formed in the Milky Way. This was also later confirmed by Helmi et al. (2018). In Fig. 19, we also plot the stars classified as low $\alpha$ and high $\alpha$ by Nissen & Schuster (2010). Seventeen of the Nissen & Schuster (2010) stars fall in the Feuillet et al. (2020) GSE selection box; of these, fourteen are low-$\alpha$ stars and three are high-$\alpha$ stars. This allows us to estimate that the Feuillet et al. (2020) GSE selection box includes about 18% of stars that are not associated with GSE. This confirms the theoretical expectation that it is impossible to dynamically select a pure sample of the accreted stars. In the case of the Nissen & Schuster (2010) stars, the chemical information allows one to disentangle the two populations; however, the two $\alpha$/Fe sequences merge at metallicities lower than about -1.0, making this chemical diagnostic unusable. Pragmatically, we decided to adopt the Feuillet et al. (2020) selection, keeping in mind that there will be a contamination, since this appears to be the least contaminated selection of GSE. We would also like to point out that there is also no guarantee that this selection of GSE provides an unbiased sample of the stars of the accreted galaxy.

The Feuillet et al. (2020) selection contains 3227 stars in our good parallax sample. Since this sample was extracted from a biased sample of stars, this MDF cannot be taken as representative of the MDF of GSE. On the other hand, we cannot apply the bias correction that has been derived for the global sample, since this sub-sample has, a priori, a different bias. We cannot apply the Feuillet et al. (2020) selection to the photometric sample to derive an appropriate bias function, since we lack the radial velocities necessary to derive $L_z$ and $J_r$. However, in Sect. 7 we derive the space dependence of the bias function. We can therefore divide the GSE sample into bins in the $\rho, |z|$ space and apply the appropriate bias function to each bin.

The bias correction results in shifting the MDF to higher metallicities, as expected; however, our MDF remains distinctly more metal-poor than that derived by Feuillet et al. (2020), as shown in Fig. 20. Both our raw and bias-corrected MDF of GSE lack the significant population between -1.0 and solar metallicity that is present in the Feuillet et al. (2020) MDF of GSE.

Our MDF is also slightly narrower than that of Feuillet et al. (2020). Although the shape of none of the MDFs is Gaussian, a Gaussian fit can provide a measure of their width. Our MDFs (both raw and bias-corrected) have a full width at half maximum of about 0.6 dex, while that of Feuillet et al. (2020) is about 0.7 dex.

A possible reason for the discrepancy is that the Casagrande et al. (2019) photometric calibration, which is used by Feuillet et al. (2020) to derive the metallicities, is validated only down to -2.0. Thus, the Feuillet et al. (2020) sample may have been selected with a bias against metal-poor stars. To gain some insight into a comparison of our metallicities and those obtained from the Casagrande et al. (2019) calibration, we cross-matched our sample with Skymapper DR 1.1 (Wolf et al. 2018) and with 2MASS (Skrutskie et al. 2006). We found a sub-sample of 20203 unique stars, and we applied the Casagrande et al. (2019, Eq. 12) calibration. This photometric metallicity correlates only weakly with our spectroscopy metallicity. The metallicity distributions of the same stars using the two different metallicity indicators do not resemble each other. The photometric metallicity is on average 0.7 dex more metal rich, but with a dispersion that is of the same order of magnitude. We do not want to push this comparison further, since our stars are dwarfs, while those used by Feuillet et al. (2020) are giants, and the systematics may well be different.

Naidu et al. (2020) also provide an MDF for GSE, and this is shown in Fig. 20. Their MDF is similar to that of Feuillet et al. (2020) around the peak, but it departs from it both in the metal-weak and in the metal-rich tail. It is clearly different from ours, that has many more metal-poor stars. We point out that the Naidu et al. (2020) selection of GSE is very different from that of Feuillet et al. (2020) and adopted by us. In fact, they write the following: “We select GSE stars by excluding previously defined structures and requiring $e > 0.7$”. Hence, it is necessary to have previously isolated all the different sub-structures identified by Naidu et al. (2020). In order to mimic the Naidu et al. (2020) selection, we selected stars with $e > 0.7$ and $-500 \text{ kpc km s}^{-1} < L_z < 500 \text{ kpc km s}^{-1}$ in our sample. This translates, in our preferred selection plane, to $18 \text{ (kpc km s}^{-1})^{1/2} < \sqrt{J_r} < 60 \text{ (kpc km s}^{-1})^{1/2}$. As stated above, this increases the contamination of stars not related to GSE, probably above the level of 50%, as estimated from applying this selection to the Nissen & Schuster (2010) sample. It is, however, not clear if these contaminating stars have an MDF different from GSE or not. We applied this selection to our sample, and the MDF we obtain is similar to what we obtain with the Feuillet et al. (2020) selection. As in the case of Feuillet et al. (2020), we may suspect that the metallicity scale of Naidu et al. (2020) is different from ours. This seems unlikely since the metal-weak tail of our bias-corrected MDF agrees very well with that of Naidu et al. (2020). A more detailed comparison will have to wait until the H3 catalogue is fully published and metallicities of known stars can be compared.

We can, however, comment on the similarity of the peak position between Feuillet et al. (2020) and Naidu et al. (2020), in spite of their disagreement in both metal-weak and metal-rich tails. Both Feuillet et al. (2020) and Naidu et al. (2020) rely on giant stars. However, neither takes into account the fact that the rate at which the stars leave the main sequence is metallicity dependent (the evolutionary flux; see Renzini & Buzzoni 1986) or the fact that the rate of evolution along the red giant branch is also metallicity dependent (Rood 1972). These effects were discussed in the case of the Galactic bulge by Zoccali et al. (2003). It may well be that the similarity in the peak of the MDFs of GSE of the two above-mentioned studies does in fact signal that both
<figure><img src="image_21.png" /><figcaption>**Fig. 21.** Observed metallicity histogram of GSE, normalised to -1.95 compared to that of the global sample. The global sample is bias corrected as described in Sect. 5.3. The GSE sample is shown both raw and corrected, as described in Sect. 7.</figcaption></figure>

these MDFs are biased by the choice of tracers, namely giant stars.<fnref n="8" />

It is remarkable that the slope of the metal-weak tail of the MDF of GSE is much steeper than what was obtained from the full sample, as shown in Fig. 21. This is true even if we consider our observed MDF, before bias correction; thus, it must be considered as a very robust finding. It could be argued that rather than comparing the GSE MDF to our mean MDF, which also relies on our 1D-corrected MDF, we should have used our 3D-corrected MDF for the full sample. There are two reasons for not doing this. First of all, in the metallicity range that is relevant for GSE, the 1D-corrected and 3D-corrected MDFs are pretty close, while they become increasingly different at lower metallicities where the 3D bias is plagued by a low number of stars per bin. Secondly, the average MDF has a much smaller error than our 1D-corrected MDF and, a fortiori, than our 3D-corrected MDF.

## 8.2. Errors on dynamical quantities and consequences on selections

In order to evaluate the uncertainty in the derived kinematical parameters, we used the *pyia* code<fnref n="9" /> to generate samples of the six parameters (coordinates, parallax, proper motions, and radial velocity) in input to *galpy* for each object. Samples are generated from a multi-variate Gaussian distribution, where the errors on the parallax, proper motions, and radial velocity; and the correlation coefficients between parallax and proper motions are only taken into account in the construction of the covariance matrix. One thousand realisations were extracted for each

object. The samples were then fed into *galpy*, and the errors were evaluated as the standard deviations of the obtained orbital parameters, after unbound orbits were rejected. For the error estimate, we only accepted objects for which more than 950 out of 1000 realisations provided meaningful results, namely 54 215 stars over 57 386. When we consider the GSE sub-sample, the number of stars retained drops from 3 227 to 2 487. We consider that stars that have too few meaningful realisations to have uncertain orbital parameters.

The question is, what is the effect of these errors on a sub-sample selected on the basis of dynamical quantities? We compared the metallicity histograms for the full GSE sample with errors (2 487 stars) and for a sub-sample selected with the supplementary condition $eL_z < 263.983 \text{ kpc km s}^{-1}$ and $eJ_r < 995.473 \text{ kpc km s}^{-1}$, where $eL_z$ and $eJ_r$ are the errors on $L_z$ and $J_r$ derived as described above, and the condition is fixed to the median of the error distribution. From this comparison, it is clear that the shape of the MDF is by and large the same if we consider the full GSE sample or only the sub-sample that has smaller errors on the orbital parameters. We thus conclude that the shape of the MDF for our GSE selection is robustly determined. In the same way, we find that the conclusions are not different if we consider the full sample of 3 227 or the sub-sample of 2 487 stars with reliable orbital parameters, according to the above definition.

## 9. Discussion

Our recommended metal-weak tail of the MDF is provided in Table 3. We deduced that the raw Hamburg/ESO MDF is essentially unbiased from its similarity with our own bias-corrected MDF and the H3 MDF. Why would the HES be unbiased, when Schörck et al. (2009) warned that it is biased in favour of the metal-poor stars and provide bias corrections? Our understanding is that globally the HES is biased; however, the follow-up on metal-poor candidates selected from the objective prism spectra has been unbiased for all stars of metallicity below -2.0. Quite likely, this was because of the low resolution of the prism spectra, which made it difficult to discriminate between a -2.0 and a -3.0 star. For this reason, the HES MDF, normalised at -1.95, is effectively unbiased.

Another question is, why does the Pristine MDF show a larger fraction of metal-poor stars at any metallicity compared to the other three MDFs? This needs to be further investigated using the Pristine data. However, we think this may be due to two factors: Youakim et al. (2020) did not use Gaia-based distance estimates to 'clean' their sample from white dwarfs, but relied on colour cuts. Our own experience, discussed in Sect. 4, shows how the sample appeared much more clear, manageable, and consistent with isochrones, once we applied a distance cut. Youakim et al. (2020) combined the Pristine filter with SDSS wide-band photometry to derive their metallicities. In Fig. 2, we show how the SDSS photometry can be discrepant with Gaia photometry. It could be that the excess in metal-poor stars is partially due to a fraction of stars with poor SDSS photometry.

Another possible issue to be considered is the lower sensitivity to metallicity of any photometric metallicity based on a narrow-band filter centred on CaII lines, for warm TO stars, with respect to cool giants. This effect is clearly illustrated in Fig. 11 of Bonifacio et al. (2019). Finally, there is a degeneracy between metallicity and surface gravity that is not accounted for in Youakim et al. (2020) and that could cause part of the discrepancy. It would certainly be useful to re-derive the Pristine MDF

<footnote n="8">⁸ After the acceptance of the paper, C. Conroy brought to our attention the fact that Naidu et al. (2000) accounted for the evolutionary flux via their selection function. However, the MDFs shown in their figures, from which we have taken the values shown in our figures, are all 'raw'. C. Conroy (priv. comm.) evaluated the metallicity dependence of an MDF based on giants, using the mock Gaia catalogue of Rybizki et al. (2018) and estimates the effect to be of the order of about 0.05 dex, which is far too small to account for the difference between our MDF of GSE and that of Naidu et al. (2020).</footnote>

<footnote n="9">⁹ https://github.com/adrn/pyia (Price-Whelan 2018)</footnote>
after cleaning the sample for stars of too-high absolute magni-
tude and uncertain wide-band photometry.

The MDF we derived for the GSE sub-structure is clearly at
odds with that of Feuillet et al. (2020) and Naidu et al. (2020).
We already pointed out how there could be an offset in the metal-
licity scales. We would like to point out that our sample and
the Feuillet et al. (2020) sample do not overlap perfectly in the
sampled galactocentric volume. The Feuillet et al. (2020) sample
preferentially contains stars with |z| < 1 kpc, while our sample
is composed mainly of stars with |z| > 1 kpc. Also, the Naidu et
al. (2020) sample covers a very different region in the R<sub>GC</sub> |Z|,
plane. In particular, their sample pushes to 50 kpc, obviously
thanks to the use of giant stars as tracers, while we only used
TO stars. If there were any gradient of the metallicity of GSE
with z or R<sub>GC</sub>, one would expect the three samples to provide
different MDFs.

10. Conclusions

We used our own analysis of SDSS spectra in order to deter-
mine the MDF of the metal-poor component of the Milky Way.
The metal-weak part of our MDF is in excellent agreement with
that determined by Allende Prieto et al. (2014) from an indepen-
dent analysis of SDSS spectra. This fact, together with the ex-
ternal comparison of a small sub-set of stars analysed at higher
resolution and/or S/N gives us confidence in our spectroscopic
metallicity, for which we estimate a precision of 0.1 dex and an
accuracy of 0.35 dex. We showed that our sample is biased, in
the sense that metal-poor stars are over-represented.

Through a comparison with a sample of stars, that we as-
sume to be unbiased with respect to metallicity, for which we de-
termined photometric metallicites, we derived the bias function.
This allowed us to provide the unbiased MDF that is given in
Table 2. When considering the metal-weak tail of the MDF, nor-
malised at -1.95, we find that our bias-corrected MDF is essen-
tially identical to that of Naidu et al. (2020) and that of Schörck
et al. (2009). We therefore advocate that the best estimate is ob-
tained by averaging these three MDFs. This average MDF and
its standard deviation is provided in Table 3, and we recommend
using this table for comparisons with theoretical predictions.

In the process of estimating the bias in our spectroscopic
sample, we introduced a new photometric metallicity estimate
based on SDSS *u* − *g* and *g* − *z* colours and on distance esti-
mates based on Gaia parallaxes. These estimates, along with ef-
fective temperatures and surface gravities, for 24 037 008 unique
TO stars is available through CDS<fnref n="10" />.

Thanks to the Gaia distances and proper motions, combined
with SDSS radial velocities, we computed Galactic orbits and
actions for a sub-set of our spectroscopic sample, characterised
by good parallaxes. From this sample, we selected likely mem-
bers of the GSE structure. We were able to derive the metallicity
bias of this GSE sample, looking at the position of these stars in
the galaxy. In this way, we derived an MDF for GSE, as shown
in Figs. 20 and 21. The most notable result is that the fraction
of metal-poor stars in GSE is clearly lower than that observed
in the Milky Way. This implies a very different formation and
evolution history of GSE with respect to the Milky Way. Fur-
thermore, although the 'progenitor' of GSE might contain more
metal-poor stars (see e.g. Jean-Baptiste et al. 2017) than what
we find in the GSE structure, our result seems to suggest that the
metal-poor tail of the Galactic halo is not made by massive dwarf
galaxies like the progenitor of GSE, but most likely by smaller

(ultra-faint) dwarfs. Indeed, according to ΛCDM models, ultra-
faint dwarfs are expected to be the basic building blocks of more
massive systems and the formation environment of extremely
metal-poor stars (see e.g. Salvadori et al. 2015, and references
therein).

Acknowledgements. This paper is dedicated to the late Roger Cayrel who has inspired the TOPOS survey as well as all of our work on metal-poor stars. He passed away on January 11th 2021, and he will be deeply missed. We are grateful to Bertrand Plez who computed all the metal-poor model atmospheres and associated synthetic spectra used in this paper. We acknowledge useful discussions with David Aguado on the contents of this paper. We gratefully acknowledge support from the French National Research Agency (ANR) funded projects "Pristine" (ANR-18-CE31-0017) and "MOD4Gaia" (ANR-15-CE31-0007). A.J.K.H., H.-G.L. and N.C. gratefully acknowledge funding by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project-ID 138713538 - SFB 881 ("The Milky Way System"), subprojects A03, A04, A05, and A11. S.S. acknowledges support from the ERC Starting Grant NEFERTITI H2020/808240 and the PRIN-MIUR2017, The quest for the first stars, prot. n. 2017T4ARJ5. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermi-Lab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is http://www.sdss3.org/. SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofísica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

References

Adelman-McCarthy, J. K., Agüeros, M. A., Allam, S. S., et al. 2008, ApJS, 175, 297-313

Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21

Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12

Allende Prieto, C., Fernández-Alvar, E., Schlesinger, K. J., et al. 2014, A&A, 568, A7

Alvarez R., Plez B., 1998, A&A 330, 1109

Aoki, W., Beers, T. C., Lee, Y. S., et al. 2013, AJ, 145, 13. doi:10.1088/0004-6256/145/1/13

Arenou, F., Luri, X., Babusiaux, C., et al. 2018, A&A, 616, A17. doi:10.1051/0004-6361/201833234

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., et al. 2018, AJ, 156, 58

Beers, T. C., Preston, G. W., & Shectman, S. A. 1985, AJ, 90, 2089. doi:10.1086/113917

Beers, T. C. 1987, Nearly Normal Galaxies. From the Planck Time to the Present, 41

<footnote n="10">¹⁰ cdurl</footnote>
Belokurov, V., Erkal, D., Evans, N. W., et al. 2018, MNRAS, 478, 611

Bessell, M. S. 1977, Proceedings of the Astronomical Society of Australia, 3, 144

Bidelman, W. P. & MacConnell, D. J. 1973, AJ, 78, 687. doi:10.1086/111475

Binney, J. 2012, MNRAS, 426, 1324. doi:10.1111/j.1365-2966.2012.21757.x

Bond, H. E. 1970, ApJS, 22, 117. doi:10.1086/190220

Bond, H. E. 1980, ApJS, 44, 517. doi:10.1086/190703

Bond, H. E. 1981, ApJ, 248, 606. doi:10.1086/159186

Bonifacio, P., Monai, S., & Beers, T. C. 2000, AJ, 120, 2065. doi:10.1086/301566

Bonifacio, P., Caffau, E., François, P., et al. 2011, Astronomische Nachrichten, 332, 251. doi:10.1002/asna.201111528

Bonifacio, P., Sbordone, L., Caffau, E., et al. 2012, A&A, 542, A87

Bonifacio, P., Caffau, E., Spite, M., et al. 2015, A&A, 579, A28

Bonifacio, P., Caffau, E., Spite, M., et al. 2018, A&A, 612, A65

Bonifacio, P., Caffau, E., Sestito, F., et al. 2019, MNRAS, 487, 3797. doi:10.1093/mnras/stz1378

Bovy, J. 2015, ApJS, 216, 29. doi:10.1088/0067-0049/216/2/29

Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127

Caffau, E., Bonifacio, P., François, P., et al. 2011, A&A, 534, A4

Caffau, E., Bonifacio, P., François, P., et al. 2012, A&A, 542, A51

Caffau, E., Bonifacio, P., François, P., et al. 2013a, A&A, 560, A15. doi:10.1051/0004-6361/201322213

Caffau, E., Bonifacio, P., Sbordone, L., et al. 2013b, A&A, 560, A71

Caffau, E., Bonifacio, P., Spite, M., et al. 2016, A&A, 595, L6

Caffau, E., Gallagher, A. J., Bonifacio, P., et al. 2018, A&A, 614, A68. doi:10.1051/0004-6361/201732475

Caffau, E., Bonifacio, P., Sbordone, L., et al. 2020, MNRAS, 493, 4677. doi:10.1093/mnras/staa589

Cargile, P. A., Conroy, C., Johnson, B. D., et al. 2020, ApJ, 900, 28. doi:10.3847/1538-4357/aba43b

Carollo, D. & Chiba, M. 2020, arXiv:2010.00235

Carney, B. W., Laird, J. B., Latham, D. W., et al. 1987, AJ, 94, 1066. doi:10.1086/114542

Casagrande, L., Wolf, C., Mackey, A. D., et al. 2019, MNRAS, 482, 2770. doi:10.1093/mnras/sty2878

Castelli, F. & Kurucz, R. L. 2003, in IAU Symposium, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, 20P, arXiv:astro-ph/0405087v1

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245. doi:10.1086/167900

Cayrel, R. 1986, A&A, 168, 81

Chen, B., Stoughton, C., Smith, J. A., et al. 2001, ApJ, 553, 184. doi:10.1086/320647

Christlieb, N., Schörck, T., Frebel, A., et al. 2008, A&A, 484, 721. doi:10.1051/0004-6361:20078748

Conroy, C., Naidu, R. P., Zaritsky, D., et al. 2019, ApJ, 887, 237. doi:10.3847/1538-4357/ab5710

Di Matteo, P., Haywood, M., Lehnert, M. D., et al. 2019, A&A, 632, A4

Feuillet, D. K., Feltzing, S., Sahlholdt, C. L., et al. 2020, MNRAS, 497, 109. doi:10.1093/mnras/staa1888

François, P., Caffau, E., Wanajo, S., et al. 2018a, A&A, 619, A10

François, P., Caffau, E., Bonifacio, P., et al. 2018b, A&A, 620, A187

François, P., Wanajo, S., Caffau, E., et al. 2020, A&A, 642, A25. doi:10.1051/0004-6361/202038028

Fuhrmann, K. 1998, A&A, 338, 161

Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1. doi:10.1051/0004-6361/201629272

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1. doi:10.1051/0004-6361/201833051

Gustafsson, B., Edvardsson, B., Eriksson, K., Graae-Jørgensen, U., Nordlund, Å., & Plez, B. 2008, A&A, 486, 951

Haywood, M., Di Matteo, P., Lehnert, M. D., et al. 2018, ApJ, 863, 11

Hartwick, F. D. A. 1976, ApJ, 209, 418. doi:10.1086/154735

Helmi, A., Babusiaux, C., Koppelman, H. H., et al. 2018, Nature, 563, 85

Hernandez, X. & Ferrara, A. 2001, MNRAS, 324, 484. doi:10.1046/j.1365-8711.2001.04346.x

Jean-Baptiste, I., Di Matteo, P., Haywood, M., et al. 2017, A&A, 604, A106. doi:10.1051/0004-6361/201629691

Johnson, H. L. & Morgan, W. W. 1953, ApJ, 117, 313. doi:10.1086/145697

Karlsson, T. 2006, ApJ, 641, L41. doi:10.1086/503604

Laird, J. B., Carney, B. W., & Latham, D. W. 1988, AJ, 95, 1843. doi:10.1086/114782

Lallement, R., Babusiaux, C., Vergely, J. L., et al. 2019, A&A, 625, A135. doi:10.1051/0004-6361/201834695

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2

Ludwig, H.-G., Bonifacio, P., Caffau, E., et al. 2008, Physica Scripta Volume T, 133, 014037. doi:10.1088/0031-8949/2008/T133/014037

Mackereth, J. T. & Bovy, J. 2018, PASP, 130, 114501. doi:10.1088/

Plez, B. 2012, Turbospectrum: Code for spectral synthesis, astrophysics Source Code Library

Naidu, R. P., Conroy, C., Bonaca, A., et al. 2020, ApJ, 901, 48. doi:10.3847/1538-4357/abaef4

Nissen, P. E. & Schuster, W. J. 2010, A&A, 511, L10. doi:10.1051/0004-6361/200913877

Oey, M. S. 2003, MNRAS, 339, 849. doi:10.1046/j.1365-8711.2003.06228.x

Placco, V. M., Frebel, A., Lee, Y. S., et al. 2015, ApJ, 809, 136. doi:10.1088/0004-637X/809/2/136

Prantzos, N. 2003, A&A, 404, 211. doi:10.1051/0004-6361:20030445

Price-Whelan, A. 2018, Zenodo

O'Donnell, J. E. 1994, ApJ, 422, 158. doi:10.1086/173713

Renzini, A. & Buzzoni, A. 1986, Spectral Evolution of Galaxies, 195. doi:10.1007/978-94-009-4598-2_19

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130. doi:10.1088/0004-637X/783/2/130

Roederer, I. U., Preston, G. W., Thompson, I. B., et al. 2014, AJ, 147, 136. doi:10.1088/0004-6256/147/6/136

Rood, R. T. 1972, ApJ, 177, 681. doi:10.1086/151747

Ryan, S. G. & Norris, J. E. 1991a, AJ, 101, 1835. doi:10.1086/115811

Ryan, S. G. & Norris, J. E. 1991b, AJ, 101, 1865. doi:10.1086/115812

Rybizki, J., Demleitner, M., Fouesneau, M., et al. 2018, PASP, 130, 074101. doi:10.1088/1538-3873/aabd70

Salvadori, S., Ferrara, A., Schneider, R., et al. 2010, MNRAS, 401, L5. doi:10.1111/j.1745-3933.2009.00772.x

Salvadori, S., Schneider, R., & Ferrara, A. 2007, MNRAS, 381, 647. doi:10.1111/j.1365-2966.2007.12133.x

Salvadori, S., Ferrara, A., Schneider, R., et al. 2010, MNRAS, 401, L5. doi:10.1111/j.1745-3933.2009.00772.x

Salvadori, S., Skúladóttir, Á., & Tolstoy, E. 2015, MNRAS, 454, 1320. doi:10.1093/mnras/stv1969

Santos, N. C., Lovis, C., Pace, G., et al. 2009, A&A, 493, 309. doi:10.1051/0004-6361:200811093

Sbordone, L., Caffau, E., Bonifacio, P., et al. 2014, A&A, 564, A109

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525. doi:10.1086/305772

Schörck, T., Christlieb, N., Cohen, J. G., et al. 2009, A&A, 507, 817

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163. doi:10.1086/498708

Soubiran, C., Le Campion, J.-F., Brouillet, N., et al. 2016, A&A, 591, A118. doi:10.1051/0004-6361/201628497

Sousa, S. G., Adibekyan, V., Delgado-Mena, E., et al. 2018, A&A, 620, A58. doi:10.1051/0004-6361/201833350

Szentgyorgyi, A., Furesz, G., Cheimets, P., et al. 2011, PASP, 123, 1188. doi:10.1086/662209

Tumlinson, J. 2006, ApJ, 641, 1. doi:10.1086/500383

van den Bergh, S. 1962, AJ, 67, 486. doi:10.1086/108757

Yong, D., Norris, J. E., Bessell, M. S., et al. 2013, ApJ, 762, 26. doi:10.1088/0004-637X/762/1/26

Yong, D., Norris, J. E., Bessell, M. S., et al. 2013, ApJ, 762, 27. doi:10.1088/0004-637X/762/1/27

York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579

Youkim, K., Starkenburg, E., Martin, N. F., et al. 2020, MNRAS, 492, 4986

Wallerstein, G. 1962, ApJS, 6, 407. doi:10.1086/190067

Wolf, C., Onken, C. A., Luvaul, L. C., et al. 2018, PASA, 35, e010. doi:10.1017/pasa.2018.5

Zoccali, M., Renzini, A., Ortolani, S., et al. 2003, A&A, 399, 931. doi:10.1051/0004-6361:20021604
## Appendix A: Details on the selection of the spectroscopic sample

In January 2016, we downloaded all the spectra from SDSS DR12 (SDSS & BOSS) with the following SQL command:

```sql
SELECT
  S.plate, S.MJD, S.fiberid, SXCobjID, S.ra, S.dec, S.type, S.fla
  flags as fl,
  S.psfMag_u as u, S.psfMag_g as g, S.psfMag_r as r, S.psfMag_i as i,
  S.psfMag_z as z,
  S.extinction_u as extu, S.extinction_g as extg, S.extinction_r as extr, S.
  extinction_i as exti,
  S.extinction_z as extz, S.survey as survey, S.programname as prog,
  sp.TEFFADOP as teff, sp.TEFFADOPUNC as teff_err, sp.LOGGADOP as glog
,
  sp.LOGGADOPUNC as glog_err, sp.FEHADOP as feh, sp.FEHADOPUNC as feh_err,
  sp.SpecObjID into mydb.dr12 from specPhotoAll as S JOIN sppParams as sp ON
  S.SpecObjID = sp.SpecObjID
WHERE
  (S.type = 6) AND
  (S.fla flags & dbo.fPhotoFlags('BRIGHT') = 0) AND
  (S.fla flags & dbo.fPhotoFlags('SATURATED') = 0) AND
  (S.fla flags & dbo.fPhotoFlags('BLENDED') = 0) AND
  (S.fla flags & dbo.fPhotoFlags('NOTCHECKED') = 0) AND
  (S.fla flags & dbo.fPhotoFlags('DEBLENDED_AS_MOVING') = 0) AND
  (S.psfMag_g-S.extinction_g-S.psfMag_z+S.extinction_z between 0.18 and
  0.70) AND
  (S.psfMag_u-S.extinction_u-S.psfMag_g+S.extinction_g > 0.7).
```

Apart from the flags used to make sure the object is classified as a star, is not too bright or saturated, and so on, this boils down to the conditions $0.18 \le (g - z)_0 \le 0.70$ and $(u - b)_0 > 0.7$. These correspond to the conditions of the TOPoS Survey, which was selected from SDSS DR9 and processed with abbo. This selection resulted in 312 009 spectra of 266 442 unique stars.

## Appendix B: Comparison samples
<table><caption>**Table B.1.** Metallicities of stars in the TOPoS sample, observed at high resolution, compared to those derived by MyGIsFOS on the SDSS spectra.</caption><thead><tr><th>Star</th><th>RA</th><th>DEC</th><th>Ref.</th><th>T<sub>eff</sub><br>K<br>Lit.</th><th>[Fe/H]<br>dex<br>Lit.</th><th>T<sub>eff</sub><br>K<br>adop.</th><th>[M/H]<br>dex<br>MyGIsFOS</th></tr><tr><th></th><th>2000</th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>SDSS J002113 - 005005</td><td>00:21:14</td><td>-00:50:05</td><td>2</td><td>6546</td><td>-3.15</td><td>6398</td><td>-2.81</td></tr><tr><td>SDSS J002558 - 101509</td><td>00:25:59</td><td>-10:15:09</td><td>4,8</td><td>6408</td><td>-3.08</td><td>6424</td><td>-2.96</td></tr><tr><td>SDSS J002749 + 140418</td><td>00:27:49</td><td>+14:04:18</td><td>2</td><td>6125</td><td>-3.37</td><td>6071</td><td>-3.59</td></tr><tr><td>SDSS J003507 - 005037</td><td>00:35:08</td><td>-00:50:38</td><td>4,8</td><td>5891</td><td>-2.95</td><td>5886</td><td>-2.83</td></tr><tr><td>SDSS J003730 + 245750</td><td>00:37:30</td><td>+24:57:51</td><td>9</td><td>6350</td><td>-3.42</td><td>6350</td><td>-3.11</td></tr><tr><td>SDSS J003954 - 001856</td><td>00:39:55</td><td>-00:18:57</td><td>4,8</td><td>6382</td><td>-2.94</td><td>6396</td><td>-3.05</td></tr><tr><td>SDSS J004029 + 160416</td><td>00:40:29</td><td>+16:04:16</td><td>2</td><td>6422</td><td>-3.27</td><td>6388</td><td>-2.87</td></tr><tr><td>SDSS J012032 - 100106</td><td>01:20:33</td><td>-10:01:07</td><td>4,8</td><td>5804</td><td>-3.50</td><td>5798</td><td>-3.72</td></tr><tr><td>SDSS J012125 - 030943</td><td>01:21:25</td><td>-03:09:44</td><td>4,8</td><td>6127</td><td>-2.91</td><td>6128</td><td>-3.46</td></tr><tr><td>SDSS J012442 - 002806</td><td>01:24:42</td><td>-00:28:07</td><td>4,8</td><td>6273</td><td>-2.79</td><td>6281</td><td>-2.00</td></tr><tr><td>SDSS J014036 + 234458</td><td>01:40:36</td><td>+23:44:58</td><td>5</td><td>5848</td><td>-4.00</td><td>5884</td><td>-3.33</td></tr><tr><td>SDSS J014721 + 021819</td><td>01:47:22</td><td>+02:18:20</td><td>4,8</td><td>5967</td><td>-3.30</td><td>5964</td><td>-3.82</td></tr><tr><td>SDSS J014828 + 150221</td><td>01:48:29</td><td>+15:02:22</td><td>4,8</td><td>6151</td><td>-3.41</td><td>6154</td><td>-1.79</td></tr><tr><td>SDSS J021238 + 013758</td><td>02:12:38</td><td>+01:37:58</td><td>5</td><td>6333</td><td>-3.59</td><td>6300</td><td>-3.17</td></tr><tr><td>SDSS J021554 + 063901</td><td>02:15:54</td><td>+06:39:01</td><td>4,8</td><td>6005</td><td>-2.75</td><td>6002</td><td>-3.41</td></tr><tr><td>SDSS J031745 + 002304</td><td>03:17:46</td><td>+00:23:04</td><td>2</td><td>5786</td><td>-3.46</td><td>5710</td><td>-3.27</td></tr><tr><td>SDSS J035925 - 063416</td><td>03:59:26</td><td>-06:34:16</td><td>4,8</td><td>6281</td><td>-3.17</td><td>6290</td><td>-3.30</td></tr><tr><td>SDSS J040114 - 051259</td><td>04:01:15</td><td>-05:12:59</td><td>4,8</td><td>5500</td><td>-3.62</td><td>5925</td><td>-3.14</td></tr><tr><td>SDSS J074748 + 264543</td><td>07:47:49</td><td>+26:45:44</td><td>4,8</td><td>6434</td><td>-3.36</td><td>6452</td><td>-3.03</td></tr><tr><td>SDSS J075338 + 190855</td><td>07:53:39</td><td>+19:08:56</td><td>4,8</td><td>6439</td><td>-2.45</td><td>6458</td><td>-2.40</td></tr><tr><td>SDSS J080336 + 053430</td><td>08:03:37</td><td>+05:34:31</td><td>4,8</td><td>6360</td><td>-2.94</td><td>6373</td><td>-2.37</td></tr><tr><td>SDSS J082118 + 181931</td><td>08:21:18</td><td>+18:19:32</td><td>2</td><td>6158</td><td>-3.80</td><td>6161</td><td>-3.28</td></tr><tr><td>SDSS J082452 + 163356</td><td>08:25:11</td><td>+16:35:00</td><td>1</td><td>5463</td><td>-3.22</td><td>5946</td><td>-0.69</td></tr><tr><td>SDSS J082521 + 040334</td><td>08:25:21</td><td>+04:03:34</td><td>2</td><td>6340</td><td>-3.46</td><td>6340</td><td>-3.60</td></tr><tr><td>SDSS J085211 + 033945</td><td>08:52:12</td><td>+03:39:45</td><td>1</td><td>6343</td><td>-3.24</td><td>6365</td><td>-2.92</td></tr><tr><td>SDSS J085232 + 112331</td><td>08:52:33</td><td>+11:23:31</td><td>4,8</td><td>6206</td><td>-3.45</td><td>6211</td><td>-3.55</td></tr><tr><td>SDSS J090533 - 020843</td><td>09:05:33</td><td>-02:08:44</td><td>4,8</td><td>5974</td><td>-3.44</td><td>5971</td><td>-2.35</td></tr><tr><td>SDSS J090733 + 024608</td><td>09:07:33</td><td>+02:46:08</td><td>2</td><td>5934</td><td>-3.44</td><td>5916</td><td>-3.19</td></tr><tr><td>SDSS J091753 + 523005</td><td>09:17:53</td><td>+52:30:05</td><td>10</td><td>5858</td><td>-3.56</td><td>5853</td><td>-3.49</td></tr><tr><td>SDSS J091913 + 232738</td><td>09:19:13</td><td>+23:27:38</td><td>4,8</td><td>5490</td><td>-3.25</td><td>5486</td><td>-3.44</td></tr><tr><td>SDSS J092912 + 023817</td><td>09:29:12</td><td>+02:38:17</td><td>5,6</td><td>5894</td><td>-4.97</td><td>5889</td><td>-1.36</td></tr><tr><td>SDSS J102915 + 172928</td><td>10:29:15</td><td>+17:29:28</td><td>3</td><td>5811</td><td>-4.73</td><td>5795</td><td>-4.29</td></tr><tr><td>SDSS J103402 + 070116</td><td>10:34:03</td><td>+07:01:17</td><td>5</td><td>6224</td><td>-4.01</td><td>6229</td><td>-3.69</td></tr><tr><td>SDSS J105002 + 242109</td><td>10:50:02</td><td>+24:21:10</td><td>10</td><td>5682</td><td>-3.90</td><td>5676</td><td>-2.70</td></tr><tr><td>SDSS J112031 - 124638</td><td>11:20:32</td><td>-12:46:39</td><td>4,8</td><td>6502</td><td>-3.27</td><td>6525</td><td>-1.70</td></tr><tr><td>SDSS J112211 - 114809</td><td>11:22:12</td><td>-11:48:09</td><td>4,8</td><td>6157</td><td>-3.00</td><td>6159</td><td>-2.24</td></tr><tr><td>SDSS J113528 + 010848</td><td>11:35:28</td><td>+01:08:49</td><td>2</td><td>6132</td><td>-3.03</td><td>6210</td><td>-2.67</td></tr><tr><td>SDSS J113723 + 255354</td><td>11:37:23</td><td>+25:53:54</td><td>5</td><td>6310</td><td>-2.70</td><td>6186</td><td>-2.34</td></tr><tr><td>SDSS J114424 - 004658</td><td>11:44:25</td><td>-00:46:59</td><td>4,8</td><td>6412</td><td>-2.80</td><td>6428</td><td>-2.21</td></tr><tr><td>SDSS J120441 + 120111</td><td>12:04:41</td><td>+12:01:12</td><td>4,8</td><td>5839</td><td>-3.37</td><td>5834</td><td>-3.55</td></tr><tr><td>SDSS J122935 + 262445</td><td>12:29:36</td><td>+26:24:46</td><td>2</td><td>6452</td><td>-3.29</td><td>6496</td><td>-3.25</td></tr><tr><td>SDSS J123055 + 000546</td><td>12:30:55</td><td>+00:05:47</td><td>4,8</td><td>6223</td><td>-3.24</td><td>6228</td><td>-3.22</td></tr><tr><td>SDSS J123404 + 134411</td><td>12:34:05</td><td>+13:44:11</td><td>4,8</td><td>5659</td><td>-3.59</td><td>5653</td><td>-3.60</td></tr><tr><td>SDSS J124304 - 081230</td><td>12:43:04</td><td>-08:12:31</td><td>10</td><td>5488</td><td>-3.48</td><td>5484</td><td>-3.96</td></tr><tr><td>SDSS J124502 - 073847</td><td>12:45:03</td><td>-07:38:47</td><td>4,8</td><td>6110</td><td>-3.21</td><td>6093</td><td>-2.91</td></tr><tr><td>SDSS J130017 + 263238</td><td>13:00:17</td><td>+26:32:39</td><td>5</td><td>6393</td><td>-3.65</td><td>6467</td><td>-3.11</td></tr><tr><td>SDSS J131249 + 001315</td><td>13:12:50</td><td>+00:13:16</td><td>4,8</td><td>6400</td><td>-2.38</td><td>6416</td><td>-3.05</td></tr><tr><td>SDSS J131456 + 113753</td><td>13:14:57</td><td>+11:37:53</td><td>4,8</td><td>6265</td><td>-3.24</td><td>6273</td><td>-3.39</td></tr><tr><td>SDSS J132508 + 222424</td><td>13:25:09</td><td>+22:24:25</td><td>4,8</td><td>6292</td><td>-2.60</td><td>6301</td><td>-2.77</td></tr><tr><td>SDSS J133718 + 074536</td><td>13:37:19</td><td>+07:45:36</td><td>1</td><td>6377</td><td>-3.48</td><td>6398</td><td>-1.76</td></tr><tr><td>SDSS J134922 + 140736</td><td>13:49:23</td><td>+14:07:37</td><td>5</td><td>6112</td><td>-3.60</td><td>6112</td><td>-2.95</td></tr><tr><td>SDSS J135331 - 032930</td><td>13:53:31</td><td>-03:29:30</td><td>4,8</td><td>6224</td><td>-3.18</td><td>6229</td><td>-3.07</td></tr><tr><td>SDSS J141249 + 013206</td><td>14:12:49</td><td>+01:32:07</td><td>4,8</td><td>5799</td><td>-2.94</td><td>5792</td><td>-3.32</td></tr><tr><td>SDSS J143632 + 091831</td><td>14:36:32</td><td>+09:18:31</td><td>2</td><td>6340</td><td>-3.40</td><td>6445</td><td>-3.73</td></tr><tr><td>SDSS J144640 + 124917</td><td>14:46:41</td><td>+12:49:17</td><td>2</td><td>6189</td><td>-3.16</td><td>6222</td><td>-2.86</td></tr><tr><td>SDSS J153747 + 281404</td><td>15:37:48</td><td>+28:14:05</td><td>4,8</td><td>6272</td><td>-3.39</td><td>6280</td><td>-3.55</td></tr><tr><td>SDSS J154246 + 054426</td><td>15:42:47</td><td>+05:44:26</td><td>2</td><td>6179</td><td>-3.48</td><td>6319</td><td>-3.16</td></tr><tr><td>SDSS J154746 + 242953</td><td>15:47:47</td><td>+24:29:53</td><td>4,8</td><td>5903</td><td>-3.16</td><td>5898</td><td>-3.19</td></tr></tbody></table>
<table><caption>**Table B.1.** continued.</caption><thead><tr><th>Star</th><th>RA</th><th>DEC</th><th>Ref.</th><th>T<sub>eff</sub><br>K<br>Lit.</th><th>[Fe/H]<br>dex<br>Lit.</th><th>T<sub>eff</sub><br>K<br>adop.</th><th>[M/H]<br>dex<br>MyGIsFOS</th></tr></thead><tbody><tr><td>SDSS J155751 + 190306</td><td>15:57:52</td><td>+19:03:06</td><td>4,8</td><td>6176</td><td>-2.94</td><td>6180</td><td>-2.59</td></tr><tr><td>SDSS J200513 - 104503</td><td>20:05:13</td><td>-10:45:03</td><td>4,8</td><td>6289</td><td>-3.41</td><td>6298</td><td>-3.63</td></tr><tr><td>SDSS J215023 + 031928</td><td>21:50:24</td><td>+03:19:28</td><td>4,8</td><td>6026</td><td>-2.84</td><td>6024</td><td>-3.40</td></tr><tr><td>SDSS J215805 + 091417</td><td>21:58:06</td><td>+09:14:17</td><td>4,8</td><td>5942</td><td>-3.41</td><td>5938</td><td>-3.53</td></tr><tr><td>SDSS J220121 + 010055</td><td>22:01:22</td><td>+01:00:56</td><td>4,8</td><td>6392</td><td>-3.03</td><td>6408</td><td>-0.47</td></tr><tr><td>SDSS J220728 + 055658</td><td>22:07:28</td><td>+05:56:59</td><td>4,8</td><td>6096</td><td>-3.26</td><td>6096</td><td>-3.58</td></tr><tr><td>SDSS J222130 + 000617</td><td>22:21:30</td><td>+00:06:17</td><td>4,8</td><td>5891</td><td>-3.14</td><td>5886</td><td>-2.93</td></tr><tr><td>SDSS J223143 - 094834</td><td>22:31:44</td><td>-09:48:34</td><td>2</td><td>6053</td><td>-3.20</td><td>6038</td><td>-3.01</td></tr><tr><td>SDSS J230814 - 085526</td><td>23:08:15</td><td>-08:55:26</td><td>2</td><td>6015</td><td>-3.01</td><td>6011</td><td>-2.32</td></tr><tr><td>SDSS J231031 + 031847</td><td>23:10:32</td><td>+03:18:48</td><td>4,8</td><td>6229</td><td>-2.91</td><td>6235</td><td>-2.47</td></tr><tr><td>SDSS J231755 + 004537</td><td>23:17:56</td><td>+00:45:38</td><td>4,8</td><td>5918</td><td>-3.54</td><td>5914</td><td>-3.41</td></tr><tr><td>SDSS J233113 - 010933</td><td>23:31:13</td><td>-01:09:33</td><td>2</td><td>6246</td><td>-3.08</td><td>6247</td><td>-3.39</td></tr></tbody></table>

**References.** (1) Caffau et al. (2011), (2) Bonifacio et al. (2012), (3) Caffau et al. (2012), (4) Caffau et al. (2013b), (5) Bonifacio et al. (2015), (6) Caffau et al. (2016), (7) Bonifacio et al. (2018), (8) François et al. (2018b), (9) François et al. (2018a), and (10) François et al. (2020).

<table><caption>**Table B.2.** Metallicities of stars in the PASTEL catalogue (Soubiran et al. 2016), compared to those derived by MyGIsFOS on the SDSS spectra.</caption><thead><tr><th>Star</th><th>RA</th><th>DEC</th><th>Ref.</th><th>T<sub>eff</sub><br>K<br>Lit.</th><th>[Fe/H]<br>dex<br>Lit.</th><th>T<sub>eff</sub><br>K<br>adop.</th><th>[M/H]<br>dex<br>MyGIsFOS</th></tr></thead><tbody><tr><td>SDSS J000219.87+292851.7</td><td>00:02:20</td><td>+29:28:52</td><td>1</td><td>6150</td><td>-3.26</td><td>6199</td><td>-2.84</td></tr><tr><td>CS 22957-0036</td><td>00:03:31</td><td>-04:44:23</td><td>2</td><td>5970</td><td>-2.28</td><td>6292</td><td>-2.07</td></tr><tr><td>SDSS J002015.45-004058.1</td><td>00:20:15</td><td>-00:40:58</td><td>1</td><td>6450</td><td>-2.81</td><td>6411</td><td>-2.70</td></tr><tr><td>SDSS J002756.76-190929.8</td><td>00:27:57</td><td>-19:09:30</td><td>1</td><td>6550</td><td>-2.71</td><td>6474</td><td>-2.88</td></tr><tr><td>2MASS J00285742-1015305</td><td>00:28:57</td><td>-10:15:31</td><td>1</td><td>6150</td><td>-2.81</td><td>6056</td><td>-2.98</td></tr><tr><td>SDSS J004150.22-095327.7</td><td>00:41:50</td><td>-09:53:28</td><td>1</td><td>6000</td><td>-2.81</td><td>5968</td><td>-2.71</td></tr><tr><td>SDSS J010026.70+004915.6</td><td>01:00:27</td><td>+00:49:16</td><td>1</td><td>5550</td><td>-3.18</td><td>5469</td><td>-3.08</td></tr><tr><td>SDSS J011150.32+144207.8</td><td>01:11:50</td><td>+14:42:08</td><td>1</td><td>6350</td><td>-2.8</td><td>6136</td><td>-2.70</td></tr><tr><td>2MASS J02091203+2120281</td><td>02:09:12</td><td>+21:20:28</td><td>1</td><td>6250</td><td>-2.89</td><td>6267</td><td>-2.85</td></tr><tr><td>SDSS J025453.33+332840.9</td><td>02:54:53</td><td>+33:28:41</td><td>1</td><td>6200</td><td>-2.82</td><td>6338</td><td>-2.84</td></tr><tr><td>SDSS J030839.27+050534.8</td><td>03:08:39</td><td>+05:05:35</td><td>1</td><td>5950</td><td>-2.19</td><td>5881</td><td>-1.94</td></tr><tr><td>SDSS J035111.26+102643.1</td><td>03:51:11</td><td>+10:26:43</td><td>1</td><td>5450</td><td>-3.18</td><td>5726</td><td>-2.45</td></tr><tr><td>SDSS J072725.15+160949.4</td><td>07:27:25</td><td>+16:09:50</td><td>1</td><td>6600</td><td>-2.92</td><td>6249</td><td>-2.99</td></tr><tr><td>2MASS J07464132+2831426</td><td>07:46:41</td><td>+28:31:43</td><td>1</td><td>6100</td><td>-2.75</td><td>6024</td><td>-2.65</td></tr><tr><td>SDSS J074859.88+175832.6</td><td>07:49:00</td><td>+17:58:33</td><td>1</td><td>6100</td><td>-2.6</td><td>6300</td><td>-2.74</td></tr><tr><td>SDSS J080428.21+515303.1</td><td>08:04:28</td><td>+51:53:03</td><td>1</td><td>5950</td><td>-3.01</td><td>5855</td><td>-2.84</td></tr><tr><td>SDSS J080917.05+090748.5</td><td>08:09:17</td><td>+09:07:49</td><td>1</td><td>6150</td><td>-3.38</td><td>6121</td><td>-3.03</td></tr><tr><td>2MASS J08145867+3337130</td><td>08:14:59</td><td>+33:37:13</td><td>1</td><td>6450</td><td>-3.28</td><td>6547</td><td>-3.23</td></tr><tr><td>2MASS J08175492+2641038</td><td>08:17:55</td><td>+26:41:04</td><td>1</td><td>6050</td><td>-2.85</td><td>6038</td><td>-2.98</td></tr><tr><td>2MASS J08273626+1052008</td><td>08:27:36</td><td>+10:52:01</td><td>1</td><td>6400</td><td>-3.17</td><td>6252</td><td>-3.39</td></tr><tr><td>SDSS J084016.15+540526.4</td><td>08:40:16</td><td>+54:05:26</td><td>1</td><td>6150</td><td>-3.25</td><td>6064</td><td>-2.86</td></tr><tr><td>SDSS J084700.49+012113.6</td><td>08:47:00</td><td>+01:21:14</td><td>1</td><td>6250</td><td>-3.2</td><td>6326</td><td>-3.09</td></tr><tr><td>Cl* NGC 2682 ES I-46</td><td>08:50:53</td><td>+11:43:40</td><td>3</td><td>5703</td><td>-0.07</td><td>5534</td><td>-0.49</td></tr><tr><td>Cl* NGC 2682 YBP 1514</td><td>08:51:01</td><td>+11:53:11</td><td>4</td><td>5076</td><td>0.09</td><td>5520</td><td>-0.41</td></tr><tr><td>SDSS J085136.67+101803.1</td><td>08:51:37</td><td>+10:18:03</td><td>1</td><td>6500</td><td>-2.89</td><td>6290</td><td>-3.11</td></tr><tr><td>SDSS J091243.72+021623.6</td><td>09:12:44</td><td>+02:16:24</td><td>1</td><td>6150</td><td>-2.68</td><td>6097</td><td>-2.60</td></tr><tr><td>SDSS J093247.29+024123.8</td><td>09:32:47</td><td>+02:41:24</td><td>1</td><td>6200</td><td>-3.14</td><td>6202</td><td>-2.66</td></tr><tr><td>SDSS J100427.70+344245.7</td><td>10:04:28</td><td>+34:42:46</td><td>1</td><td>6100</td><td>-3.09</td><td>6091</td><td>-2.97</td></tr><tr><td>2MASS J10330141+4001035</td><td>10:33:01</td><td>+40:01:04</td><td>1</td><td>6600</td><td>-3.06</td><td>6391</td><td>-2.64</td></tr><tr><td>SDSS J1036+1212</td><td>10:36:50</td><td>+12:12:20</td><td>1</td><td>5850</td><td>-3.47</td><td>5833</td><td>-3.08</td></tr><tr><td>SDSS J110610.48+034321.9</td><td>11:06:10</td><td>+03:43:22</td><td>1</td><td>6300</td><td>-2.88</td><td>6203</td><td>-2.82</td></tr><tr><td>2MASS J11082167+1747468</td><td>11:08:22</td><td>+17:47:47</td><td>1</td><td>6050</td><td>-3.17</td><td>5971</td><td>-2.86</td></tr><tr><td>SDSS J112051.73+302724.4</td><td>11:20:52</td><td>+30:27:25</td><td>1</td><td>5750</td><td>-3.14</td><td>5752</td><td>-2.97</td></tr><tr><td>SDSS J112813.56+384148.9</td><td>11:28:14</td><td>+38:41:49</td><td>5</td><td>6449</td><td>-3.53</td><td>6616</td><td>-2.96</td></tr><tr><td>SDSS J114723.53+151044.7</td><td>11:47:24</td><td>+15:10:45</td><td>1</td><td>6500</td><td>-2.96</td><td>6528</td><td>-2.93</td></tr><tr><td>SDSS J115906.18+542512.6</td><td>11:59:06</td><td>+54:25:13</td><td>1</td><td>5700</td><td>-3.26</td><td>5724</td><td>-3.14</td></tr><tr><td>SDSS J120441.38+120111.5</td><td>12:04:41</td><td>+12:01:12</td><td>6</td><td>5467</td><td>-4.34</td><td>5834</td><td>-3.55</td></tr><tr><td>SDSS J121307.22+445040.9</td><td>12:13:07</td><td>+44:50:41</td><td>1</td><td>6650</td><td>-3.2</td><td>6411</td><td>-3.79</td></tr></tbody></table>
<table><caption>**Table B.2.** continued.</caption><thead><tr><th>Star</th><th>RA</th><th>DEC</th><th>Ref.</th><th>T<sub>eff</sub><br>K</th><th>[Fe/H]<br>dex</th><th>T<sub>eff</sub><br>K</th><th>[M/H]<br>dex</th></tr><tr><td></td><td></td><td>2000</td><td></td><td>Lit.</td><td>Lit.</td><td>adop.</td><td>MyGIsFOS</td></tr></thead><tbody><tr><td>SDSS J123055.25+000546.9</td><td>12:30:55</td><td>+00:05:47</td><td>1</td><td>6150</td><td>-3.34</td><td>6228</td><td>-3.22</td></tr><tr><td>2MASS J12330008+3407583</td><td>12:33:00</td><td>+34:07:58</td><td>1</td><td>6300</td><td>-2.87</td><td>6175</td><td>-2.62</td></tr><tr><td>SDSS J124502.68-073847.1</td><td>12:45:03</td><td>-07:38:47</td><td>1</td><td>6100</td><td>-3.17</td><td>6093</td><td>-2.91</td></tr><tr><td>2MASS J13001719+2632385</td><td>13:00:17</td><td>+26:32:39</td><td>1</td><td>6450</td><td>-3.53</td><td>6467</td><td>-3.40</td></tr><tr><td>SDSS J130402.25+323909.1</td><td>13:04:02</td><td>+32:39:09</td><td>1</td><td>6050</td><td>-2.9</td><td>6075</td><td>-2.32</td></tr><tr><td>SDSS J131201.47+245006.9</td><td>13:12:01</td><td>+24:50:07</td><td>1</td><td>6250</td><td>-2.72</td><td>6193</td><td>-2.86</td></tr><tr><td>SDSS J133453.44+002238.6</td><td>13:34:53</td><td>+00:22:38</td><td>1</td><td>5650</td><td>-3.03</td><td>5636</td><td>-3.00</td></tr><tr><td>SDSS J133841.16+120415.2</td><td>13:38:41</td><td>+12:04:15</td><td>1</td><td>6300</td><td>-2.86</td><td>6284</td><td>-2.80</td></tr><tr><td>SDSS J1349-0229</td><td>13:49:14</td><td>-02:29:43</td><td>1</td><td>6200</td><td>-3.24</td><td>6138</td><td>-1.25</td></tr><tr><td>SDSS J140035.31+075317.7</td><td>14:00:35</td><td>+07:53:18</td><td>1</td><td>6250</td><td>-2.98</td><td>6269</td><td>-3.01</td></tr><tr><td>SDSS J140813.88+623942.1</td><td>14:08:14</td><td>+62:39:42</td><td>1</td><td>6300</td><td>-2.97</td><td>6267</td><td>-2.66</td></tr><tr><td>2MASS J14100174+5350180</td><td>14:10:02</td><td>+53:50:18</td><td>1</td><td>6100</td><td>-3.42</td><td>6109</td><td>-2.93</td></tr><tr><td>SDSS J141207.32+560931.9</td><td>14:12:07</td><td>+56:09:32</td><td>1</td><td>6600</td><td>-3.19</td><td>6413</td><td>-3.21</td></tr><tr><td>SDSS J142441.88+561535.0</td><td>14:24:42</td><td>+56:15:35</td><td>1</td><td>6350</td><td>-2.97</td><td>6234</td><td>-2.61</td></tr><tr><td>SDSS J142518.09+113713.9</td><td>14:25:18</td><td>+11:37:14</td><td>1</td><td>6300</td><td>-3.08</td><td>6253</td><td>-3.01</td></tr><tr><td>SDSS J142541.33+574207.5</td><td>14:25:41</td><td>+57:42:08</td><td>1</td><td>6450</td><td>-3.29</td><td>6343</td><td>-3.39</td></tr><tr><td>SDSS J143451.02+103626.4</td><td>14:34:51</td><td>+10:36:26</td><td>1</td><td>6400</td><td>-3.21</td><td>6319</td><td>-3.40</td></tr><tr><td>SDSS J143654.45+030143.2</td><td>14:36:54</td><td>+03:01:43</td><td>1</td><td>6150</td><td>-3.6</td><td>6180</td><td>-3.10</td></tr><tr><td>SDSS J143708.91+523146.7</td><td>14:37:09</td><td>+52:31:47</td><td>1</td><td>6300</td><td>-2.9</td><td>6227</td><td>-2.81</td></tr><tr><td>SDSS J143759.06+583723.6</td><td>14:37:59</td><td>+58:37:24</td><td>1</td><td>6450</td><td>-3.02</td><td>6384</td><td>-2.98</td></tr><tr><td>SDSS J150217.16+311316.5</td><td>15:02:17</td><td>+31:13:16</td><td>1</td><td>6350</td><td>-2.86</td><td>6181</td><td>-2.89</td></tr><tr><td>SDSS J151534.44+450317.7</td><td>15:15:34</td><td>+45:03:18</td><td>1</td><td>6300</td><td>-3.34</td><td>6227</td><td>-3.06</td></tr><tr><td>SDSS J151646.69+433331.6</td><td>15:16:47</td><td>+43:33:32</td><td>1</td><td>6500</td><td>-2.91</td><td>6403</td><td>-2.98</td></tr><tr><td>2MASS J15215863+3437292</td><td>15:21:59</td><td>+34:37:29</td><td>1</td><td>6200</td><td>-2.79</td><td>6127</td><td>-2.97</td></tr><tr><td>SDSS J152301.85+494210.6</td><td>15:23:02</td><td>+49:42:11</td><td>1</td><td>5850</td><td>-3.06</td><td>5764</td><td>-2.83</td></tr><tr><td>SDSS J152810.51+491526.8</td><td>15:28:11</td><td>+49:15:27</td><td>1</td><td>6450</td><td>-2.99</td><td>6467</td><td>-2.84</td></tr><tr><td>2MASS J15531082+2511404</td><td>15:53:11</td><td>+25:11:40</td><td>1</td><td>5850</td><td>-3.3</td><td>5852</td><td>-2.88</td></tr><tr><td>SDSS J160303.73+291709.5</td><td>16:03:04</td><td>+29:17:10</td><td>1</td><td>6000</td><td>-3.36</td><td>6071</td><td>-2.94</td></tr><tr><td>SDSS J162311.84+391319.6</td><td>16:23:12</td><td>+39:13:20</td><td>1</td><td>6350</td><td>-3.19</td><td>6253</td><td>-2.78</td></tr><tr><td>SDSS J162603.61+145844.3</td><td>16:26:04</td><td>+14:58:44</td><td>1</td><td>6400</td><td>-2.99</td><td>6327</td><td>-2.18</td></tr><tr><td>SDSS J163331.43+390742.6</td><td>16:33:31</td><td>+39:07:43</td><td>1</td><td>6300</td><td>-2.88</td><td>6237</td><td>-3.02</td></tr><tr><td>SDSS J164005.30+370907.8</td><td>16:40:05</td><td>+37:09:08</td><td>1</td><td>6450</td><td>-3.39</td><td>6382</td><td>-3.24</td></tr><tr><td>SDSS J164610.19+282422.2</td><td>16:46:10</td><td>+28:24:22</td><td>1</td><td>6100</td><td>-3.05</td><td>6163</td><td>-2.28</td></tr><tr><td>SDSS J165016.66+224213.9</td><td>16:50:17</td><td>+22:42:14</td><td>1</td><td>6600</td><td>-2.56</td><td>6381</td><td>-3.03</td></tr><tr><td>2MASS J16593474+3515545</td><td>16:59:35</td><td>+35:15:54</td><td>1</td><td>6050</td><td>-3.24</td><td>6182</td><td>-2.74</td></tr><tr><td>SDSS J172846.88+065701.9</td><td>17:28:47</td><td>+06:57:02</td><td>1</td><td>6350</td><td>-2.85</td><td>6578</td><td>-3.31</td></tr><tr><td>SDSS J183045.75+414126.8</td><td>18:30:46</td><td>+41:41:27</td><td>1</td><td>6550</td><td>-3.01</td><td>6546</td><td>-3.77</td></tr><tr><td>SDSS J183414.28+202335.5</td><td>18:34:14</td><td>+20:23:35</td><td>1</td><td>6550</td><td>-1.98</td><td>6353</td><td>-2.03</td></tr><tr><td>2MASS J20051348-1045030</td><td>20:05:13</td><td>-10:45:03</td><td>1</td><td>6600</td><td>-3.34</td><td>6298</td><td>-3.63</td></tr><tr><td>SDSS J205252.68+010939.3</td><td>20:52:53</td><td>+01:09:39</td><td>1</td><td>6050</td><td>-2.8</td><td>6032</td><td>-2.68</td></tr><tr><td>SDSS J211850.12-064055.8</td><td>21:18:50</td><td>-06:40:56</td><td>1</td><td>6650</td><td>-2.91</td><td>6436</td><td>-2.73</td></tr><tr><td>2MASS J21231081-0820393</td><td>21:23:11</td><td>-08:20:39</td><td>1</td><td>6350</td><td>-2.88</td><td>6172</td><td>-3.01</td></tr><tr><td>SDSS J212841.25-075629.3</td><td>21:28:41</td><td>-07:56:29</td><td>1</td><td>6150</td><td>-2.76</td><td>6067</td><td>-2.85</td></tr><tr><td>SDSS J220743.35+201752.3</td><td>22:07:43</td><td>+20:17:52</td><td>1</td><td>6200</td><td>-2.42</td><td>6433</td><td>-1.86</td></tr><tr><td>SDSS J220845.57+061341.3</td><td>22:08:46</td><td>+06:13:41</td><td>1</td><td>6450</td><td>-2.85</td><td>6300</td><td>-2.68</td></tr><tr><td>V* FV Aqr</td><td>22:12:03</td><td>-08:45:46</td><td>2</td><td>6200</td><td>-2.58</td><td>6568</td><td>-2.21</td></tr><tr><td>CS 22886-0012</td><td>22:13:26</td><td>-08:43:43</td><td>2</td><td>5650</td><td>-2.88</td><td>6056</td><td>-2.55</td></tr><tr><td>SDSS J230026.35+055956.2</td><td>23:00:26</td><td>+05:59:56</td><td>1</td><td>6450</td><td>-2.94</td><td>6343</td><td>-2.93</td></tr><tr><td>SDSS J233403.22+153829.3</td><td>23:34:03</td><td>+15:38:29</td><td>1</td><td>6550</td><td>-2.91</td><td>6196</td><td>-2.81</td></tr><tr><td>SDSS J234939.71+383217.8</td><td>23:49:40</td><td>+38:32:18</td><td>1</td><td>6250</td><td>-3.2</td><td>6284</td><td>-2.67</td></tr></tbody></table>

**References.** (1) Aoki et al. (2013), (2) Roederer et al. (2014), (3) Santos et al. (2009), (4) Sousa et al. (2018), (5) Yong et al. (2013a), (6) Placco et al. (2015)

**Appendix C: Metallicity estimates from the *p* index**

The procedure is as follows:

1. For each star, the effective temperature is estimated from the $(g-z)_0$ colour using formula given in footnote 3, and the *p* index is computed.
2. The values of theoretical bolometric corrections in the G band and *p* indices are read for a set of eight effective temperatures and six surface gravities, and nine metallicities are read.

3. For each star, the parallax (obtained as $1/r_{est}$, where $r_{est}$ is the geometrical distance estimate from Bailer-Jones et al. 2018), the effective temperature, and the G0 magnitude is read.

4. For an initial estimate of log g = 4.2, a first estimate of metallicity is derived from the *p* index, as well as a first estimate of the bolometric correction, by a spline interpolation in the theoretical table.

5. A new gravity estimate is derived using the current bolometric correction estimate.

6. A new estimate of metallicity and bolometric correction is derived with the current surface gravity.

7. Steps 5 and 6 are iterated until the log g difference from the previous step is less than 0.05 dex.