# A search for distant, pulsating red giants in the southern halo

Nicolas Mauron, Kamo S. Gigoyan, Tim R. Kendall, Kelly M. Hambleton

► To cite this version:

Nicolas Mauron, Kamo S. Gigoyan, Tim R. Kendall, Kelly M. Hambleton. A search for distant, pulsating red giants in the southern halo. Astronomy and Astrophysics - A&A, 2021, 650, pp.A146. 10.1051/0004-6361/201937005 . hal-03217629

HAL Id: hal-03217629

https://hal.science/hal-03217629

Submitted on 22 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A search for distant, pulsating red giants in the southern halo*

N. Mauron<fnref n="1" />, K. S. Gigoyan<fnref n="2" />, T. R. Kendall<fnref n="3" />, and K. M. Hambleton<fnref n="4" />

Received 28 October 2019 / Accepted 16 December 2020

## ABSTRACT

To investigate the asymptotic giant branch (AGB) population in the Galactic halo, we search for pulsating AGB stars at a heliocentric distance *D* > 50 kpc. Our research is based on the Catalina Southern Survey (CSS) catalogue of variables, comprising 1286 long-period variables (LPVs) with $\delta$ < -20°. We first focus on the 77 stars in the cap $|b| > 30°$ for which spectral M-type or C-type classification can be derived from Hamburg-ESO objective prism spectra. Most of these are oxygen-rich (M-type) and very few are carbon rich. The periods are in the range 100–500 days, and CSS amplitudes are up to 3 mag. In this small sample, no halo AGB star is fainter than $K_{s0}$ = 12.5. This may be due to the scarcity of AGBs in the outer halo, or insufficient instrumental depth. Leaving aside spectral information, we then searched for even fainter pulsators ($K_s$ > 12.5) in the entire CSS catalogue. Gaia astrometry makes it possible to identify some contaminants. Our final result is the identification of ten candidate distant LPVs. If these ten stars obey the fundamental mode K-band period–luminosity relation used for Miras and small-amplitude Miras, their distances are between 50 and 120 kpc from the Sun. In a diagram showing distance versus Gaia tangential velocity, these ten stars have positions consistent with that of other objects in the halo, such as globular clusters and dwarf galaxies. We detect some underluminous AGBs that deserve further study. Finally, the halo LPVs ressemble the slow redder variable of globular clusters when colour and periods are compared.

**Key words.** stars: late-type – Galaxy: halo

## 1. Introduction

The stellar populations located in the Galactic halo contain information about the general properties of the Galaxy, such as its mass, formation and evolution (Helmi 2008; Bland-Hawthorn & Gerhard 2016; Deason et al. 2021). They are also essential for comparing the diverse halos of similar galaxies (Monachesi et al. 2016; Harmsen et al. 2016). Large and recent surveys use RR Lyr variable stars (Hernitschek et al. 2016; Sesar et al. 2017; Stringer et al. 2021), blue horizontal branch stars (Thomas et al. 2018; Starkenburg et al. 2019), K-type giants (Xue et al. 2015; Janesh et al. 2016; Yang et al. 2019), or M-type giants (Bochanski et al. 2014). In this paper we focus on the population of evolved long-period variables (LPV). We follow previous works by Huxor & Grebel (2015) centred on carbon stars, by Grady et al. (2019) who used LPVs as age indicators in the Galactic components and as tracers of merging debris, and Mauron et al. (2019a,b) who studied M-type LPVs in the northern halo and Sagittarius tidal arms.

Long-period variables evolve along the asymptotic giant branch (AGB) and deserve attention because they are more luminous than some of the tracers cited above. Miras or small-amplitude Miras are rarer than giants by more than an order of magnitude, but, like RR Lyr, follow the period–luminosity (PL)

relation. There are old LPVs in many globular clusters (GCs) with periods in the range of ~30–250 days. If the nascent Galaxy contained a large number of now completely dissolved GCs (Gnedin & Ostriker 1997; Odenkirchen et al. 2001; Lee et al. 2004), one might expect a fossil population of cluster LPVs in the halo. Similar and more numerous LPVs can also originate from dislocated dwarf galaxies.

The old AGB populations in GCs are M-type. These cluster stars are in most cases semi-regulars (SRa or SRb type), or, more rarely, Miras possessing a large flux amplitude $\Delta V > 2.5$ mag (Clement 2001, 2017). The halo AGB population also comprises cool N-type carbon (C) stars that have intermediate ages; from the census of halo C stars of Huxor & Grebel (2015), there are about 40 stars at more than 50 kpc from the Sun, with 9 at more than 100 kpc. Roughly two-thirds of cool N-type C stars belong to the Sgr stream. Concerning the M stars, Mauron et al. (2019b) presented a catalogue covering the section of sky from $\delta$ > -20°. These latter authors found 57 stars at *D* > 50 kpc and 6 at > 100 kpc; again, about two-thirds of them are in the Sgr stream.

Here, our goal is to increase the known sample of LPVs with *D* > 50 kpc by examining the complementary region with $\delta$ < -20°. We use the Catalina Southern Survey (CSS) LPV catalogue of Drake et al. (2017), and separate the M-type from C-type stars by exploiting objective prism plates. Attention was paid to contaminants: in Mauron et al. (2019a), we discovered that numerous faint and distant stars are due to instrumental artefacts. Thanks to the Gaia EDR3 data, some supplementary contaminating stars can be discarded. After analysing the data in Sect. 2, we discuss our results in Sect. 3. We finally conclude

<footnote n="1">¹ Université de Montpellier, Laboratoire Univers et Particules de Montpellier CNRS-IN2P3, Place Bataillon, 34095 Montpellier, France
e-mail: nicolas.mauron@umontpellier.fr</footnote>

<footnote n="2">² NAS RA V.A. Ambartsumian Byurakan Astrophysical Observatory (BAO), Byurakan, Armenia</footnote>

<footnote n="3">³ 18 Market Square, Northampton NN1 2DL, UK</footnote>

<footnote n="4">⁴ Department of Astrophysics and Planetary Science, Villanova University, 800 Lancaster Avenue, Villanova, PA 19085, USA</footnote>

<footnote>* A detailed catalogue of the 77 high-latitude M or C stars (Full Table 5) is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/650/A146</footnote>
that very few LPVs are as distant as 50 kpc in the CSS survey, and only ten distant candidates are found.

## 2. Analysis

The CSS catalogue contains 1286 LPVs distributed over an area of ~10 000 square degrees at a declination of $\delta < -20^{\circ}$. In particular, it provides periods *P*, time-averaged magnitudes $m_{CSS}$, and amplitudes $\Delta m_{CSS}$, based on observations made from 2005 to 2013. Its limiting magnitude is $m_{CSS} \sim 17$. The CSS instrument has no filter, and for red stars similar to those of this work, $m_{CSS}$ is close to a I-band magnitude. A majority of these 1286 LPVs are in the Galactic disc, with very few at high Galactic latitude, as shown in Fig. 1 of Mauron et al. (2019a). Also, a fraction of the halo of the Large Magellanic Cloud is covered, and a number of LPVs are detected in the Fornax and Sculptor dwarf galaxies located at distances of ~140 and 90 kpc, respectively (van den Bergh 2000). This shows that the CSS catalogue is a suitable resource for detecting very distant LPVs ($D > 50$ kpc).

