Free Access
Issue
A&A
Volume 654, October 2021
Article Number A107
Number of page(s) 14
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202040074
Published online 19 October 2021

© ESO 2021

1. Introduction

Blue horizontal-branch stars are old objects found on the blue side of RR Lyra at the bottom of the instability strip. They represent a late stage in the evolution of low-mass stars with initial masses in the range of ∼0.8 M to ∼2.3 M. A star will leave the main-sequence once core hydrogen burning ceases and the inert helium core contracts. This contraction heats the hydrogen shell around the core, thereby allowing fusion to begin, and the star evolves onto the red giant branch (RGB), where a large amount of the outer envelope is shed by the strong solar wind combined with the tenuous gravitational binding forces at large radii. The inert helium core will continue to contract and becomes electron-degenerate. Finally the temperature in the helium core rises enough to allow core helium burning to begin and the helium flash occurs. This increases the temperature further and the core ceases to be degenerate. At this point, stars move to the blue horizontal-branch with masses in the range of ∼0.5−1.0 M (Montenegro et al. 2019) and display helium burning in their core as well as hydrogen shell burning (Moehler 2001; Ruhland et al. 2011; Paunzen et al. 2019). The mass of the hydrogen shell may vary from 0.02 M to over 0.2 M. The mass of the helium-burning core is ∼0.5 M. Blue horizontal-branch stars with a higher mass of hydrogen shell will appear cooler and, hence, redder (Moehler 2001; Paunzen et al. 2019).

The existence of horizontal-branch stars was first proven in photometric studies of the M 92 and M 3 globular clusters in the 1950s (Arp et al. 1952; Sandage 1953). These observations explained the differences between the colour-magnitude diagrams from open and globular clusters first observed by Harlow Shapley in his study of the structure of the Galaxy between 1914 and 1919. Catalogues of blue horizontal-branch stars were first published in the 1980s (Pier 1982; Beers et al. 1988) and have continued to develop as larger and higher-quality datasets have become available (Beers et al. 1996, 2007; Montenegro et al. 2019).

The small scatter in their absolute magnitude and their relative abundance makes blue horizontal-branch stars a particularly useful standard candle in determining distances (Clewley & Jarvis 2006). This, and the fact that they are metal-poor population II stars makes blue horizontal-branch stars ideal for studying Galactic structure in older parts of the Galaxy (Starkenburg et al. 2019; Newberg et al. 2003; Monaco et al. 2003; Clewley & Jarvis 2006; Yanny et al. 2009; Niederste-Ostholt et al. 2010) or for studying the enclosed masses of the Milky Way within different radii (Sommer Larsen & Christensen 1986; Sommer-Larsen et al. 1989; Norris & Hawkins 1991; Xue et al. 2008).

Observationally, most blue horizontal-branch stars are expected to have an A-type (Xue et al. 2008) or very late B-type spectra. The A-type spectra seen in blue horizontal-branch stars can be distinguished from population I main-sequence A-type stars by their stronger Balmer jump, stronger and deeper Balmer lines, and the lack of metal spectral lines (Smith et al. 2010). Those with an effective temperature of Teff ≥ 10 000 K have spectra that are similar to the much younger population I B-type main-sequence stars but with weaker Helium lines, which can be distinguished spectroscopically (see Fig. 1). The main challenge in producing a catalogue of low contamination is differentiating between blue horizontal-branch stars and younger, population I main-sequence A-type and B-type stars as well as blue stragglers, which are older, population II main-sequence stars that have been rejuvenated through acquisition of additional mass (Rain et al. 2021) and that display A-type spectra. The additional mass may be gained by Roche Lobe overflow in a close binary system (McCrea 1964), through stellar merger (Hills & Day 1976; Davies et al. 1994), or both. Blue stragglers, which have similar effective temperatures but higher gravities than blue horizontal-branch stars, are also fainter (Starkenburg et al. 2019) and are expected to be found in older parts of the Galaxy (Santucci et al. 2015), where younger main-sequence stars of similar effective temperatures no longer exist.

thumbnail Fig. 1.

Typical B-type spectra displaying the difference between main-sequence and blue horizontal-branch stars. Upper panel: example of a typical SDSS B-type spectrum from a blue- horizontal-branch star (SDSS J131916.15–011404.9 from Geier et al. 2015). Lower panel: example of a typical LAMOST B-type spectrum from a main-sequence star. Some reference lines for hydrogen and helium are indicated by vertical lines. SDSS and LAMOST spectra are comparable due to their similar resolutions and spectral coverage.

The goal of this work is to use the second data release (DR2) and early third data release (EDR3) of Gaia to produce a catalogue of blue horizontal-branch candidate stars with a low and quantifiable blue straggler and population I main-sequence contamination. This method identifies likely blue horizontal-branch stars without relying on spectra, thus providing important input information for upcoming large spectroscopic surveys. The blue horizontal-branch star catalogues of Xue et al. (2008) and Behr (2003) were used as reference catalogues to calibrate the colour and absolute magnitude selection criteria in the Gaia DR2 colour magnitude space.

2. Selection of known blue horizontal-branch stars

The blue horizontal-branch dataset of Xue et al. (2008) was considered the reference for determining the Gaia DR2 colour and absolute magnitude selection criteria. Xue et al. (2008) used the seventh data release of the Sloan Digital Sky Survey (SDSS DR7, Abazajian et al. 2009) data to obtain 2558 blue horizontal-branch stars claiming contamination of less than 10%. The contaminants were considered to be main-sequence stars, particularly blue stragglers. This reference dataset was created using several steps (detailed in Xue et al. 2008, and references therein). Firstly, by applying colour cuts in the (u − g)×(g − r) parameter space, then by considering the Balmer line widths at 20% below the local continuum plotted against the depth of the Balmer line compared to the local continuum, and finally by using the ‘scale width versus shape’ method (Clewley et al. 2002) where the Balmer lines are fit to a Sérsic profile (Sérsic 1968).

The Behr (2003) dataset was used together with the Xue et al. (2008) dataset for the purposes of comparing reference datasets with the Gaia DR2 selections outlined in subsequent sections. This was done to confirm that the selection criteria developed were valid. The Behr (2003) dataset comprises 74 blue horizontal-branch stars selected from 6 metal-poor globular clusters in the Globular Cluster Catalogue (Harris 1996).

3. Constructing the blue horizontal-branch star catalogue from Gaia DR2

The European Space Agency’s Gaia mission was initiated with the aim of charting a three-dimensional (3D) map of the Milky Way by surveying over one billion stars in our Galaxy. This is about 1% of the stars in our Galaxy. The procedure involved creating a database of the position, parallax, brightness, colour, and proper motions of all objects visible to the satellite (i.e. down to an apparent GaiaG magnitude of ∼20.7 mag, Gaia Collaboration 2016). Each Gaia data release provides a dataset that is more complete and with reduced uncertainties than its predecessor.

