Open Access
Issue
A&A
Volume 698, May 2025
Article Number A109
Number of page(s) 14
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202453125
Published online 03 June 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

The third data release (DR3) from the Gaia mission includes astrometry and broad-band photometry for a total of 1.8 billion objects based on 34 months of satellite operations (Gaia Collaboration 2023). A subset of ∼1.7 million stars have been identified as long-period variable (LPV) star candidates with amplitudes > 0.1 mag in the Gaia photometric G-band (Lebzelter et al. 2023). This subset includes stars on the asymptotic giant branch (AGB), which are in general variable stars and are therefore a major component of the LPV population. The measurement of parallaxes by Gaia holds the promise that sizeable samples of AGB stars with reliable distances will become available in the Milky Way and especially in the solar neighbourhood. They can be used for an observational validation of AGB evolutionary models. However, the selection of AGB stars from LPV catalogues needs to remove stars of different nature, for example young stellar objects (YSOs), which can mimic the variability behaviour of AGB LPV stars (Lucas et al. 2024). Validation samples of well-curated AGB stars are therefore valuable to test the selection criteria for AGB stars in LPV catalogues.

The majority of low- and intermediate-mass stars (∼0.8 − 8 M) evolve through the AGB phase, where they experience high mass-loss rates (≳10−7 M yr−1). When the mass-loss rate surpasses 10−6 M yr−1, the dust shell eventually becomes opaque to visible light (Habing 1996), hiding stars in the advanced AGB phase from detection by Gaia. The dust of the circumstellar envelope (CSE) absorbs the radiation of the central star and thermally re-emits it mainly in the infrared (IR). At the same time, the gas of the CSE can be the source of strong maser emissions, for example by the OH molecule at 1612, 1665, and 1667 MHz. In this paper, we use the term ‘OH/IR stars’ to refer to all AGB stars with OH maser emission. This includes the classical dust-enshrouded OH/IR stars as well as optically bright LPVs, for example Mira variables.

According to evolutionary models, at the end of AGB evolution a different chemical composition is expected for the stellar surface and in succession for the CSE, depending on their progenitor mass. At solar metallicities, the stars end carbon rich (C/O > 1) for progenitor masses 1.5 M ≲ M ≲ 3 M, while below and above this mass range, the stars terminate the AGB evolution with an oxygen-rich chemistry (Di Criscienzo et al. 2016; Ventura et al. 2018). Depending on the metallicity, these boundaries can vary by several tenths of a solar mass (Tosi et al. 2022).

Although these models are commonly accepted and applied to explain the AGB population in the Magellanic Clouds (Dell’Agli et al. 2015a; Dell’Agli et al. 2015b; Pastorelli et al. 2019, 2020), their application to Galactic AGB stars need to be validated, as the Milky Way has a higher metallicity. This validation is pending because the lack of reliable distances makes determinations of stellar luminosities and the subsequent matches of individual stars with evolutionary tracks uncertain. For a small sample of optically obscured oxygen-rich stars close to the end of AGB evolution (‘extreme OH/IR stars’) that are assumed to be located in the Galactic bulge, Jiménez-Esteban & Engels (2015) found stellar luminosities within the mass range in which the stars should already have been converted into carbon stars. In contrast to this result, Blommaert et al. (2018) questioned the use of a common distance to these OH/IR stars and argued all OH/IR stars in the bulge have low-mass stars (∼1.5 M) as progenitors. A local AGB sample with reliable Gaia parallaxes may allow us to test the models anew.

In addition to distances, the determination of reliable stellar luminosities in the Galaxy requires well-covered spectral energy distributions (SEDs) to avoid imprecise bolometric corrections. The SEDs can today be built with unprecedented coverage from a wealth of photometric surveys that are easily accessible through the Virtual Observatory (VO). These catalogues cover wavelength ranges from the UV to the far-IR. Based on theoretical models of SEDs representative of the AGB evolutionary phase, gaps in the observational coverage of the SEDs can be filled, and bolometric corrections can be avoided completely.

This work is the first step to create a Gaia Catalogue of Galactic AGB stars with reliable distances from Gaia-DR3, focusing on O-rich AGB stars identified by their OH maser emission. In forthcoming works, we will add a sample of C-rich AGB stars and a sample of O-rich AGB stars without detections of OH maser emission.

The paper is structured as follows. In Sect. 2 we describe the identification of OH maser emitting AGB stars in Gaia-DR3. In Sect. 3 we construct their SEDs using VO tools, obtain their bolometric fluxes from SED fitting, and determine stellar luminosities when trustworthy Gaia parallaxes are available. The final Gaia OH/IR star sample is defined in Sect. 4. This is followed by the description of the properties of the sample in Sect. 5. Finally, we present our conclusions in Sect. 6.

2. Identifying OH maser sources in Gaia DR3

As starting point for the construction of the Gaia OH/IR star sample, we used the database of circumstellar OH masers1 by Engels & Bunzel (2015) in its version 2.5 (Engels & López-Martí 2024). This database (hereafter, the OH maser database) is a compilation of masers that were predominantly detected at 1612 MHz until 2022, and it is restricted to stars in the Milky Way. Overall, the database contains 16 879 individual observations with detections of 2905 different objects. This number was decreased to 2857 sources after we removed some entries that were erroneous or duplicated. We call this sample of sources from the OH maser database ‘the OH maser sample’.

The majority of sources in this sample are AGB stars. Based on the definition given in the Introduction, they are considered OH/IR stars in this paper.

The source coordinates in the OH maser database were obtained from many different bibliographic resources, and as a result, their accuracy largely varies, from milli-arcseconds to several arcminutes in the worst cases. Moreover, a direct cross-match of a catalogue of OH masers with an optical catalogue such as Gaia might yield too many spurious counterparts because counterparts of the OH masers are often strongly obscured optically. Therefore, matches with field stars, especially in regions with a high source density such as the Galactic plane, cannot be excluded. For these reasons, our cross-matching procedure followed intermediate steps to improve the source coordinates with the help of infrared catalogues and to ensure that the OH maser counterparts were detectable at shorter wavelengths. This procedure is explained below.

2.1. Cross-match with infrared catalogues

Because of the substantial mass-loss rates, AGB stars are strong infrared sources and are generally detected by surveys such as the Infrared Astronomical Satellite (IRAS), the Two-micron All Sky Survey (2MASS), and/or the Wide-field Infrared Survey Explorer (WISE). This allows us to improve the coordinates of the OH maser database. In addition, a steep rise in the spectral energy distribution in the near-infrared (NIR) indicates that a source is deeply obscured by circumstellar dust and is therefore expected to have a faint or even undetectable optical (e.g. Gaia) counterpart. Furthermore, the infrared photometry might help us to distinguish between correct and spurious cross-matches, and to some extent, it might also be useful to identify the different types of objects included in the OH maser database.

The catalogues we used in this first step of the cross-match procedure were the AllWISE Source catalogue (3.4–22 μm; Cutri et al. 2013) and the 2MASS Point Source catalogue (1.2–2.2 μm; Cutri et al. 2003). Both catalogues have a similar astrometric accuracy: In AllWISE, positions for non-saturated sources with a high signal-to-noise ratio (S/N; 8 < W1 < 12 mag) are accurate to about 50 mas. In the case of 2MASS, the source positions are accurate to 70–80 mas in the range 9 < Ks < 14 mag and to 120 mas for brighter sources. The astrometric accuracy of fainter sources decreases monotonically with brightness. However, 2MASS has a better spatial resolution than AllWISE (1″/pixel versus 1.375″/pixel), and it may therefore resolve sources that are blended in the latter catalogue. This allowed us to improve the identification of the OH maser counterpart in the infrared in case of confusion in the AllWISE data or when the mid-infrared source was very bright, which sometimes made it difficult to correctly place the photocentre of the star on the WISE images. The higher resolution of 2MASS also helped us to identify the most likely counterparts in the Gaia DR3 data, whose spatial resolution is even better.

