We have started an ambitious program to determine if the full diversity of extinction laws is real or if some of it is due to calibration or methodological issues. Here we start by analyzing the information on NIR extinction in a 2MASS stellar sample with good quality photometry and very red colours. We calculate the extinction at 1 $\mu$m, $A_1$, and the power-law exponent, $\alpha$ ($A_\lambda = A_1 \lambda^{-\alpha}$), for the 2MASS stars located in the extinction trajectory in the $H-K$ vs. $J-H$ plane expected for red giants with $A_1 > 5$ mag. We test the validity of the assumption about the nature of those stars, whether a single or multiple values of $\alpha$ are needed, and the spatial variations of the results. Most ($\sim$83%) of those stars can indeed be explained by high-extinction red giants and the rest is composed of extinguished AGB stars (mostly O-rich), blended sources, and smaller numbers of other objects, a contaminant fraction that can be reduced with the help of Gaia} DR2 data. Galactic red giants experience a NIR extinction with $\alpha\sim 2.27$ and an uncertainty of a few hundredths of a magnitude. There is no significant spread in $\alpha$ even though our sample is widely distributed and has a broad range of extinctions. Differences with previous results are ascribed to the treatment of non-linear photometric effects and/or the contaminant correction. Future research should concentrate in finding the correct functional form for the NIR extinction law. In the appendix we detail the treatment of non-linear photometric effects in the 2MASS bands.

This paper describes the characterization of the GALANTE photometric system, a seven intermediate- and narrow-band filter system with a wavelength coverage from 3000 $\r{A}$ to 9000 $\r{A}$ . We describe the photometric system presenting the full sensitivity curve as a product of the filter sensitivity, CCD, telescope mirror, and atmospheric transmission curves, as well as some first- and second-order moments of this sensitivity function. The GALANTE photometric system is composed of four filters from the J-PLUS photometric system, a twelve broad-to-narrow filter system, and three exclusive filters, specifically designed to measure the physical parameters of stars such as effective temperature $T_{\rm eff}$, $\log(g)$, metallicity, colour excess $E(4405-5495)$, and extinction type $R_{5495}$. Two libraries, the Next Generation Spectral Library (NGSL) and the one presented in Ma\'iz Apell\'aniz & Weiler (2018), have been used to determine the transformation equations between the Sloan Digital Sky Survey ($\textit{SDSS}$) $\textit{ugriz}$ photometry and the GALANTE photometric system. We will use this transformation to calibrate the zero points of GALANTE images. To this end, a preliminary photometric calibration of GALANTE has been made based on two different $\textit{griz}$ libraries ($\textit{SDSS}$ DR12 and ATLAS All-Sky Stellar Reference Catalog, hereinafter $\textit{RefCat2}$). A comparison between both zero points is performed leading us to the choice of $\textit{RefCat2}$ as the base catalogue for this calibration, and applied to a field in the Cyg OB2 association.

This is the second paper of a series making use of Herschel/PACS spectroscopy of evolved stars in the THROES catalogue to study the inner regions of their circumstellar envelopes (CSEs). We analyze the CO emission spectra, including a large number of high-$J$ CO lines (from $J$=14-13 to $J$=45-44), as a proxy for the warm molecular gas in the CSEs of a sample of bright carbon-rich stars spanning different evolutionary stages from the Asymptotic Giant Branch (AGB) to the young planetary nebulae (PNe) phase. We use the rotational diagram (RD) to derive rotational temperatures ($T_{\rm rot}$) and masses ($M_{\rm H_2}$) of the envelope layers where the CO transitions arise. We also obtain a first order estimate of the mass-loss rates and assess the impact of the opacity correction for a range of characteristic envelope radii. We use multi-epoch spectra for the well studied C-rich envelope IRC+10216 to investigate the impact of CO flux variability on the values of $T_{\rm rot}$ and $M_{\rm H_2}$. PACS sensitivity allowed the study of higher rotational numbers than before, indicating the presence of a significant amount of warmer gas ($\sim$200-900 K) not traceable with lower-$J$ CO observations at sub-mm/mm wavelengths. The masses are in the range $\sim10^{-2}-10^{-5}\,\rm M_{\odot}$, anti-correlated with temperature. For some strong CO emitters we infer a double temperature (warm $T_{\rm rot}\sim$400 K and hot $T_{\rm rot}\sim$820 K) component. From the analysis of IRC+10216, we corroborate that the effect of line variability is perceptible on the $T_{\rm rot}$ of the hot component only, and certainly insignificant on $M_{\rm H_2}$ and, hence, the mass-loss rate. Therefore, the parameters derived from the RD are robust even when strong line flux variability occurs, with the major source of uncertainty in the estimate of the mass-loss rate being the size of the CO-emitting volume.