Assessments of the accuracy of the astrometric and photometric measurements in the second data release (Gaia Collaboration 2018a) were published by Lindegren et al. (2018) and Evans et al. (2018), respectively. Assessments of the accuracy of the astrometry and photometry in the early third data release (Gaia Collaboration 2021) were published by Lindegren et al. (2020) and Riello et al. (2020), respectively. We used the TOPCAT application (Taylor 2005) to cross-match the Xue et al. (2008) blue horizontal-branch stars with the Gaia DR2 dataset to provide an initial calibration point for finding such stars in the Gaia DR2 parameter space.

The methodology used was based on the procedure used by Gentile Fusillo et al. (2018) to create a white dwarf catalogue, as well as by Pelisoli & Vos (2019) to create a catalogue of extremely low-mass white dwarf candidates, and by Geier et al. (2019) to create their catalogue of hot subluminous stars. Some deviation from these methodologies was required to address the challenge of main-sequence contamination as the blue horizontal-branch crosses the main-sequence leading to similar parameter values in the Gaia DR2 and EDR3 colour magnitude parameter space.

3.1. Gaia DR2 colour and absolute magnitude criteria selection

The Xue et al. (2008) reference dataset was cross-matched with Gaia DR2. The absolute G magnitude was then calculated as:

As the absolute G magnitude was calculated using the Gaia DR2 parallax measurement for distance, it was necessary to restrict the reference dataset to only those stars with a reliable parallax measurement. Using the parallax quality criterion of parallax error < 20% (parallax_over_error > 5), as in Geier et al. (2019) and Pelisoli & Vos (2019), left only 39 of the original 2560 objects from the Xue et al. (2008) reference dataset, as shown in Fig. 2. These 39 stars were found to lie in the ranges of:

thumbnail Fig. 2.

Result of cross-matching the Xue et al. (2008) reference dataset with Gaia DR2. Upper panel: Xue et al. (2008) objects plotted as turquoise dots on the Gaia DR2 colour-magnitude diagram exhibiting a large scatter in absolute magnitude due to incorrect parallax measurements. Those with parallax error < 20% are plotted as red large dots. Lower panel: Lindegren et al. (2018) ‘clean subset selection C’ (black) compared with Xue et al. (2008) objects with parallax error < 20% (red) on the same scales as the upper plot in the Gaia DR2 colour-magnitude diagram.

in the Gaia DR2 colour absolute G magnitude parameter space.

These colour and absolute magnitude cut-off criteria (henceforth, referred to as the initial CMD criteria) were applied to the Gaia DR2 dataset. We also applied quality criteria that were based on the astrometry and photometry quality criteria of the ‘clean subset Selection C’ by Lindegren et al. (2018). We modified these criteria to be less stringent by allowing parallax errors up to 20% (as opposed to 10%) and by relaxing the parallax lower limit to the requirement of parallax being non-negative (as opposed to being greater than 10 mas). Using the unmodified Lindegren et al. (2018) criteria would have limited the selection to distances below ∼100 pc, where predominantly main-sequence Galactic disc objects are found. Furthermore, we did not apply the astrometric_excess_noise < 1 mas criterion as we wanted to keep the blue horizontal-branch candidates that may exist in binary systems. As can be seen in Belokurov et al. (2020), in particular, Fig. 7, there is an over-density of sources with a re-normalised unit weight error (RUWE) > 1 in the blue horizontal-branch region of the Gaia DR2 colour-magnitude diagram, indicating that many blue horizontal-branch candidates may exist in unresolved binary systems. Lindegren et al. (2018) show that the astrometric excess noise to be a proxy for RUWE.

We retained the Lindegren et al. (2018) quality criteria for the GBP and GRP photometry, which require that the errors on the GBP and GRP flux measurements are < 10% and that the flux excess factor (phot_bp_rp_excess_factor), namely:

lies in the range defined by the following cut-offs:

Applying the initial CMD criteria and the modified Lindegren et al. (2018) criteria for astrometric and photometric quality to Gaia DR2 resulted in a set of over 9 million objects. These stars are expected, based on their position on the colour-magnitude diagram (CMD), to be blue horizontal-branch stars, blue-stragglers, and population I main-sequence A-type and B-type stars. The measured parallaxes of these objects lead us to believe that they exist in either the Galactic disc or the Galactic halo.

It is known that the relatively young population I main-sequence A-type and B-type stars are not expected in the much older Galactic halo. We expected that focusing on the Galactic halo will reduce the population I main-sequence contamination of the dataset. To do this we applied the stringent halo selection criteria (VT > 200 km s−1, |b|> 60°) to the 9 million objects found so far. Using a tangential velocity cut-off to differentiate between disk and halo objects is an accepted method (Gaia Collaboration 2018b; Koppelman & Helmi 2020) that uses only the two tangential components of velocity. It should be noted that any halo objects with low tangential velocities but with high radial velocities along the line of sight will not be selected. This is more likely to be the case in the direction of the motion of the heliocentric system within the Galactic reference frame. We calculated the tangential velocity as:

This resulted in a set of 43 040 objects that, when plotted on the Gaia DR2 CMD that showed a clear clustering in the region where blue horizontal-branch stars would be expected. This was confirmed by overlaying the reference datasets from Xue et al. (2008) and Behr (2003) (see Fig. 3).

thumbnail Fig. 3.

Gaia DR2 colour absolute magnitude diagram showing the blue horizontal-branch star region. Top panel: Gaia DR2 objects that satisfy the modified (Lindegren et al. 2018) criteria (gray points) with those with the stringent halo criteria applied (small blue dots). Xue et al. (2008) objects (red circles) and Behr (2003) objects (yellow triangles) are superimposed. The green trapezium denotes the Gaia DR2 colour-magnitude coordinate space (final CMD criteria) considered as the blue horizontal-branch region. The main-sequence turn-off can be seen at around (GBP − GRP) = 0.6 and Gabs = 3. Blue stragglers are found extending the main-sequence beyond that. Our selection was designed to minimise the blue-straggler contamination and focus on the blue horizontal-branch stars. Bottom panel: for comparison, the Lindegren et al. (2018) ‘clean subset selection C’ comprised of disc stars within 100 pc which are primarily main-sequence stars (black) with red trapezium denoting the Gaia DR2 colour-magnitude coordinate space considered as the blue horizontal-branch region.

We were then able to define modified CMD selection criteria that would allow us to focus more closely on the blue horizontal-branch region in Gaia DR2 colour-magnitude space. We defined these modified criteria as:

henceforth the final Gaia CMD criteria. Applying the final CMD selection criteria to the 9 million initial CMD objects reduced the selection to 305 622 objects comprising of blue horizontal-branch stars, blue-stragglers, and population I main-sequence A-type and B-type stars. It is known that regions of high main-sequence stellar population density and hence, high potential population I main-sequence contamination exist in the Galactic plane and the Magellanic Clouds. These high apparent stellar density regions are also prone to having astrometric solutions that are not accurate, but still precise, and therefore not easily removed by a quality assessment alone (e.g. Gentile Fusillo et al. 2018). Objects in these regions were removed by excluding regions in the Galactic plane and the Magellanic Clouds with a local apparent stellar population density > 50 000 per deg2:

the Large Magellanic Cloud:

and the Small Magellanic Cloud:

This further reduced the number of objects in our selection to 104 763. Fainter objects in crowded regions are known to be susceptible to photometric errors in the GBP and GRP bands (Lindegren et al. 2018). Thus, additional filtering was added whereby stars with apparent neighbours within 5 arcsec were rejected (Pelisoli & Vos 2019; Gentile Fusillo et al. 2018; Boubert et al. 2018). This close apparent neighbour filtering criterion had, together with the Galactic plane and Magellanic Cloud filtering, the additional benefit of removing the majority of objects where reddening from extinction in dust clouds is not negligible. We did not apply extinction corrections to the Gaia DR2 colours for these objects due to the low resolution of existing dust maps (Pelisoli & Vos 2019). Removing any remaining objects with close apparent neighbours further reduced the number of selected objects to 91 719 which are expected to comprise of blue horizontal-branch stars, blue stragglers, and a reduced proportion of population I main-sequence A-type and B-type stars (see Fig. 4).

thumbnail Fig. 4.

Sky distribution of parallax selection candidates. Regions of high stellar number density (> 50 000 stars per square degree) have been excluded.

3.1.1. Galactic halo selection

The stringent halo selection criteria were relaxed in small iterative increments in order to find the least restrictive cut-offs before a disproportionate increase in the number of candidate objects was observed indicating population I main-sequence and non-halo object contamination. A cut-off for Galactic latitude was not applied and a value of tangential velocity of VT > 145 km s−1 was found for the parallax candidate objects. This was considered the Galactic halo selection criterion for the parallax error < 20% dataset (see Table 2) and was applied to the parallax candidates to produce a cleaner set of halo, blue horizontal-branch stars and, thus, a cleaner selection of parallax candidates.

3.1.2. Verification of the Galactic halo selection using LAMOST spectra

The 91 719 objects found prior to applying any halo selection criteria were cross-matched with the SDSS data release 16 (DR16, Ahumada et al. 2020) and LAMOST fifth data release (DR5; Yao et al. 2019) surveys to find spectra for the blue horizontal-branch candidate objects. The relatively close proximity of the parallax candidates leads to their being brighter than more distant candidates. The low apparent magnitudes (lower than the SDSS DR16 photometric saturation limit) of the objects being cross-matched with the SDSS DR16 spectra meant that less than 100 matches were found which we considered to be not statistically significant and thus, may not be representative of the dataset. LAMOST DR5 spectra were, however, found for ∼24 000 stars. This represented roughly one third of the objects found had a similar magnitude distribution and were thus considered to be a statistically significant and representative population.

The LAMOST stellar parameter pipeline (Xiang et al. 2015) found that ∼21 500 stars had A-type spectra, as would be expected for most blue horizontal-branch stars, for blue stragglers, and for population I main-sequence A-type stars (Xue et al. 2008). Furthermore, about 880 were found to have B-type spectra, as would be expected for hotter blue horizontal-branch stars (see the upper panel of Fig. 1 for a typical example) and population I main-sequence B-type stars (see the lower panel of Fig. 1 for a typical example). The remaining ∼1600 spectra were considered to be too noisy for a spectral classification. Visual inspection of these spectra confirmed the LAMOST classification for the A-type and B-type spectra and less than 1% of the LAMOST objects were manually re-categorised.

The spectra of objects categorised by the LAMOST stellar parameter pipeline as being B-type stars were found to have distinct helium lines as would be expected for population I main-sequence B stars. Blue horizontal-branch B-type stars display much subtler helium lines. Thus, we used the presence of large numbers of LAMOST B-type spectra as a proxy for the presence of population I main-sequence stars and, hence, as a proxy for non-halo objects.

The application of the halo selection criterion of VT > 145 km s−1 reduced the number candidate objects by around 90% while reducing the number of LAMOST cross-matched objects with B-type spectra from ∼880 to 6, a reduction of over 99%. This represents an effective filter with which we can discern between disc and halo objects, thereby facilitating the removal of a large potential source of population I main-sequence contamination.

3.2. Gaia reduced proper motion criteria selection

Selecting only the Gaia DR2 stars with parallax errors of < 20% results in stars that are more than a few hundred pc away being rejected. New selection criteria that do not use the Gaia parallax measurement needed to be found in order to find candidate stars in this large population of more distant objects. The parallax was used in calculating the absolute magnitude and the tangential velocity in the parallax_error < 20% dataset.

Reduced proper motion (H) is calculated from total proper motion (μ) and can be used as a proxy for absolute magnitude (Gentile Fusillo et al. 2015; Geier et al. 2019):

Figure 5 shows the distribution of reduced proper motions seen in the objects from the parallax selection. This was used to define the selection criteria for the blue horizontal-branch candidates from the parallax error ≥ 20% objects (henceforth, called the proper motion candidates).

thumbnail Fig. 5.

Gaia DR2 reduced proper motion versus colour diagram and reduced proper motion frequency plot used to define the selection criteria for the objects with a parallax error ≥ 20%. Left panel: reduced proper motion (H) versus colour (GBP − GRP) diagram for the blue horizontal-branch candidates with parallax error < 20%. Right panel: distribution of reduced proper motion (H) seen in the blue horizontal-branch candidates with a parallax error < 20%.

The total_proper_motion_over_error was also calculated in order to provide a quality criterion to prevent stars with large proper motion uncertainties being used (Gentile Fusillo et al. 2015). This is analogous to using parallax_over_error as a quality criterion in the parallax error < 20% dataset. A quality cut-off of total proper motion error < 20% was used in line with the parallax error < 20% cut-off used previously.

The plot of Gaia DR2 colour versus reduced proper motion for the parallax candidates dataset (Fig. 5) gave a distribution of reduced proper motions with , σH = 0.81. A range of ±3σ about the mean reduced proper motion gave a 99% certainty of capturing all blue horizontal-branch candidates using reduced proper motion as a proxy for absolute magnitude.

The same photometric quality criteria as used in selecting the parallax candidate dataset were also applied for the proper motion candidate dataset. The astrometric quality criteria used for the parallax candidate dataset were removed but the criterion parallax_over_error ≤ 5 was applied to ensure that no stars were to be found in both datasets.