Of the 1286 CSS LPVs, 93 have been spectroscopically classified by exploiting the Hamburg-ESO objective prism survey (HES, Christlieb et al. 2004). Of these 93 stars, 16 are false variables due to a polluting Mira (their list is in Sect. 2). In Sect. 2.1, we focus on the remaining 77 stars and analyse their properties. In Sect. 2.2, we proceed by considering objects too faint to be classified with HES. Finally, in Sect. 2.3, we extend our search over the whole CSS LPV catalogue. The astrometric data that we use are mostly from the Gaia EDR3 catalogue.

### 2.1. M-type and C-type stars from the CSS catalogue and the Hamburg-ESO survey

Figure 1 shows a colour–magnitude diagram for stars that could be classified M or C on HES plates. The interstellar $E_{B-V}$ colour excess from Schlafly & Finkbeiner (2011) is used to take into account Galactic extinction<fnref n="1" />. More precisely, we adopt $A_V = 3.10 E_{B-V}$, $A_J = 0.87 E_{B-V}$, and $A_{K_s} = 0.35 E_{B-V}$. We consider that all extincting dust is located between the Sun and the sample star, because most of the LPVs are much higher above the Galactic plane than the thickness (~100 pc) of the Galactic dust layer. It can be seen that M-type stars largely dominate our sample. Stars belonging to the Large Magellanic Cloud (LMC), Fornax, or Sculptor are indicated. For one object, KG-14 in the LMC, it is not clear whether the spectrum is C or M-type, but its colour favours a carbon type that we adopt here. One can see that a small number of carbon-rich LPV stars are detected over the imprint of the HES plates. Table 1 lists their 2MASS identifications and main properties.

In Fig. 1, faint stars ($K_s > 8$) appear slightly bluer than brighter stars, indicating that our sample is not homogeneous. This inhomogeneity is reinforced by the diagram of Fig. 2, where CSS amplitude is plotted versus period. There is first a main group with period from ~100 to 300 days, and amplitude $\Delta m_{CSS}$ rising from 0.6 to ~3 mag. Most of the C stars are slightly to the right of this group. A second family of 10 LPVs exist with $P > 220$ days and relatively smaller amplitude ($\Delta m_{CSS} < 0.8$ mag) compared to objects of the main group. We show later that this large-period small-amplitude group is also separated with a colour based on Gaia and 2MASS photometry.

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Near-infrared colour–magnitude diagram of spectrally classified stars, showing the dereddened $K_{s0}$ magnitude as a function of $(J - K_s)_0$ from 2MASS photometry. Uncertainties are comparable to symbol sizes. Red objects are to the right, faint objects at the top. Circles: M-type stars. Filled symbols: carbon stars. The labels correspond to the KG number: KG-2, 21, 62, 67 are in the halo. Objects labelled ‘3 Scl’ and ‘4 Scl’ are in Sculptor. The spectrum of KG-14 (filled triangle) is unclear, but its colour favors a carbon type.</figcaption></figure>

<table><caption>**Table 1.** Properties of the stars labelled in Fig. 1.</caption><thead><tr><th>KG</th><th>2MASS name</th><th><i>m</i><sub>CSS</sub></th><th><i>P</i></th><th><i>K</i><sub>s</sub></th><th>Sp.</th><th></th></tr></thead><tbody><tr><td>2</td><td>00165576-4400406</td><td>14.71</td><td>440</td><td>7.07</td><td>C</td><td>halo</td></tr><tr><td>3</td><td>00595893-3328351</td><td>16.23</td><td>196</td><td>12.27</td><td>C</td><td>Scl</td></tr><tr><td>4</td><td>01012084-3353047</td><td>16.18</td><td>91:</td><td>13.14</td><td>M</td><td>Scl</td></tr><tr><td>5</td><td>01015347-6518233</td><td>15.71</td><td>162</td><td>12.09</td><td>C</td><td>LMC</td></tr><tr><td>8</td><td>02391532-3415083</td><td>17.71</td><td>196</td><td>13.62</td><td>C</td><td>For</td></tr><tr><td>14</td><td>04360427-6347499</td><td>17.76</td><td>324</td><td>11.21</td><td>C</td><td>LMC</td></tr><tr><td>21</td><td>04563115-3129327</td><td>14.26</td><td>147</td><td>10.65</td><td>C</td><td>halo</td></tr><tr><td>62</td><td>21271642-3051573</td><td>14.29</td><td>233</td><td>9.49</td><td>C</td><td>halo</td></tr><tr><td>67</td><td>21434114-3414311</td><td>14.78</td><td>158</td><td>10.98</td><td>C</td><td>halo</td></tr></tbody></table>

**Notes.** We list the KG number, 2MASS identification, Catalina CSS-band mean magnitude, period *P* (in days), $K_s$ magnitude, spectral classification and membership. The spectral type of KG-14 is uncertain on HES plates, but its colour favors a C type, adopted here. The last column gives membership. The coordinates of KG-62 are erroneous in Mauron et al. (2019a).

Below, we first comment on the C stars, and discuss later on the more numerous M stars.

### 2.1.1. C stars

The brighter C-type star, KG-2, has a 2MASS single-epoch magnitude $K_s = 7.07$, a period of $P = 440$ days, and peak-to-peak $m_{CSS}$-band amplitude $\Delta m_{CSS} = 2.5$ mag. This star is HG-2 in the catalogue of Huxor & Grebel (2015). Neither the origin nor the metallicity is known for this star located at high Galactic latitude ($b = -71^{\circ}$). If this star is similar to solar-metallicity, cool C stars of the solar neighbourhood, with a Gaia-based distance, then $M_{K_s} \approx -8.5 \pm 0.8$ from Fig. 11 of Abia et al. (2020). If it is similar to carbon stars of the more metal-poor Small Magellanic

<footnote n="1">¹ The colour excess is provided by the NASA/IPAC Galactic Dust Reddening and Extinction service at https://irsa.ipac.caltech.edu/applications/DUST</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Catalina amplitude $\Delta m_{CSS}$ versus period. Circles: M-type stars. Filled circles: carbon stars. The two stars in the upper right-hand corner are the halo carbon star KG-2 ($K_s = 7.1, D \sim 14$ kpc, discussed in the text) and the M-type star WOH S24 (alias KG-15; $K_s = 10.0$ at 51 kpc), located in the Large Magellanic Cloud.</figcaption></figure>

Cloud (SMC), we obtain a very similar estimate, $M_{K_s} \approx -8.60$ by applying the period-luminosity relation<fnref n="2" />.

Consequently, we find the distance of KG-2 from the Sun and from the Galactic plane to be $D=14$ kpc. Concerning the uncertainty on $D$, we have $D = 10^{**}0.2(K_s - M_{K_s} + 5)$ and errors on $M_{K_s}$ and $K_s$ have to be considered. The first is typically $\sigma_{M_{K_s}} = 0.15$ mag (i.e. the dispersion of the PL relation; Whitelock et al. 2009, their Table 2). The second originates from the difference between the required but unknown time-average $K_s$-band flux and the available single-epoch 2MASS photometry. This difference is $\delta K_s = 0.2\Delta m_{CSS} + 0.08$ mag (see details in Mauron et al. 2019b). Adding these errors in quadrature leads to a relative error on $D$ of 9% for $\Delta m_{CSS} = 0.2$, to 30% for $\Delta m_{CSS} = 2.5$. Therefore, the distance of KG-2 is $14 \pm 4$ kpc. This distance is larger than the one derived by Huxor & Grebel (11$\pm$4) because they use the Wesenheit indices and consequently take into account circumstellar extinction better for this very red star. KG-2 is angularly located far from the lane of C stars in the Sgr arms and its origin is unclear. One possibility is that it may have been ejected from the Galactic disc, as found by Grady et al. (2019) for some Miras.