GALANTE is an optical photometric survey with seven intermediate/narrow filters that has been covering the Galactic Plane since 2016 using the Javalambre T80 and Cerro Tololo T80S telescopes. The P.I.s of the northern part (GALANTE NORTE) are Emilio J. Alfaro and Jes\'us Ma\'iz Apell\'aniz. and the P.I. of the southern part (GALANTE SUR) is Rodolfo H. Barb\'a. The detector has a continuous 1.4 degr x 1.4 degr field of view with a sampling of 0.55"/pixel and the seven filters are optimized to detect obscured early-type stars. The survey includes long, intermediate, short, and ultrashort exposure times to reach a dynamical range close to 20 magnitudes, something never achieved for such an optical project before. The characteristics of GALANTE allow for a new type of calibration scheme using external Gaia, Tycho-2, and 2MASS data that has already led to a reanalysis of the sensitivity of the Gaia G filter. We describe the project and present some early results. GALANTE will identify the majority of the early-type massive stars within several kpc of the Sun and measure their amount and type of extinction. It will also map the Halpha nebular emission, identify emission-line stars, and do other studies of low- and intermediate-mass stars.

We investigated the infrared (IR) dust emission from PN IC 418, using a detailed model controlled by a previous determination of the stellar properties and the characteristics of the photoionized nebula, keeping as free parameters the dust types, amounts and distributions relative to the distance of the central star. The model includes the ionized region and the neutral region beyond the recombination front (Photodissociation region, or PDR), where the [OI] and [CII] IR lines are formed. We succeeded in reproducing the observed infrared emission from 2 to 200~\mm. The global energy budget is fitted by summing up contributions from big grains of amorphous carbon located in the neutral region and small graphite grains located in the ionized region (closer to the central star). Two emission features seen at 11.5 and 30~\mm are also reproduced by assuming them to be due to silicon carbide (SiC) and magnesium and iron sulfides (Mg$_x$Fe$_{1-x}$S), respectively. For this, we needed to consider ellipsoidal shapes for the grains to reproduce the wavelength distribution of the features. Some elements are depleted in the gaseous phase: Mg, Si, and S have sub-solar abundances (-0.5 dex below solar by mass), while the abundance of C+N+O+Ne by mass is close to solar. Adding the abundances of the elements present in the dusty and gaseous forms leads to values closer to but not higher than solar, confirming that the identification of the feature carriers is plausible. Iron is strongly depleted (3 dex below solar) and the small amount present in dust in our model is far from being enough to recover the solar value. A remaining feature is found as a residue of the fitting process, between 12 and 25~\mm, for which we do not have identification.