Application of the total proper motion over error cut-off, the reduced proper motion and colour criteria, the photometric quality criteria, crowded region criterion, the Galactic latitude criterion, and the parallax error cut-off to the Gaia DR2 dataset resulted in the proper motion selection dataset.

Figure 6 shows the Gaia DR2 objects that conform to the total proper motion over error cut-off, the photometric quality and the crowded region criteria, the subset that conforms to the blue horizontal-branch and halo criteria with the reference objects (Xue et al. 2008) overlain.

thumbnail Fig. 6.

Gaia DR2 colour and reduced proper motion plot: objects conforming to the crowded region criterion as well as the total proper motion and photometric quality criteria (grey dots), objects also conforming to the blue horizontal-branch and halo criteria (blue crosses), and Xue et al. (2008) objects (red squares).

3.2.1. Galactic halo selection

As the parallax is used in calculating the tangential velocity this could no longer be reliably used as the halo selection criterion. We used the Galactic latitude for this purpose. As done previously, we applied the stringent halo selection criterion (absolute Galactic latitude |b|> 60°) and we relaxed this criterion in small iterative increments in order to find the least stringent cut-offs before a disproportionate increase in the number of candidate objects was observed indicating population I main-sequence and non-halo object contamination. The halo selection criterion found was |b|> 50°.

3.2.2. Verification of the Galactic halo selection using LAMOST spectra

As previously done for the parallax selection, the halo selection criterion was verified by comparing the number of LAMOST B-spectra removed to the total number of objects removed. Application of the halo selection criterion of |b|> 50° to the ∼24 000 LAMOST cross-matched stars, presented in Sect. 3.1.2, reduced the number candidate objects by around 92%, while reducing the number of LAMOST cross-matched objects with B-type spectra from ∼880 to 2, a reduction of over 99.8%. This represents an effective filter for discerning between disc and halo objects, thereby facilitating the removal of a large potential source of population I main-sequence contamination.

4. Application of the Gaia GR2 selection criteria to Gaia EDR3

Upon the release of Gaia EDR3, the processing steps outlined in the sections above were performed using Gaia EDR3 and the results were compared to those achieved with Gaia DR2. Firstly, we did this to verify that the methodology and criteria developed for the Gaia DR2 dataset was still applicable in the latest release. Secondly, this was done to create a catalogue of blue horizontal-branch stars using Gaia EDR3.

The 2558 blue horizontal-branch reference objects from Xue et al. (2008) were found in Gaia EDR3 absolute magnitude-colour space. The improved accuracy of the parallax measurements compared to Gaia DR2 meant that 121 objects were found to have a parallax error of less than 20% in Gaia EDR3 (compared to 39 in Gaia DR2) see Fig. 7. The position of these reference objects in the Gaia EDR3 colour and absolute magnitude space verified that using the same colour and absolute magnitude criteria that were found using Gaia DR2 was valid.

thumbnail Fig. 7.

Xue et al. (2008) blue horizontal-branch reference objects (red dots) displayed in DR2 (upper panel) and EDR3 (lower panel). Those objects with a parallax error of less than 20% are shown as blue squares.

Thus, the same selection criteria were applied to Gaia EDR3 and the same stringent halo selection criteria were applied (see Fig. 8). We inspected the results and found that the final DR2 CMD criteria as defined in Sect. 3.1

thumbnail Fig. 8.

Top panel: objects that satisfy the modified (Lindegren et al. 2018) criteria (gray points) with those with the stringent halo criteria and the crowded region criterion applied (small light blue dots). Xue et al. (2008) objects (red circles) and Behr (2003) objects (yellow triangles) are superimposed. The green trapezium denotes the Gaia DR2 colour-magnitude coordinate space (final CMD criteria) considered as the (blue horizontal-branch region) displayed in DR2 (upper panel) and EDR3 (lower panel).

were also applicable to EDR3.

The observed reduced proper motion distribution from Gaia EDR3 parallax selection was then used in defining the proper motion candidate selection criteria. We compared the reduced proper motion distribution of these objects from the Gaia DR2 and Gaia EDR3 datasets and found a close agreement despite the larger number of objects found in Gaia EDR3 (see Fig. 9).

thumbnail Fig. 9.

Distribution of reduced proper motions for the parallax selection in Gaia DR2 (red) and Gaia EDR3 (blue).

5. Comparison of Gaia DR2 and EDR3 results

We found that 7307 of the 8556 objects in the Gaia DR2 parallax selection were also present in the Gaia EDR3 parallax selection. The 1249 parallax candidate objects that are in the DR2 dataset but no longer in the EDR3 dataset are found on the Gaia CMD with increasing concentration towards the −1 < Gabs < (2.5 − 2.77(GBP − GRP)) cut-off criterion. This criterion removes millions of objects, some of which will have measured values that deviate from their true value by multiple standard deviations causing an increase in contamination towards this highly populated region in the Gaia colour-magnitude space.

By comparing the distance (calculated as the reciprocal of parallax) distributions of the parallax selections from Gaia DR2 and Gaia EDR3, we can demonstrate that we have been able to identify more distant candidate objects associated with higher apparent magnitudes (see Fig. 10).

thumbnail Fig. 10.

Distance distribution of parallax selection objects in Gaia DR2 (red) and Gaia EDR3 (blue).

Making such a comparison for objects with high uncertainties on their parallaxes may be misleading. As blue horizontal-branch stars are considered to have nearly constant absolute magnitudes we can consider apparent magnitude as a proxy for distance. The comparison of Gaia DR2 and Gaia EDR3 catalogues can be meaningfully made in both the parallax selection and the proper motions selection in Fig. 11.

thumbnail Fig. 11.

Distribution of apparent magnitudes for the parallax selection (upper panel) for Gaia DR2 (red) and Gaia EDR3 (blue). Distribution of apparent magnitudes for proper motion selection (lower panel) in Gaia DR2 (red) and Gaia EDR3 (blue).

We compared the proper motions candidate objects found in Gaia DR2 and the parallax candidate objects found in Gaia EDR3. We observed that 1737 objects that were found in the Gaia DR2 proper motions catalogue were now to be found in the Gaia EDR3 parallax dataset. Our examination of the apparent magnitudes of these 1737 objects showed them to be at the fainter (average apparent G magnitude 14.38 mag) and, hence, at the more distant end of the dataset. Thus, the increased accuracy of the parallax measurements made in Gaia EDR3 enables the identification of halo blue horizontal-branch objects over the full sky out to a greater distance than was previously possible using these methods in Gaia DR2.

6. Estimation of catalogue contamination and completeness