In crowded areas or in areas with extended diffuse infrared emission, we were unable to reliably identify the counterparts in one or both of the infrared catalogues. This issue was especially severe towards the Galactic centre. Therefore, OH masers within 2 degrees of the Galactic centre were removed from the cross-match procedure. We removed 279 sources from the OH maser sample for this reason, and the cross-match was then made for the remaining 2578 sources.

We used a radius of 15″ for the cross-match between the OH maser database and AllWISE, and we used a radius of of 3″ between AllWISE and 2MASS. All matches were visually inspected to assess their reliability. In some cases, an incorrect counterpart had been assigned because it happened to be the one nearest to the input coordinates. Sometimes, the coordinates from the OH maser database were so inaccurate that the actual IR counterpart lay well beyond the search radius. If an object was found in the visual inspection that seemed to be the correct counterpart of the OH maser source, then this AllWISE or 2MASS counterpart was assigned to the source.

Because of the shorter wavelengths, it was expected that 2MASS contained fewer counterparts to the OH maser catalogue than AllWISE: Some sources would be deeply obscured and not detectable in this range. On the other hand, all 2MASS counterparts of OH maser sources were expected to have counterparts in the AllWISE Source Catalog. However, in practice, we found a number of sources without an AllWISE counterpart that had a counterpart in 2MASS; these sources are visible in the WISE images, but lack an entry in the catalogue, probably due to issues in the source extraction (they are mostly highly saturated sources).

Unique IR counterparts were identified for 2427 OH masers, but for 151 detected masers listed in the OH maser database, no counterpart was found or it was not possible to assign a reliable counterpart in either of the infrared catalogues. The reasons were manifold. Part of the OH maser sources may have no (unambiguous) infrared counterpart because the coordinates are grossly inaccurate (> 1′) in the original literature or because the maser was misclassified and has an interstellar origin. In other cases, the identification of the correct IR counterpart was difficult due to crowding or blending with nearby sources or surrounding diffuse infrared emission.

2.2. Cross-match with Gaia DR3

The next step was to cross-match the sample of OH maser sources with unique IR counterparts with the Gaia DR3 catalogue (Gaia Collaboration 2023) using the 2MASS coordinates as input. When the source had no 2MASS counterpart, the AllWISE coordinates were used. The search radius in Gaia DR3 was 1.5″.

As with the infrared cross-matches, we assessed the reliability of the identified counterparts by visual inspection. Because Gaia does not provide any images, we used the images from the Digital Sky Survey 2 (DSS-2; Lasker et al. 1996) and the Pan-STARRS survey (Chambers & Pan-STARRS Team 2017) for this assessment. They cover a similar brightness range as Gaia, but the spatial resolution is far lower. We unexpectedly found a few sources with AllWISE counterparts that had a counterpart in Gaia DR3 even though they lacked a 2MASS counterpart. As for the sources with 2MASS counterparts but no AllWISE counterparts, the near-infrared source in these cases was visible in the 2MASS images, but was not included in the point-source catalogue.

When there were several possible counterparts, we checked the Gaia colours and favoured the counterpart with the reddest colours or, if the colours of the possible counterparts were similar, the closest source to the input position. In some cases, the reliability of the possible counterparts could not be assessed, and the source remained without a Gaia identification.

The cross-match with the Gaia DR3 catalogue resulted in 1487 OH maser sources with unique Gaia counterparts and 940 sources with no (unique) optical identification.

2.3. Parallaxes and distances

Although 88% of the OH maser sources with unique Gaia counterparts have parallax measurements, these are not reliable in many cases, either because they are negative and thus spurious (21% of the sample), or because the astrometric quality is insufficient.

Recent studies of Gaia distances for AGB stars have been made: Andriantsaralaza et al. (2022) analysed the uncertainties of the Gaia parallaxes on these stars by comparing them with parallaxes measured with maser observations with very long baseline interferometry (VLBI). More recently, Bhattacharya et al. (2024) presented a novel method for estimating statistical distances to AGB stars using the SED shape and a brightness estimate. They then compared the obtained distances with the Gaia distances. Both studies came to similar conclusions: It is difficult for Gaia to provide reliable parallaxes for sources with relative parallax errors higher than 20%, and the Gaia distances for distances beyond 2 kpc are underestimated.

Based on these results, we defined parallaxes as “good” when the relative parallax error was lower than 20% and the astrometric excess noise was lower than 1.5 mas. This latter parameter measures the disagreement, expressed as an angle, between the observations of a source and the best-fitting standard astrometric model (Lindegren et al. 2021). Following these criteria, the parallaxes of only 277 sources (about 19% of the sample) were regarded as “good”.

For all sources with a good parallax, we extracted the geometric distances from the catalogue of Bailer-Jones et al. (2021). These distances, which are calculated from Gaia DR3 data using a probabilistic approach, are expected to be more reliable than a simple inversion of the parallax. We did not consider the photogeometric distances that these authors also provided because they were computed using colour information based on a galaxy model, and it is highly unlikely that this prior is applicable to our extremely red sources.

2.4. Overall cross-match results

Based on the cross-match with the AllWISE, 2MASS, and Gaia DR3 catalogues, we identified 1487 sources from the OH maser sample with unique Gaia DR3 counterparts (hereafter, the Gaia sample), 277 of which (19% of the sample) had distance estimates from Bailer-Jones et al. (2021).

The majority of sources in the Gaia sample (1424, or 96%) have counterparts in the near-infrared (2MASS) and the mid-infrared (AllWISE). Eight sources only have an AllWISE counterpart, and 55 sources (4% of the sample) only have a 2MASS counterpart.

In addition, 940 sources from the OH maser database without Gaia counterparts (about 33% of the OH maser sample) had detections in AllWISE and/or 2MASS. Of these, 729 sources (78% of this infrared sample) had counterparts in both catalogues, 178 (19%) were included only in AllWISE, and 33 (about 4%) were included only in 2MASS. The remaining 151 sources in the OH maser sample were not detected in either of these catalogues, had ambiguous counterparts, or were removed because they are located within 2 degrees of the Galactic centre. The sources without Gaia counterparts are not further considered here. Table 1 summarises our cross-matching results.

Table 1.

Results of the cross-match of detected OH masers in the database of circumstellar OH masers (Engels & López-Martí 2024) with the AllWISE, 2MASS, and Gaia DR3 catalogues.

We recall that although the majority of sources in the OH maser database are expected to be oxygen-rich AGB stars, other types of sources such as post-AGB stars, red supergiants (RSG), or unrecognised YSOs are included as well. Therefore, additional steps were required to remove them from the Gaia sample. These steps included the construction and analysis of the spectral energy distributions (SEDs) of the sources. They are explained in the next sections.

3. Spectral energy distributions

The SEDs of the sources from the Gaia sample were built and analysed with the help of the Virtual Observatory SED Analyzer2 (VOSA; Bayo et al. 2008). This is a web-based VO-compliant tool developed by the Spanish Virtual Observatory designed to query accessible photometric catalogues and to obtain stellar physical parameters (e.g. temperature, bolometric flux, and mass-loss rate) of the studied sources by fitting theoretical models to the observational SEDs.

3.1. Photometry compilation and SED construction

First of all, we compiled available photometric data from the ultraviolet to the far-infrared in order to obtain a well-covered SED. To do this, we uploaded the AllWISE, 2MASS, and Gaia photometry collected in the previous step, along with the distance estimations, if available, to VOSA, and queried the VO for additional photometric information from catalogues available through the application.