The C stars KG-21, KG-62, and KG-67 are warm, that is $(m_{CSS} - K_s)_0 \sim 3.5-4.8$, compared to KG-2 with $(m_{CSS} - K_s)_0 \sim 7.6$. They are short-period variables (see Table 1) and have $K_s \approx 9-11$, which is fainter than KG-2. While KG-62 is named HG-108 in the census of Huxor & Grebel (2015), KG-21 and KG-67 are not in the list of these latter authors. If they are members of the Sgr tidal arms, their [Fe/H] is probably between $-1.2$ to $-0.5$, which is the range of metallicity in the Sgr arms (see Fig. 2 of Carlin et al. 2018). Because the SMC has a mean average abundance [Fe/H] $\sim -0.73$ (van den Bergh 2000) or $\sim -1.0$ (McConnachie 2012), our best choice is to adopt again the

<table><caption>**Table 2.** Sixteen false LPV objects.</caption><thead><tr><th>Star</th><th>α</th><th>δ</th><th>Mira</th><th>θ</th><th>P<sub>KG</sub></th><th>P<sub>Mira</sub></th></tr></thead><tbody><tr><td>KG-22</td><td>85.43000</td><td>-38.9572</td><td>SS Col</td><td>26</td><td>364</td><td>358</td></tr><tr><td>KG-36</td><td>192.93978</td><td>-26.7755</td><td>EP Hya</td><td>23</td><td>168</td><td>168</td></tr><tr><td>KG-42</td><td>207.26013</td><td>-28.5105</td><td>W Hya</td><td>514</td><td>400</td><td>389</td></tr><tr><td>KG-41</td><td>207.25596</td><td>-28.4919</td><td>W Hya</td><td>447</td><td>397</td><td>389</td></tr><tr><td>KG-44</td><td>209.18330</td><td>-25.5373</td><td>FT Hya</td><td>32</td><td>214</td><td>216</td></tr><tr><td>KG-45</td><td>209.18796</td><td>-25.5651</td><td>FT Hya</td><td>72</td><td>214</td><td>216</td></tr><tr><td>KG-46</td><td>209.43125</td><td>-31.1054</td><td>TW Cen</td><td>128</td><td>267</td><td>273</td></tr><tr><td>KG-47</td><td>209.43687</td><td>-31.0558</td><td>TW Cen</td><td>54</td><td>267</td><td>273</td></tr><tr><td>KG-48</td><td>209.77937</td><td>-25.8551</td><td>FQ Hya</td><td>28</td><td>179</td><td>177</td></tr><tr><td>KG-52</td><td>219.36920</td><td>-20.3208</td><td>LY Lib</td><td>28</td><td>289</td><td>287</td></tr><tr><td>KG-55</td><td>302.31581</td><td>-57.3415</td><td>BR Pav</td><td>36</td><td>251</td><td>250</td></tr><tr><td>KG-56</td><td>302.31626</td><td>-57.3173</td><td>BR Pav</td><td>52</td><td>252</td><td>250</td></tr><tr><td>KG-63</td><td>322.63793</td><td>-53.9635</td><td>X Ind</td><td>51</td><td>225</td><td>224</td></tr><tr><td>KG-83</td><td>340.79609</td><td>-41.5202</td><td>DS Gru</td><td>46</td><td>255</td><td>260</td></tr><tr><td>KG-89</td><td>352.03795</td><td>-47.4441</td><td>RU Phe</td><td>49</td><td>291</td><td>288</td></tr><tr><td>KG-91</td><td>359.12315</td><td>-49.7759</td><td>R Phe</td><td>44</td><td>267</td><td>263</td></tr></tbody></table>

**Notes.** The columns are the KG name, 2MASS coordinates (J2000) in degrees, polluting Mira star, separation angle in arcsec between the CSS object and polluant, period of the KG object in days (from CSS), period P of the Mira (from GCVS; for SS Col, period is from the AAVSO VSX catalog).

SMC C-type period luminosity relation (see footnote above). We find that KG-21, KG-62, and KG-67 are at distances of $28 \pm 9, 24 \pm 5$, and $35 \pm 4$ kpc, respectively, with $Z = 17, 17$, and $26$ kpc, respectively. The uncertainties on $D$ essentially originate from $\Delta m_{CSS} = 2.7, 1.6$, and $0.5$ mag, for KG-21, KG-62, and KG-67, respectively.

The star KG-8 is located in Fornax and we classified the object as C-type. This star is listed as not-periodic in the near-infrared monitoring program by Whitelock et al. (2009). In their Table 3, this star is F25006. However, a period of 196 days is clearly seen in the better quality CSS data<fnref n="3" />. Its amplitude is $\Delta m_{CSS} = 1.3$ mag. The SMC PL relation for C stars leads to a distance of $140 \pm 25$ kpc, which agrees with $D(\text{Fornax}) = 147 \pm 12$ kpc given by McConnachie (2012). Agreement is also found for KG-3 located in Sculptor and classified C-type. With $K_s = 12.27, P = 196$ days, and $\Delta m_{CSS} = 1.2$, our method yields $D = 75 \pm 13$ kpc, which is compatible with $D(\text{Sculptor}) = 86 \pm 6$ kpc from McConnachie (2012).

These agreements on external galaxies supports the view that the distances of KG-2, KG-21, KG-62, and KG-67 estimated above are plausible. This means that they are largely out of the Galactic disc, and thus we consider them to be in the halo. Our distances imply heights from the Galactic plane of between 10 and 26 kpc, much larger than 1 kpc seen for cool C stars of the disc (Abia et al. 2020).

## 2.1.2. M stars

The list of the 16 LPVs that could be classified M-type, but are in fact false variables, is given in Table 2. They are in general field dwarfs with significant proper motion, angularly located close to a polluting bright Mira.

Among the remaining 77 stars, there are five M-type stars that are angularly close or in the LMC. We checked that their periods and $K_s$ magnitudes provide a correct distance when the Soszyński et al. PL relation for M-type LMC AGBs is applied.

<footnote n="2">² Soszyński et al. (2007) provides for SMC C-stars of sequence C the following relation $M_{K_s} = -7.15 - 4.22 \log_{10}(P/200)$, where P is in days, $K_s$ is the time-averaged magnitude determined using I-band light curves (Sect. 3 of their paper), and where a SMC distance modulus of 18.99 is used. We note that, throughout this paper, we assume all stars (M-type or C-type) to lie on the pulsation sequence C. Their distances would increase by a factor 1.8 if they were to lie on the more luminous sequence C. For more details, see Table 2 and Figs. 3 and 4 of Soszyński et al. (2007).</footnote>

<footnote n="3">³ CSS light curves can be seen on the Catalina site http://nunuku.caltech.edu/cgi-bin/getcssconedb_releaseImg.cgi</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** APOGEE radial velocities of stars located within one degree from the centre of the Sculptor dwarf galaxy. The LPV KG-4 is indicated. Each point represents the average of 11 individual APOGEE visits. This average is plotted along the *X* axis, and the scatter around this average is on *Y* axis. The median error of an individual visit is only 0.18 km s⁻¹. The high scatter for KG-4 is consistent with its variability.</figcaption></figure>

