We revisit the subject of Be star candidates towards the Magellanic Clouds, previously studied by the authors using SPM4 proper motions. We obtain GAIA DR2 parallaxes and proper motions for 2357 and 994 LMC and SMC Be candidates, respectively. Parallaxes and proper motions vs. color V-I easily reveal the presence of the redder galactic contaminant foreground, as concluded in our previous work, but this time we do find a few red Be stars candidates consistent with being true Magellanic objects. A membership assessment to each Magellanic Cloud is done for each Be candidate, based on the distribution of their parallaxes and proper motions. From a compilation of published catalogues of spectroscopically confirmed Be stars, we found that 40 (LMC) and 64 (SMC) of these Be candidates, are in fact Be stars. Near infrared IRSF JHKs magnitudes were obtained for about 70% the Be stars candidates with GAIA DR2 astrometric data. Mid-infrared SAGE IRAC 3.6, 4.5, 5.8 and 8.0 $\mu$m magnitudes were obtained for about 85% as well. 6 LMC and 7 SMC confirmed Be stars show optical, near- or mid-infrared colours redder than what has been typically measured for Classical Be stars. Several of the Be candidates follow those redder-than-expected colours distributions suggesting the existence of more red Magellanic Be stars.

Statistical pattern recognition methods have provided competitive solutions for variable star classification at a relatively low computational cost. In order to perform supervised classification, a set of features is proposed and used to train an automatic classification system. Quantities related to the magnitude density of the light curves and their Fourier coefficients have been chosen as features in previous studies. However, some of these features are not robust to the presence of outliers and the calculation of Fourier coefficients is computationally expensive for large data sets. We propose and evaluate the performance of a new robust set of features using supervised classifiers in order to look for new Be star candidates in the OGLE-IV Gaia south ecliptic pole field. We calculated the proposed set of features on six types of variable stars and on a set of Be star candidates reported in the literature. We evaluated the performance of these features using classification trees and random forests along with K-nearest neighbours, support vector machines, and gradient boosted trees methods. We tuned the classifiers with a 10-fold cross-validation and grid search. We validated the performance of the best classifier on a set of OGLE-IV light curves and applied this to find new Be star candidates. The random forest classifier outperformed the others. By using the random forest classifier and colour criteria we found 50 Be star candidates in the direction of the Gaia south ecliptic pole field, four of which have infrared colours consistent with Herbig Ae/Be stars. Supervised methods are very useful in order to obtain preliminary samples of variable stars extracted from large databases. As usual, the stars classified as Be stars candidates must be checked for the colours and spectroscopic characteristics expected for them.

We present a proper motion investigation of a sample of Be star candidates towards the Magellanic Clouds, which has resulted in the identification of separate populations, in the Galactic foreground and in the Magellanic background. Be stars are broadly speaking B-type stars that have shown emission lines in their spectra. In this work, we studied a sample of 2446 and 1019 Be star candidates towards the LMC and SMC respectively, taken from the literature and proposed as possible Be stars due to their variability behaviour in the OGLE-II I band. JHKs magnitudes from the IRSF catalog and proper motions from the SPM4 catalog, were obtained for 1188 and 619 LMC and SMC Be stars candidates, respectively. Color-color and vector-point diagrams were used to identify different populations among the Be star candidates. In the LMC sample, two populations with distinctive infrared colours and kinematics were found, the bluer sample is consistent with being in the LMC and the redder one with belonging to the Milky Way disk. This settles the nature of the redder sample which had been described in previous publications as a possible unknown subclass of stars among the Be candidates in the LMC. In the SMC sample, a similar but less evident result was obtained, since this apparent unknown subclass was not seen in this galaxy. We confirm that in the selection of Be stars by their variability, although generally successful, there is a higher risk of contamination by Milky Way objects towards redder B$-$V and V$-$I colors.

The investigation of the non-linearity of the Leavitt law is a topic that began more than seven decades ago, when some of the studies in this field found that the Leavitt law has a break at about ten days. The goal of this work is to investigate a possible statistical cause of this non-linearity. By applying linear regressions to OGLE-II and OGLE-IV data, we find that, in order to obtain the Leavitt law by using linear regression, robust techniques to deal with influential points and/or outliers are needed instead of the ordinary least-squares regression traditionally used. In particular, by using $M$- and $MM$-regressions we establish firmly and without doubts the linearity of the Leavitt law in the Large Magellanic Cloud, without rejecting or excluding Cepheid data from the analysis. This implies that light curves of Cepheids suggesting blending, bumps, eclipses or period changes, do not affect the Leavitt law for this galaxy. For the SMC, including this kind of Cepheids, it is not possible to find an adequate model, probably due to the geometry of the galaxy. In that case, a possible influence of these stars could exist.