We retrieved infrared photometry from the following catalogues: the IRAS catalogue of Point Sources (version 2.0) and the IRAS Faint Source catalogue (Beichman et al. 1988); the AKARI/FIS All-Sky Point Source Catalogue (Yamamura et al. 2010) and the AKARI/IRC Mid-Infrared All-Sky Survey (Ishihara et al. 2010); and the MSX6C Infrared Point Source catalogue (version 2.3; Egan et al. 2003). In the latter case, we discarded bands B1 and B2 because the quality of the photometry is low.

We also queried the Deep Near-Infrared Survey (DENIS; Epchtein et al. 1997), as well as the VISTA Variables in the Via Lactea Survey (VVV; Saito et al. 2012) and the VISTA Hemisphere Survey (VHS; McMahon et al. 2013). To avoid fitting issues because our sources are variable in the near-infrared, we only used the VISTA Z and Y bands and the DENIS I band, and we omited the J, H, and Ks bands in all these cases because they overlap with the 2MASS photometry.

In addition to the optical bands from DENIS and the VISTA surveys, we also retrieved photometry from the PanSTARRS 1 Survey (Chambers & Pan-STARRS Team 2017), which provides photometry in the g, r, i, z, and y bands. For a limited number of sources, ultraviolet photometry from the Galaxy Evolution Explorer (GALEX) was also available (Bianchi & GALEX Team 2000).

Bad data and upper limits are already flagged by VOSA based on the information provided by the catalogues. As a general rule, we accepted this information as it was provided by the tool. However, in the few cases in which the visual inspection of the SEDs suggested that a data point that was not considered bad was problematic, we manually flagged it as such.

3.2. SED shapes

The SEDs constructed with the help of VOSA were visually inspected and classified into four main types based on their appearance in a plot of λFλ versus wavelength (see Fig. 1 for some examples). We describe them below.

  • ‘Single-peaked SEDs’ display a clear maximum in the infrared domain. These SEDs correspond to the shape expected for typical AGB stars. The observed distribution is a superposition of emission from the stellar photosphere and the dust in the CSE.

  • ‘Double-peaked SEDs’ have two local maxima, one maximum in the mid-infrared, and the other in the optical to near-infrared domain. These two components are identified with the circumstellar shell and the central object, respectively. This type of SED is characteristic of the post-AGB phase, when the shell is being cleared and the central object is already emerging.

  • ‘Flat SEDs’ present a rather flat distribution in the mid- to far-infrared domain, typically with a steep rise at shorter wavelengths (near-infrared and optical), and no flux drop at the longest wavelengths. They might be a borderline case between the two types described above, or the flat appearance might just be caused by a lack of coverage at the longest wavelengths, where the drop is expected to occur.

  • ‘Peculiar SEDs’ have shapes that do not match any of the previous descriptions. Some of these sources seem to be interlopers (sources with OH masers that are not evolved stars, mainly YSOs), while in other cases, the peculiar shape may be the result of associating one or more incorrect counterpart in the photometric catalogues. It may also be caused by poor-quality data.

The SEDs of about 88% of our sources are classified as single-peaked or flat, and the objects were thus considered AGB star candidates. The next step was to compare them with the predictions from theoretical models.

thumbnail Fig. 1.

Examples of the different types of SEDs found in the Gaia sample. The photometry we used is listed in the insets of the panels, and the solid line connects the data points to guide the eye. From left to right and from top to bottom: Single-peaked SED, double-peaked SED, flat SED, and peculiar SED. According to the literature, the plotted sources are an OH/IR star, a post-AGB star, another OH/IR star, and a YSO, respectively. An SED graph for any source in the final catalogue can be retrieved through the online service (see Appendix A).

3.3. Model fitting

The goals of this step were to confirm the AGB nature of the sources and to estimate their bolometric fluxes, and when a distance value was available, also their absolute luminosities.

For each parameter combination in the chosen model grid, VOSA computes the synthetic photometry for the different filters used to construct the SED, building a synthetic SED that is compared with the observed SED. The best-fitting model is determined by the smallest reduced χ2, although users can choose a different model if they consider it more appropriate. In the calculation of the reduced χ2 for each model of the grid, VOSA weights the photometric data with the inverse of the square of the photometric error.

In our case, the theoretical models used for the fits were the Grid of Red supergiant and Asymptotic giant branch star ModelS (GRAMS) for oxygen-rich stars (Sargent et al. 2011). This model grid explores the following parameters: stellar effective temperature from 2100 K to 4700 K; dust shell inner radii of 3, 7, 11, and 15 R*; and 10.0 μm optical depth from 10−4 to 26. We set no constraints in this parameter space because were only interested in deriving bolometric fluxes and luminosities.

Interstellar extinction corrections were applied to sources with a good parallax determination. The value of the interstellar extinction was obtained from the 3D map of Galactic interstellar dust from Lallement et al. (2022). These authors used Gaia EDR3 photometric data combined with 2MASS measurements to derive the extinction density at λ = 5500 Å towards any direction in the sky in a volume around the Sun of 6 × 6 × 0.8 kpc3. Knowing the distances to the sources, it is thus possible to use these data to estimate the individual integrated extinction in the line of sight for each of them.

About 30% of the sources with a good parallax fall outside the volume that is covered by the extinction map. All these sources are far from the Galactic plane, with Galactic heights larger than ∼350 pc. In these cases, we considered the maximum line-of-sight integrated extinction provided by the map as the interstellar extinction value for the source. Since most of the interstellar dust is located close to the Galactic plane, the cumulative extinction above or below 350–400 pc from the Galactic plane is probably low, and the effect on the overall SED is negligible.

The obtained extinction values were then used in VOSA to deredden the observed photometry before model fitting and bolometric flux estimation. VOSA uses the extinction law by Fitzpatrick (1999) improved by Indebetouw et al. (2005) in the infrared3.

The model fitting of an observational SED built from a variety of single-epoch catalogues was a challenging task. First, the agreement between the model and the observations was hampered by the high variability of the AGB star candidates, especially in the optical domain, where the variability amplitude reaches several magnitudes (see e.g. Fig. 9 of Jiménez-Esteban et al. 2021). In addition, the longest-wavelength bands (beyond 60 μm) were in many cases affected by contamination. To avoid that the fitting results were dominated too much by these effects, we decided to run the fits only in the range between 1 and 60 μm, so that the data points outside this range were not included in any calculations.

On the other hand, the photometry data provided by the different catalogues also have very different absolute errors. The smallest errors usually corresponded to the shortest wavelengths, and therefore, these photometric points had the heaviest weights in the fitting process. The consequence was that in many cases, the NIR data are very well fitted, while the mid- and far-infrared domains are completely discrepant from the best-fitting model. Thus, after several tests, we decided to upload the photometry onto VOSA with zero errors for all bands. When this is done, VOSA considers an error of 10% of the corresponding flux in all cases. This is still lower than the typical variability amplitude of our sources, but better than the usually much smaller absolute photometric errors provided by the catalogues, especially in the near-infrared range.

The fits run with these conditions (zero input errors, and only in the wavelength range between 1 and 60 μm) agreed better overall in the central part of the SED (which contains the peak of the distribution in the case of a single-peaked SED) than the fits performed over the entire range and/or considering absolute photometric errors. Because VOSA computes the bolometric flux using the observed photometry whenever possible and it only interpolates with the model when no observations are available, this ensured that the errors introduced in this computation through variability or contamination were minimised.

All fits were visually inspected to ensure their reliability. When the SED clearly deviated from the model, we considered it a poor fit. Based on this inspection, the majority of SEDs classified as single-peaked (92%) turned out to have relatively good fitting results. A fraction of the flat SEDs also displayed good fits; for some of these sources, the models suggested that a flux drop indeed occurred beyond the range covered by the observations. As expected, all double-peaked and peculiar SEDs had poor fits.