These stars are KG-13, KG-15, KG-16, KG-17, and KG-18, and we obtain 50 ± 14, 58 ± 19, 53 ± 7, 51 ± 13, and 54 ± 10 kpc, respectively, in good agreement with the LMC distance of 51 kpc (adopting a distance modulus of 18.52±0.09, from McConnachie 2012).

The star KG-4 is located at 0.345 deg from the centre of the Sculptor dwarf galaxy, which is 1.8 times the half-light radius of this galaxy ($r_h = 0.19$ deg, from McConnachie 2012). This star is M-type according to our examination of HES plates, and the Catalina database suggests $P = 91$ d, $\Delta m_{CSS} = 0.5$ mag. However, the ASAS-SN survey (Shappee et al. 2014; Jayasinghe et al. 2018) offers a slightly better light curve and a comparable time span, and finds KG-4 to be aperiodic. If we apply the SMC PL relation for M-type stars as explained below, we find $D = 65 \pm 7$ kpc, smaller than, but marginally compatible with the distance of Sculptor 86 ± 6 (McConnachie 2012). Its radial velocity determined by the APOGEE experiment (Jönsson et al. 2020) shows that KG-4 is member of Sculptor, and the important velocity scatter of individual visits is consistent with its pulsating variability (Fig. 3).

The next finding of this section is the identification of five M-type stars (KG-19, KG-23, KG-72, KG-88, and KG-93) that are faint ($m_{CSS} \sim 15$), with *P* in the range 140–240 days and $\Delta m_{CSS} > 0.8$ mag. Their CSS names are CRTS J044714.7-622118, CRTS J054202.5-553828, CRTS J220821.6-291754, CRTS J230418.1-405245, and CRTS J235855.3-430054, respectively. We obtain their distances by adopting the SMC period–luminosity relation for M-type stars from Soszyński et al. (2007):

$$M_{K_s} = -7.34 - 4.14 \log_{10}(P/200). \qquad (1)$$

Similarly to C stars above, this SMC relation is chosen because the mean SMC metallicity ([Fe/H] = −0.75 to −1.0) coincides with the average metallicity distribution of the Sgr tidal arms and because many of our stars can be members of these arms. Therefore, we apply relation (1), which implies distances of 50–64 kpc. If we had adopted a LMC PL relation, the

distances would be systematically (but only) 10% smaller. The main characteristics of these five stars are in Table 3.

The third result of this section is that among the LPVs that could be spectrally classified, there are none fainter than $K_s$ 0 = 12.5. We did find fainter stars that could be seen and classified M, but their variability is spurious (see Mauron et al. 2019a for details), or they are members of dwarf galaxies as indicated. This lack of faint halo LPVs could be due to an instrumental limitation. For example, identifying faint LPVs may have been difficult in the process of making the CSS catalogue of Drake et al. (2017) that forms the basis of our study. Also, spectral classification is difficult or impossible for faint stars. Alternatively, it could be due to a strong decrease in the LPV number with increasing distance beyond ~30 kpc. Consequently, in the following section, we focus on the CSS LPVs that lie in the HES area, but could not be classified C or M.

## 2.2. Stars without M/C classification from the HES

There are 35 LPVs that could not be classified because no spectrum is visible on the Hamburg plates or because of overlapping spectra. Their Catalina $m_{CSS}$ magnitudes range from 15 to 18. Of these 35, 7 are in Fornax, 19 can be associated with the Large Magellanic Cloud, and 6 are false faint Miras (because their light is contaminated by a true nearby, bright Mira).

Among the remaining 3 unclassified LPVs, the first (hereafter named A) has remained unnoticed in the literature so far, but is interestingly faint ($K_s = 14.2$). The second (called B) was mentioned by Grady et al. (2019) among their six most distant stars, and could be at ~100 kpc. These two stars are discussed below. The CSS names of A and B are CRTS J050449.9-401314 and CRTS J135224.2-314452. Their properties are given in Table 3. Finally, the third star is a confirmed halo C-type star (2MASS J22065366-2506282) with $m_{CSS} = 14.07$, $K_s = 8.95$, $P = 328$ days, and named HG-111 by Huxor & Grebel (2015). It is possible that this star escaped our classification because a deep minimum occurred when the HES plate was exposed.

To conclude, when considering the 35 LPVs lying in the area of the Hamburg survey but not classified M or C, we find only two interesting LPVs, stars A and B listed in Table 3.

## 2.3. All other stars from the CSS catalogue

In the preceding sections, the stars A and B are faint ($K_s = 13.6$ and 14.2) and presumably distant, at about 100 kpc if the PL relation applies. Here, we present the extension of our search to all CSS LPVs, and not only those in the HES imprint. We first noticed one peculiar object, 2MASS J15204636-2533017, with $l = 340^\circ$, $b = +26^\circ$, $m_{CSS} = 14.93$, $K_s = 11.70$, $J - K_s = 1.01$. Its period is 323 days, and $D \sim 90$ kpc if the $K_s$-band PL relation is valid ($M_{K_s} = -8.2$). However, a strong indication that it is closer to us is its proper motion of $1.4 \pm 0.1$ mas yr⁻¹, which implies a relatively high tangential velocity ($V_t$) of 640 km s⁻¹ if the large distance is correct⁴. In addition, the light-curve amplitude is rather small, $\Delta m_{CSS} = 0.32$ mag. Therefore, this star is likely not a pulsating AGB. It might be an ellipsoidal red giant variable similar to those of the LMC, which are ~2.5 mag less luminous than AGB stars (see Fig. 2 of Soszyński et al. 2007).

To find very distant LPVs, we selected objects with $K_{s0} > 12.5$ and obeying ($m_{CSS} - K_s) > 2.0$ as our 77 sample stars do. We found 24 stars: eight are false variable stars, one is in Sculptor, and nine are in Fornax. Of the six remaining