We made an estimation of the population I main-sequence and blue-straggler contamination levels using the SDSS SEGUE Stellar Parameter Pipeline (SSPP) DR10 (Lee et al. 2008; Smolinski et al. 2011), based on its effective temperature (Teff) and surface gravity (log g) parameters. The method used and the cut-offs applied were taken from Santucci et al. (2015 and references therein), as noted in Sect. 3.2.1.

Cross-matching the parallax candidate objects with the SDSS SEGUE database found 556 stars within the final Gaia CMD criteria (Table 1) and 276 stars within the Galactic halo selection criterion (Table 2). Applying the cut-offs used in Santucci et al. (2015), we found that 70% plotted in the blue horizontal-branch region, 10% plotted in the population I main-sequence region, and 20% plotted in the blue straggler region (Fig. 12).

thumbnail Fig. 12.

Differentiating between blue horizontal-branch stars and blue stragglers. Upper panel: Gaia EDR3 objects that satisfy the modified (Lindegren et al. 2018) and the final CMD criteria (grey dots) with those with the stringent halo criteria and the crowded region criterion applied (blue circles); SDSS SEGUE objects plotting in the blue horizontal-branch region (yellow squares); and SDSS SEGUE objects plotting in the blue straggler region (magenta triangles) in TEFF-LOGG space. Lower panel: SDSS SEGUE TEFF × LOGG plot with blue horizontal-branch (yellow square) and blue straggler stars (magenta triangles).

Table 1.

Selection criteria applied to Gaia DR2 to define blue horizontal-branch star candidate selection from colour and absolute G magnitude (parallax selection).

Table 2.

Comparison of the Gaia DR2 and EDR3 selection criteria and results.

We found that modifying the final Gaia CMD cut-off from:

to the more restrictive

gave a smaller catalogue but with a lower contamination level, with 86% of objects in the blue horizontal-branch region, 8% of objects in the population I main-sequence region, and 6% in the blue-straggler region. This more restrictive selection was not used in the final catalogue, but it demonstrates that a higher level of purity can be achieved at the expense of completeness.

Cross-matching the proper motion candidate objects with the SDSS SEGUE database found 488 stars within the Gaia colour and reduced proper motion criteria (see Sect. 3.2) and 107 stars within the Galactic halo selection criterion (b < −50° or b > 50°). Applying the same cut-offs, we found that 68% plotted in the blue horizontal-branch region, 8% plotted in the population I main-sequence region, and 24% plotted in the blue-straggler region.

The catalogue of blue horizontal-branch candidate objects is thus estimated to have a purity level of ∼70%. The main contaminants are population I main-sequence stars with a contamination level of ∼10%, and a blue-straggler contamination of ∼20%. Using more restrictive CMD selection criteria can reduce these levels of contamination for the candidate objects selected using parallax measurements, but at the cost of reduced completeness of the catalogue.

It should be noted that cross-matching the catalogues of blue horizontal-branch candidates with the SDSS SEGUE database only finds a relatively small number of matches, which are not evenly distributed over the sky. Thus, the calculated levels of purity and contamination can only be considered to be estimates indicating a general level achieved.

The catalogue does not include those blue horizontal-branch stars that do not conform to the various astrometric and photometric quality criteria applied, or those with close apparent neighbours, or those that exist in the direction of crowded regions. The actual number of such objects has not been estimated.

The distance distribution of parallax selection objects peaks at around ∼5.5 kpc and goes up to ∼10 kpc. The proportion of parallax candidates that are found at |b|> 50° is 22%. This apparent overpopulation at higher Galactic latitudes (only ∼8% of a sphere’s volume is at |b|> 50°) indicates that the crowded region and apparent neighbour filters, which dominate at lower Galactic latitudes, have a beneficial effect on the catalogue’s purity at the expense of the catalogue’s completeness. It also implies that neither the catalogue’s completeness nor its purity is uniform in all directions.

There are 30 343 of the 336 911 final Gaia EDR3 CMD criteria objects that conform to the VT > 145 km s−1 criterion. Applying the crowded region criterion and the near apparent neighbour criterion to these 30 343 objects reduces the number to 16 794. This shows that these criteria have removed roughly half of the potential candidates and will thus have a significant effect on the completeness of the catalogue.

A further adverse effect on the completeness is due to the high line-of-sight radial velocity of some halo objects, as noted in Sect. 3.2. Using the halo stellar velocity distribution functions given in Sect. 3.1.3 in Anguiano et al. (2020), we estimate that ∼5% of halo objects will not pass the VT > 145 km s−1 criterion and will, thus, be missing in the catalogue. The number of thin-disk stars with an apparent tangential velocity greater than 145 km s−1 is considered to be negligible. When we make the assumption all blue horizontal-branch stars have an absolute magnitude of 1, then we can calculate an approximate distance to the reduced proper motion selection objects using:

Using this equation, we find that the distance distribution of these objects peaks at ∼12 kpc with 50% of the objects closer than ∼22 kpc but with some candidates found out to ∼60 kpc.

7. Catalogue of blue horizontal-branch objects

The ADQL query for the parallax candidate dataset given in Appendix A returns 30 343 Gaia EDR3 sources, of which 21 751 have no apparent neighbours within 5 arcsec in the Gaia database (Fig. 13). The dataset of parallax candidates did not have a cut-off applied for Galactic latitude, but the criterion that there should be no apparent neighbour within 5 arcsec does, however, remove the many candidate objects in the Galactic plane. This means that our selection places a greater emphasis on having an uncontaminated catalogue rather than on a complete full-sky catalogue. The criterion that the parallax error should be < 20% results in a catalogue of objects that reside within a limited distance from Earth. Using the distance calculated from the parallax, we found that the parallax candidates are found at a minimum distance of 200 pc, a mean distance of 4657 pc, and a maximum of 11 888 pc (Fig. 10 upper). The Gaia EDR3 apparent G magnitude of the parallax selection objects ranged from 6.0 mag to 16.1 mag with a mean value of 14.2 (Fig. 11, upper panel).

thumbnail Fig. 13.

Catalogue of parallax selection objects. Upper panel: Gaia DR2 colour-magnitude diagram showing the blue horizontal-branch parallax candidates (blue) overlain on a Lindegren et al. (2018) ‘clean subset selection C’ (black). Lower panel: sky distribution of blue horizontal-branch parallax candidates.

We compared the distance calculated from parallax to the theoretical distance calculated from the apparent magnitude while assuming that blue horizontal-branch stars have a nearly constant absolute magnitude (Sirko et al. 2003; Deason et al. 2011). For this purpose, we assumed a constant absolute Gaia EDR3 G magnitude of 1 mag (taken from Fig. 3, upper panel). This approximation is not a state-of-the-art approach, but it demonstrates the distance ranges in question and shows that the increasing scatter with increasing parallax error is as expected (Fig. 14).

thumbnail Fig. 14.