For nine sources in the sample, the use of the restricted wavelength range had the consequence that the number of available photometric points (excluding upper limits) was too low to perform the fit. All of them displayed single-peaked SEDs, except for one, which had a peculiar SED.

4. The Gaia OH/IR star sample

After the SED inspection, the fits for 1190 sources were considered good. This is nearly 80% of the Gaia sample obtained by the cross-match process described in Sect. 2. Although these sources are likely OH/IR stars, a good fit does not guarantee that other types of sources are completely removed from the sample. (To reiterate: We use the designation ‘OH/IR star’ as a general term to refer to AGB stars with OH maser emission). Thus, in order to check the reliability of the identification of OH/IR stars and to remove interlopers from the sample, we cross-matched the whole Gaia sample with the SIMBAD database.

SIMBAD had counterparts for 99.5% of the Gaia sample (1481 sources), the classification of most of which was compatible with an OH/IR star. The largest number of interlopers were post-AGB stars (49 sources), followed by 19 red supergiant stars (RSGs), 13 planetary nebulae, 13 YSOs, two interstellar medium sources, and three multiple stellar systems. In total, the classification of 99 sources (6.7% of the sources with SIMBAD counterparts) was not compatible with an OH/IR star. Of these, only 15 had a good model fit, namely all the RSGs and one post-AGB star. This is ∼1.3% of the sources in the Gaia sample with a good fit. Thus, the expected contamination of the sources with good model fits is very low (∼1–2%).

Consequently, we defined the ‘Gaia OH/IR star sample’ as sources with a good model fit and with a SIMBAD classification compatible with an AGB star. In total, 1172 sources (79%) out of the 1487 members of the Gaia sample were assigned to this sample. Of these, 222 had Bailer-Jones geometric distances and thus luminosity estimates.

We present these 1172 sources in an online table (Table 2), along with their Gaia-DR3, SIMBAD, AllWISE, and 2MASS identifications, their Gaia coordinates, brightnesses, and parallaxes, their geometric distances from Bailer-Jones et al. (2021), their derived bolometric fluxes and, when distances are available, their luminosities from the SED fitting. In the remainder of the paper, we focus on this sample of stars. This table, along with the observational photometry data we used to build and analyse the SEDs, is electronically available in an online archive (see Appendix A) and at the CDS VizieR catalogue service4. Additionally, the sources that were excluded from the Gaia OH/IR star sample are presented in a second online table (Table 3), along with their Gaia-DR3, SIMBAD, WISE, and 2MASS identifications, Gaia coordinates, their SED type, the reliability of the model fit (good or poor) as assessed by eye, the SIMBAD classification when available, and the reason for rejection.

Table 2.

Gaia sample of OH/IR stars.

Table 3.

Sources with Gaia-DR3 counterparts that were excluded from the Gaia sample of OH/IR stars.

The Gaia OH/IR star sample should consist mostly of long-period variable stars and consequently is expected to be contained in the second Gaia catalogue of LPV candidates (Lebzelter et al. 2023). However, only 475 (40%) OH/IR stars could be recovered in the LPV catalogue. The recovery rate is strongly dependent on the optical brightness. Figure 2 shows the number of OH/IR stars that we recovered as function of G-band brightness. The number of OH/IR stars included in the LPV catalogue decreases strongly for objects with G ≥ 17 mag. Optically bright OH/IR stars G ≤ 13 mag are also frequently missed. Lebzelter et al. (2023) obtained recovery rates ≥80% relative to the LPV published by ASAS-SN and OGLE-III independent of Gaia brightnesses. We found high recovery rates like this relative to the Gaia OH/IR star sample only in the brightness range 13 ≤ G ≤ 17.

thumbnail Fig. 2.

Objects of the Gaia OH/IR star sample that are included and missed (blue and red histogram, respectively) in the second Gaia catalogue of LPV candidates (Lebzelter et al. 2023).

We also cross-matched our sample with the sample of nearly 66 000 Mira variables included in the OGLE Collection of Variable Stars (Iwanek et al. 2022). We found 345 sources in common with this catalogue, which is nearly 30% of the Gaia OH/IR sample. The majority of the remaining sources in our sample are located outside the overlapping area between the two catalogues. Most Gaia OH/IR sources within the overlap region that lack a match in the OGLE catalogue are brighter than G ∼ 12 mag, and they lie outside the brightness range of the OGLE survey.

5. Properties of the Gaia sample of OH/IR stars

In this section, we discuss the main properties of the Gaia OH/IR star sample. A more detailed scientific analysis is deferred until the remaining samples (carbon-rich AGB stars and oxygen-rich AGB stars without OH masers) are presented in forthcoming papers.

5.1. Spatial distribution

Figure 3 shows the distribution in the sky of the Gaia OH/IR star sample. As explained in Sect. 2.2, an area of two degrees around the Galactic centre was removed from the sample due to crowding issues.

thumbnail Fig. 3.

Sky distribution of the Gaia sample of OH/IR stars. The sources with good parallaxes (relative parallax error < 20% and astrometric excess noise < 1.5 mas) are shown in red, and the rest are plotted in grey.

The sources are spread throughout the sky, although there is a concentration toward the inner Galactic plane. This was previously reported for known Galactic O-rich AGB stars (Ishihara et al. 2011, and references therein). However, it is important to recall that this sky distribution may be biased because it originates from a compilation of OH maser sources from the literature, which is dominated by some surveys that were made along the Galactic plane (Engels & Bunzel 2015).

Figure 4 shows the distribution of Bailer-Jones et al. (2021) geometric distances for the 222 OH/IR stars with good parallaxes. These stars are located within 4.0 kpc, and only ∼10% of them are located beyond 2 kpc. The mean and median distance of this sample are 1.29 and 1.17 kpc, respectively.

thumbnail Fig. 4.

Histogram showing the geometric distances from Bailer-Jones et al. (2021) for the OH/IR stars with good parallax measurements.

These sources are also shown in Fig. 3 in red. They are less concentrated towards the inner Galactic plane than the whole OH/IR sample. This suggests that this is a local sample of stars.

Figure 5 shows the relative position with respect to the Sun of the OH/IR stars with good parallax and Bailer-Jones et al. (2021) distances. The sample is more numerous toward the inner Galaxy, and it is confined to distances ≤1.1 kpc in the outer Galaxy. There also seems to be some east-west asymmetry. This projected spatial distribution is very similar to the distribution found by Ishihara et al. (2011) for O-rich AGB stars in the Galactic plane for the distance range covered by our sample. Thus, although we cannot completely rule out with the current data that these features are caused by selection effects, they seem to be real.

thumbnail Fig. 5.

Spatial distribution of the OH/IR stars with good parallax projected onto the Galactic plane. The Sun is located at the centre (0,0), and the direction towards the Galactic centre is upwards along the vertical dashed line.

Finally, we show in Fig. 6 the distribution of height above the Galactic plane of the OH/IR sources with Bailer-Jones et al. (2021) distances. In this case, the distribution is very symmetric with respect to the Galactic plane, and the number of sources decreases with height.

thumbnail Fig. 6.

Galactic height distribution of OH/IR sources with geometric distances from Bailer-Jones et al. (2021).

To summarise, the group of OH/IR stars with available distances trace a more or less local sample of OH/IR stars that is somewhat more numerous in the direction of the inner part of the Galaxy and is distributed roughly symmetrically with respect to the Galactic plane. It is expected that the OH/IR stars with a poor parallax determination are mainly located in the direction of the inner galaxy and at farther distances.