$$^4 V_t(\text{km s}^{-1}) = 4.75D(\text{kpc}) \times \mu(\text{mas yr}^{-1}).$$
<table><caption>**Table 3.** Properties of our candidate, distant CSS LPV stars (see text for their selection).</caption><thead><tr><th>Name</th><th>α</th><th>δ</th><th>l</th><th>b</th><th>m<sub>CSS</sub></th><th>P</th><th>Δm<sub>CSS</sub></th><th>K<sub>s</sub></th><th>J - K<sub>s</sub></th><th>A<sub>V</sub></th><th>D<sub>PL</sub></th><th>δD<sub>PL</sub></th><th>μ</th><th>V<sub>t</sub></th></tr></thead><tbody><tr><td>KG-19</td><td>71.811385</td><td>-62.355530</td><td>272</td><td>-38</td><td>15.13</td><td>139</td><td>0.78</td><td>11.82</td><td>1.03</td><td>0.10</td><td>50</td><td>7</td><td>1.85</td><td>440</td></tr><tr><td>KG-23</td><td>85.510835</td><td>-55.641182</td><td>263</td><td>-32</td><td>15.65</td><td>167</td><td>1.28</td><td>11.65</td><td>1.08</td><td>0.15</td><td>54</td><td>10</td><td>1.51</td><td>390</td></tr><tr><td>KG-72</td><td>332.090176</td><td>-29.298389</td><td>19</td><td>-54</td><td>15.34</td><td>160</td><td>1.20</td><td>11.79</td><td>0.99</td><td>0.05</td><td>56</td><td>10</td><td>1.00</td><td>240</td></tr><tr><td>KG-88</td><td>346.075719</td><td>-40.879162</td><td>354</td><td>-64</td><td>14.95</td><td>141</td><td>1.27</td><td>12.29</td><td>0.77</td><td>0.03</td><td>63</td><td>12</td><td>1.26</td><td>380</td></tr><tr><td>KG-93</td><td>359.730982</td><td>-43.014912</td><td>333</td><td>-71</td><td>15.15</td><td>235</td><td>1.56</td><td>11.41</td><td>0.87</td><td>0.03</td><td>64</td><td>14</td><td>0.75</td><td>230</td></tr><tr><td>A</td><td>76.208748</td><td>-40.221203</td><td>244</td><td>-37</td><td>17.44</td><td>108</td><td>0.55</td><td>14.16</td><td>1.06</td><td>0.08</td><td>119</td><td>14</td><td>0.26</td><td>150</td></tr><tr><td>C</td><td>141.992615</td><td>-34.560715</td><td>262</td><td>+12</td><td>17.33</td><td>139</td><td>1.30</td><td>13.61</td><td>1.20</td><td>0.53</td><td>112</td><td>21</td><td>0.07:</td><td>40:</td></tr><tr><td>D</td><td>154.438690</td><td>-29.276445</td><td>267</td><td>+23</td><td>16.36</td><td>184</td><td>0.45</td><td>12.68</td><td>1.11</td><td>0.20</td><td>93</td><td>10</td><td>0.29</td><td>130</td></tr><tr><td>B</td><td>208.101207</td><td>-31.747896</td><td>318</td><td>+29</td><td>16.92</td><td>128</td><td>1.50</td><td>13.64</td><td>1.02</td><td>0.18</td><td>108</td><td>22</td><td>0.55</td><td>280</td></tr><tr><td>E</td><td>217.840607</td><td>-36.831627</td><td>324</td><td>+22</td><td>15.76</td><td>110</td><td>0.80</td><td>12.81</td><td>0.88</td><td>0.26</td><td>64</td><td>9</td><td>0.72</td><td>220</td></tr></tbody></table>

**Notes.** $m_{CSS}$ is the CSS time-average magnitude. $P$ is the period in days. $\Delta m_{CSS}$ is the CSS amplitude (peak to peak). $K_s$ and $J - K_s$ are from 2MASS. $A_V$ is the Galactic V-band extinction in mag, and $A_K/A_V = 0.11$. $D_{PL}$ is the heliocentric distance in kpc derived with the PL sequence C relation for M-type (Eq. (1)). $\delta D$, in kpc, is the distance uncertainty (see text for details). $\mu$ is the Gaia EDR3 proper motion (see also Table 4). $V_t$ is the tangential velocity in km s$^{-1}$ derived from $D_{PL}$ and $\mu$. The KG stars are M-type. The spectral types of stars ABCDE are unknown.

<table><caption>**Table 4.** Gaia EDR3 parallaxes ($\pi$ in mas) and proper motions (in mas yr$^{-1}$) for our ten candidate distant stars.</caption><thead><tr><th>Star</th><th>π</th><th>μ<sub>RA</sub></th><th>μ<sub>Dec</sub></th></tr></thead><tbody><tr><td>KG-19</td><td>+0.005 ± 0.021</td><td>+1.814 ± 0.024</td><td>-0.367 ± 0.031</td></tr><tr><td>KG-23</td><td>+0.023 ± 0.027</td><td>+1.507 ± 0.033</td><td>-0.042 ± 0.036</td></tr><tr><td>KG-72</td><td>-0.001 ± 0.033</td><td>+0.115 ± 0.035</td><td>-0.994 ± 0.030</td></tr><tr><td>KG-88</td><td>-0.040 ± 0.041</td><td>-0.320 ± 0.037</td><td>-1.211 ± 0.029</td></tr><tr><td>KG-93</td><td>-0.018 ± 0.036</td><td>+0.461 ± 0.037</td><td>-0.590 ± 0.034</td></tr><tr><td>A</td><td>-0.060 ± 0.072</td><td>+0.231 ± 0.076</td><td>-0.118 ± 0.081</td></tr><tr><td>C</td><td>+0.010 ± 0.068</td><td>-0.031 ± 0.055</td><td>+0.066 ± 0.056</td></tr><tr><td>D</td><td>-0.036 ± 0.057</td><td>-0.162 ± 0.058</td><td>-0.238 ± 0.048</td></tr><tr><td>B</td><td>-0.049 ± 0.070</td><td>-0.419 ± 0.079</td><td>-0.363 ± 0.090</td></tr><tr><td>E</td><td>-0.002 ± 0.041</td><td>-0.514 ± 0.043</td><td>-0.500 ± 0.039</td></tr></tbody></table>

**Notes.** No correction was added to $\pi$ because of the 0.030 mas zero-point of parallaxes, as recommended in the case of individual stars (Arenou et al. 2018; their Sect. 4.5).

stars, one has a proper motion of 4.6 mas yr$^{-1}$ (discussed below). The remaining five include stars A and B, as expected. The three others are referred to as C, D, and E below. The CSS names of stars C, D, and E are CRTS J092758.1-343339, CRTS J101745.2-291635, and CRTS J143121.7-364954, respectively. All five stars have proper motions of less than 0.8 mas yr$^{-1}$. Tables 3 and 4 give the properties of these five stars (sorted by Right Ascension).

The star with $\mu = 4.6$ mas yr$^{-1}$ is 2MASS J06091200-4717207 and serves as another demonstration of the fact that the PL relation is not always valid. This star has an amplitude of $\Delta m_{CSS} = 1.4$ mag. Its magnitudes are $m_{CSS} = 18.14$ and $K_s = 14.3$. The light curve is periodic with $P = 116$ days. If it was a classical cool semi-regular, its period would put it at 125 kpc, but its proper motion would imply an excessively large tangential velocity of ~2700 km s$^{-1}$. Therefore, it is more probable that its distance is about ten times smaller. This star is underluminous with respect to the main (sequence C, i.e. fundamental mode) classical PL relation.

# 3. Discussion

## 3.1. Catalogue

We provide a catalogue of the 77 LPV stars that were specifically classified on the HES plates. They obey $\delta < -20^\circ$, and

almost all obey $|b| > 30^\circ$. This latitude limit is due to the HES imprint, although KG-25, KG-26, KG-40, and KG-54 are between $b = -26^\circ$ and $b = -29^\circ$, and KG-30 is at $b = +22^\circ$. An abridged version of the first five lines is given in Table 5. We list: J2000 coordinates, Catalina mean magnitude, period, CSS-band amplitude, $K_s$(2MASS), flag, distance in kpc, uncertainty on distance, star denomination KG-nn, M or C type, and a comment. The flag is set to 1 for normal halo LPV cases. The flag is >1 for members of dwarf galaxies or other peculiarities. The distances of the halo LPVs are in the range 3–65 kpc, and their heights above or below the Galactic plane are more than 5 kpc for 57 of them.

Our list can be compared to that of Grady et al. (2019). After selecting stars with $\delta < -20^\circ$ and $|b| > 30^\circ$ in the Grady list, we obtain 34 objects. Inspecting this sample shows that 8 are in the LMC, 2 are in Fornax, 1 is in Sculptor, 2 are false variables$^5$. There are 14 are in common with our sample, and finally, 7 remain that are not in our list. These 7 are not within the HES plates or impossible to classify with certainty. Inversely, of our initial 93 stars in the same sky imprint, we eliminate 16 false variable, 7 in the LMC, 2 in Sculptor, 1 in Fornax. Consequently, this cleaning procedure results in samples of 21 and 67 stars from the Grady list and our selection, respectively.