The universality and linearity of the Leavitt law are hypotheses commonly adopted in studies of galaxy distances using Cepheid variables as standard candles. In order to test these hypotheses, we obtain slopes of the Leavitt law using linear regressions of fundamental-mode Cepheids observed by the Optical Gravitational Lensing Experiment project in the Magellanic Clouds. We find that the slopes in $VI$-bands and in the Wesenheit index behave exponentially, indicating non-linearity. We also find that the slopes obtained using long-period Cepheids can be considered as universal in the $VI$-bands, but not in the Wesenheit index.

Mennickent et al.and Sabogal et al.identified a large number of Classical Be (CBe) candidates in the L&SMC based on their photometric variability using the OGLEII database. They classified these stars into four different groups based on the appearance of their variability. We studied the infrared properties of the sample as well as the spectroscopic properties of a subsample. We cross-correlated the optical sample with the IRSF catalog to obtain the J, H, Ks magnitudes of all the four types of stars in the L&SMC. Spectra of 120 stars belonging to the types 1, 2 and 3 were analysed to study their spectral properties. Among the four types, the type 4 stars is the dominant group. The NIR colour-colour diagrams suggest that the type 4 stars in the LMC have a subclass, which is not found in our Galaxy or in the SMC. The main type 4 sample which is \sim 49% of the total sample has NIR properties similar to the Galactic CBe stars and the SMC type 4 stars. Though the new subclass of type 4 stars have high E(B - V) \sim 0.75, they are not located close to regions with high reddening. The type 3 stars (\sim 6% & 7.3% in the L&SMC) are found to have large H{\alpha} EW in the SMC and some are found to have large NIR excess. This small fraction of stars are unlikely to be CBe stars. The type 2 stars are found in larger fraction in the SMC, when compared to the LMC. The spectroscopic and the NIR properties suggest that these could be CBe stars. The spectroscopic sample of type 1 stars which show H{\alpha} in emission and confirmed as CBe stars are more abundant in the SMC by a factor of 2.6. If the effect of metallicity is to cause more CBe stars in the SMC, when compared to the LMC, then type 1, type 2 and type 4 stars follow this rule, with an enhancement of 2.6, 2.4 and 1.3 respectively.

We present new L-band spectra of 13 outbursting Be stars obtained with ISAAC at the ESO Paranal observatory. These stars can be classified in three groups depending on the presence/absence of emission lines and the strength of Br$\alpha$ and Pf$\gamma$ emission lines relative to those of Humphreys lines from transitions 6--14 to the end of the series. These groups are representative of circumstellar envelopes with different optical depths. For the group showing Br$\alpha$ and Pf$\gamma$ lines stronger than Humphreys lines, the Humphreys decrement roughly follow the Menzel case-B for optically thin conditions. For the group showing comparable Br$\alpha$, Pf$\gamma$ and Humphreys emission line strengths, the Humphreys decrements moves from an optically thin to an optically thick regime at a transition wavelength which is characteristic for each star, but typically is located around 3.65--3.75 $\mu$m (transitions 6--19 and 6--17). Higher order Humphreys lines probe optically thin inner regions even in the optically thicker envelopes. We find evidence of larger broadening in the infrared emission lines compared with optical lines, probably reflecting larger vertical velocity fields near the star. The existence of the aforementioned groups is in principle consistent with the proposed description by de Wit et al. (2006) for Be star outbursts in terms of the ejection of an optically thick disk that expands and becomes optically thin before dissipation into the interstellar medium. Time resolved L-band spectroscopy sampling the outburst cycle promises to be an unique tool for testing Be star disk evolution.

Detailed studies of Be stars in environments with different metallicities like the Magellanic Clouds or the Galactic bulge are necessary to understand the formation and evolution mechanisms of the circumstellar disks. However, a detailed study of Be stars in the direction of the bulge of our own galaxy has not been performed until now. We report the first systematic search for Be star candidates in the direction of the Galactic Bulge. We present the catalogue, give a brief description of the stellar variability seen, and show some light curve examples. We searched for stars matching specific criteria of magnitude, color and variability in the I band. Our search was conducted on the 48 OGLE II fields of the Galactic Bulge.This search has resulted in 29053 Be star candidates, 198 of them showing periodic light variations. Nearly 1500 stars in this final sample are almost certainly Be stars, providing an ideal sample for spectroscopic multiobject follow-up studies.