5.2. Magnitudes and colours

We now discuss the magnitude and colour distributions of the OH/IR stars with Gaia counterparts. Only sources with good photometry are considered, that is, a photometry quality flag A, B, or C for WISE photometry; a photometry quality flag A, B, C, or D for 2MASS photometry; and a relative error lower than 10% for the Gaia fluxes. Furthermore, all sources with BP > 20.5 mag were removed from the magnitude and colour distributions using this filter because of a known bias that produces systematically unreal blue colours for faint red sources (Fabricius et al. 2021).

Figure 7 shows the magnitude distribution of the whole OH/IR sample (semi-filled histograms) and of the sources with good parallaxes (filled histograms). Four bands are displayed, namely the Gaia G-band, 2MASS Ks, and the WISE W1 and W4 bands. In the G band, the number of detected sources increases towards higher magnitudes. The drop at the highest observed magnitudes is most likely due to incompleteness. Most sources with G <  10 mag also have good parallaxes. This is expected because brighter sources should have better quality detections and better Gaia astrometric parameters. On the other hand, the W4 magnitude distribution of the whole sample and of the good-parallax sample are similar. This suggests that there is no intrinsic difference between the sources with a good parallax and the rest, except that they are probably closer and less affected by interstellar extinction.

thumbnail Fig. 7.

Magnitude distributions of the Gaia OH/IR sample (semi-filled histograms) and of the subsample of sources with good astrometry (filled histograms). From left to right and from top to bottom: Gaia G band, 2MASS Ks band, and the WISE W1 and W4 bands. The percentages refer to the fraction of sources with respect to the total Gaia OH/IR sample of 1172 sources.

Figure 8 shows the distributions of two Gaia colours (BP–RP and G–RP), as well as two infrared colours (J − Ks and W1–W4), for the likely OH/IR stars with Gaia counterparts (semi-filled histograms) and of the subsample of sources with a good astrometry (filled histograms).

thumbnail Fig. 8.

Colour distributions of of the Gaia OH/IR sample (semi-filled histograms) and of the subsample of sources with a good astrometry (filled histograms). From left to right and from top to bottom: Gaia BP–RP and G–RP, 2MASS J − Ks, and WISE W1–W4 colours. The percentages refer to the fraction of sources with respect to the total Gaia OH/IR sample of 1172 sources.

As expected, the sources tend to have quite red optical colours, as is evident from the BP–RP colour distribution. This colour is available for only 11% of the Gaia detected OH/IR stars because many faint sources with red colours were only detected by Gaia in the RP band, and as explained above, sources with BP magnitudes > 20.5 mag were removed. The plot appears to be truncated at BP–RP ∼ 8.3 simply because of the Gaia detection capabilities; the OH/IR that are left undetected are likely to extend the colour range to even higher values.

The BP–RP distribution of the Gaia OH/IR sample peaks at rather red colours BP–RP ∼ 6 mag. This is different from the colour distributions found for Gaia LPV in general. Mowlavi et al. (2018) for DR2 and Lebzelter et al. (2023) for DR3 considered LPV with G-band amplitudes QR5(G) > 1 mag predominantly consisting of Mira variables, which make up 10–20% of all the LPVs that were identified with QR5(G) > 0.2 mag. In amplitude-colour diagrams, this group of Mira variable candidates stands out and has colours BP–RP >  4 mag (Mowlavi et al. 2018). The fraction of stars with BP–RP >  4 in our sample is ∼85%, indicating that the small amplitude QR5(G) < 1 mag variables are underrepresented compared to the general LPV population seen by Gaia.

The good-parallax sources are not at all different in this respect. However, they are among the bluest sources in the near-infrared, as shown in the J − Ks colour distribution. This again suggests a lower extinction than for the whole sample.

5.3. Bolometric fluxes and luminosities

The left panel of Fig. 9 shows the distribution of bolometric fluxes obtained from the SED fitting for our OH/IR sample. The stars with good parallaxes are among the brightest objects, as expected.

thumbnail Fig. 9.

Left panel: Distribution of bolometric fluxes (left) of the Gaia OH/IR sample. The grey and red histograms correspond to the whole sample and to stars with good parallaxes (good ω), respectively. Right panel: Luminosity distribution for the subsample with good parallaxes.

The luminosity of the sources in this latter subsample was computed by VOSA from the bolometric flux and the Bailer-Jones geometric distances. We recall that the bolometric flux was obtained after correcting the observed photometry for interstellar extinction (see Sect. 3.3). This distribution is shown in the right panel of Fig. 9. Most luminosity values are as expected for AGB stars, typically between 1000 and 10 000 L, with an extended tail up to ∼25 000 L. The peak of the distribution is around 4000 L.

The luminosities of five sources are lower than 1000 L: IRAS 17528+1144, IRAS 18042–2517, IRAS 20234–1357, SY Sco, and R Vir. All but IRAS 20234–1357 are in the Gaia-DR3 catalogue of LPV stars (Lebzelter et al. 2023), and their Gaia light curves show clear Mira-like variability, with periods ranging from ∼150 to ∼600 d. The remaining source is included in the ASAS-SN catalogue of variable stars (Jayasinghe et al. 2018). It also shows Mira-like variability with a period of ∼550 d. Hence, all five sources therefore appear to be genuine AGB stars. Most likely, their distances are underestimated, and so are their luminosities.

5.4. Gaia Hertzsprung-Russell diagram

Using the Bailer-Jones et al. (2021) distances and the dereddened Gaia magnitudes, we constructed the Gaia Hertzsprung-Russell (HR) diagram for the OH/IR sample. This diagram displays the G-band absolute magnitude versus one Gaia colour.

Two versions of this diagram are shown in Fig. 10: In the left panel, we show the commonly used representation with the BP–RP colour. Due to the quality criteria for the photometry imposed in Sect. 5.2 and especially because we omitted the Gaia BP photometry for faint red sources, only one-third of the sources with distance estimations can be plotted in this diagram (65 sources). We therefore also show a representation with the G–RP colour in the right panel of Fig. 10, where a larger number of sources (N = 150) from our OH/IR sample can be displayed. In both diagrams, we also plot a sample of Gaia nearby stars within 100 pc following the same selection criteria as used by Jiménez-Esteban et al. (2023). As the interstellar extinction for these stars is negligible, no dereddening was applied to their Gaia photometry.

thumbnail Fig. 10.

Gaia DR3 HR diagrams using the extinction-corrected BP–RP and G–RP colours (left and right panel, respectively) for the subsample of sources with good parallaxes (red circles), compared to the location of the 100 pc stellar population (in grey). In each diagram, only sources with good photometry according to the criteria outlined in Sect. 5.2 are displayed.

In the HR diagram using the BP–RP colour, the sources are spread over a wide colour range 2.3 ≲ (BP − RP)≲7.3 mag. The absolute G-Band magnitude Gabs decreases from Gabs ∼ −4.6 to Gabs ∼ +1.2 with increasing BP–RP colour. The locus of the Gaia OH/IR sample in this diagram is similar to the locus of the LPV candidates in the first Gaia catalogue (Mowlavi et al. 2018), and this is the expected location for AGB stars (Gaia Collaboration 2019).

In Gaia-selected LPV samples, the bulk of objects has colours 2 ≲ BP−RP ≲ 4 (Lebzelter et al. 2023) and shows semi-regular variability (Gavras et al. 2023). The Gaia OH/IR star sample, in contrast, has predominantly redder Gaia colours, which are frequent among large-amplitude periodically pulsating Mira variables (Mowlavi et al. 2018; Gavras et al. 2023). This reflects the larger fraction of OH maser emitters among the AGB stars with higher mass-loss rates and hence redder colours (see the discussion in Lewis et al. 1990 for IRAS colours).