The main reason for which the Grady sample is smaller than ours is that these latter authors apply an amplitude cut, favouring high-amplitude variables that preferentially obey the sequence-C Mira PL relation. In contrast, the stars we seek and find are generally periodic, low-amplitude SRa/SRb variables. Grady et al. adopt the definition of Gaia amplitude introduced by Deason et al. (2017): $A_{Gaia} = \sqrt{(N_{obs}) \times (\sigma_{FG}/FG)}$ where $N_{obs}$ is the number of observations and FG is the G-band flux. Also, when building their sample, Grady et al. require $\log(A_{Gaia}) > -0.55$. However, this cut discards numerous interesting stars for us. To illustrate this point, we first consider all Miras and periodic SRa and SRb of the General Catalogue of Variable Stars (Samus et al. 2017) and we plot the histogram of $\log(A_{Gaia})$ (Fig. 4). It can be seen that the large majority (73%) of periodic SRa/SRb variables is missed with the Grady amplitude cut.

<footnote>$^5$ CRTS J213033.1-535749 at $\alpha = 322.6379$ deg., $\delta = -53.9635$ deg., $m_{CSS} = 17.65$, $K_s = 15.6$, P(CSS) = 223 days, polluted by X Ind located at 51" with P(GCVS) = 225 d; and CRTS J224311.1-413113 at $\alpha = 340.7961$, $\delta = -41.5202$, $m_{CSS} = 16.85$, $K_s = 15.2$, P(CSS) = 255 d, polluted by DS Gru at 46" with P(GCVS) = 260 d.</footnote>
<table><caption>**Table 5.** Catalogue of Catalina southern LPVs with spectral classification and distances.</caption><thead><tr><th>α (deg)</th><th>δ (deg)</th><th>m<sub>CSS</sub></th><th>P(d)</th><th>Δm<sub>CSS</sub></th><th>K<sub>s</sub></th><th>Flag</th><th>D(kpc)</th><th>δD</th><th>Name</th><th>Spect.</th><th>Note</th></tr></thead><tbody><tr><td>3.1427</td><td>-22.9213</td><td>12.93</td><td>171.1</td><td>1.93</td><td>9.22</td><td>1</td><td>16</td><td>4</td><td>KG-01</td><td>M-type</td><td></td></tr><tr><td>4.2324</td><td>-44.0113</td><td>14.71</td><td>440.1</td><td>2.40</td><td>7.07</td><td>1</td><td>13</td><td>4</td><td>KG-02</td><td>C-type</td><td></td></tr><tr><td>14.9956</td><td>-33.4764</td><td>16.23</td><td>196.0</td><td>1.07</td><td>12.27</td><td>6</td><td>75</td><td>12</td><td>KG-03</td><td>C-type</td><td>in Sculpt</td></tr><tr><td>15.3368</td><td>-33.8847</td><td>16.18</td><td>91.5</td><td>0.75</td><td>13.14</td><td>6</td><td>59</td><td>8</td><td>KG-04</td><td>M-type</td><td>in Sculpt ?</td></tr><tr><td>15.4728</td><td>-65.3065</td><td>15.71</td><td>162.3</td><td>0.50</td><td>12.09</td><td>3</td><td>59</td><td>7</td><td>KG-05</td><td>C-type</td><td>C in LMC</td></tr></tbody></table>

**Notes.** The full table is available at the CDS.

The same conclusion is reached when we consider 135 M-type pulsators of the Sagittarius leading arm<fnref n="6" />. These stars satisfactorily obey the main PL relation. Again, the Grady cut would eliminate more than half of these Sgr arm M-type pulsators. Consequently, we chose not to apply this *Gaia* amplitude selection here, but we use *Gaia* parallaxes and proper motions (via tangential velocities) wherever possible to eliminate stars for which the PL distance is questionable.

## 3.2. Colour comparison with AGBs of globular clusters

The colours of our halo LPVs can be compared to those of globular clusters. The *Gaia*-2MASS colour $G_{RP} - K_s$ is adopted here because it offers sufficient wavelength baseline, and because $G_{RP}$ is the *Gaia* filter for which amplitudes are the lowest, compared to G or $G_{BP}$. In addition *Gaia* provides homogeneous photometry for both our LPV sample and the globular cluster stars. The interstellar extinction for the $G_{RP}$ passband is 1.83 $E(B - V)$ (Wang & Chen 2019). We first consider the four metal-poor clusters M2, M3, M5, M13 at metallicities [Fe/H] of -1.65, -1.50, -1.29 and -1.53, respectively. For these clusters, García-Hernández et al. (2015) found a total of 44 AGB stars of first or second generation. The colour histograms of Fig. 5 show that the LPVs studied here are almost all redder than those of these four clusters.

A complementary view is obtained by plotting colours versus periods when periods are available (Fig. 6). We consider two globular cluster possessing a substantial number of identified LPVs: ω Cen with [Fe/H] = -1.60 and with 15 known red, slow variable stars (Lebzelter & Wood 2016), and 47 Tuc ([Fe/H] = -0.70) with 38 variables (Lebzelter & Wood 2005). Figure 6 shows that the majority of the cluster LPVs have $P < 100$ days, with a few located on a rising band with larger periods and a redder colour. We note that most of our halo LPVs, including the ten distant stars (except one), are located on this band, suggesting that many cluster-like LPVs with periods in the range ~30 to 100 days remain to be identified in the halo.

Figure 6 also shows that the large-period small-amplitude LPVs seen in Fig. 2 and plotted with triangles in Fig. 6 are clearly separated from the other ‘normal’ LPVs. Because they are relatively blue for their periods, we suspect that those stars do not obey the P-L relation, in a way similar to the variable V13 of 47 Tuc in its long-period mode (Fig. 14 of Lebzelter & Wood 2005). Therefore, in our catalogue, distances for these stars are not given.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Histograms of the *Gaia* DR2 amplitude for Mira and SRa/SRb periodic variables of the General Catalogue of Variable Stars (Samus et al. 2017). The *Gaia* amplitude has been defined by Deason et al. (2017) and combines the *Gaia* G flux, its fluctuations, and the number of measurements (see Sect. 3.1, 2nd paragraph). The *Gaia* amplitude histogram is plotted in log here.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Histograms of the colour ($G_{RP} - K_s$), corrected for interstellar extinction, for the M-type halo LPVs (continuous line), and for the LPV of metal-poor globular clusters M2, M3, M5, M13, whose metallicities are between -1.3 to -1.65 (dotted line).</figcaption></figure>

## 3.3. Distant stars

In Tables 3 and 4, we identified ten cases, with amplitudes $Δm_{CSS} > 0.45$ mag and $K_s$ between 11 and 14. These objects are considered as best candidate distant AGBs, the major finding of this study. However, the applicability of the PL relation in this case remains debatable. If this is the case, distances between 55 and 120 kpc are derived. Table 4 shows that their *Gaia* DR2 proper motions are typically smaller than ~1 mas yr⁻¹. This translates to plausible tangential velocities of less than ~450 km s⁻¹. In Fig. 7, we plot these ten tangential velocities as a function of distance from the Sun (filled circles).