We present seeing-limited narrow-band mid-IR GTC/CanariCam images of the spatially extended fullerene-containing planetary nebula (PN) IC 418. The narrow-band images cover the C60 fullerene band at 17.4 {\mu}m, the polycyclic aromatic hydrocarbon like (PAH-like) feature at 11.3 {\mu}m, the broad 9-13 {\mu}m feature, and their adjacent continua at 9.8 and 20.5 {\mu}m. We study the relative spatial distribution of these complex species, all detected in the Spitzer and Infrared Space Observatory (ISO) spectra of IC 418, with the aim of getting observational constraints to the formation process of fullerenes in H-rich circumstellar environments. A similar ring-like extended structure is seen in all narrow-band filters, except in the dust continuum emission at 9.8 {\mu}m, which peaks closer to the central star. The continuum-subtracted images display a clear ring-like extended structure for the carrier of the broad 9-13 {\mu}m emission, while the spatial distribution of the (PAH-like) 11.3 {\mu}m emission is not so well defined. Interestingly, a residual C60 17.4 {\mu}m emission (at about 4-{\sigma} from the sky background) is seen when subtracting the dust continuum emission at 20.5 {\mu}m. This residual C60 emission, if real, might have several interpretations; the most exciting being perhaps that other fullerene-based species like hydrogenated fullerenes with very low H-content may contribute to the observed 17.4 {\mu}m emission. We conclude that higher sensitivity mid-IR images and spatially resolved spectroscopic observations (especially in the Q-band) are necessary to get some clues about fullerene formation in PNe.