It must be noted that the BP–RP colours are averages, and their uncertainty is influenced by the coverage of the light curves by Gaia observations. The colours can vary by 2–3 mag depending on the phase of the pulsation cycle, and together with the brightness variations, they cause loops in the HR diagram (Gaia Collaboration 2019; Mowlavi et al. 2018). Because the light curve will be better covered, the location in the diagram for a given star may be different in future GAIA releases, while the area covered by AGB stars in the HR diagram is not expected to change globally considered.

The increasing intrinsic faintness in the G band with BP–RP colour can be attributed to the increasing circumstellar extinction. Quite significant bolometric corrections are needed to infer the stellar luminosities from the G-band absolute magnitudes. This is shown in Fig. 11, where the bolometric corrections BC(G) = Mbol − Gabs for the subsample of sources with good parallaxes are given. The bolometric magnitudes Mbol were obtained from the luminosities determined in Sect. 5.3. We fitted the data with a third-degree polynomial, giving

BC ( G ) = 1.662 1.701 ( BP RP ) + 0.317 ( BP RP ) 2 0.031 ( BP RP ) 3 . $$ \begin{aligned} \mathrm{BC(G)}&= 1.662 - 1.701\, (\mathrm{BP} - \mathrm{RP})\nonumber \\&\quad + 0.317\, (\mathrm{BP} - \mathrm{RP})^2 - 0.031\, (\mathrm{BP} - \mathrm{RP} )^3. \end{aligned} $$(1)

thumbnail Fig. 11.

Bolometric corrections BC(G) = Mbol − Gabs for the subsample of sources from the Gaia OH/IR stars with good parallaxes. The outlier in the diagram is IRAS 18319-1001, as discussed in the text. The solid curve shows the fit to the data as given in Sect. 5.4, and the broken curve is the bolometric correction as a function of BP–RP colour proposed by Lebzelter et al. (2019) without interstellar extinction.

The fit is shown as the solid curve in Fig. 11. This fit was previously introduced by Lebzelter et al. (2019) using the original Gaia photometry. Their results are shown in the figure as a broken curve. The bolometric correction curve provides a means to obtain stellar luminosities using (dereddened) Gaia photometry alone. From the scatter around the curve, we estimate an uncertainty of ΔMbol = 0.5(1σ) mag in the colour range 4 ≤ BP − RP ≤ 7, which is reasonably well populated by objects.

A few objects are detached from the main locus of our sample in the Gaia HR diagram. The most obvious case is IRAS 18319-1001, which is placed close to the tip of the red giant branch in the HR diagram using the BP–RP colour (Gabs ≈ –0.4 mag, BP–RP ≈ 1.6 mag, and G–RP ≈ 0.5 mag). In this region, type II Cepheid variable stars are also expected (Gaia Collaboration 2019; Rimoldini et al. 2019). However, IRAS 18319-1001 is the counterpart of the OH maser G22.136-0.389 that was detected by the recent THOR OH maser survey (Beuther et al. 2019), which displays a convincing double-peaked profile typical of an OH/IR star. Inspection of the SED shows that the infrared part is very well fitted with a GRAMS model, in accordance with its classification as an OH/IR star. However, the optical data are brighter than the model, and the brightness toward shorter wavelengths decreases not as strongly. Although the separation between the 2MASS counterpart and the selected Gaia counterpart is about 1″, the association is almost certainly erroneous, as judged from a detailed comparison of the counterpart positions in Pan-STARRS and 2MASS images. The Gaia optical counterpart is a field star in the line of sight of the apparently optically obscured OH/IR star, which dominates the IR emission. Thus, the peculiar position of IRAS 18319-1001 in Fig. 10 (and Fig. 11) is quite likely due to a misidentification. This case shows that although great care was taken to verify the correct association of the Gaia counterparts, a few cases of contamination with field stars may remain.

In the HR diagram that is based on the G–RP colour (Fig. 10, right panel), the brightness of a group of four sources with Gabs ≥ 1.0 and G–RP < 1.5 appears to be too low for their G–RP colour. These are IRAS 11525-5057, IRAS 17020-5254 (both with Gabs ∼ 1.4 mag and G–RP ∼ 1.1 mag), and IRC +30292, IRAS 14591-4438 (both with Gabs ∼ 2.4 mag and G–RP ∼ 1.3 mag). Gaia-DR3 provides epoch photometry for these four sources. All show Mira-like variability in the Gaia G-band light curves. Their periods range from ∼490 to 590 days and the peak-to-peak amplitudes are ∼3 mag in G band, which means that they are genuine long-period large-amplitude variables. The variability information supports an AGB nature of these sources. In contrast to the case of IRAS 18319-1001, their detached positions in the G–RP HR diagram with respect to the bulk of the Gaia OH/IR sample are therefore not caused by incorrect identifications of their Gaia counterpart, but most probably by problems with the Gaia parallax and/or brightness determination.

In addition to IRAS 18319-1001, a group of six sources (IRAS 18413-4616, MSB 38, IRAS 18006-3213, V840 Aql, V718 Cyg, and V CVn) is detached from the main OH/IR star locus in the left panel of Fig. 10 due to their bluer BP–RP colours (Gabs ∼ –4.6 to –2.2 mag, and BP–RP ∼ 2.3–3.4 mag). Gaia G-band light curves are available for all these sources. Four of them (IRAS 18413-4616, MSB 38, V840 Aql, and IRAS 18006-3213) show an irregular variability with a low peak-to-peak amplitude (< 0.3 mag). V CVn is a semi-regular variable with a period of ∼200 d (Cadmus 2024) and a photometric peak-to-peak variability of 1 mag in the G band. The remaining two sources, V718 Cyg and IRC +00363, show Mira-like variability with long periods (∼514 d and ≲950 d, respectively) and a G-band photometric peak-to-peak variability of ∼1 mag.

To summarise, the sources with the bluest colours in the HR diagrams are likely a mix of AGB stars with low-amplitude irregular and large-amplitude semi-regular and regular variability. A general study of the variability properties of the sources in the Gaia OH/IR sample will be presented in a forthcoming paper after the Gaia catalogue of Galactic AGB stars is complete.

6. Conclusions

We presented the first step towards the construction of the Gaia catalogue of Galactic AGB stars by providing a sample of OH/IR stars detected by Gaia DR3. We cross-matched the list of sources with detected OH maser emission (OH/IR stars) from the most recent version of the database of circumstellar OH masers of Engels & López-Martí (2024) with AllWISE, 2MASS, Gaia DR3, and other catalogues in order to retrieve astrometric and photometric information that we used to construct the SEDs of the objects from the near-ultraviolet to the far-infrared. We identified ∼1500 OH maser sources with Gaia detections, while more than 900 OH maser sources in the database with infrared identifications are too obscured to be detected by Gaia. Whenever distances based on trustworthy parallaxes were available, we obtained the interstellar extinction to each source and used it to deredden the observed SED. We then used the GRAMS model grid to fit the SEDs and obtain bolometric fluxes. Stellar luminosities were also determined for sources with reliable distances.

The AGB stars were separated from other OH maser sources based on their SED shape, the fitting results, and information from the literature. The final ‘Gaia OH/IR star sample’ contains 1172 sources, 222 of which have distance and luminosity estimations. These latter sources represent a bright and local subsample, given that most lie at distances within 2 kpc. As discussed by Lebzelter et al. (2023), a radius of 2 kpc qualifies as a reasonable border for the solar neighbourhood, in which most of the objects with good Gaia parallaxes are found.