<footnote n="6">⁶ We select them with $B_{Sgr}$ between -15° and +15°, and $Λ_{Sgr} < 180°$. They can be seen in Fig. 4 of Mauron et al. (2019b) (centre panel) along the sinusoid line. The definition of longitude and latitude with respect to the Sgr stream plane is given by Belokurov et al. (2014).</footnote>
<figure><img src="image_7.png" /></figure>

**Fig. 6.** Colour ($G_{RP} - K_s$), corrected for interstellar extinction, plotted as a function of pulsation period. Circles represent the M-type LPVs of our sample. Triangles are the large-period small-amplitude group in our sample seen in Fig. 2. Filled circles are the discovered distant LPVs listed in Table 3. The LPVs of the globular clusters ω Cen ($[\text{Fe/H}] = -1.6$) and 47 Tuc ($[\text{Fe/H}] = -0.70$) are shown for comparison, and are indicated by plus signs and crosses, respectively. Note that a few GC LPVs with $log(P) > 2.1$ are located along the rising red main group. Uncertainties are comparable to symbol size.

For comparison, we also plot two other families of objects: all Galactic GCs, and the few dwarf galaxies for which Gaia DR2 motions have been measured. This figure shows that our candidates have reasonable $V_t$ velocities. Most GCs have $D < 20$ kpc. For these, the average of their velocity ($\sim 250$ km s⁻¹) corresponds to the reflex circular velocity of the Sun, whereas the dispersion ($\pm 250$ km s⁻¹) reflects their intrinsic motions. The halo clusters at $D > 30$ kpc have $V_t$ between 0 and 400 km s⁻¹. The dwarf galaxies have a similar range. In conclusion, given the large uncertainties on the small proper motions, the positions of our ten distant candidates shows them to be reasonably similar to other halo populations (clusters and dwarf galaxies), supporting their halo membership.

## 3.4. An underluminous AGB

We see above that 2MASS J06091200–4717207 is significantly underluminous, by approximately 5 mag. The period and amplitude of this star, 116 days and 1.4 mag, respectively, are close to or larger than those of our ten candidates. To investigate the rarity of such stars, we sought similar stars in the solar neighbourhood by considering the General Catalogue of Variable Stars (Samus et al. 2017). We required Gaia parallaxes of excellent quality, that is $\pi/\delta\pi > 10$, and we found several underluminous cases. One example is CO Sge, an M-type Mira with $P = 190$ d and $\Delta V = 3.6$ mag, which is confirmed by the light curve of the ASAS-SN survey (Shappee et al. 2014; Kochanek et al. 2017). Its parallax is $\pi = 3.45 \pm 0.18$ mas, and $D = 290$ pc. The flux $K_s = 4.70 \pm 0.02$ is not saturated in 2MASS. Because of variability, the mean $K_s$-band signal is within 1.4 mag of this magnitude<fnref n="7" />. We derive from Gaia a mean $K_s$-band absolute magnitude

<figure><img src="image_8.png" /></figure>

**Fig. 7.** Tangential velocities versus heliocentric distances. No correction for the Sun Galactic orbital motion was applied. Our ten stars in Table 3 are plotted with big filled circles. Squares are globular clusters from Harris (2010), with distance relative uncertainty of 5%, and proper motions from Vasiliev (2019). Diamonds are dwarf galaxies: at the top, from left to right: Sgr, Boo I, Sextans, Carina, and Fornax; at the bottom, Ursa Minor, Draco, and Sculptor. Distances are from van den Bergh (2000) and proper motions from Gaia DR2 (Gaia Collaboration 2018). Errors for dwarf galaxies are comparable to the diamond size.

$M_{K_s} = -2.6 \pm 1.4$, while the SMC PL relation gives $-6.4$. Consequently, CO Sge is a striking example of a large-amplitude pulsator with a luminosity deficit of 3–4 mag. If our stars are similarly underluminous, their distance $D$ and velocity $V_t$ are smaller than proposed here.

## 3.5. Contaminants

Here, we address qualitatively the issue of contamination by other halo variable stars in our sample of 77 candidate halo AGB stars. For that purpose, the very recent paper by Chen et al. (2020) (published when our work was being achieved) provides very useful information. Their work is based on the Zwicky Transient Facility (Bellm et al. 2019; Graham et al. 2019) scanning the entire northern sky in three nights, in two passbands ($g, r$), and to magnitude $r \sim 20.5$ ($5\sigma$); its Data Release 2 spans over 470 days. Chen et al. (2020) analyse the DR2 light curves and provide a catalogue of ~800 000 periodic variables. They classify Mira and SR stars with periods $P > 80$ days for Miras and $P > 20$ days for SRs, together with cuts in amplitude, with Amp($g$) > 2.4 for Miras and Amp($g$) < 2.4 for SRs. These latter authors find as many as 12 000 Miras and 120 000 SRs, in part because ZTF covers the Galactic plane.

It is clear from their Fig. 6, which shows luminosity versus period, that a small number of binaries (classified “EW” by them) approach the locus of SRs. Most of these binaries have periods of 0.1–1days, but some are found with $P \sim 100$ d. Their light shapes obey characteristic relations linking the Fourier coefficients, and these relations are used to identify them (Table 1 of Chen et al.). Given that Catalina photometric quality is poorer than that of ZTF, it is very plausible that binaries contaminate the halo SR populations in which we have searched for candidates. For example, in Fig. 11, panel b, of Chen et al. (2020) the top two binary “EW” ZTF curves are quite similar

<footnote n="7">⁷ We found a $K_s$ to $m_{CSS}$ amplitude ratio of $\Delta K/\Delta m_{CSS} \sim 0.4$ for this kind of star (Mauron et al. 2019b, their Fig. 8). A smaller ratio, ~0.15, is given for the Johnson V-band for cool AGBs in Sect. 2.2 of Groenewegen (2006).</footnote>
to the two ZTF SRs shown at the top of their Fig. 12, panel a. If observed with Catalina, the curves would appear almost identical.

In conclusion, our halo LPVs may be contaminated by some binary variables, but this should improve in the future by exploit-
ing the modern ZTF survey that surpasses Catalina in many aspects that are useful to separate true halo AGBs from these contaminants, particularly for small amplitudes. We also note that, for the region $\delta < -20^\circ$ not covered by ZTF, the Catalina survey has presently accumulated unique data over more than 16 years and is ongoing, and this will offer even better light curves.

4. Summary and conclusions

To search for AGB stars in the outer halo (>50 kpc), we started with the CSS catalogue of LPV of Drake et al. (2017) (with $\delta < -20^\circ$, 10 000 deg$^2$, limiting magnitude ~17–18, time span 8 years). By focussing first on the CSS LPVs that are located within the imprint of the southern Hamburg-ESO objective prism plates ($|b| > 30^\circ$), we were able to classify (M or C) 77 of them, after discarding 16 artefacts. The large majority are M-type. Only two new, uncatalogued C stars were found. After excluding members of known dwarf galaxies, all are found to be brighter than $K_s = 12.5$. This could be due to the scarcity of the sought distant objects, or to some instrumental limitation.

