Study of extremely reddened AGB stars in the Galactic bulge
Key Words.:Stars: AGB and post-AGB – Stars: circumstellar matter – Stars: variable: general – Stars: evolution – Infrared: stars – Galaxy: bulge
Context:Extremely reddened asymptotic giant branch stars (AGB) lose mass at high rates of 10 M yr. This is the very last stage of AGB evolution, in which stars in the mass range 2.0 – 4.0 M (for solar metallicity) should have been converted to C stars already. The extremely reddened AGB stars in the Galactic bulge are however predominantly O-rich, implying that they might be either low-mass stars or stars at the upper end of the AGB mass range.
Aims:Our goal is to determine the mass range of the most reddened AGB stars in the Galactic bulge.
Methods:Using Virtual Observatory tools, we constructed spectral energy distributions of a sample of 37 evolved stars in the Galactic bulge with extremely red IRAS colours. We fitted DUSTY models to the observational data to infer the bolometric fluxes. Applying individual corrections for interstellar extinction and adopting a common distance, we determined luminosities and mass-loss rates, and inferred the progenitor mass range from comparisons with AGB evolutionary models.
Results:The observed spectral energy distributions are consistent with a classification as reddened AGB stars, except for two stars, which are proto-planetary nebula candidates. For the AGB stars, we found luminosities in the range 3000 – 30,000 L and mass-loss rates 10 – M yr. The corresponding mass range is 1.1 – 6.0 M assuming solar metallicity.
Conclusions:Contrary to the predictions of the evolutionary models, the luminosity distribution is continuous, with many O-rich AGB stars in the mass range in which they should have been converted into C stars already. We suspect that bulge AGB stars have higher than solar metallicity and therefore may avoid the conversion to C-rich. The presence of low-mass stars in the sample shows that their termination of the AGB evolution also occurs during a final phase of very high mass-loss rate, leading to optically thick circumstellar shells.
Low- and intermediate-mass stars pass along the asymptotic giant branch (AGB) towards the end of their evolution. Stellar evolution models predict that stars approaching the tip of the AGB experience thermal pulses. In this so-called ‘thermal pulses AGB’ (TP-AGB) phase, stars steadily increase in luminosities and mass-loss rates. This increase is interrupted on timescales of several ten thousand years by short-term variations due to the pulses. The heavy mass loss leads to the formation of circumstellar envelopes (CSE) of gas and dust. If the mass-loss rate surpasses 10 M yr, the dust shell eventually becomes opaque to visible light (Habing, 1996). In addition, the stars may change their photospheric chemistry from O- to C-rich as a consequence of the so-called third dredge-up following a thermal pulse. The range of initial masses, in which the star will end as C-rich, depends on the metallicity. This range, defined on the main sequence, is 2 M 4 M for solar metallicities (Marigo & Girardi, 2007). At lower metallicities, as for example in the Magellanic Clouds, this mass range is extended to lower masses, significant increasing the fraction of C stars in the AGB population at the end of the TP-AGB phase (Blum et al., 2006). The stars ending the TP-AGB O-rich must have developed either from low-mass stars with main-sequence masses M 2 M (Bertelli et al., 2008), for which the TP-AGB phase is terminated before the conversion is completed, or from intermediate-mass stars with M 4 M, in which hot bottom burning (HBB) is active at the base of the convective envelope, preventing the formation of carbon.
Low-mass stars are expected to experience only moderate mass-loss rates during most of their AGB lifetime, leading to the formation of optically thin CSE. Low-mass AGB stars have been discovered predominantly in the optical and near-infrared wavelength range. Classical Mira variables belong to these kind of stars. AGB stars descending from progenitors with larger main-sequence masses are expected to be more luminous and to have higher mass-loss rates. They probably spend a larger fraction of their TP-AGB life in an optically obscured state (Mouhcine, 2002), which constrains their detection to the infrared or radio wavelength range. OH/IR stars, originally discovered by their radio OH maser emission and later identified in the infrared are an example of these kind of stars. OH maser emission is however also widespread among Mira variables, so that the term ‘OH/IR star’ encompasses low- and intermediate-mass stars.
If low-mass stars develop high mass-loss rates at the end of their evolution on the TP-AGB, they will form optically thick CSEs at lower mass-loss rates than intermediate-mass stars because they have more compact shells with smaller outflow velocities. Also, because of their frequency they should outnumber the intermediate-mass stars even among the obscured stars. However, in the Galactic disk, Chen et al. (2001) and Likkel (1989) found that the reddest (e.g. with highest mass-loss rates) OH/IR stars are dominated by progenies of main-sequence intermediate-mass stars. A dearth of low-mass stars among infrared bright sources in clusters of the Magellanic Clouds was also noted by van Loon et al. (2005). Therefore, low-mass stars may never evolve into the high mass-loss regime observed for the very red OH/IR stars, or it might be a very brief phase because the envelope mass above the core is rapidly depleted.
These stars should have relatively low luminosities, L 8000 L for M 2 M (Bertelli et al., 2008). Because of the uncertain distances in the Galactic disk and the strong variability, mean luminosities for AGB stars are difficult to determine and objects with relatively low luminosities are difficult to single out. A way to overcome this pitfall is to study AGB stars in the Galactic bulge, where a common distance can be used (van der Veen & Habing, 1990). AGB stars are luminous enough to be detectable there, especially in the mid-infrared. During their heavily obscured stage, they emit most of their radiation in the 10 – 30 m range, which is less affected by the significant interstellar extinction towards many parts of the Galactic bulge.
The Galactic bulge is the central part of the Milky Way galaxy, consisting of a dominant old population of stars (t 8 Gyr; Zoccali et al. 2003; Vanhollebeke et al. 2009) and a smaller population of younger stars with no well constrained ages of 200 Myr – 7 Gyr (van Loon et al., 2003) or 1 – 3 Gyr (Groenewegen & Blommaert, 2005), to which the AGB stars in the Galactic bulge belong. In Baade’s window two metallicity populations have been found, of which the metal-poor population is associated with the old stellar population and the metal-rich population is consistent with a younger stellar population supporting a bar (Hill et al., 2011). A metallicity gradient was noted along the minor axis of the bulge (Zoccali et al., 2008) with the metallicity decreasing towards higher latitudes. No luminous C-rich AGB stars are known in the Galactic bulge (Blanco & Terndrup, 1989), suggesting that the younger stellar population in the Galactic bulge may be restricted to low-mass stars with AGB evolutionary timescales too short to experience enough dredge-ups to convert them into C stars (Buell, 2013). Then, the Galactic bulge OH/IR stars should also exclusively be descendants from low-mass stars, including the extremely red stars. This was implied by the findings of van der Veen & Habing (1990), who derived low progenitor masses of 1.0 – 1.4 M for a sample of these red sources from the IRAS survey.
Although subsequent studies of AGB stars in the Galactic bulge also concluded that their samples originate from low-mass stars, luminosity distributions are not necessarily compatible with an exclusive origin from low-mass stars. For example, Glass et al. (2001) studied the variability properties of long-period, large-amplitude variable stars near the Galactic centre, and found a period-luminosity (P-L) relation that is compatible with the existence of a young population with initial masses in the range 2.5 – 3.0 M. In a later study on Mira variables in the optical gravitational lensing experiment (OGLE) bulge fields, Groenewegen & Blommaert (2005) found a different P-L relation, which they explained with a population with low initial masses (1.5 – 2.0 M). However, possibly more massive stars with P 700 days, such as obscured OH/IR stars, might have been missed in the I-band of the OGLE survey. Blommaert et al. (2006) studied a small sample of Galactic bulge AGB stars with near-infrared spectroscopy, and found a luminosity range from 1700 to 7700 L with a mean of 4100 L, as expected for a population of AGB stars with low initial masses. Again, their study was restricted to objects with relatively low mass-loss rates, and may therefore have excluded more massive stars. This can be inferred from Ojha et al. (2007), who focussed their study on high mass-loss AGB stars, and found a wider luminosity distribution with a peak at 8000 L and an appreciable number of sources with a luminosity exceeding 20,000 L.
A number of studies have been specifically focussed on OH/IR stars. Blommaert et al. (1998) and Wood et al. (1998) found that the OH/IR stars of the Galactic centre do not fit the classical Mira P-L relations (Feast et al., 1989; Whitelock et al., 1991). They systematically have lower luminosity and higher expansion velocity for a given period. OH/IR stars with periods as long as 1000 days, expected to have higher initial masses, were found to have low luminosities corresponding to Miras with masses below 1.5 M. These findings were explained by Wood et al. (1998) with a higher metallicity in the Galactic centre compared to the rest of the bulge. This fits into the general view that at low latitudes the metal-rich stellar population of the disk merges with the generally more metal-poor population of the bulge prevailing at larger latitudes. The origin from low-mass stars was confirmed by Ortiz et al. (2002) for other OH/IR stars observed by the ISOGAL survey. They concluded that the small number of objects more luminous than M = –5.5 (L = 12,000 L) indicates that the initial masses of stars in the sample are rarely higher than 3 M, but instead the masses are about 1 – 2 M, on average.
All AGB star samples studied in the Galactic bulge were dominated by objects with blue or moderate reddened colours. This may have biased the studies towards low-mass objects, which would explain the scarcity of high-luminosity sources. In this paper, we study a sample of AGB stars of the Galactic bulge with extremely reddened colours. These stars are expected to be at the very end of the TP-AGB evolution. If low-mass stars evolve into the high mass-loss regime, the sample should include objects with rather low luminosities. If not, they must have originated from intermediate mass stars with correspondingly higher luminosities.
Currently, with new infrared surveys available, the spectral energy distributions (SED) of heavily obscured objects are covered with unprecedented coverage. Together with theoretical models of the CSEs of AGB stars, reliable bolometric fluxes can be derived without invoking uncertain bolometric corrections. Finally, it is also possible to obtain individual extinctions for each source leading to more reliable luminosities.
The sample of extremely reddened AGB stars is presented in Sect. 2. Section 3 describes the construction of the spectral energy distributions (SED), their fitting with DUSTY models, and the estimation of luminosities and mass-loss rates. The nature of the sources in the sample is discussed in Sect. 4, and the conclusions are presented in Sect. 5.
2 The sample
Our sample of AGB stars of the Galactic bulge (hereafter: GB sample) was selected from the sample of sources in Jiménez-Esteban et al. (2006). They studied around 100 IRAS sources distributed throughout the Milky Way with very red IRAS colours ( 0.75 mag), located in a characteristic region of the IRAS two-colour diagram along the sequence of colours predicted for O-rich AGB stars with increasing mass loss (Bedijn, 1987), and having a high (VAR 50) IRAS variability index (Beichman et al., 1988). The high value of this index implies that all sources should be large-amplitude variables. By comparison with the colours of the ‘Arecibo sample’ of O-rich AGB stars (Jiménez-Esteban et al., 2005), they concluded that their sample of extremely reddened IRAS sources was also mainly made up of O-rich AGB stars highly obscured by their CSEs. To define the GB sample presented here, we selected the extremely reddened IRAS sources from Jiménez-Esteban et al. (2006) with angular distances 10 from the Galactic centre.
Of the 39 sources selected in this region of the sky, two (IRAS 17317–3331 and IRAS 17411–3154) do not belong to the bulge. They are known to be foreground sources, located at a distance of 3.2 and 1.2 kpc, respectively (Cohen et al., 2005; Yuasa et al., 1999). The remaining 37 sources, presented in Table 1, define our sample of extremely reddened AGB stars in the Galactic bulge.
Although the sample is small, it is probably rather complete within its selection boundaries. Its flux density distribution at 25 m peaks at 7 Jy, while the weakest source has 4 Jy (Table 2). AGB stars in the bulge matching the colour selection criterion 0.75 mag and 7 Jy are probably all contained in the sample, except those that are missing because of incompleteness of the IRAS catalogue. Surpassing the flux limit requires a minimum luminosity, which depends on the optical thickness of the CSE. This minimum luminosity is the lowest for sources in which the bulk of the emission is emitted in the 15–35 m wavelength range. This minimum luminosity is also not biased by extinction, which is low at these wavelengths. As we show below, the minimum luminosity required to surpass the flux limit is low enough for high optical thickness, so that low-mass stars will also be present in the sample, if they develop the corresponding high mass-loss rates.
In accordance with the two foreground disk sources identified, we also expect contamination of the sample by background disk sources. The number of background disk sources can be roughly estimated as follows. Considering the foreground disk sources as “typical” for contaminating disk sources and using 4 Jy as the flux limit, all such disk sources within the solar orbit around the centre of the Milky Way should have been detected by IRAS. Assuming the same disk source density in front and behind the Galactic bulge, the number of foreground and background disk sources is coarsely proportional to the surface area of the Galactic disk segments in front and behind the bulge. Using a distance to the Galactic centre of 8.0 ( 0.5) kpc (Reid, 1993), the adopted 10 angular size of the Galactic bulge corresponds to a radius of 1.4 kpc. Based on this consideration, we expect 4 times more background sources than foreground sources. This yields about 8 sources out of 37, which could be background sources and whose luminosities would be underestimated in the following.
We searched in the literature for information on the predominant chemistry in the CSE of our objects. Six sources have IRAS Low Resolution Spectra (LRS) spectral classifications by Kwok et al. (1997). IRAS 17292–2727, IRAS 17316–3523, IRAS 17367–2722, IRAS 17418–2713, and IRAS 17495–2534 were classified as LRS class A, meaning that their spectra present the 9.7m silicate dust feature in absorption. IRAS 18092–2347 was classified as LRS class I, meaning that the LRS was noisy or incomplete. For the rest of the sources, we searched for LRS spectra in the Calgary database111http://www.iras.ucalgary.ca/saldatabase/index.html. Only the LRS of IRAS 17367–3633 was found and visually classified. It shows a clear 9.7m silicate dust feature in absorption. We also visually inspected the LRS of IRAS 18092–2347, which also showed the silicate dust feature in absorption. In addition, Golriz et al. (2014) have recently observed three of our sources that lack LRS classifications (IRAS 17251–2821, IRAS 17276–2846, and IRAS 17323–2424) with the Spitzer-Infrared Spectrograph. All three Spitzer spectra show the 9.7m silicate in absorption. We classified these five sources as LRS class A.
|IRAS||WISE coordinates||IRAS||Radio data||v||A||Comment|
|name||RA(J2000)||Dec(J2000)||Gal Long ()||Gal Lat ()||-||-||LRS||Var||OH||HO||SiO||[km s]||mag|
IRAS colours as defined in Jiménez-Esteban et al. (2005); IRAS Low Resolution Spectra; Classification by this work (see Sections 2 and 3.1); IRAS variability index. Radio data. Y: detection, N: no detection, Sp: single peak; Expansion velocity from the OH maser velocity range, or from CO data for IRAS 17495–2534; SED peculiarities (see Section 3.1). DP: double peaked, NIR-Exc: near-infrared excess, FIR-Exc: far-infrared excess, PPN: proto-planetary nebula candidate.
In addition, all sources of our sample except IRAS 17128–3528 were observed for OH maser emission (te Lintel Hekkert et al., 1991; David et al., 1993; Sevenster et al., 1997). Of these, eight were not detected, two (IRAS 17350–2413 and IRAS 18019–3121) presented single-peak OH emission, and the rest showed the typical double-peak spectrum commonly shown by OH/IR stars. Expansion velocities of the CSEs were derived for our objects from the velocity separations between the two OH peaks. For IRAS 17495–2534, with no OH maser emission detected, the expansion velocity was obtained from CO data (Loup et al., 1993). Four sources (IRAS 17207–3632, IRAS 17392–3020, IRAS 17418–2713, and IRAS 18092–2347) were observed for HO maser emission (Deacon et al., 2007; Suárez et al., 2007). Only IRAS 17392–3020 was detected by Deacon et al. (2007) in April 2003, although it was not detected two years later, in September 2005, by Suárez et al. (2007). Finally, two sources (IRAS 17292–2727 and IRAS 17341–3529) were detected as SiO masers (Nyman et al., 1998; Nakashima & Deguchi, 2003). Evidence for the O-rich nature of 30 sources is therefore provided by either detection of a maser and/or the presence of the silicate dust feature.
In Table 1 we present the observational parameters for the sources in our sample and summarise the above information. The IRAS object name is listed in the first column. Columns 2 to 5 give the equatorial and Galactic WISE coordinates (see Sect 3.1). Columns 6 and 7 give the IRAS colours. The IRAS LRS spectral classification and the IRAS variability index are given in the eight and ninth columns, respectively. In columns 10 to 12, we summarise the detection status for OH, HO, and SiO maser emission, respectively. The expansion velocity, if available, is listed in column 13, and the estimated individual extinction in the line of sight used to de-redden the observational photometry (see Sect. 3.2) is in column 14. Some comments on individual SEDs, discussed below, are given in the last column.
3.1 Spectral energy distributions
To create the SEDs of the sources of our sample we collected, when available, photometric information from astronomical catalogues. To do that, we took advantage of the Virtual Observatory222http://www.ivoa.net (VO). In particular, we made use of the Multiple Cone Search utility of Topcat333http://www.star.bris.ac.uk/mbt/topcat/, an interactive graphical viewer and editor for tabular data that allows the user to examine, analyse, combine, and edit astronomical tables. With Topcat, we cross-matched our sample with several available astronomical catalogues within the VO, using different search radii depending on the particular catalogue.
This way, we populated the SEDs from the near- to the far-infrared, wavelength range where these objects radiate the majority of their energy. We collected data from the following catalogues: Two Micron All Sky Survey (2MASS; Cutri et al. 2003) and VISTA variables in the Vía Láctea catalogue (VVV-DR1; Saito et al. 2012) at 1.25, 1.65, and 2.2 m (filter J, H, and K, respectively); GLIMPSE (Churchwell et al., 2009) at 3.56, 4.51, 5.76, and 7.96 m (filters I1, I2, I3, and I4, respectively); WISE catalogue (Wright et al., 2010) at 3.4, 4.6, 11.6, and 22.1 m (filters W1, W2, W3, and W4, respectively); MSX6C Infrared Point Source Catalogue (Egan et al., 2003) at 8.28, 12.13, 14.65, and 21.3 m (filters A, C, D, and E, respectively); AKARI Point Source Catalogue (Ishihara et al., 2010; Yamamura et al., 2010) at 9, 18, 65, 90, 140, and 160 m (filters S09, S18, S65, S90, S140, and S160, respectively); and IRAS Point Sources Catalogue (Beichman et al., 1988) at 12, 25, 60 and 100 m (filters 12, 25, 60, and 100, respectively). Since WISE data were not available through the VO, they were downloaded directly from the NASA/IPAC Infrared Science Archive.
We plotted the SED for each source and visually inspected them to check the consistency of the photometric data. When deviating photometric data were found, we used Aladin444http://aladin.u-strasbg.fr/ (Bonnarel et al., 2000), another VO tool, to visualise the astronomical images of the field and to check the reliability of the selected counterpart. Aladin allows users to visualise and analyse digitised astronomical images, and superimpose entries from astronomical catalogues or databases available from the VO services.
We found few inconsistencies, basically because of chance coincidences with field stars within our search radii. In all these cases the offending counterparts were removed. In the case of missing counterparts both in 2MASS and VISTA catalogues, we used a nearby (30″) faint star to define an upper flux limit, preferably in the VISTA catalogue.
The photometric data are presented in Table 2, and the SEDs are plotted in Fig. 7. This figure shows the observational data after correcting for interstellar extinction and the DUSTY model that best fits the observations (see Sect. 3.2). Most of the sources present the typical SEDs of an AGB star heavily obscured by an optically thick CSE, with the maximum of their emission in the mid-IR. Note that most of the sources were not detected in the near-IR, or only barely in the K band.
Two sources (IRAS 17276–2846 and IRAS 17504–3312) show a peculiar double-peaked SED (Fig. 7). In the case of IRAS 17276–2846, the photospheric emission extends into the visual so that we added VVVDR1 data for filters Z (0.88 m) and Y (1.02 m) to the SED. These SEDs appear as soon as the mass-loss rates have greatly decreased and the obscuration has substantially diminished. The bluer peak at 2 m is then attributed to the reddened stellar photosphere, and the redder peak corresponds to the emission at longer wavelengths by the cool circumstellar dust (e.g. see Fig. 9 in Sánchez Contreras et al., 2008). This kind of drop in mass-loss rate can be expected in the aftermath of a thermal pulse or when the stellar envelope has nearly been lost at the end of AGB evolution. Stars with this characteristic SED are generally considered to be proto-planetary nebula (PPN) (Oudmaijer, 1996), and therefore we assume that these stars have left the AGB. Since post-AGB evolution is outside the scope of the present paper, the two PPN candidates are not included in the following analyses.
|Flux in Jy|
Data without errors are upper limits; 2MASS data; VVV-DR1 data; No confirmed detection; VVV-DR1 fluxes in Jy 10: F=7612 and F=30174.
|Filter ( (m))||D||S18||E||W4||25||60||S65||S90||100||S140||S160|
|Flux in Jy|
Data without errors are upper limits; 2MASS data; VVV-DR1 data; No confirmed detection; VVV-DR1 fluxes in Jy 10: F=7612 and F=30174.
3.2 Observed bolometric flux
The bulge sample of AGB stars is nowadays covered over a large part of their SED by a multitude of infrared surveys. This allows us to determine the bolometric flux with much higher accuracy compared to earlier times, when (uncertain) bolometric corrections or (unsuitable) black-body fitting had to be applied to cover the unobserved parts of the SED (e.g. van der Veen & Breukers, 1989; Ortiz et al., 2002). However, the photometric data available is not complete for all sources so that a determination of the bolometric flux by integration of an interpolated curve between the available data points will be unsatisfactory. Also bolometric corrections (although small ones) for the spectral ranges outside the infrared are still required. We therefore decided to fit the observed SEDs with models of stars obscured by a dusty radial symmetric CSE and to determine the bolometric fluxes from the models. Because of the large distance of the Galactic bulge, the stars may be observed through large columns of intervening dust, such that extinction may have a noticeable influence on the shape of the SEDs and on the observed bolometric flux.
To de-redden the photometric data, we corrected the flux at 2.2 m ( band) using the extinction maps towards the Galactic bulge given by Gonzalez et al. (2012). The absolute extinctions towards the sources of our sample are 0.13 A 2.45 mag (see Table 1). Two sources (IRAS 16582–3059 and IRAS 17030–3053) are not covered by the maps and we assigned = 0.15 mag to them. About half of the sources are also covered by the extinction maps of Nidever et al. (2012). Their absolute extinctions are 0.14 0.04 mag higher for 1.0 mag than those of Gonzalez et al. (2012). Because the uncertainties given in both papers are 0.1 mag, we considered this difference as not significant. For higher extinctions, the scatter between the maps is larger with = 0.3 mag.
To correct the other photometric bands, we determined the extinction coefficients A/A from the extinction curves given by Chen et al. (2013) up to = 8 m and from Gao et al. (2009) up to 24 m. The extinction coefficients beyond = 24 m were determined by extrapolation using a power law with = –1.7.
The de-reddened SEDs were modelled using the radiation transport code DUSTY (Ivezić et al., 1999). Because the observed photometry taken at different epochs is likely affected by variability, the photometric data are expected to scatter around the mean SED. The relatively sparse coverage of the SEDs with photometry and the variability preclude the determination of accurate models. We therefore decided to use a standard model, in which most of the parameters are fixed, and only the optical depth of the shell was left as variable parameter.
Combining the information from the IRAS LRS classifications and the maser detections, O-rich chemistry can be inferred for 28 objects. Carbon stars are extremely rare in the Galactic bulge (Blanco & Terndrup, 1989). Hence, we decided to assume O-rich chemistry for the remaining seven objects, following Lewis (1992), who found that infrared sources without OH maser detection but with IRAS colours similar to those of OH/IR stars are most probably O-rich variable AGB stars. Consequently, we only used models with O-rich chemistry.
The standard model assumes a central star with an effective temperature T = 2500 K, and a dust condensation temperature T = 1000 K. We used the optical constants for amorphous cold silicates from Ossenkopf et al. (1992) and took the standard MRN (Mathis et al., 1977) dust size distribution with , where is the number density and is the size of the grains. The grain sizes were limited to 0.005 0.25 m. The inner radius of the CSE is determined by DUSTY from the dust condensation temperature assumed. The outer radius of the shell was set to = 100 . Instead of describing the density distribution in the shell by an analytic profile, the option to calculate the wind structure from hydrodynamics was chosen, which allows us to determine mass-loss rates from the models (Ivezic & Elitzur, 1995). Models were calculated for optical depths 0.2 50 and the best model was found by minimising the deviations from the de-reddened flux densities. In several cases two models were fitted representing the star at different variability states. Because of the requirement of a very red IRAS colour, most stars possess a CSE with a high optical depth 10, and only four have 5 10. The bolometric fluxes determined vary by a factor 20 and are in the range = – Watt m.
The results are given in Table 3. It lists the object name in the first column. The best-fit optical depth , and the bolometric flux delivered by DUSTY are given in the two following columns. The last two columns give luminosities and mass-loss rates discussed in the next Section.
The uncertainties of are approximately 0.2 dex, which we attribute mostly to variability. For the cases where models for different variability states are given, the difference is 0.3 dex. We conclude therefore that the peak-to-peak luminosity variation of the large-amplitude variable objects is 2. This number is probably a lower limit, as Suh (2004) determined for two extremely obscured OH/IR stars in the Galactic disk luminosity variations by a factor 3 – 3.6. The bolometric fluxes determined here cannot be individually corrected for variability, but in general we expect them to be close to the mean bolometric fluxes. The deviations of the observed bolometric fluxes from the mean fluxes is at most half of the peak-to-peak variations, and very likely even lower in those cases where average bolometric fluxes are available from fits at different variability states.
The de-reddened SEDs and the models are shown in Fig. 7. The observed SEDs are the result of the emission from the circumstellar dust with high optical depths, which in addition are reddened further by interstellar extinction. The removal of the interstellar extinction makes the SEDs bluer and brighter than the observed SEDs. As an example, the influence of extinction on photometry and model fit is shown in Fig. 1 for IRAS 173923020. We found that the increase in logarithmic bolometric flux surpasses 0.2 dex only if A 1.0 mag, which applies for four sources in the sample. Therefore, the largest contributions to the uncertainties of the determination is in general from the variability and not from interstellar extinction.
For several sources the DUSTY models predict less flux than observed in the near infrared at 3 m and in the far infrared at wavelengths 60 160 m. These sources are marked in Table 1 as ‘NIR-Exc’ and ‘FIR-Exc’. The observed excess emission makes only a minor contribution to the overall luminosity and has been neglected for the luminosity determination.
The near-infrared excess presented by IRAS 17128–3528 and IRAS 18092–2508 (see Fig. 7) is evident from the unexpectedly high brightness in the J, H, and K filters measured by the VVV-DR1 survey. We verified that this is not due to missidentification or contamination of neighbouring objects. These stars may have developed deviations of their dust distribution from spherical symmetry or their circumstellar shell may have been diluted after the end of the strong AGB mass loss, so that the light of the central star could escape to form the excess. They might currently be starting the post-AGB evolution and may later develop double-peaked SEDs as seen in the PPN candidates IRAS 17276–2846 and IRAS 17504–3312 (cf. Sect. 3.1). Alternatively, they might belong to binary systems, where the excess emission comes from the companion. In any case, additional photometry is required to verify the near-IR excesses for these sources.
The SEDs of IRAS 17292–2727, IRAS 17367–2722, IRAS 17367–3633, and IRAS 17583–3346 (Fig. 7) present a far-infrared excess. This far-infrared excess emission can not be attributed to photometric errors or variability since consistent flux densities have been measured at different times by more than one survey. Such excess emission has been seen before in a number of objects, in which the excess can be followed over a large wavelength interval extending up to the submillimeter and millimeter ranges. To account for the excess flux an additional dust component of larger cool grains is required to model the SEDs of these stars (Sánchez Contreras et al., 2007).
|17251–2821||23.3 – 28.2||11.44 – 11.80||7.2 – 3.1||1.7 – 3.4|
|17341–3529||5.1 – 6.2||11.05 – 11.30||17.7 – 9.9||2.1 – 3.4|
|17350–2413||9.0 – 10.9||11.50 – 11.95||6.3 – 2.2||0.7 – 1.7|
|17367–3633||10.9 – 13.2||10.80 – 11.15||31.4 – 14.0||4.7 – 9.3|
|17418–2713||50.0 – 50.0||10.67 – 11.02||42.4 – 18.9||18.0 – 40.4|
|17495–2534||15.9 – 28.2||10.70 – 11.00||39.6 – 19.8||11.1 – 14.0|
|17521–2938||41.3 – 50.0||11.60 – 11.88||4.9 – 2.6||2.3 – 3.7|
|18040–2726||34.2 – 34.2||11.35 – 11.65||8.9 – 4.4||2.8 – 5.7|
|18040–2953||19.3 – 23.3||11.30 – 11.65||9.9 – 4.4||2.4 – 4.6|
|18091–2437||15.9 – 19.3||11.50 – 11.85||6.3 – 2.8||1.2 – 2.3|
in Watt m; Luminosity estimated using the distance of 8.0 kpc to the Galactic centre. v = 14 km s was assumed. Otherwise v is derived from either OH or CO radio data.
We estimated the luminosity from the bolometric fluxes , assuming a common distance to all of the sources in the sample equivalent to the distance to the Galactic centre (8.0 0.5 kpc). The luminosity values obtained are tabulated in the fourth column of Table 3. In case of two fits, representing a brighter and a fainter phase, both luminosities are given in the table.
To quantify the uncertainties of the luminosity determination, we take the uncertainties of the distances and of the bolometric flux determination into account. Taking the Galactic centre distance uncertainty and the estimated size of the bulge ( 2.8 kpc) into account, the distance dependent uncertainty of the luminosity is 35%. This is of the same order as the uncertainty of the bolometric flux due to variability and extinction (see Sect. 3.2). Therefore, luminosities of individual objects are only accurate to within a factor of 2.
Figure 2 shows a histogram of the luminosities in our Galactic bulge sample. In the case of two model fits, we used the average luminosity. The range of luminosities found goes from 3000 to 30,000 L, peaking at about 4500 L.
As mentioned in Section 2, the minimum luminosity required by AGB stars to be present in the current sample (e.g. 7 Jy) depends on the optical depth of the shells. Stars at a distance of 8.0 kpc and a SED represented by one of the DUSTY models would need for 20 merely a minimum luminosity of 2000 – 3000 L. For lower optical depths 5 20 this minimum luminosity rises to 8500 L. Luminosities of the order of a few thousand solar luminosities in Table 3 are therefore found only for high optical depths . Stars with CSEs of lower optical depths could enter the sample only with larger luminosities.
3.4 Mass-loss rates
The calculations of mass-loss rates by DUSTY assume that the mass loss is radiatively driven. The model rates and terminal outflow velocities are given for stars with a luminosity L, fixed gas-to-dust ratio and dust grain bulk density g cm. Using the observed luminosities (Table 3) and the expansion velocities v listed in Table 1 as proxies for terminal outflow velocities, the mass-loss rates were adjusted according to the scaling relations given in the DUSTY manual (Ivezić et al., 1999). For stars without OH maser or CO observations v = 14 km s was adopted. The mass-loss rates obtained, are in the range – 3 M yr. The rates for individual stars are listed in Table 3. Because of uncertainties inherent to the DUSTY code of 30% and the uncertainties of the luminosities, the error of individual mass-loss rates is also a factor 2.
The scaling relations adjust the product to accommodate the observed expansion velocities. With the DUSTY values for the gas-to-dust ratio and dust grain bulk density , the model outflow velocities are generally only about half as large as the observed velocities. To increase the model velocities, the product needs to be decreased, as v. As g cm is already considered as a conservatively low density for astronomical silicates (Draine & Lee, 1984), the gas-to-dust ratio is probably smaller. For the 27 sources with observed expansion velocities (cf. Table 1), is required to bring the model and observed expansion velocities into agreement. This value corresponds to the lower end of the range determined for OH/IR stars in the Galactic disk by Justtanont et al. (2006).
The mass-loss rates in general are rather high. To check their reliability, we compared them with those of Groenewegen (2006), who calculated mid- and far-infrared colours for mass losing AGB and post-AGB stars based on own models. From the different dust compositions Groenewegen (2006) employed, we used his pure silicate models for comparison. The range of mid-IR colours he considered covers only about 40% (the blue part) of the colours of the GB sample, but among these the mass-loss rates are in reasonable agreement. The mean ratio between Groenewegen’s and our mass-loss rates is 0.9 0.4.
Figure 3 plots the mass-loss rates against the luminosities of the sample. Because of scaling relations the two quantities are not independent from each other, and in the classical case a relation (van Loon et al., 1999; Elitzur & Ivezić, 2001) is expected. The scatter of the mass-loss rates for a given luminosity is caused mainly by the scatter of optical depths required to model the SEDs and to a lesser degree by the variations of the expansion velocities. In general, the luminosities and mass-loss rates cover the range between the classical limit in the -plane and the empirical limit found by van Loon et al. (1999) in the Large Magellanic Cloud (LMC). The classical limit is determined by the assumption that the photons transfer momentum to the dust by a single scattering event. This limit is superseded however at high optical depths, where photons experience multiple scattering events, before escaping from the dust shell. Compliance with this empirical limit was found also for the AGB stars in the central regions of M33, which have a near solar metallicity comparable with the Galactic bulge (Javadi et al., 2013).
Few stars exceed the empirical limit. For example, IRAS 174182713 has the largest mass-loss rate in the sample of M yr, and exceeds the empirical limit by a factor . The very steep rise of the SED in the near-infrared requires a high optical depth for the model dust shell, which leads to a high mass-loss rate in combination with the inferred high luminosity of 30,000 L. However, given the uncertainties for the luminosities and mass-loss rates, a systematic difference of the upper limit of the -relation of AGB stars in the bulge and the LMC cannot be inferred.
4 The nature of the sample
4.1 The distribution of luminosities and mass-loss rates
Our sample is made up of the most reddened AGB stars of the Galactic bulge known so far. We found a large range of luminosities with a peak in the luminosity distribution at L 4500 L consistent with low-mass progenitors, and with a tail extending to high luminosities.
This result is in good agreement with other studies of less reddened AGB stars in the Galactic bulge. The luminosity distribution of Blommaert et al. (1998) for a sample of OH/IR stars located close to the Galactic centre spans a similar range and has a median value of 4600 L. The median luminosity of Galactic centre OH/IR stars monitored by Wood et al. (1998) for variability is 6400 L, and their range of luminosities extends to L 30,000 L consistent with an upper limit for the mass range M. Both studies relied on measurements in the near-infrared (K, L-band), making the colours of their stars bluer than in our sample.
Intermediate in colours between the Galactic centre OH/IR stars of the just above mentioned studies and the colours of our sample are the OH/IR stars observed by the ISOGAL survey (Ortiz et al., 2002) and of the MSX-selected selected sample of AGB stars (Ojha et al., 2007). Their luminosity distributions peak around 8000 L, and objects with higher luminosities are found in both samples. The red IRAS sources studied by van der Veen & Habing (1990) are very similar to our sample. Their sources cover the colour range 0.0 1.3, which overlaps in its red part with the colour range ( 0.75 mag) of the present sample. Most of their sources are bluer than ours, with only seven sources in common. They find the luminosities to be strongly peaking at 5000–5500 L and observed a tail of high luminosity sources extending to well above 20,000 L.
We therefore conclude that the luminosity distributions of the different samples of AGB stars observed in the Galactic bulge are consistent with our result in showing that the majority of the stars have luminosities well below 10,000 L, but that stars with luminosities up to the AGB limit are present as well. There is no evidence for differences in the luminosity distribution with colour, which would indicate the presence of mass segregation.
Mass-loss rates in the range from to M yr for Galactic bulge AGB stars were derived by Ojha et al. (2007), using the models of Groenewegen (2006). There are nine sources in common with our sample. The deviations in the mass-loss rates are a factor 0.6 – 1.5, except for IRAS 174182713, in which the rates differ by a factor of 10. Ojha et al. (2007) underestimate the rate in this case, because of their adoption of a mean luminosity for all objects to scale the model mass-loss rates. IRAS 174182713 is actually a factor 4 more luminous, which explains at least part of the difference in the mass-loss rates.
4.2 Predictions from AGB evolutionary models
High mass-loss rates of the order M yr or higher (dubbed superwind hereafter) imply that the stars lose most of their mass in a relatively short time. Actually the mass-loss rate determines the stellar lifetime on the AGB. According to the models, during the evolution on the TP-AGB, lasting 0.25 – 2.2 million years, luminosities and mass-loss rates steadily increase, until the last 60,000 – 120,000 years when the stars enter the superwind phase (Vassiliadis & Wood, 1993). These time ranges vary slightly depending on the metallicity assumed and how the mass-loss process is implemented in model calculations (Blöcker, 1995; Marigo & Girardi, 2007; Weiss & Ferguson, 2009). The stars of our sample have mass-loss rates in the range from ( 1 – 30) M yr, characteristic of the superwind phase, and we assume that they will leave the AGB within the next years. Hence, a comparison with stellar evolutionary models has to assign these stars to the low-temperature, high-luminosity ends of the evolutionary tracks on the AGB for low- and intermediate mass stars.
The final phases of these tracks on the AGB are characterised by similar effective temperatures of 2500 – 3000 K, and models for different zero-age main-sequence (ZAMS) masses are then basically distinguished by the final luminosity reached and in part by their chemical composition (Marigo & Girardi, 2007; Weiss & Ferguson, 2009). While Vassiliadis & Wood (1993) and Blöcker (1995) do not discuss the conversion of O-rich to C-rich chemistry in their AGB models due to the third dredge-up, Marigo & Girardi (2007) predict that stars in the mass range 2.0 – 4.0 M with solar metallicity (Z = 0.019), or 1.5 – 4.0 M with typical metallicity for the LMC, (Z = 0.008), end as C stars. In Weiss & Ferguson (2009), the corresponding mass range for solar metallicity is 1.8 – 3.0 M and, for Z = 0.008, no model with M 5.0 M ends O-rich. For most stars of our sample (80%), however, the O-rich chemistry is evident from the presence of masers and/or a signature of silicates in the infrared. Hence, if these stars are indeed close to the end of their TP-AGB evolution, they should have evolved either from intermediate-mass stars with M 4.0 M experiencing HBB, or from low-mass stars of at least solar metallicity. This conclusion would be avoided by the suppression of carbon star formation for the rest of the masses, which may occur for higher than solar metallicity or for increased helium content (Karakas, 2014).
Figure 4 compares the locations of the GB sample stars in the Hertzsprung-Russell diagram with the models of Bertelli et al. (2008, 2009), describing the TP-AGB stellar evolution for ZAMS masses 0.7 M 6.0 M, assuming solar composition (Z = 0.017, Y = 0.26)555Available at http://stev.oapd.inaf.it/YZVAR/. The effective temperatures of the stars are not known, but we assume that they are 3000 K. To avoid crowding in the figure, we assigned random effective temperatures to the stars in the range 2500 2800 K. Contrary to the model predictions, the mass range 1.1 M 6.0 M is continuously covered by the sample, including the range in which the stars should end as C-rich. Most of the stars, in which the chemistry of the envelope is unknown, have relatively low luminosity in line with our assumption that they are all oxygen rich.
4.3 The high-luminosity group
The luminosity distribution of the GB sample (Fig. 2) has a steep increase at low luminosities, peaks at 4500 L, and falls off less steeply towards higher luminosities. For the following discussion, we name sources with L 10,000 L the ‘high-luminosity group’ and those with L 7000 L the ‘low-luminosity group’.
In the ‘high-luminosity group’ two sources (IRAS 174182713 and IRAS 17495–2534) have the highest luminosities ( 30,000 10,000 L; Table 3). Apart from the possibility that they could be foreground sources and hence are actually less luminous, their luminosities are compatible with intermediate-mass stars (M 4 – 6 M) experiencing HBB, preventing them from becoming C stars. Both show the absorption feature from amorphous silicates at 10 m (Table 1), confirming their O-rich chemical composition. Furthermore, IRAS 17418–2713 shows crystalline silicate emission features in the infrared (García-Hernández et al., 2007), while IRAS 17495–2534 is the only source known so far showing these features in absorption (Speck et al., 2008). Speck et al. argue that the rarity of such a distinctive feature is due to their emergence only in massive AGB stars, which is consistent with our findings.
The rest of the group, making up about one quarter of the whole GB sample, has luminosities between 10,000 – 25,000 L. These luminosities are predicted for the mass range, where for solar composition, stars end as C-rich. For lower metallicities, this group will even increase, as conversion to C-rich chemistry also occurs for smaller masses. However, masers and silicate features in the infrared are observed in most stars from this group, ruling out a C-rich chemistry. It is unlikely that a C-rich central star coexists with OH maser emission in all these stars, because such a configuration can last only for 1000 years, which is the time needed to replace the O-rich material in the OH maser shell at 5 cm, assuming typical outfow velocities of 15 km s. As the effective temperatures are unknown, these stars in principle could have temperatures 3000 K, assigning them to an earlier evolutionary phase, where the chemical composition is still O-rich, and the luminosities are almost the same (cf. Fig. 4). However, none of the model calculations (Vassiliadis & Wood, 1993; Blöcker, 1995; Marigo & Girardi, 2007; Weiss & Ferguson, 2009) predict the observed high mass-loss rates at less advanced phases of TP-AGB evolution.
The presence of O-rich objects at the tip of the AGB with luminosities, corresponding to masses where the stars should be already converted to C stars, may have two possibility. Either phases of very high mass-rates can occur, before the O C transition is experienced, or for some of the stars from this mass range the transition might be delayed or avoided. The first possibility is unlikely except in the case of higher mass stars, because it is likely that very little mass is left over after a superwind phase to continue the evolution with further thermal pulses, allowing an O C transition before the stars leave the AGB. Therefore, the second possibility may be more rewarding to follow. Weiss & Ferguson (2009) find that models with super-solar metallicities can avoid the conversion to C/O 1 before the end of the AGB simply because they start with smaller C/O ratios, and therefore more thermal pulses are needed to reach C/O 1. Also, Karakas (2014) finds that at a super-solar metallicity Z=0.03, no carbon stars are formed except in a narrow mass range of 3.25–4 M. Therefore, in environments where a mix of metallicities is present, a mix of stars originating from similar ZAMS masses, but with different chemistries, might also be present. It would imply that most of the objects from the GB sample were born in bulge regions with higher than solar metallicities. A similar line of argument was forwarded by Jura et al. (1993) to explain the differences in the period distribution of Mira variables in the solar neighbourhood and the Galactic centre. They suggested that Mira variables have higher metallicity in the Galactic centre and are therefore oxygen rich, while stars of the same mass in the solar neighbourhood are less metal rich and already converted into carbon stars.
The metallicity argument can be tested for consistency in comparison with the spatial distribution of the GB sample with the global metallicity distribution in the Galactic bulge. Fig. 5 shows the metallicity map of Gonzalez et al. (2013), which clearly confirms the vertical metallicity gradient with metal-rich stars dominating the inner bulge in regions closer to the Galactic plane (). The GB sample members are overplotted. They are completely confined to the more metal rich part of the bulge. Although the mean metallicities in the bins of the metallicity map do not reach super-solar metallicities, higher metallicities may occur on smaller spatial scales, where the GB sample members may be located. Therefore, at least this comparison does not contradict the model predictions that the chemistry of intermediate-mass stars at the termination of TP-AGB evolution depends on the metal content of the environment in which they were born.
4.4 The low-luminosity group
Almost half of the sample has luminosities L 7000 L (log L 3.85) and mass-loss rates (1 – 5) M yr. These kinds of combinations of relatively low luminosities and very high mass-loss rates are not predicted by the models either, at least as long the stars are considered to be on the TP-AGB and in between two thermal pulses (the interpulse phase). Assuming that the stars are in the final phase of their AGB evolution (the last couple of interpulse periods) and making the requirement that the luminosities do not surpass 7000 L, these stars must have descended from 1.1 M 1.8 M main-sequence stars (Fig. 4).
In models of TP-AGB evolution, mass-loss rates exceeding M are not assumed to occur before large amplitude pulsations of the stellar envelope have been established. Most members of the GB sample have not been monitored for variability, except the few sources in common with van der Veen & Habing (1990), for which they report large amplitude variability with periods P days. The other sources usually have high IRAS variability index (, see Table 1), indicating that they are also, in the majority, large amplitude variable stars. We therefore exclude the remote possibility that the members of the low-luminosity group are intermediate-mass stars (M ¿ 3.0 M) on the early AGB preceding the TP-AGB, although the observed luminosities (L 7000 L) are achieved by models during this phase (Bertelli et al., 2008).
During the TP-AGB, the mass-loss rates steadily increase, modulated by the intervening thermal pulses, which affect the surface luminosity and radius of the stars. During the interpulse phase, the stellar luminosity reaches its maximum shortly before the next thermal pulse happens and is referred here to as quiescent luminosity. The models of Vassiliadis & Wood (1993) provide the mass-loss rates ( M yr) observed for M 2.0 M during the last couple of interpulse periods, but the correspondent model luminosities are too high. Blöcker (1995) presents models with M = 1, 3, 4, 5, and 7 M. While the 1 M model does not reach the TP-AGB phase, the M = 3 M model has quiescent luminosities in the required range, but does not reach mass-loss rates surpassing M yr. Models for larger masses do reach the required mass-loss rates, but are again too luminous. Marigo & Girardi (2007) and Weiss & Ferguson (2009) also arrive at the same result: the high observed mass-loss rates are not predicted for the mass ranges corresponding to the observed quiescent luminosities. Given the short times the thermal pulses last ( several hundred years, Vassiliadis & Wood 1993), it is unlikely that all these stars are currently experiencing a boost of the mass-loss rate following the brief increases of the luminosity.
A possible explanation for the discrepancies between the observed properties and the model predictions may originate in the assumptions concerning the relation between mass-loss rates and other properties in the TP-AGB models. As no complete mass-loss theory for AGB stars is available, mass-loss rates are empirically coupled to pulsational period. The periods are theoretically determined from relations involving mass and radius, which are provided by the models. While these relations provide resonable values in general, there might be brief phases where the predictions are insufficient. At least for the low luminosity part of our sample, which had main-sequence masses M 2 M, the models predict mass-loss rates that are too low in the brief final evolutionary phase.
4.5 Comparison with the AGB population in the LMC
In the LMC, the AGB population has been studied extensively by the SAGE project (Blum et al., 2006). The chemistry of the population was determined from photometry using basically the J–[3.6] vs. [3.6] colour-magnitude diagram (CMD). In Fig. 6 the boundaries for different chemistries are delineated in the [3.6]–[8.0] vs. [8.0] CMD, which is more appropriate for very obscured stars. While C stars are in general found to be redder than O-rich stars, Blum et al. also define a group of ‘extreme stars’ with even redder colours ( 3.1 mag), for which the chemistry could not be determined. The few known stars of this group, which have been observed spectroscopically, are predominantly C-rich. Along with the model predictions that the lowest effective temperatures and hence the largest mass-loss rates are achieved when the stars become C-rich, the ‘extreme stars’ are uniformly classified as predominantly C-rich (Srinivasan et al., 2011; Riebel et al., 2012). However, a minority population of O-rich stars is not excluded and contamination with red supergiants and young stellar objects is probably also present (Srinivasan et al., 2009; Blum et al., 2014).
We determined synthetic apparant magnitudes in the [3.6] and [8.0] bands of the GLIMPSE photometric system (Churchwell et al., 2009) for all stars of the GB sample using the DUSTY model fits. Comparison with the actually observed magnitudes by GLIMPSE gives an uncertainty of 1 mag. To compare the location of the sample in the CMD of Blum et al. (2006), a distance modulus of 4.0 mag was used to account for the distance difference between the LMC and the Galactic bulge. As shown in Fig. 6, the GB sample stars are located in a sparsely populated region extending the branch populated by the ‘extreme stars’ towards redder colours 2.5 mag. This region of the CMD was also found by Gruendl & Chu (2009) to be populated by ‘evolved stars’. The reddest of them (dubbed EROS for ‘Extremely Red ObjectS’) were spectroscopically observed by Gruendl et al. (2008), showing C-rich chemistry. It follows that while the SAGE ‘extreme stars’ and the even redder ‘evolved stars’ of Gruendl & Chu (2009) in the LMC are predominantly C-rich, the corresponding population in the Galactic bulge represented by the GB sample seem to be O-rich.
Riebel et al. (2012) determined dust mass-loss rates for the SAGE stars using model SEDs. Almost all ‘extreme stars’ were matched with C-rich models. Most dust mass-loss rates for O-rich stars are in the range – M yr, and for C-rich stars M yr. With a gas-to-dust ratio of 200 their total mass-loss rates are well below the rates achieved during the superwind. The separation of O-rich and C-rich sources in the SAGE CMDs is very likely valid only as long as mass-loss rates are low, the CSEs are optically thin and therefore the stellar temperature is able to influence the infrared colours. As soon as the mass-loss rates increase, and especially during superwind phases, the stars migrate in CMDs towards redder colours, independent of their current chemistry. Gruendl et al. (2008) find a mean luminosity of 7100 L for their EROS sources, and conclude that these stars must have descended from 1.5 – 2.5 M main-sequence stars. The total mass-loss rates derived are in the range of (4 – 23) M yr, much larger than derived for the SAGE AGB stars. The EROS luminosities, main-sequence masses, and mass-loss rates are in perfect agreement with those derived here for the GB sample. The main difference is the chemical composition. As we previously concluded for the GB sample, Gruendl et al. (2008) noted that the mass-loss rates of the EROS sources typically exceed the maximum expected for both O-rich and C-rich AGB stars with luminosities in the observed range.
The reddest bright 25 m sources towards the Galactic bulge with IRAS colours expected for AGB stars were found to span a luminosity range from 3000 to 30,000 L and mass-loss rates between – 3 M yr. For the lower luminosity group with L 7000 L the main-sequence progenitor masses are in the range 1.1 – 1.8 M, assuming solar metallicity. With mass-loss rates (1 – 5) M yr their stellar envelopes are lost in mere 10,000 – 200,000 years and AGB evolution is terminated. One concludes that the heavily obscured AGB stars do not necessarily descend from massive AGB stars (M 4 M) only. At higher luminosities (10,000 – 25,000 L), AGB evolutionary models predict that for metallicities similar to the Sun or less these high-mass loss stars should already be converted to C-rich chemistry. On the contrary, in the bulge these stars are mostly O-rich. The models predict that the conversion to C-rich chemistry can be delayed or inhibited for higher than solar metallicities. Our sources might have been born in regions of the bulge, where several generations of stars have leftover processed matter with higher than solar metallicity. Their spatial distribution favourably coincides with regions of higher mean metallicities in the Galactic bulge.
Acknowledgements.This paper has benefitted from discussions in depth with the referee J. Th. van Loon. This work was partially funded by the Spanish MICINN under the Consolider-Ingenio 2010 Program grant CSD2006-00070 First Science with the GTC666http://www.iac.es/consolider-ingenio-gtc, and through grants AyA2011-24052. This work has been supported by the CoSADIE Coordination Action (FP7, Call INFRA-2012-3.3 Research Infrastructures, project 312559). FJE acknowledges financial support from the ARCHES project (7th Framework of the European Union, n 313146). Support was also granted by Deutsche Forschungsgemeinschaft through grant EN 176/30.
- Bedijn (1987) Bedijn, P. J. 1987, A&A, 186, 136
- Beichman et al. (1988) 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
- Bertelli et al. (2008) Bertelli, G., Girardi, L., Marigo, P., & Nasi, E. 2008, A&A, 484, 815
- Bertelli et al. (2009) Bertelli, G., Nasi, E., Girardi, L., & Marigo, P. 2009, A&A, 508, 355
- Blöcker (1995) Blöcker, T. 1995, A&A, 297, 727
- Blanco & Terndrup (1989) Blanco, V. M. & Terndrup, D. M. 1989, AJ, 98, 843
- Blommaert et al. (2006) Blommaert, J. A. D. L., Groenewegen, M. A. T., Okumura, K., et al. 2006, A&A, 460, 555
- Blommaert et al. (1998) Blommaert, J. A. D. L., van der Veen, W. E. C. J., van Langevelde, H. J., Habing, H. J., & Sjouwerman, L. O. 1998, A&A, 329, 991
- Blum et al. (2006) Blum, R. D., Mould, J. R., Olsen, K. A., et al. 2006, AJ, 132, 2034
- Blum et al. (2014) Blum, R. D., Srinivasan, S., Kemper, F., Ling, B., & Volk, K. 2014, AJ, 148, 86
- Bonnarel et al. (2000) Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33
- Buell (2013) Buell, J. F. 2013, MNRAS, 428, 2577
- Chen et al. (2013) Chen, B. Q., Schultheis, M., Jiang, B. W., et al. 2013, A&A, 550, A42
- Chen et al. (2001) Chen, P. S., Szczerba, R., Kwok, S., & Volk, K. 2001, A&A, 368, 1006
- Churchwell et al. (2009) Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213
- Cohen et al. (2005) Cohen, M., Parker, Q. A., & Chapman, J. 2005, MNRAS, 357, 1189
- Cutri et al. (2003) Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, in Two Micron All Sky Survey (2MASS) All-Sky Catalog of Point Sources
- David et al. (1993) David, P., Lesqueren, A. M., & Sivagnanam, P. 1993, A&A, 277, 453
- Deacon et al. (2007) Deacon, R. M., Chapman, J. M., Green, A. J., & Sevenster, M. N. 2007, ApJ, 658, 1096
- Draine & Lee (1984) Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89
- Egan et al. (2003) Egan, M. P., Price, S. D., & Kraemer, K. E. 2003, in The Midcourse Space Experiment (MSX) Point Source Catalog Version 2.3
- Elitzur & Ivezić (2001) Elitzur, M. & Ivezić, Ž. 2001, MNRAS, 327, 403
- Feast et al. (1989) Feast, M. W., Glass, I. S., Whitelock, P. A., & Catchpole, R. M. 1989, MNRAS, 241, 375
- Gao et al. (2009) Gao, J., Jiang, B. W., & Li, A. 2009, ApJ, 707, 89
- García-Hernández et al. (2007) García-Hernández, D. A., Perea-Calderón, J. V., Bobrowsky, M., & García-Lario, P. 2007, ApJ, 666, L33
- Glass et al. (2001) Glass, I. S., Matsumoto, S., Carter, B. S., & Sekiguchi, K. 2001, MNRAS, 321, 77
- Golriz et al. (2014) Golriz, S. S., Blommaert, J. A. D. L., Vanhollebeke, E., et al. 2014, MNRAS, 443, 3402
- Gonzalez et al. (2013) Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2013, A&A, 552, A110
- Gonzalez et al. (2012) Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2012, A&A, 543, A13
- Groenewegen (2006) Groenewegen, M. A. T. 2006, A&A, 448, 181
- Groenewegen & Blommaert (2005) Groenewegen, M. A. T. & Blommaert, J. A. D. L. 2005, A&A, 443, 143
- Gruendl & Chu (2009) Gruendl, R. A. & Chu, Y.-H. 2009, ApJS, 184, 172
- Gruendl et al. (2008) Gruendl, R. A., Chu, Y.-H., Seale, J. P., et al. 2008, ApJ, 688, L9
- Habing (1996) Habing, H. J. 1996, A&A Rev., 7, 97
- Hill et al. (2011) Hill, V., Lecureur, A., Gómez, A., et al. 2011, A&A, 534, A80
- Ishihara et al. (2010) Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1
- Ivezić et al. (1999) Ivezić, u., Nenkova, M., & Elitzur, M. 1999, User Manual for DUSTY. http://www.pa.uky.edu/moshe/dusty
- Ivezic & Elitzur (1995) Ivezic, Z. & Elitzur, M. 1995, ApJ, 445, 415
- Javadi et al. (2013) Javadi, A., van Loon, J. T., Khosroshahi, H., & Mirtorabi, M. T. 2013, MNRAS, 432, 2824
- Jiménez-Esteban et al. (2005) Jiménez-Esteban, F. M., Agudo-Mérida, L., Engels, D., & García-Lario, P. 2005, A&A, 431, 779
- Jiménez-Esteban et al. (2006) Jiménez-Esteban, F. M., García-Lario, P., Engels, D., & Perea Calderón, J. V. 2006, A&A, 446, 773
- Jura et al. (1993) Jura, M., Yamamoto, A., & Kleinmann, S. G. 1993, ApJ, 413, 298
- Justtanont et al. (2006) Justtanont, K., Olofsson, G., Dijkstra, C., & Meyer, A. W. 2006, A&A, 450, 1051
- Karakas (2014) Karakas, A. I. 2014, MNRAS, 445, 347
- Kwok et al. (1997) Kwok, S., Volk, K., & Bidelman, W. P. 1997, ApJS, 112, 557
- Lewis (1992) Lewis, B. M. 1992, ApJ, 396, 251
- Likkel (1989) Likkel, L. 1989, ApJ, 344, 350
- Loup et al. (1993) Loup, C., Forveille, T., Omont, A., & Paul, J. F. 1993, A&AS, 99, 291
- Marigo & Girardi (2007) Marigo, P. & Girardi, L. 2007, A&A, 469, 239
- Mathis et al. (1977) Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425
- Mouhcine (2002) Mouhcine, M. 2002, A&A, 394, 125
- Nakashima & Deguchi (2003) Nakashima, J.-I. & Deguchi, S. 2003, PASJ, 55, 229
- Nidever et al. (2012) Nidever, D. L., Zasowski, G., & Majewski, S. R. 2012, ApJS, 201, 35
- Nyman et al. (1998) Nyman, L.-Å., Hall, P. J., & Olofsson, H. 1998, A&AS, 127, 185
- Ojha et al. (2007) Ojha, D. K., Tej, A., Schultheis, M., Omont, A., & Schuller, F. 2007, MNRAS, 381, 1219
- Ortiz et al. (2002) Ortiz, R., Blommaert, J. A. D. L., Copet, E., et al. 2002, A&A, 388, 279
- Ossenkopf et al. (1992) Ossenkopf, V., Henning, T., & Mathis, J. S. 1992, A&A, 261, 567
- Oudmaijer (1996) Oudmaijer, R. D. 1996, A&A, 306, 823
- Reid (1993) Reid, M. J. 1993, ARA&A, 31, 345
- Riebel et al. (2012) Riebel, D., Srinivasan, S., Sargent, B., & Meixner, M. 2012, ApJ, 753, 71
- Saito et al. (2012) Saito, R. K., Hempel, M., Minniti, D., et al. 2012, A&A, 537, A107
- Sánchez Contreras et al. (2007) Sánchez Contreras, C., Le Mignant, D., Sahai, R., Gil de Paz, A., & Morris, M. 2007, ApJ, 656, 1150
- Sánchez Contreras et al. (2008) Sánchez Contreras, C., Sahai, R., Gil de Paz, A., & Goodrich, R. 2008, ApJS, 179, 166
- Sevenster et al. (1997) Sevenster, M. N., Chapman, J. M., Habing, H. J., Killeen, N. E. B., & Lindqvist, M. 1997, A&AS, 122, 79
- Speck et al. (2008) Speck, A. K., Whittington, A. G., & Tartar, J. B. 2008, ApJ, 687, L91
- Srinivasan et al. (2009) Srinivasan, S., Meixner, M., Leitherer, C., et al. 2009, AJ, 137, 4810
- Srinivasan et al. (2011) Srinivasan, S., Sargent, B. A., & Meixner, M. 2011, A&A, 532, A54
- Suárez et al. (2007) Suárez, O., Gómez, J. F., & Morata, O. 2007, A&A, 467, 1085
- Suh (2004) Suh, K.-W. 2004, ApJ, 615, 485
- te Lintel Hekkert et al. (1991) te Lintel Hekkert, P., Caswell, J. L., Habing, H. J., et al. 1991, A&AS, 90, 327
- van der Veen & Breukers (1989) van der Veen, W. E. C. J. & Breukers, R. J. L. H. 1989, A&A, 213, 133
- van der Veen & Habing (1990) van der Veen, W. E. C. J. & Habing, H. J. 1990, A&A, 231, 404
- van Loon et al. (2003) van Loon, J. T., Gilmore, G. F., Omont, A., et al. 2003, MNRAS, 338, 857
- van Loon et al. (1999) van Loon, J. T., Groenewegen, M. A. T., de Koter, A., et al. 1999, A&A, 351, 559
- van Loon et al. (2005) van Loon, J. T., Marshall, J. R., & Zijlstra, A. A. 2005, A&A, 442, 597
- Vanhollebeke et al. (2009) Vanhollebeke, E., Groenewegen, M. A. T., & Girardi, L. 2009, A&A, 498, 95
- Vassiliadis & Wood (1993) Vassiliadis, E. & Wood, P. R. 1993, ApJ, 413, 641
- Weiss & Ferguson (2009) Weiss, A. & Ferguson, J. W. 2009, A&A, 508, 1343
- Whitelock et al. (1991) Whitelock, P., Feast, M., & Catchpole, R. 1991, MNRAS, 248, 276
- Wood et al. (1998) Wood, P. R., Habing, H. J., & McGregor, P. J. 1998, A&A, 336, 925
- Wright et al. (2010) Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868
- Yamamura et al. (2010) Yamamura, I., Makiuti, S., Ikeda, N., et al. 2010, VizieR Online Data Catalog, 2298, 0
- Yuasa et al. (1999) Yuasa, M., Unno, W., & Magono, S. 1999, PASJ, 51, 197
- Zoccali et al. (2008) Zoccali, M., Hill, V., Lecureur, A., et al. 2008, A&A, 486, 177
- Zoccali et al. (2003) Zoccali, M., Renzini, A., Ortolani, S., et al. 2003, A&A, 399, 931