This is the first of a series of papers presenting the THROES (A caTalogue of HeRschel Observations of Evolved Stars) project, intended to provide a comprehensive overview of the spectroscopic results obtained in the far-infrared (55-670 microns) with the Her- schel space observatory on low-to-intermediate mass evolved stars in our Galaxy. Here we introduce the catalogue of interactively reprocessed PACS (Photoconductor Array Camera and Spectrometer) spectra covering the 55-200 microns range for 114 stars in this category for which PACS range spectroscopic data is available in the Herschel Science Archive (HSA). Our sample includes objects spanning a range of evolutionary stages, from the asymptotic giant branch to the planetary nebula phase, displaying a wide variety of chemical and physical properties. The THROES/PACS catalogue is accessible via a dedicated web-based inter- face (https://throes.cab.inta-csic.es/) and includes not only the science-ready Herschel spectroscopic data for each source, but also complementary photometric and spectroscopic data from other infrared observatories, namely IRAS (Infrared Astronomical Satellite), ISO (Infrared Space Observatory) or AKARI, at overlapping wavelengths. Our goal is to create a legacy-value Herschel dataset that can be used by the scientific community in the future to deepen our knowledge and understanding of these latest stages of the evolution of low-to-intermediate mass stars.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

The investigation of post-AGB objects (proto-planetary nebulae) is very important from the standpoint of physical and chemical changes occurring during the late stages of stellar evolution. The Toru\'n catalogue of Galactic post-AGB and related objects is an evolutive catalogue containing astrometric, photometric and spectroscopic data as well as HST images for all known post-AGB objects and candidates in our Galaxy. This free-access catalogue can serve as an ideal tool to study different groups of post-AGB objects, especially due to the fact that all information is gathered in one place. The second release of our catalogue introduces a simple classification scheme of post-AGB objects and includes a significant number of new objects, photometric data, spectra and images. Here, using objects from the catalogue we consider the problem of the termination of the AGB phase.

We present an analysis of a high resolution (R~50,000) optical spectrum of the central region of the proto-planetary nebula CRL 2688. This object is thought to have recently moved off the AGB, and display abundance patterns of CNO and heavy elements that can provide us with important clues to understand the nucleosynthesis, dredge-up and mixing experienced by the envelope of the central star during its AGB stage of evolution. The analysis of the molecular features, presumably originated from the circumstellar matter provides further constraints on the chemistry and velocity of the expanding shell, expelled as a consequence of the strong mass loss experienced by the central star. We confirm that the central star shows a spectrum typical of an F-type supergiant with Teff=7250 K, log g=0.5 and [Fe/H]=-0.3 dex. We find that the abundance pattern of this object is characterized by enhancements of Carbon ([C/Fe]=0.6), Nitrogen ([N/Fe]=1.0) and Na ([Na/Fe]=0.7), similar to other previously known carbon-rich post-AGB stars. Yttrium is also enhanced while the [Ba/Y] ratio is very low (-1.0), indicating that only the light s-process elements are enhanced. The Zinc abundance is found to be normal, [Zn/Fe]=0.0, suggesting that there is no depletion of refractory elements. The Halpha, Na I and K I resonance lines show prominent emission components, whose helio-centric radial velocities are offsetted by -41 km/s relative to the photospheric metal-absorption lines. The molecular features of C_2 and CN also show emission components, whose velocities are consistent with the emission components of the Halpha, Na I, and K I lines. On the other hand, their absorption components are more highly blue shifted than the corresponding emission components, which suggests that the regions where the emission and absorption components arise are expanding at different velocities.

The mission planning software developed for the Herschel Project is presented: The Herschel Inspector and mid/Longterm Scheduler (HILTS) and the short-term scientificMission Planning system (SMPS).

Visualization, querying, statistical analysis and mid-long-term scheduling are common concerns for any observatory. HILTS is a Java tool developed for the Herschel project to address all these issues in a unified way.

We present an analysis of elemental abundances of He, N, O, Ne, S, and Ar in Magellanic Cloud planetary nebulae (PNe), and focus initially on 14 PNe in the Small Magellanic Cloud (SMC). We derived the abundances from a combination of deep, high dispersion optical spectra, as well as mid-infrared (IR) spectra from the Spitzer Space Telescope. A detailed comparison with prior SMC PN studies shows that significant variations among authors of relative emission line flux determinations lead to systematic discrepancies in derived elemental abundances between studies that are >~0.15 dex, in spite of similar analysis methods. We used ionic abundances derived from IR emission lines, including those from ionization stages not observable in the optical, to examine the accuracy of some commonly used recipes for ionization correction factors (ICFs). These ICFs, which were developed for ions observed in the optical and ultraviolet, relate ionic abundances to total elemental abundances. We find that most of these ICFs work very well even in the limit of substantially sub-Solar metallicities, except for PNe with very high ionization. Our abundance analysis shows enhancements of He and N that are predicted from prior dredge-up processes of the progenitors on the AGB, as well as the well known correlations among O, Ne, S, and Ar that are little affected by nucleosynthesis in this mass range. We identified MG_8 as an interesting limiting case of a PN central star with a ~3.5 M_sun progenitor in which hot-bottom burning did not occur in its prior AGB evolution. We find no evidence for O depletion in the progenitor AGB stars via the O-N cycle, which is consistent with predictions for lower-mass stars. We also find low S/O ratios relative to SMC H_II regions, with a deficit comparable to what has been found for Galactic PNe.

We present AKARI/IRC observations of a sample of six extremely red IRAS sources, of which three are variable OH/IR stars and the rest are early post-AGB stars. The OH/IR stars show a red continuum with the expected strong 10 micron silicate absorption feature, while the post-AGB stars show an even redder continuum accompanied with a comparably weak silicate absorption. We modelled the spectral energy distributions with DUSTY. While for the OH/IR stars a reasonable fit can be obtained with almost pure silicate dust, the post-AGB stars require a mixture of silicate and carbon-rich dust. We assume that in the latter objects the inner dust shell is carbon-rich, while the outer shells are still oxygen-rich.

We present AKARI/IRC and Spitzer/IRS observations of a selected sample of galactic IRAS sources considered to be heavily obscured AGB/post-AGB stars based on their characteristic IRAS colours. All of them are completely invisible in the optical range but extremely bright in the infrared. Based on AKARI and Spitzer spectroscopy and using DUSTY we are able to determine the dominant chemistry of their circumstellar shells as well as the properties of the dust grains contained in these shells. Most of the sources are found to be C-rich (being the reddest C-rich stars observed so far). We find only molecular absorptions (and no PAH features) such as acetylene (C2H2) at 13.7 micron, indicative of an early post-AGB stage. We shortly discuss our findings in the context of stellar evolution during the hidden "transition phase" from AGB stars to Planetary Nebulae.

We present the results of our preliminary study of all known Galactic PNe (included in the Kerber 2003 catalog) which are detected by the AKARI/FIS All-Sky Survey as identified in the AKARI/FIS Bright Source Catalog (BSC) Version Beta-1.

We present first preliminary results from AKARI/FIS pointed observations of post-AGB stars and planetary nebulae (PNe). A first analysis of the radial (azimuthally averaged) profile of the observed sources shows no evidence for excess emission due to the presence of circumstellar dust. No (detached) circumstellar faint dust-shells are seen in the image maps. Also, we present here first results of aperture flux photometry at wavelengths of 65, 90, 140 and 160 micron. Results are compared with IRAS flux densities as well as the beta release of the FIS Bright Source Catalog. Finally, spectral energy distributions are given, by way of an example, for two individual targets in our sample.

The AKARI and Spitzer satellites provided an unique opportunity to observe a variety of stars, which are considered as departing from the Asymptotic Giant Branch (AGB) and have started their post-AGB evolution recently. Most of these stars are absent optically and are bright in the mid-IR wavelength range. Spectra of close to 200 objects have been obtained. For all of them the 1-60 micron spectral energy distribution has been constructed using photometric data from various surveys. We report here on the results of Spitzer observations of 88 IRAS selected post-AGB candidates and discuss them in comparison to the results of the AKARI observations of post-AGB candidates reported elsewhere in these proceedings. The dust compositions can be divided broadly in oxygen- and carbon-rich types, but a variety of intermediate types have been found. Among the oxygen-rich stars amorphous dust prevails, but a few sources show emission features from crystalline dust. The spectra from carbon-rich shells may be completely featureless, may show emission features from PAHs or a molecular absorption line from C2H2. We found also sources with a neon emission line at 12.8 micron. More than a third of all sources show a near-infrared excess at < 5 micron and almost all of them show evidence of C-rich dust in their shells. We postulate that the emerging post-AGB wind after the end of AGB evolution contains always carbon-rich dust irrespective of the chemistry of the former AGB star.

We performed a near-IR/optical monitoring programme from 1999 to 2005 in order to study the variability properties of the 'Arecibo sample of OH/IR stars' (periods, amplitudes, and colour variations). Here we describe this multi-wavelength long-term monitoring programme. Data analysis is still in process. Our ultimate goal is to study in particular the oxygen-rich AGB stars with masses larger than 2 Solar Mass, which are probably rare among AGB stars discovered optically and/or in the near-IR, but are common in samples discovered in the mid- to far-IR.

We report the detection of nebular emission lines in the optical and mid-infrared spectra of IRAS 17347-3139, a heavily obscured OH/IR star which may be rapidly evolving from the AGB to the PN stage. The presence of emission lines is interpreted as a clear indication that the ionization of its circumstellar envelope has already started. This source belongs to the rare class of objects known as `OHPNe' displaying both OH maser and radio continuum emission. However, unlike the rest of stars in this class, prominent C-rich dust features are detected in its mid-infrared spectrum, which makes the analysis of this star particularly interesting.

We present the results of a near-infrared monitoring programme of a selected sample of stars, initially suspected to be Mira variables and OH/IR stars, covering more than a decade of observations. The objects monitored cover the typical range of IRAS colours shown by O-rich stars on the Asymptotic Giant Branch and show a surprisingly large diversity of variability properties. 16 objects are confirmed as large-amplitude variables. Periods between 360 and 1800 days and typical amplitudes from 1 to 2 magnitudes could be determined for nine of them. In three light curves we find a systematic decrease of the mean brightness, two light curves show pronounced asymmetry. One source, IRAS 07222-2005, shows infrared colours typical of Mira variables but pulsates with a much longer period (approx. 1200 days) than a normal Mira. Two objects are ither close to (IRAS 03293+6010) or probably in (IRAS 18299-1705) the post-AGB phase. In IRAS 16029-3041 we found a systematic increase of the H-K colour of approximately 1 magnitude, which we interpret as evidence of a recent episode of enhanced mass loss. IRAS 18576+0341, a heavily obscured Luminous Blue Variable was also monitored. The star showed a continued decrease of brightness over a period of 7 years (1995 - 2002).

In this paper we present optical and near-infrared finding charts, accurate astrometry (~1") and single-epoch near-infrared photometry for 371 IRAS sources, 96% of those included in the so-called Arecibo sample of OH/IR stars (Eder et al. 1988; Lewis et al. 1990a; Chengalur et al. 1993). The main photometric properties of the stars in the sample are presented and discussed as well as the problems found during the process of identification of the optical/near-infrared counterparts. In addition, we also identify suitable reference stars in each field to be used for differential photometry purposes in the future. We find that 39% of the sources (144 in number) have no optical counterpart, 8 of them being invisible even at near infrared wavelengths. The relative distribution of sources with and without optical counterpart in the IRAS two-colour diagram and their characteristic near infrared colours are interpreted as the consequence of the increasing thickness of their circumstellar shells. Among the objects not detected at near infrared wavelengths four non-variable sources are proposed to be heavily obscured post-AGB stars which have just very recently left the AGB. Eight additional objects with unusually bright and/or blue near-infrared colours are identified as candidate post-AGB stars and/or proto-planetary nebulae.

We present optical and near-infrared finding charts taken from the DSS and 2MASS surveys of 94 IRAS sources selected from the GLMP catalogue (Garc\'{\i}a-Lario 1992), and accurate astrometry (~0.2") for most of them. Selection criteria were very red IRAS colours representative for OH/IR stars with optically thick circumstellar shells and the presence of variability according to the IRAS variability index (VAR>50). The main photometric properties of the stars in this `GLMP sample' are presented, discussed and compared with the correspondent properties of the `Arecibo sample' of OH/IR stars studied in Jim\'enez-Esteban et al. (2005a). We find that 37% of the sample (N=34) have no counterpart in the 2MASS, implying extremely high optical depths of their shells. Most of the sources identified in the 2MASS are faint (K>~8) and of very red colour in the near-infrared, as expected. The brightest 2MASS counterpart (K=5.3mag) was found for IRAS 18299--1705. Its blue colour H--K=1.3 suggests that IRAS 8299--1705 is a post-AGB star. A couple of GLMP sources have faint but relatively blue counterparts. They might be misidentifed field stars or stars hich experienced recently a drop of their mass loss rates. The `GLMP sample' in general is made of oxygen-rich AGB stars, which are highly obscured by their circumstellar shells. They belong to the same population as the reddest OH/IR stars in the `Arecibo sample'.

The `O-rich AGB sequence' is a sequence of colours describing the location of O-rich AGB stars in the IRAS two-colour diagram [12]--[25] vs [25]--[60]. We propose an evolutionary scenario for this sequence in which all stars, independent of their progenitor mass, start the AGB phase in the blue part of the `O-rich AGB sequence' and then evolve toward redder colors, although only the more massive stars would reach the very end of the `O-rich AGB sequence'. The sources located in the blue part of the sequence are mainly Mira variables, whose mean period is increasing with the IRAS colours. Most of them will evolve into O-rich Type II (and III) Planetary Nebulae. Part of the stars located in the red part of the sequence will change their chemical composition from O-rich to C-rich during their evolution in the AGB phase, and might evolve into C-rich Type II Planetary Nebulae. Hot bottom burning may prevent the conversion to carbon stars of the rest of sources located in the red part of the sequence and they will end up as N-rich Type I Planetary Nebulae.