To search for fainter, even more distant LPVs, we consid-
ered CSS stars in the whole Drake et al. catalogue, and required
$K_s > 12.5$. After discarding again false variable sources and
LPVs located in the LMC, Fornax, and Sculptor dwarf galax-
ies, we find five candidate distant stars. In our small sample of
ten objects (Table 3), amplitudes $\Delta m_{CSS}$ are in the range 0.45–
1.5 mag, and periods are between 100 and 230 days. We assume
that the fundamental mode PL relation can provide their dis-
tance, and that they are similar to semi-regulars of the SMC.
This system was chosen for its supposed low metal abundance
comparable to that of our halo sample. The derived distances
are found to be between ~50 and 120 kpc; their *Gaia* parallaxes
are not significant, and their tiny proper motions are compatible
with large distances and Galactic halo membership. These ten
candidates are the result of this research. In addition, a catalogue
of the 77 M-type and C-type stars will be made available at the
CDS. During our analysis, we also found several cases of under-
luminous LPVs that deserve attention. Finally, the halo M-type
LPVs studied here resemble the redder slow variables of glob-
ular clusters with $P > 100$ d, but the more numerous LPVs with
smaller periods ($P \sim 30–100$ d) that are found in clusters remain
to be identified in the halo.

The superb light curves and depth of the ZTF for northern stars, the data points acquired over seven additional years by the Catalina experiment for southern stars, and the future light curves of the Large Synoptic Survey Telescope can be considered to develop this line of research. In addition, Gaia proper motions will also improve in the near future. When follow-up spectroscopy providing radial velocities is achieved, we will hopefully have greater knowledge of the population of AGB stars in the outer halo of our Galaxy.

Acknowledgements. The authors thank the anonymous referee for remarks that greatly improved the manuscript. K. H. acknowledges support for this work from the NASA ADAP program under NASA grant 18-ADAP18-228. N. M. thanks Sebastian Otero (AAVSO) for very useful comments, Henri Reboul and Denis Puy for encouragements, and the staff of LUPM for kind assistance, especially

Stéphane Nou and Lydie Le Clainche. N. M. is grateful to the Université de Montpellier and the Institut National de Physique Nucléaire et de Physique des Particules (IN2P3; Centre National de Recherche Scientifique) for very appreciated support. This work makes use of the Catalina database (California Institute of Technology, and NASA), the Two Micron All Sky Survey (U. of Massachusetts and IPAC/Caltech funded by NASA and NSF). We very much appreciated to use the on-line archive of the Hamburg-ESO (Hamburg Observatory and European Southern Observatory) and express thanks to D. Engels. We also exploited the unique data of the European Space Agency mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/gaia. The Gaia archive website is https://archives.esac.esa.int/gaia. Finally, we used the SIMBAD/Vizier facilities offered by and operated at CDS, Strasbourg, France.

**References**

Abia, C., de Laverny, P., Cristallo, S., et al. 2020, A&A, 633, A135

Arenou, F., Luri, X., Babusiaux, C., et al. 2018, A&A, 616, A17

Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002

Belokurov, A., Koposov, S. E., Evans, N. W., et al. 2014, MNRAS, 437, 116

Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529

Bochanski, J., Willman, B., West, A., et al. 2014, AJ, 147, 76

Carlin, J. L., Sheffield, A. A., Cunha, K., & Smith, V. V. 2018, ApJ, 859, L10

Chen, X., Wang, S., Deng, L., et al. 2020, ApJS, 249, 18

Christlieb, N., Reimers, D., & Wisotzki, L. 2004, Messenger, 117, 40

Clement, C. 2001, AJ, 122, 2587

Clement, C. 2017, Wide-Field Variability Surveys: A 21<sub>st</sub> Century Perpective; EPJ Web of Conferences

Deason, A. J., Belokurov, V., Erkal, D., et al. 2017, MNRAS, 467, 2636

Deason, A. J., Erkal, D., Belokurov, V., et al. 2021, MNRAS, 501, 5964

Drake, A., Djorgovski, S. G., Catelan, M., et al. 2017, MNRAS, 469, 3688

Gaia Collaboration (Helmi, A., et al.) 2018, A&A, 616, A12

García-Hernández, D. A., Meszaros, Sz., Monelli, M., et al. 2015, ApJ, 815, L4

Gnedin, O. Y., & Ostriker, J. P. 1997, ApJ, 474, 223

Grady, J., Belokurov, V., & Evans, N. W. 2019, MNRAS, 483, 3022

Graham, M. J., Kulkarni, S. R., Bellm, E. C., et al. 2019, PASP, 131, 78001

Groenewegen, M. A. T. 2006, in Why Galaxies Care About AGB Stars: Their Importance as Actors and Probes, eds. F. Kerschbaum, C. Charbonnel, & R. F. Wing, ASP Conf. Ser., 378, 433

Harmsen, B., Monachesi, A., Bell, E. F., et al. 2016, MNRAS, 466, 1491

Harris, W. E. 2010, ArXiv e-prints [arXiv:1012.3224]

Helmi, A. 2008, A&ARv, 15, 145

Hernitschek, N., Schlafly, E. F., Sesar, B., et al. 2016, ApJ, 817, 73

Huxor, A. P., & Grebel, E. K. 2015, MNRAS, 453, 2653

Janesh, W., Morrison, H. L., Ma, Z., et al. 2016, ApJ, 816, 80

Jayasinghe, T., Kochanek, C. S., Stanek, K. Z., et al. 2018, MNRAS, 477, 3145

Jönsson, H., Holtzman, J. A., Allende Prieto, C., et al. 2020, AJ, 160, 120

Kochanek, C. S., Shappee, B. J., Stanek, K. Z., et al. 2017, PASP, 129, 104502

Lebzelter, T., & Wood, P. R. 2005, A&A, 441, 1117

Lebzelter, T., & Wood, P. R. 2016, A&A, 585, A111

Lee, K. H., Lee, H. M., Fahlman, G. G., & Sung, H. 2004, ApJ, 128, 2838

Mauron, N., Gigoyan, K. S., Gigoyan, K. K., et al. 2019a, Astrophysics, 62, 202

Mauron, N., Maurin, L. P. A., & Kendall, T. R. 2019b, A&A, 626, A112

McConnachie, A. W. 2012, AJ, 144, 4

Monachesi, A., Bell, E. F., Radburn-Smith, D. J., et al. 2016, MNRAS, 457, 1419

Odenkirchen, M., Grebel, E. K., Rockosi, C. M., et al. 2001, ApJ, 548, L165

Samus, N. N., Kazarovets, E. V., Durlevitch, O. V., et al. 2017, Astron. Rep., 61, 80

Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103

Sesar, B., Hernischek, N., Mitrović, S., et al. 2017, AJ, 153, 204

Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48

Soszyński, I., Dziembowski, W. A., Udalski, A., et al. 2007, Acta Astron., 57, 201

Starkenburg, E., Youakim, K., Martin, N., et al. 2019, MNRAS, 490, 5757

Stringer, K. M., Drlica-Wagner, A., Macri, L., et al. 2021, ApJ, 911, 109

Thomas, G. F., McConnachie, A. W., Ibata, R., et al. 2018, MNRAS, 481, 5223

van den Bergh, S. 2000, The Galaxies of the Local Group (Cambridge University Press)

Vasiliev, E. 2019, MNRAS, 484, 2832

Wang, S., & Chen, X. 2019, ApJ, 877, 116

Whitelock, P. A., Menzies, J. W., Feast, M. W., et al. 2009, MNRAS, 394, 795

Xue, X.-X., Rix, H.-W., Ma, Z., et al. 2015, ApJ, 809, 144

Yang, C., Xue, X., Li, J., et al. 2019, ApJ, 880, 65