Parallax candidates plotted comparing the distance calculated from parallax to the theoretical distance assuming that all blue horizontal-branch objects have the same absolute magnitude (1 mag). The parallax error is plotted on the colour axis.

The ADQL query for the proper motion data set given in Appendix A returns 42 091 Gaia EDR3 sources, of which 40 538 have no apparent neighbours within 5 arcsec in the Gaia database (Fig. 15). We applied a Galactic latitude cut-off to the proper motion candidates. This was done to ensure that contamination was minimised while consciously accepting that this will not be a complete catalogue. The distribution of Gaia EDR3 apparent G magnitudes ranges from a minimum of 11.3 mag to a maximum of 20.1 mag and a mean of 17.6 mag (Fig. 16).

thumbnail Fig. 15.

Catalogue of proper motion selection objects. Upper panel: Gaia DR2 colour-reduced proper motion diagram showing the blue horizontal-branch proper motion candidates (blue) overlain on all Gaia DR2 objects with total proper motion errors of < 1% (gray). Lower panel: sky distribution of blue horizontal-branch proper motion candidates.

thumbnail Fig. 16.

Distribution of Gaia EDR3 apparent G magnitudes of blue horizontal-branch candidate objects selected using parallax (blue) and using proper motions (turquoise).

8. Summary and conclusions

In this work, we compiled a catalogue (see Table 3) comprising of two subsets:

  1. 16 794 blue horizontal-branch parallax candidate stars (with parallax error < 20%);

  2. 40 583 blue horizontal-branch proper motion candidate stars (with parallax error ≥ 20%).

Table 3.

Structure of the online version of the blue horizontal-branch candidate star catalogue.

The purity of the catalogue is estimated to be ∼70% with a population I main-sequence contamination level of less than ∼10% seen in both subsets of the final catalogue and a blue straggler contamination levels of ∼20% in both subsets of the catalogue. These purity and contamination estimates have been made using the SDSS SEGUE results.

We used reference datasets (Xue et al. 2008; Behr 2003) to define the colour and absolute magnitude region of the Gaia DR2 colour-magnitude diagram where blue horizontal-branch stars would be expected. Additional criteria were used based on Lindegren et al. (2018) to ensure that the objects being considered had data that were fit for the purpose. Modifications and additions that were made to the Lindegren et al. (2018) astrometry and photometry quality criteria were based upon similar publications, providing Gaia DR2 based catalogues of extremely low-mass white dwarfs (Pelisoli & Vos 2019), hot sub-dwarf stars (Geier et al. 2019), and white dwarfs (Gentile Fusillo et al. 2018).

The population I main-sequence contamination was reduced by focusing on objects in the Galactic halo. The great age of the Milky Way’s halo is such that comparatively young main-sequence stars are not expected to be present in large numbers.

These steps were subsequently repeated using Gaia EDR3, having first confirmed that the criteria selected using Gaia DR2 were applicable to Gaia EDR3. A comparison of the Gaia DR2 and Gaia EDR3 based catalogues showed that a richer catalogue has been made possible with the improved accuracy of the astrometric measurements.

Our catalogue is a first step in creating a more complete, full-sky catalogue of blue horizontal-branch stars once more reliable distance measurements as well as more reliable and unambiguous photometry data become available. It is known that the Herzsprung-Russel diagram on its own is not sufficient to disentangle the blue horizontal-branch from the main-sequence but this catalogue may provide the selection criteria upon which further surveys can be based. In particular, upcoming large spectroscopic surveys (e.g. 4MOST, WEAVE, DESI) can construct their target lists using catalogues such as ours, which identify objects without relying on existing spectra.

Acknowledgments