For more than 80% of the Gaia OH/IR stars, no parallaxes in the Gaia DR3 release match our adopted quality criteria because of circumstellar and interstellar extinction and, quite likely, because they are located beyond the solar neighbourhood. The OH maser database is dominated by a few OH maser surveys that were conducted at low Galactic latitudes and reached the Galactic centre at ∼8 kpc, which is reflected by the concentration of the Gaia OH/IR star sample towards the inner Galactic plane. In contrast, and as expected, the local subsample is evenly distributed on the sky.

The causes for the lack of good-quality Gaia parallaxes of the Gaia OH/IR stars is also evident from the brightness and colour distributions. The local subsample is bright (mostly G ≲ 12), while in general, the Gaia OH/IR stars are optically faint (12 ≲ G ≲ 22). The G-band brightness distribution continuously increases up to G ∼ 19.5, and this would probably continue if the sizeable fraction of OH maser sources that is not detected by Gaia were not excluded by the survey detection limit. When circumstellar and interstellar extinction becomes significant, the SEDs become redder. While the local subsample has mostly colours J − Ks ≲ 2, the remaining Gaia OH/IR stars have typical colours 2 ≲ J − Ks ≲ 7.

The Gaia OH/IR star sample is a well-curated sample of O-rich AGB stars that can be used as a benchmark sample for further studies of LPV with Gaia. It is the first sample in the Gaia Catalogue of Galactic AGB Stars we are assembling. In the next steps, a sample of carbon-rich AGB stars and another sample of oxygen-rich AGB stars without OH maser detections will be added.

Data availability

Full Tables 2 and 3 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/698/A109.

The Gaia Catalogue of Galactic OH/IR Stars can also be queried via an online service available at the Spanish Virtual Observatory (SVO) website (see Appendix A).


Acknowledgments