IP was partially funded by the Deutsche Forschungsgemeinschaft under grant GE2506/12-1 and by the UK’s Science and Technology Facilities Council (STFC), grant ST/T000406/1. 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. This research has also made extensive use of TOPCAT (http://www.starlink.ac.uk/topcat/, Taylor 2005), and of NASA’s Astrophysics Data System (http://adsabs.harvard.edu/).

References

  1. Ahumada, R., Allende Prieto, C., Almeida, A., et al. 2020, ApJS, 249, 3 [NASA ADS] [CrossRef] [Google Scholar]
  2. Anguiano, B., Majewski, S. R., Hayes, C. R., et al. 2020, AJ, 160, 43 [NASA ADS] [CrossRef] [Google Scholar]
  3. Arp, H. C., Baum, W. A., & Sandage, A. R. 1952, AJ, 57, 4A [CrossRef] [Google Scholar]
  4. Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [Google Scholar]
  5. Beers, T. C., Preston, G. W., & Shectman, S. A. 1988, ApJS, 67, 461 [NASA ADS] [CrossRef] [Google Scholar]
  6. Beers, T. C., Wilhelm, R. J., Doinidis, S. P., & Mattson, C. 1996, ApJS, 103, 433 [NASA ADS] [CrossRef] [Google Scholar]
  7. Beers, T. C., Rossi, S., Wilhelm, R. J., et al. 2007, ApJS, 168, 277 [NASA ADS] [CrossRef] [Google Scholar]
  8. Behr, B. B. 2003, ApJS, 149, 67 [NASA ADS] [CrossRef] [Google Scholar]
  9. Belokurov, V., Penoyre, Z., Oh, S., et al. 2020, MNRAS, 496, 1922 [Google Scholar]
  10. Boubert, D., Strader, J., Aguado, D., et al. 2018, MNRAS, 486, 2618 [Google Scholar]
  11. Clewley, L., & Jarvis, M. J. 2006, MNRAS, 368, 310 [NASA ADS] [Google Scholar]
  12. Clewley, L., Warren, S. J., Hewett, P. C., et al. 2002, MNRAS, 337, 87 [NASA ADS] [CrossRef] [Google Scholar]
  13. Davies, M. B., Benz, W., & Hills, J. G. 1994, ApJ, 424, 870 [Google Scholar]
  14. Deason, A. J., Belokurov, V., & Evans, N. W. 2011, MNRAS, 416, 2903 [NASA ADS] [CrossRef] [Google Scholar]
  15. Evans, D. W., Riello, M., De Angeli, F., et al. 2018, A&A, 616, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Gaia Collaboration (Brown, A. G. A., et al.) 2018a, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Gaia Collaboration (Babusiaux, C., et al.) 2018b, A&A, 616, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 650, C3 [EDP Sciences] [Google Scholar]
  20. Geier, S., Kupfer, T., Heber, U., et al. 2015, A&A, 577, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Geier, S., Raddi, R., Gentile Fusillo, N. P., & Marsh, T. R. 2019, A&A, 621, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Gentile Fusillo, N. P., Gänsicke, B. T., & Greiss, S. 2015, MNRAS, 448, 2274 [Google Scholar]
  23. Gentile Fusillo, N. P., Tremblay, P. E., Gänsicke, B. T., et al. 2018, MNRAS, 482, 4570 [Google Scholar]
  24. Harris, W. E. 1996, AJ, 112, 1487 [Google Scholar]
  25. Hills, J. G., & Day, C. A. 1976, Astrophys. Lett., 17, 87 [NASA ADS] [Google Scholar]
  26. Koppelman, H. H., & Helmi, A. 2020, A&A, 645, A69 [Google Scholar]
  27. Lee, Y. S., Beers, T. C., Sivarani, T., et al. 2008, AJ, 136, 2022 [Google Scholar]
  28. Lindegren, L., Hernandez, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Lindegren, L., Bastian, U., Biermann, M., et al. 2020, A&A, 649, A4 [Google Scholar]
  30. McCrea, W. H. 1964, MNRAS, 128, 147 [Google Scholar]
  31. Moehler, S. 2001, PASP, 113, 788 [Google Scholar]
  32. Monaco, L., Bellazzini, M., Ferraro, F. R., & Pancino, E. 2003, ApJ, 597, L25 [NASA ADS] [CrossRef] [Google Scholar]
  33. Montenegro, K., Minniti, D., Alonso-García, J., et al. 2019, ApJ, 872, 206 [NASA ADS] [CrossRef] [Google Scholar]
  34. Newberg, H. J., Yanny, B., Grebel, E. K., et al. 2003, ApJ, 596, L191 [Google Scholar]
  35. Niederste-Ostholt, M., Belokurov, V., Evans, N. W., & Peñarrubia, J. 2010, ApJ, 712, 516 [NASA ADS] [CrossRef] [Google Scholar]
  36. Norris, J. E., & Hawkins, M. R. S. 1991, ApJ, 380, 104 [CrossRef] [Google Scholar]
  37. Paunzen, E., Bernhard, K., Hümmerich, S., et al. 2019, A&A, 622, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Pelisoli, I., & Vos, J. 2019, MNRAS, 488, 2892 [Google Scholar]
  39. Pier, J. R. 1982, AJ, 87, 1515 [NASA ADS] [CrossRef] [Google Scholar]
  40. Rain, M. J., Ahumada, J. A., & Carraro, G. 2021, A&A, 650, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Riello, M., De Angeli, F., Evans, D. W., et al. 2020, A&A, 649, A3 [Google Scholar]
  42. Ruhland, C., Bell, E. F., Rix, H. W., & Xue, X. X. 2011, ApJ, 731, 119 [NASA ADS] [CrossRef] [Google Scholar]
  43. Sandage, A. R. 1953, AJ, 58, 61S [CrossRef] [Google Scholar]
  44. Santucci, R. M., Placco, V. M., Rossi, S., et al. 2015, ApJ, 801, 116 [Google Scholar]
  45. Sérsic, J. L. 1968, Atlas de Galaxias Australes (Cordoba: Obs. Astron.) [Google Scholar]
  46. Sirko, E., Goodman, J., Knapp, G. R., et al. 2003, AJ, 127, 899 [Google Scholar]
  47. Smith, K. W., Bailer-Jones, C. A., Klement, R. J., & Xue, X. X. 2010, A&A, 522, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Smolinski, J. P., Lee, Y. S., Beers, T. C., et al. 2011, AJ, 141, 89S [CrossRef] [Google Scholar]
  49. Sommer Larsen, J., & Christensen, P. R. 1986, MNRAS, 219, 537 [NASA ADS] [CrossRef] [Google Scholar]
  50. Sommer-Larsen, J., Christensen, P. R., & Carter, D. 1989, MNRAS, 238, 225 [NASA ADS] [CrossRef] [Google Scholar]
  51. Starkenburg, E., Youakim, K., Martin, N., et al. 2019, MNRAS, 490, 5757 [Google Scholar]
  52. Taylor, M. B. 2005, ASP Conf. Ser., 347, 29 [Google Scholar]
  53. Xiang, M. S., Liu, X.-W., Yuan, H.-B., et al. 2015, MNRAS, 448, 822 [NASA ADS] [CrossRef] [Google Scholar]
  54. Xue, X. X., Rix, H. W., Zhao, G., et al. 2008, ApJ, 684, 1143 [Google Scholar]
  55. Yanny, B., Newberg, H. J., Johnson, J., et al. 2009, ApJ, 700, 1282 [NASA ADS] [CrossRef] [Google Scholar]
  56. Yao, S., Wu, X.-B., Ai, Y. L., et al. 2019, ApJS, 240, 6 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: ADQL Queries

For the parallax candidates (parallax error < 20% dataset):

SELECT *

FROM gaiaedr3.gaia_source

WHERE parallax_over_error > 5

AND parallax > 0

AND phot_bp_mean_flux_over_error > 10

AND phot_rp_mean_flux_over_error > 10

AND phot_bp_rp_excess_factor < 1.3+0.06*power(bp_rp,2)

AND phot_bp_rp_excess_factor > 1.0+0.015*power(bp_rp,2)

AND phot_g_mean_mag+5+5*(log10(parallax/1000))> -1

AND phot_g_mean_mag+5+5*(log10(parallax/1000))< (2.5-(2.77*bp_rp))

AND bp_rp> -0.1

AND bp_rp< 0.53

AND (sqrt(power(pmra,2)+power(pmdec,2)))*(4.74/parallax)> 145

For the proper motion candidates (parallax error ≥20% dataset):

SELECT *

FROM gaiaedr3.gaia_source

WHERE parallax_over_error < = 5

AND phot_bp_mean_flux_over_error > 10

AND phot_rp_mean_flux_over_error > 10

AND phot_bp_rp_excess_factor < 1.3+0.06*power(bp_rp,2)

AND phot_bp_rp_excess_factor > 1.0+0.015*power(bp_rp,2)

AND phot_g_mean_mag+5+5*log10((sqrt((pmra*pmra)+(pmdec*pmdec)))/1000)< 11.88

AND phot_g_mean_mag+5+5*log10((sqrt((pmra*pmra)+(pmdec*pmdec)))/1000)> 6.80

AND (sqrt((pmra*pmra)+(pmdec*pmdec)))/(sqrt((pmdec_error*pmdec_error)+(pmra_error*pmra_error)))> 5

AND bp_rp> -0.1

AND bp_rp< 0.53

AND abs(b)>=50

The following ADQL call was used to find the number of objects within Gaia DR2 that are within 5 arcsec (0.001388889 deg) of a coordinate:

SELECT mine.source_id, count(*) as within_5_arcsec

from user_username.filename as mine

join gaiadr2.gaia_source as gaia

on 1 = contains(

point(’ICRS’, mine.ra, mine.dec),

circle(’ICRS’, gaia.ra, gaia.dec, 0.00138889)

)

group by mine.source_id

This returns a file containing the Gaia DR2 source_id and the number of objects to be found within 5 arcsec. If there are no apparent neighbours within 5 arcsec then a value of 1 (the star itself) will be returned.

All Tables

Table 1.

Selection criteria applied to Gaia DR2 to define blue horizontal-branch star candidate selection from colour and absolute G magnitude (parallax selection).

Table 2.

Comparison of the Gaia DR2 and EDR3 selection criteria and results.

Table 3.

Structure of the online version of the blue horizontal-branch candidate star catalogue.

All Figures

thumbnail Fig. 1.

Typical B-type spectra displaying the difference between main-sequence and blue horizontal-branch stars. Upper panel: example of a typical SDSS B-type spectrum from a blue- horizontal-branch star (SDSS J131916.15–011404.9 from Geier et al. 2015). Lower panel: example of a typical LAMOST B-type spectrum from a main-sequence star. Some reference lines for hydrogen and helium are indicated by vertical lines. SDSS and LAMOST spectra are comparable due to their similar resolutions and spectral coverage.

In the text
thumbnail Fig. 2.

Result of cross-matching the Xue et al. (2008) reference dataset with Gaia DR2. Upper panel: Xue et al. (2008) objects plotted as turquoise dots on the Gaia DR2 colour-magnitude diagram exhibiting a large scatter in absolute magnitude due to incorrect parallax measurements. Those with parallax error < 20% are plotted as red large dots. Lower panel: Lindegren et al. (2018) ‘clean subset selection C’ (black) compared with Xue et al. (2008) objects with parallax error < 20% (red) on the same scales as the upper plot in the Gaia DR2 colour-magnitude diagram.

In the text
thumbnail Fig. 3.

Gaia DR2 colour absolute magnitude diagram showing the blue horizontal-branch star region. Top panel: Gaia DR2 objects that satisfy the modified (Lindegren et al. 2018) criteria (gray points) with those with the stringent halo criteria applied (small blue dots). Xue et al. (2008) objects (red circles) and Behr (2003) objects (yellow triangles) are superimposed. The green trapezium denotes the Gaia DR2 colour-magnitude coordinate space (final CMD criteria) considered as the blue horizontal-branch region. The main-sequence turn-off can be seen at around (GBP − GRP) = 0.6 and Gabs = 3. Blue stragglers are found extending the main-sequence beyond that. Our selection was designed to minimise the blue-straggler contamination and focus on the blue horizontal-branch stars. Bottom panel: for comparison, the Lindegren et al. (2018) ‘clean subset selection C’ comprised of disc stars within 100 pc which are primarily main-sequence stars (black) with red trapezium denoting the Gaia DR2 colour-magnitude coordinate space considered as the blue horizontal-branch region.

In the text
thumbnail Fig. 4.

Sky distribution of parallax selection candidates. Regions of high stellar number density (> 50 000 stars per square degree) have been excluded.

In the text
thumbnail Fig. 5.

Gaia DR2 reduced proper motion versus colour diagram and reduced proper motion frequency plot used to define the selection criteria for the objects with a parallax error ≥ 20%. Left panel: reduced proper motion (H) versus colour (GBP − GRP) diagram for the blue horizontal-branch candidates with parallax error < 20%. Right panel: distribution of reduced proper motion (H) seen in the blue horizontal-branch candidates with a parallax error < 20%.

In the text
thumbnail Fig. 6.

Gaia DR2 colour and reduced proper motion plot: objects conforming to the crowded region criterion as well as the total proper motion and photometric quality criteria (grey dots), objects also conforming to the blue horizontal-branch and halo criteria (blue crosses), and Xue et al. (2008) objects (red squares).

In the text
thumbnail Fig. 7.

Xue et al. (2008) blue horizontal-branch reference objects (red dots) displayed in DR2 (upper panel) and EDR3 (lower panel). Those objects with a parallax error of less than 20% are shown as blue squares.

In the text
thumbnail Fig. 8.

Top panel: objects that satisfy the modified (Lindegren et al. 2018) criteria (gray points) with those with the stringent halo criteria and the crowded region criterion applied (small light blue dots). Xue et al. (2008) objects (red circles) and Behr (2003) objects (yellow triangles) are superimposed. The green trapezium denotes the Gaia DR2 colour-magnitude coordinate space (final CMD criteria) considered as the (blue horizontal-branch region) displayed in DR2 (upper panel) and EDR3 (lower panel).

In the text
thumbnail Fig. 9.

Distribution of reduced proper motions for the parallax selection in Gaia DR2 (red) and Gaia EDR3 (blue).

In the text
thumbnail Fig. 10.

Distance distribution of parallax selection objects in Gaia DR2 (red) and Gaia EDR3 (blue).

In the text
thumbnail Fig. 11.

Distribution of apparent magnitudes for the parallax selection (upper panel) for Gaia DR2 (red) and Gaia EDR3 (blue). Distribution of apparent magnitudes for proper motion selection (lower panel) in Gaia DR2 (red) and Gaia EDR3 (blue).

In the text
thumbnail Fig. 12.

Differentiating between blue horizontal-branch stars and blue stragglers. Upper panel: Gaia EDR3 objects that satisfy the modified (Lindegren et al. 2018) and the final CMD criteria (grey dots) with those with the stringent halo criteria and the crowded region criterion applied (blue circles); SDSS SEGUE objects plotting in the blue horizontal-branch region (yellow squares); and SDSS SEGUE objects plotting in the blue straggler region (magenta triangles) in TEFF-LOGG space. Lower panel: SDSS SEGUE TEFF × LOGG plot with blue horizontal-branch (yellow square) and blue straggler stars (magenta triangles).

In the text
thumbnail Fig. 13.

Catalogue of parallax selection objects. Upper panel: Gaia DR2 colour-magnitude diagram showing the blue horizontal-branch parallax candidates (blue) overlain on a Lindegren et al. (2018) ‘clean subset selection C’ (black). Lower panel: sky distribution of blue horizontal-branch parallax candidates.

In the text
thumbnail Fig. 14.

Parallax candidates plotted comparing the distance calculated from parallax to the theoretical distance assuming that all blue horizontal-branch objects have the same absolute magnitude (1 mag). The parallax error is plotted on the colour axis.

In the text
thumbnail Fig. 15.

Catalogue of proper motion selection objects. Upper panel: Gaia DR2 colour-reduced proper motion diagram showing the blue horizontal-branch proper motion candidates (blue) overlain on all Gaia DR2 objects with total proper motion errors of < 1% (gray). Lower panel: sky distribution of blue horizontal-branch proper motion candidates.

In the text
thumbnail Fig. 16.

Distribution of Gaia EDR3 apparent G magnitudes of blue horizontal-branch candidate objects selected using parallax (blue) and using proper motions (turquoise).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.