Authors thank to John Fabio Aguilar Sánchez for the invaluable and useful discussion on distance uncertainties, and to Raquel Murillo-Ojeda for the help in the preparation of the SED graphs for the archive. This work is funded by ESA through the Faculty of the European Space Astronomy Centre (ESAC) - Funding reference 4000139151/22/ES/CM, and by grant No. PID2020-112949GB-I00 by the Spanish Ministry of Science, Innovation/State Agency of Research MCIN/AEI/10.13039/501100011033. This publication 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. It also made use of data products from the Wide-field Infrared Survey Explorer (WISE), which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration; and from the Two Micron All Sky Survey (2MASS), which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This work used as well data from the following catalogues and surveys: ESO’s Visible and Infrared Survey Telescope for Astronomy (VISTA); the Digitized Sky Survey 2 (DSS-2), produced at the Space Telescope Science Institute under US Government grant NAG W-2166; the Galaxy Evolution Explorer GALEX-GR5 catalogues, a NASA mission led by the California Institute of Technology; the Deep Near-infrared southern Sky survey (DENIS); the AKARI/IRC and the AKARI/FIS mid-infrared surveys (ISAS/JAXA); the MSX6C Infrared Point Source Catalog; and the IRAS Point Source Catalog (NASA). We also made use of data and images from the Pan-STARRS1 Surveys (PS1), which have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. We used the following VO-compliant tools: Topcat (Taylor 2005); Aladin, developed at the Centre de Donées astronomiques de Strasbourg (CDS), France; ESASky (Baines et al. 2017), developed at ESA/ESAC, Madrid, Spain, by the ESAC Science Data Centre (ESDC); and VOSA (Bayo et al. 2008), developed under the Spanish Virtual Observatory (https://svo.cab.inta-csic.es) project funded by MCIN/AEI/10.13039/501100011033/ through grant PID2020-112949GB-I00. VOSA has been partially updated by using funding from the European Union’s Horizon 2020 Research and Innovation Programme, under Grant Agreement n° 776403 (EXOPLANETS-A). This research also made use of the VizieR catalogue access tool and the SIMBAD database, both developed and maintained at the CDS, and of the ADS bibliographic services.

References

  1. Andriantsaralaza, M., Ramstedt, S., Vlemmings, W. H. T., & De Beck, E. 2022, A&A, 667, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147 [Google Scholar]
  3. Baines, D., Giordano, F., Racero, E., et al. 2017, PASP, 129, 028001 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bayo, A., Rodrigo, C., Barrado y Navascués, D., et al. 2008, A&A, 492, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Beichman, C. A., Neugebauer, G., Habing, H. J., Clegg, P. E., & Chester, T. J. 1988, in Infrared astronomical satellite (IRAS) catalogs and atlases. volume 1: Explanatory supplement [Google Scholar]
  6. Beuther, H., Walsh, A., Wang, Y., et al. 2019, A&A, 628, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bhattacharya, R., Medina, B. M., Pihlström, Y. M., et al. 2024, ApJ, 969, 109 [Google Scholar]
  8. Bianchi, L., & GALEX Team 2000, Mem. Soc. Astron. It., 71, 1123 [NASA ADS] [Google Scholar]
  9. Blommaert, J. A. D. L., Groenewegen, M. A. T., Justtanont, K., & Decin, L. 2018, MNRAS, 479, 3545 [CrossRef] [Google Scholar]
  10. Cadmus, R. R. 2024, AJ, 167, 200 [Google Scholar]
  11. Chambers, K. C., & Pan-STARRS Team 2017, in American Astronomical Society Meeting Abstracts, 229, 223.03 [Google Scholar]
  12. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources [Google Scholar]
  13. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2013, Explanatory Supplement to the AllWISE Data Release Products [Google Scholar]
  14. Dell’Agli, F., García-Hernández, D. A., Ventura, P., et al. 2015a, MNRAS, 454, 4235 [Google Scholar]
  15. Dell’Agli, F., Ventura, P., Schneider, R., et al. 2015b, MNRAS, 447, 2992 [Google Scholar]
  16. Di Criscienzo, M., Ventura, P., García-Hernández, D. A., et al. 2016, MNRAS, 462, 395 [NASA ADS] [CrossRef] [Google Scholar]
  17. Egan, M. P., Price, S. D., & Kraemer, K. E. 2003, in American Astronomical Society Meeting, 203, 57.08 [Google Scholar]
  18. Engels, D., & Bunzel, F. 2015, A&A, 582, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Engels, D., & López-Martí, B. 2024, in Cosmic Masers: Proper Motion Toward the Next-Generation Large Projects, eds. T. Hirota, H. Imai, K. Menten, & Y. Pihlström, IAU Symposium, 380, 368 [Google Scholar]
  20. Epchtein, N., de Batz, B., Capoani, L., et al. 1997, Messenger, 87, 27 [Google Scholar]
  21. Fabricius, C., Luri, X., Arenou, F., et al. 2021, A&A, 649, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Fitzpatrick, E. L. 1999, PASP, 111, 63 [Google Scholar]
  23. Gaia Collaboration (Eyer, L., et al.) 2019, A&A, 623, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Gavras, P., Rimoldini, L., Nienartowicz, K., et al. 2023, A&A, 674, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Habing, H. J. 1996, A&ARv, 7, 97 [NASA ADS] [CrossRef] [Google Scholar]
  27. Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931 [NASA ADS] [CrossRef] [Google Scholar]
  28. Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Ishihara, D., Kaneda, H., Onaka, T., et al. 2011, A&A, 534, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Iwanek, P., Soszyński, I., Kozłowski, S., et al. 2022, ApJS, 260, 46 [NASA ADS] [CrossRef] [Google Scholar]
  31. Jayasinghe, T., Kochanek, C. S., Stanek, K. Z., et al. 2018, MNRAS, 477, 3145 [Google Scholar]
  32. Jiménez-Esteban, F. M., & Engels, D. 2015, A&A, 579, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Jiménez-Esteban, F. M., Engels, D., Aguado, D. S., González, J. B., & García-Lario, P. 2021, MNRAS, 505, 6051 [CrossRef] [Google Scholar]
  34. Jiménez-Esteban, F. M., Torres, S., Rebassa-Mansergas, A., et al. 2023, MNRAS, 518, 5106 [Google Scholar]
  35. Lallement, R., Vergely, J. L., Babusiaux, C., & Cox, N. L. J. 2022, A&A, 661, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Lasker, B. M., Doggett, J., McLean, B., et al. 1996, in Astronomical Data Analysis Software and Systems V, eds. G. H. Jacoby, & J. Barnes, Astronomical Society of the Pacific Conference Series, 101, 88 [NASA ADS] [Google Scholar]
  37. Lebzelter, T., Trabucchi, M., Mowlavi, N., et al. 2019, A&A, 631, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Lebzelter, T., Mowlavi, N., Lecoeur-Taibi, I., et al. 2023, A&A, 674, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Lewis, B. M., Eder, J., & Terzian, Y. 1990, ApJ, 362, 634 [Google Scholar]
  40. Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 [EDP Sciences] [Google Scholar]
  41. Lucas, P. W., Smith, L. C., Guo, Z., et al. 2024, MNRAS, 528, 1789 [NASA ADS] [CrossRef] [Google Scholar]
  42. McMahon, R. G., Banerji, M., Gonzalez, E., et al. 2013, Messenger, 154, 35 [Google Scholar]
  43. Mowlavi, N., Lecoeur-Taïbi, I., Lebzelter, T., et al. 2018, A&A, 618, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Pastorelli, G., Marigo, P., Girardi, L., et al. 2019, MNRAS, 485, 5666 [Google Scholar]
  45. Pastorelli, G., Marigo, P., Girardi, L., et al. 2020, MNRAS, 498, 3283 [Google Scholar]
  46. Rimoldini, L., Holl, B., Audard, M., et al. 2019, A&A, 625, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Saito, R. K., Hempel, M., Minniti, D., et al. 2012, A&A, 537, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Sargent, B. A., Srinivasan, S., & Meixner, M. 2011, ApJ, 728, 93 [NASA ADS] [CrossRef] [Google Scholar]
  49. Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, Astronomical Society of the Pacific Conference Series, 347, 29 [Google Scholar]
  50. Tosi, S., Dell’Agli, F., Huerta-Martinez, E., & Ventura, P. 2022, Universe, 8, 270 [NASA ADS] [CrossRef] [Google Scholar]
  51. Ventura, P., Karakas, A., Dell’Agli, F., García-Hernández, D. A., & Guzman-Ramirez, L. 2018, MNRAS, 475, 2282 [NASA ADS] [Google Scholar]
  52. Yamamura, I., Makiuti, S., Ikeda, N., et al. 2010, VizieR On-line Data Catalog: II/298 [Google Scholar]

Appendix A: Online catalogue service

Table A.1.

Description of the catalogue.

To help the astronomical community use the Gaia Catalogue of Galactic AGB Stars, we have developed an online catalogue service that can be accessed from a webpage5 or through a Virtual Observatory ConeSearch6. Table A.1 lists the data accessible through this service.

The service implements a simple search interface that permits queries by coordinates and radius as well as by other parameters of interest. The user can also select the maximum number of sources to return (with values from 10 to unlimited).

The result of the query is a HTML table with all the sources found in the catalogue fulfilling the search criteria, which can be downloaded as a VOTable or as a CSV file. Detailed information on the output fields can be obtained placing the mouse over the question mark (‘?’) located close to the name of the column. The results of a query can be easily transferred to other VO applications (e.g. Topcat) thanks to the Simple Application Messaging Protocol (SAMP)7, which allows Virtual Observatory applications to communicate with each other in a seamless and transparent manner.

The service has been designed as a live archive. Its initial version uses the data published in this paper. In following versions, it will be extended with more AGB star samples that will be published in forthcoming papers. Additional information, such as source variability properties obtained from Gaia DR3 light curves, will also be eventually added.

All Tables

Table 1.

Results of the cross-match of detected OH masers in the database of circumstellar OH masers (Engels & López-Martí 2024) with the AllWISE, 2MASS, and Gaia DR3 catalogues.

Table 2.

Gaia sample of OH/IR stars.

Table 3.

Sources with Gaia-DR3 counterparts that were excluded from the Gaia sample of OH/IR stars.

Table A.1.

Description of the catalogue.

All Figures

thumbnail Fig. 1.

Examples of the different types of SEDs found in the Gaia sample. The photometry we used is listed in the insets of the panels, and the solid line connects the data points to guide the eye. From left to right and from top to bottom: Single-peaked SED, double-peaked SED, flat SED, and peculiar SED. According to the literature, the plotted sources are an OH/IR star, a post-AGB star, another OH/IR star, and a YSO, respectively. An SED graph for any source in the final catalogue can be retrieved through the online service (see Appendix A).

In the text
thumbnail Fig. 2.

Objects of the Gaia OH/IR star sample that are included and missed (blue and red histogram, respectively) in the second Gaia catalogue of LPV candidates (Lebzelter et al. 2023).

In the text
thumbnail Fig. 3.

Sky distribution of the Gaia sample of OH/IR stars. The sources with good parallaxes (relative parallax error < 20% and astrometric excess noise < 1.5 mas) are shown in red, and the rest are plotted in grey.

In the text
thumbnail Fig. 4.

Histogram showing the geometric distances from Bailer-Jones et al. (2021) for the OH/IR stars with good parallax measurements.

In the text
thumbnail Fig. 5.

Spatial distribution of the OH/IR stars with good parallax projected onto the Galactic plane. The Sun is located at the centre (0,0), and the direction towards the Galactic centre is upwards along the vertical dashed line.

In the text
thumbnail Fig. 6.

Galactic height distribution of OH/IR sources with geometric distances from Bailer-Jones et al. (2021).

In the text
thumbnail Fig. 7.

Magnitude distributions of the Gaia OH/IR sample (semi-filled histograms) and of the subsample of sources with good astrometry (filled histograms). From left to right and from top to bottom: Gaia G band, 2MASS Ks band, and the WISE W1 and W4 bands. The percentages refer to the fraction of sources with respect to the total Gaia OH/IR sample of 1172 sources.

In the text
thumbnail Fig. 8.

Colour distributions of of the Gaia OH/IR sample (semi-filled histograms) and of the subsample of sources with a good astrometry (filled histograms). From left to right and from top to bottom: Gaia BP–RP and G–RP, 2MASS J − Ks, and WISE W1–W4 colours. The percentages refer to the fraction of sources with respect to the total Gaia OH/IR sample of 1172 sources.

In the text
thumbnail Fig. 9.

Left panel: Distribution of bolometric fluxes (left) of the Gaia OH/IR sample. The grey and red histograms correspond to the whole sample and to stars with good parallaxes (good ω), respectively. Right panel: Luminosity distribution for the subsample with good parallaxes.

In the text
thumbnail Fig. 10.

Gaia DR3 HR diagrams using the extinction-corrected BP–RP and G–RP colours (left and right panel, respectively) for the subsample of sources with good parallaxes (red circles), compared to the location of the 100 pc stellar population (in grey). In each diagram, only sources with good photometry according to the criteria outlined in Sect. 5.2 are displayed.

In the text
thumbnail Fig. 11.

Bolometric corrections BC(G) = Mbol − Gabs for the subsample of sources from the Gaia OH/IR stars with good parallaxes. The outlier in the diagram is IRAS 18319-1001, as discussed in the text. The solid curve shows the fit to the data as given in Sect. 5.4, and the broken curve is the bolometric correction as a function of BP–RP colour proposed by Lebzelter et al. (2019) without interstellar extinction.

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.