The Galactic emission-line object MWC 349A is one of the brightest radio stars in the sky. The central object is embedded in an almost edge-on oriented Keplerian rotating thick disc that seems to drive a rotating bipolar wind. The dense disc is also the site of hot molecular emission such as the CO bands with its prominent band heads in the near-infrared spectral range. Despite numerous studies, the nature of MWC 349A is still controversial with classifications ranging from a pre-main sequence object to an evolved supergiant. We collected new high-resolution near-infrared spectra in the K- and L-bands using the GNIRS spectrograph at Gemini-North to study the molecular disc of MWC 349A, and in particular to search for other molecular species such as SiO and the isotope 13CO. The amount of 13CO, obtained from the 12CO/13CO ratio, is recognized as an excellent tool to discriminate between pre-main sequence and evolved massive stars. We find no signatures of SiO band emission, but detect CO band emission with considerably lower intensity and CO gas temperature compared to previous observations. Moreover, from detailed modelling of the emission spectrum, we derive an isotope ratio of 12CO/13CO = 4 +- 1. Based on this significant enrichment of the circumstellar environment in 13CO we conclude that MWC 349A belongs to the group of B[e] supergiants, and we discuss possible reasons for the drop in CO intensity.

A significant number of Be stars show a second Balmer discontinuity (sBD) attributed to an extended circumstellar envelope (CE). The fast rotational velocity of Be stars undoubtedly plays a significant role in the formation of the CE. However, Bn stars, which are also B-type rapidly rotating stars, do not all present clear evidence of being surrounded by circumstellar material. We aim to characterize the populations of Be and Bn stars, and discuss the appearance of the sBD as a function of the stellar parameters. We expect to find new indices characterizing the properties of CEs in Be stars and properties relating Be and Bn stars. Correlations of the aspect and intensity of the sBD and the emission in the H$\alpha$ line with the stellar parameters and the $V\!\sin i$ are presented. Some Bn stars exhibit the sBD in absorption, which may indicate the presence of rather dense CEs. Six Bn stars show emission in the H$\alpha$ line, so they are reclassified as Be stars. The sBD in emission appears in Be stars with $V\!\sin i \lesssim 250$ km\,s$^{-1}$, and in absorption in both Be and Bn stars with \mbox{$V\!\sin i \gtrsim 50$ km\,s$^{-1}$}. Low-mass Be and Bn stars share the same region in the Hertzsprung-Russell diagram. The distributions of rotational to critical velocity ratios of Be and Bn stars corresponding to the current stellar evolutionary stage are similar, while distributions inferred for the zero-age main sequence have different skewness. We found emission in the H$\alpha$ line and signs of a CE in some Bn stars, which motivated us to think that Bn and Be stars probably belong to the same population. It should be noted that some of the most massive Bn stars could display the Be phenomenon at any time. The similarities found among Be and Bn stars deserve to be more deeply pursued.

Massive stars present strong stellar that which are described by the radiation driven wind theory. Accurate mass-loss rates are necessary to properly describe the stellar evolution across the Hertzsprung--Russel Diagram. We present a self-consistent procedure that coupled the hydrodynamics with calculations of the line-force, giving as results the line-force parameters, the velocity field, and the mass-loss rate. Our calculations contemplate the contribution to the line-force multiplier from more than $\sim 900,000$ atomic transitions, an NLTE radiation flux from the photosphere and a quasi-LTE approximation for the occupational numbers. A full set of line-force parameters for $T_\text{eff}\ge 32,000$ K, surface gravities higher than 3.4 dex for two different metallicities are presented, with their corresponding wind parameters (terminal velocities and mass-loss rates). The already known dependence of line-force parameters on effective temperature is enhanced by the dependence on $\log g$. The terminal velocities present a stepper scaling relation with respect to the escape velocity, this might explain the scatter values observed in the hot side of the bistability jump. Moreover, a comparison of self-consistent mass-loss rates with empirical values shows a good agreement. Self-consistent wind solutions are used as input in FASTWIND to calculate synthetic spectra. We show, comparing with the observed spectra for three stars, that varying the clumping factor, the synthetic spectra rapidly converge into the neighbourhood region of the solution. It is important to stress that our self-consistent procedure significantly reduces the number of free parameters needed to obtain a synthetic spectrum.

Variable B supergiants (BSGs) constitute a heterogeneous group of stars with complex photometric and spectroscopic behaviours. They exhibit mass-loss variations and experience different types of oscillation modes, and there is growing evidence that variable stellar winds and photospheric pulsations are closely related. To discuss the wind properties and variability of evolved B-type stars, we derive new stellar and wind parameters for a sample of 19 Galactic BSGs by fitting theoretical line profiles of H, He, and Si to the observed ones and compare them with previous determinations. The mass-loss rate of three stars has been obtained for the first time. All stars follow the known empirical wind momentum-luminosity relationships, and the late BSGs show the trend of the mid BSGs. We find an empirical relationship that associates the amplitude of mass-loss variations with photometric/spectroscopic variability on timescales of tens of days. We find that stars located on the cool side of the bi-stability jump show a decrease in the terminal and escape velocities ratio, while their corresponding mass-loss rates are similar to or lower than the values found for stars on the hot side. Particularly, for those variable stars a decrease in the terminal and escape velocities ratio is accompanied by a decrease in the mass-loss rate. Our results also suggest that radial pulsation modes with periods longer than 6 days might be responsible for the wind variability in the mid/late-type BSGs. These radial modes might be identified with strange modes, which are known to facilitate (enhanced) mass loss. On the other hand, we propose that the wind behaviour of stars on the cool side of the bi-stability jump could fit with predictions of the delta-slow hydrodynamics solution for radiation-driven winds with highly variable ionization.

Be stars are rapid rotators surrounded by a gaseous disk envelope whose origin is still under debate. This envelope is responsible for observed emission lines and large infrared excess. To progress in the understanding of the physical processes involved in the disk formation, we estimate the disk parameters for a sample of Be stars and search for correlations between these parameters and stellar properties. We performed spectro-interferometric observations of 26 Be stars in the region of the Br$\gamma$ line to study the kinematical properties of their disks through the Doppler effect. Observations were performed at the Paranal observatory with the VLTI/AMBER interferometer. This instrument provides high spectral and high spatial resolutions. We modeled 18 Be stars with emission in the Br$\gamma$ line. The disk kinematic is described by a quasi-Keplerian rotation law, with the exception of HD28497 that presents a one-arm density-wave structure. Using a combined sample, we derived a mean value for the velocity ratio V/Vc=0.75, and found that rotation axes are probably randomly distributed in the sky. Disk sizes in the line component model are in the range of 2-13 stellar radii and do not correlate with the effective temperature or spectral type. However, we found that the maximum size of a stable disk correlates with the rotation velocity at the inner part of the disk and the stellar mass. We found that, on average, the Be stars of our combined sample do not rotate at their critical velocity. However, the centrifugal force and mass of the star defines an upper limit size for a stable disk configuration. For a given rotation, high-mass Be stars tend to have more compact disks than their low-mass counterparts. It would be interesting to follow up the evolution of the disk size in variable stars to better understand the formation and dissipation processes of their circumstellar disks.

B[e] Supergiants are luminous evolved massive stars. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dusty regions usually found in rings or disk-like structures. For a better comprehension of the mechanisms behind the formation of these rings, detailed knowledge about their structure and dynamics is essential. To address that, we obtained high-resolution optical and near-infrared spectra for 8 selected Galactic B[e] Supergiants, for which CO emission has been detected. Assuming Keplerian rotation for the disk, we combine the kinematics obtained from the CO bands in the near-IR with those obtained by fitting the forbidden emission [OI] $\lambda$5577, [OI] $\lambda\lambda$6300,6363, and [CaII] $\lambda\lambda$7291,7323 lines in the optical to probe the disk structure. We find that the emission originates from multiple ring structures around all B[e] Supergiants, with each one of them displaying a unique combination of rings regardless of whether the object is part of a binary system. The confirmed binaries display spectroscopic variations of their line intensities and profiles as well as photometric variability, whereas the ring structures around the single stars are stable.

The advanced stages of several high-mass stars are characterized by episodic mass loss shed during phases of instability. Key for assigning these stars a proper evolutionary state is to assess the composition and geometry of their ejecta alongside the stellar properties. We selected five hot LBV candidates in M33 to refine their classification, investigate their circumstellar environments and explore their evolutionary properties. Being accessible targets in the near-infrared, we conducted medium-resolution spectroscopy with GNIRS/GEMINI in the $K-$band to investigate their molecular circumstellar environments. Two stars were found to display CO emission, which was modeled to emerge from a circumstellar or circumbinary Keplerian disk/ring. The identification of the carbon isotope $^{13}$C and, for one of the two stars, a significantly low $^{12}$CO/$^{13}$CO ratio, implies an evolved stellar state. As both CO emission stars are highly luminous and hence do not undergo a red supergiant phase, we suggest that stripping processes and equatorial high-density ejecta due to fast rotation are responsible for the enrichment of the stellar surface with processed material from the core. A candidate B[e]SG displays an absorption CO profile, which may be attributed to a jet or stellar pulsations. The featureless infrared spectra of two stars suggest a low-density gas shell or dissipation of the molecule due to the ionizing temperature of the star. We propose spectroscopic monitoring of our targets to evaluate the stability of the CO molecule and assess the time-dependent dynamics of the circumstellar gas structures.

We find that the emission line object OGLEJ005039.05-725751.4, a member of the cluster OGLE-CL SMC 64, exhibits a peculiar light curve pattern repeating with a recurrence time of 141.45 days. The light curve resembles periodic outbursts with a duty cycle of 20%. A second long-cycle of 2500 days is also detected in the photometric dataset. Two X-SHOOTER spectra obtained at minimum and maximum reveal a Be star dominating at minimum light resembling the Classical Be star 48 Lib. The larger H$\alpha$ emission, the stronger NaD absorption and the appearance of emission in the infrared Ca II triplet at maximum, might indicate periodic mass transfer in a complex binary system.

Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $\Omega$-slow solution characterized by a dense and slow wind. Here, in this work we study the transition region of the solutions where the fast solution changes to the $\Omega$-slow. Using both time-steady and time-dependent numerical codes, we study this transition region for different equatorial models of B-type stars. In all the cases, at certain range of rotational speeds, we found a region where the fast and $\Omega$-slow solution can co-exist. We find that the type of solution obtained in this co-existence region depends heavily on the initial conditions of our models. We also test the stability of the solutions within the co-existence region by performing base density perturbations in the wind. We find that under certain conditions, the fast solution can switch to a $\Omega$-slow solution, or vice versa. Such solution switching may be a possible contributor of material injected into the circumstellar environment of Be stars, without requiring rotational speeds near critical values.

B[e] supergiants (SGs) are massive post-main-sequence stars, surrounded by a complex circumstellar (CS) environment. The aim of this work is to investigate the structure and kinematics of the CS disc of the B[e] SG LHA 120-S 35. We used high-resolution optical spectra obtained in different years to model the forbidden emission lines and determine the kinematical properties of their line-forming regions, assuming Keplerian rotation. We also used low-resolution near-infrared (IR) spectra to explore the variability of molecular emission. LHA 120-S 35 displays spectral variability in both optical and IR regions. The P-Cygni line profiles of H I, as well as those of Fe II and O I, suggest the presence of a strong bipolar clumped wind. We distinguish density enhancements in the P-Cygni absorption component of the first Balmer lines, which show variations in both velocity and strength. The P-Cygni profile emission component is double-peaked, indicating the presence of a rotating CS disc. We also observe line-profile variations in the permitted and forbidden features of Fe II and O I. In the IR, we detect variations in the intensity of the H I emission lines as well as in the emission of the CO band-heads. Moreover, we find that the profiles of each [Ca II] and [O I] emission lines contain contributions from spatially different (complete or partial) rings. Globally, we find evidence of detached multi-ring structures, revealing density variations along the disc. We suggest that LHA 120-S 35 has passed through the red-supergiant (RSG) phase and evolves back bluewards in the Hertzsprung-Russell diagram. The formation of the complex CS structure could be the result of the wind-wind interactions of the post-RSG wind with the previously ejected material from the RSG. However, the presence of a binary companion can not be excluded. Finally, we find that LHA 120-S 35 belongs to a young stellar cluster.

Context. Stellar physical properties of star clusters are poorly known. Aims. Our goals are to perform a spectrophotometric study of the B star population in open clusters to derive accurate stellar parameters, search for the presence of circumstellar envelopes, and discuss the characteristics of these stars. Methods. The BCD spectrophotometric system is a powerful tool to obtain fundamental parameters and infer the main properties of open clusters: distance modulus, color excess, and age. We inspected the Balmer discontinuity to seek circumstellar disks and identify Be-star candidates. High resolution spectra in the H$\alpha$ region are used to confirm the Be nature. Results. We provide Teff , log g, Mv , Mbol and spectral types for a sample of 68 stars in the field of the open clusters NGC 6087, NGC 6250, NGC 6383, and NGC 6530, as well as the cluster distances, ages and reddening. Then, based on a sample of 230 B stars in the direction of the 11 open clusters studied along this series of three papers, we report 6 new Be stars, 4 blue straggler candidates, and 15 B-type stars (called Bdd) with a double Balmer discontinuity. We also find that the majority of the Be stars are dwarfs and present a maximum at the spectral type B2-B4 in young and intermediate-age open clusters. Another maximum of Be stars is observed at the spectral type B6-B8 in open clusters older than 40 Myr, where the population of Bdd stars also becomes relevant. Conclusions. Our results support previous statements that the Be phenomenon is present along the whole main sequence band and occurs in very different evolutionary states. We find clear evidence of an increase of stars with circumstellar envelopes with cluster age. The Be phenomenon reaches its maximum in clusters of intermediate age and the number of B stars with circumstellar envelopes (Be+Bdd stars) is also high for the older clusters.

The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue.

As the disk formation mechanism(s) in Be stars is(are) as yet unknown, we investigate the role of rapidly rotating radiation-driven winds in this process. We implemented the effects of high stellar rotation on m-CAK models accounting for: the shape of the star, the oblate finite disk correction factor, and gravity darkening. For a fast rotating star, we obtain a two-component wind model, i.e., a fast, thin wind in the polar latitudes and an $\Omega$-slow, dense wind in the equatorial regions. We use the equatorial mass densities to explore H$\alpha$ emission profiles for the following scenarios: 1) a spherically symmetric star, 2) an oblate shaped star with constant temperature, and 3) an oblate star with gravity darkening. One result of this work is that we have developed a novel method for solving the gravity darkened, oblated m-CAK equation of motion. Furthermore, from our modeling we find a) the oblate finite disk correction factor, for the scenario considering the gravity darkening, can vary by at least a factor of two between the equatorial and polar directions, influencing the velocity profile and mass-loss rate accordingly, b) the H$\alpha$ profiles predicted by our model are in agreement with those predicted by a standard power-law model for following values of the line-force parameters: $1.5 \lesssim k \lesssim 3$, $ \, \alpha \sim 0.6$ and $\, \delta \gtrsim 0.1$, and c) the contribution of the fast wind component to the H$\alpha$ emission line profile is negligible; therefore, the line profiles arise mainly from the equatorial disks of Be stars.

We assume that stars may undergo surface differential rotation to study its impact on the interpretation of $V\!\sin i$ and on the observed distribution $\Phi(u)$ of ratios of true rotational velocities $u=V/V_\rm c$ ($V_\rm c$ is the equatorial critical velocity). We discuss some phenomena affecting the formation of spectral lines and their broadening, which can obliterate the information carried by $V\!\sin i$ concerning the actual stellar rotation. We studied the line broadening produced by several differential rotational laws, but adopted Maunder's expression $\Omega(\theta)=\Omega_o(1+\alpha\cos^2\theta)$ as an attempt to account for all of these laws with the lowest possible number of free parameters. We studied the effect of the differential rotation parameter $\alpha$ on the measured $V\!\sin i$ parameter and on the distribution $\Phi(u)$ of ratios $u=V/V_\rm c$. We conclude that the inferred $V\!\sin i$ is smaller than implied by the actual equatorial linear rotation velocity $V_\rm eq$ if the stars rotate with $\alpha<0$, but is larger if the stars have $\alpha>0$. For a given $|\alpha|$ the deviations of $V\!\sin i$ are larger when $\alpha<0$. If the studied Be stars have on average $\alpha<0$, the number of rotators with $V_\rm eq\simeq0.9V_\rm c$ is larger than expected from the observed distribution $\Phi(u)$; if these stars have on average $\alpha>0$, this number is lower than expected. We discuss seven phenomena that contribute either to narrow or broaden spectral lines, which blur the information on the rotation carried by $V\!\sin i$ and, in particular, to decide whether the Be phenomenon mostly rely on the critical rotation. We show that two-dimensional radiation transfer calculations are needed in rapid rotators to diagnose the stellar rotation more reliably.

The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star, while a cold CO envelope is encircling the borders of the optical nebula from the south to the west.

B[e] Supergiants are a phase in the evolution of some massive stars for which we have observational evidence but no predictions by any stellar evolution model. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dust regions usually found in rings or disk-like structures. However, the detailed structure and the formation of the circumstellar environment are not well-understood, requiring further investigation. To address that we initiated an observing campaign to obtain a homogeneous set of high-resolution spectra in both the optical and NIR (using MPG-ESO/FEROS, GEMINI/Phoenix and VLT/CRIRES, respectively). We monitor a number of Galactic B[e] Supergiants, for which we examined the [OI] and [CaII] emission lines and the bandheads of the CO and SiO molecules to probe the structure and the kinematics of their formation regions. We find that the emission from each tracer forms either in a single or in multiple equatorial rings.

We aim to improve our knowledge on the structure and dynamics of the circumstellar disk of the LMC B[e] supergiant LHA 120-S 73. High-resolution optical and near-IR spectroscopic data were obtained over a period of 16 and 7 years, respectively. The spectra cover the diagnostic emission lines from [CaII] and [OI], as well as the CO bands. These features trace the disk at different distances from the star. We analyzed the kinematics of the individual emission regions by modeling their emission profiles. A low-resolution mid-infrared spectrum was obtained as well, which provides information on the composition of the dusty disk. All diagnostic emission features display double-peaked line profiles, which we interpret as due to Keplerian rotation. We find that LHA 120-S 73 is surrounded by at least four individual rings of material with alternating densities (or by a disk with strongly non-monotonic radial density distribution). Moreover, we find that the molecular ring must have gaps or at least strong density inhomogeneities, or in other words, a clumpy structure. The mid-infrared spectrum displays features of oxygen- and carbon-rich grain species, which indicates a long-lived, stable dusty disk. We cannot confirm the previously reported high value for the stellar rotation velocity. The line profile of HeI 5876 A is strongly variable in both width and shape and resembles of those seen in non-radially pulsating stars. A proper determination of the real underlying stellar rotation velocity is hence not possible. The existence of multiple stable and clumpy rings of alternating density recalls ring structures around planets. Although there is currently insufficient observational evidence, it is tempting to propose a scenario with one (or more) minor bodies or planets revolving around LHA 120-S 73 and stabilizing the ring system, in analogy to the shepherd moons in planetary systems.

Blue supergiant stars are known to display photometric and spectroscopic variability that is suggested to be linked to stellar pulsations. Pulsational activity in massive stars strongly depends on the star's evolutionary stage and is assumed to be connected with mass-loss episodes, the appearance of macroturbulent line broadening, and the formation of clumps in the wind. To investigate a possible interplay between pulsations and mass-loss, we carried out an observational campaign of the supergiant 55 Cyg over a period of five years to search for photospheric activity and cyclic mass-loss variability in the stellar wind. We modeled the H, He I, Si II and Si III lines using the nonlocal thermal equilibrium atmosphere code FASTWIND and derived the photospheric and wind parameters. In addition, we searched for variability in the intensity and radial velocity of photospheric lines and performed a moment analysis of the line profiles to derive frequencies and amplitudes of the variations. The Halpha line varies with time in both intensity and shape, displaying various types of profiles: P Cygni, pure emission, almost complete absence, and double or multiple peaked. The star undergoes episodes of variable mass-loss rates that change by a factor of 1.7-2 on different timescales. We also observe changes in the ionization rate of Si II and determine a multiperiodic oscillation in the He I absorption lines, with periods ranging from a few hours to 22.5 days. We interpret the photospheric line variations in terms of oscillations in p-, g-, and strange modes. We suggest that these pulsations can lead to phases of enhanced mass loss. Furthermore, they can mislead the determination of the stellar rotation. We classify the star as a post-red supergiant, belonging to the group of alpha Cyg variables.

B[e] supergiants (B[e]SGs) are evolved massive stars in a short-lived transition phase. During this phase, these objects eject large amounts of material, which accumulates in a circumstellar disk-like structure. The expelled material is typically dense and cool, providing the cradle for molecule and dust condensation and for a rich, ongoing chemistry. Very little is known about the chemical composition of these disks, beyond the emission from dust and CO revolving around the star on Keplerian orbits. As massive stars preserve an oxygen-rich surface composition throughout their life, other oxygen-based molecules can be expected to form. As SiO is the second most stable oxygen compound, we initiated an observing campaign to search for first-overtone SiO emission bands. We obtained high-resolution near-infrared L-band spectra for a sample of Galactic B[e]SGs with reported CO band emission. We clearly detect emission from the SiO first-overtone bands in CPD-52 9243 and indications for faint emission in HD 62623, HD 327083, and CPD-57 2874. From model fits, we find that in all these stars the SiO bands are rotationally broadened with a velocity lower than observed in the CO band forming regions, suggesting that SiO forms at larger distances from the star. Hence, searching for and analyzing these bands is crucial for studying the structure and kinematics of circumstellar disks, because they trace complementary regions to the CO band formation zone. Moreover, since SiO molecules are the building blocks for silicate dust, their study might provide insight in the early stage of dust formation.

Accurate massloss rate estimates are crucial keys to study wind properties of massive stars and test different evolutionary scenarios. From a theoretical point of view, this implies to solve a complex set of differential equations in which the radiation field and the hydrodynamics are strong coupled. The use of analytical expression to represent the radiation force and the solution of the equation of motion have many advantages over numerical integrations. Therefore, in this work, we present an analytical expression as solution of the equation of motion for radiation driven winds, in terms of the force multipliers parameters. This analytical expression is obtained by employing the line acceleration expression given by Villata (1992) and the methodology proposed by Mueller and Vink (2008). On the other hand, we find useful relationships to determine the parameters for the line acceleration given by Mueller and Vink (2008) in terms of the force multiplier parameters.

The evolutionary phase of B[e] stars is difficult to establish due to the uncertainties in their fundamental parameters. For instance, possible classifications for the Galactic B[e] star MWC 137 include pre-main-sequence and post-main-sequence phases, with a large range in luminosity. Our goal is to clarify the evolutionary stage of this peculiar object, and to study the CO molecular component of its circumstellar medium. To this purpose, we modeled the CO molecular bands using high-resolution K-band spectra. We find that MWC 137 is surrounded by a detached cool (T = 1900 +-100 K) and dense (N = (3 +- 1) x 10^21 cm^(-2)) ring of CO gas orbiting the star with a rotational velocity, projected to the line of sight, of 84 +- 2 km/s. We also find that the molecular gas is enriched in the isotope 13C, excluding the classification of the star as a Herbig Be. The observed isotopic abundance ratio (12C / 13C = 25 +- 2) derived from our modeling is compatible with a proto-PN, main-sequence or supergiant evolutionary phase. However, based on some observable characteristics of MWC 137, we propose that the supergiant scenario seems to be the most plausible. Hence, we suggest that MWC 137 could be in an extremely short-lived phase, evolving from a B[e] supergiant to a blue supergiant with a bipolar ring nebula.

We present a mini-survey of Galactic B[e] stars mainly undertaken with the Large Binocular Telescope (LBT). B[e] stars show morphological features with hydrogen emission lines and an infrared excess, attributed to warm circumstellar dust. In general, these features are assumed to arise from dense, non-spherical, disk-forming circumstellar material in which molecules and dust can condensate. Due to the lack of reliable luminosities, the class of Galactic B[e] stars contains stars at very different stellar evolutionary phases like Herbig AeBe, supergiants or planetary nebulae. We took near-infrared long-slit K-band spectra for a sample of Galactic B[e] stars with the LBT-Luci I. Prominent spectral features, such as the Brackett gamma line and CO band heads are identified in the spectra. The analysis shows that the stars can be characterized as evolved objects. Among others we find one LBV candidate (MWC314), one supergiant B[e] candidate with 13CO (MWC137) and in two cases (MWC623 and AS 381) indications for the existence of a late-type binary companion, complementary to previous studies. For MWC84, IR spectra were taken at different epochs with LBT-Luci I and the GNIRS spectrograph at the Gemini North telescope. The new data show the disappearance of the circumstellar CO emission around this star, previously detectable over decades. Also no signs of a recent prominent eruption leading to the formation of new CO disk emission are found during 2010 and 2013.

B[e] supergiants (B[e]SGs) are transitional objects in the post-main sequence evolution of massive stars. The small number of B[e]SGs known so far in the Galaxy and the Magellanic Clouds indicates that this evolutionary phase is short. Nevertheless, the strong aspherical mass loss occurring during this phase, which leads to the formation of rings or disk-like structures, and the similarity to possible progenitors of SN1987A emphasize the importance of B[e]SGs for the dynamics of the interstellar medium as well as stellar and galactic chemical evolution. The number of objects and their mass loss behavior at different metallicities are essential ingredients for accurate predictions from stellar and galactic evolution calculations. However, B[e]SGs are not easily identified, as they share many characteristics with luminous blue variables (LBVs) in their quiescent (hot) phase. We present medium-resolution near-infrared K-band spectra for four stars in M 31, which have been assigned a hot LBV (candidate) status. Applying diagnostics that were recently developed to distinguish B[e]SGs from hot LBVs, we classify two of the objects as bonafide LBVs; one of them currently in outburst. In addition, we firmly classify the two stars 2MASS J00441709+4119273 and 2MASS J00452257+4150346 as the first B[e]SGs in M 31 based on strong CO band emission detected in their spectra, and infrared colors typical for this class of stars.

During the post-main sequence evolution massive stars pass through several short-lived phases, in which they experience enhanced mass loss in the form of clumped winds and mass ejection events of unclear origin. The discovery that stars populating the blue, luminous part of the Hertzsprung-Russell diagram can pulsate hence suggests that stellar pulsations might influence or trigger enhanced mass loss and eruptions. We present recent results for two objects in different phases: a B[e] star at the end of the main sequence and a B-type supergiant.

Massive evolved stars in transition phases, such as Luminous Blue Variables (LBVs), B[e] Supergiants (B[e]SGs), and Yellow Hypergiants (YHGs), are not well understood, and yet crucial steps in determining accurate stellar and galactic evolution models. The circumstellar environments of these stars reveal their mass-loss history, identifying clues to both their individual evolutionary status and the connection between objects of different phases. Here we present a survey of 25 such evolved massive stars (16 B[e]SGs, 6 LBVs, 2 YHGs, and 1 Peculiar Oe star), observed in the K-band with the Spectrograph for INtegral Field Observation in the Near-Infrared (SINFONI; R = 4500) on the ESO VLT UT4 8 m telescope. The sample can be split into two categories based on spectral morphology: one group includes all of the B[e]SGs, the Peculiar Oe star, and two of the LBVs, while the other includes the YHGs and the rest of the LBVs. The difference in LBV spectral appearance is due to some objects being in a quiescent phase and some objects being in an active or outburst phase. CO emission features are found in 13 of our targets, with first time detections for MWC 137, LHA 120-S 35, and LHA 115-S 65. From model fits to the CO band heads, the emitting regions appear to be detached from the stellar surface. Each star with ^12CO features also shows ^13CO emission, signaling an evolved nature. Based on the level of ^13C enrichment, we conclude that many of the B[e]SGs are likely in a pre-Red Supergiant phase of their evolution. There appears to be a lower luminosity limit of log L/L_solar = 5.0 below which CO is not detected. The lack of CO features in several high luminosity B[e]SGs and variability in others suggests that they may in fact be LBV candidates, strengthening the connection between these two very similar transition phases.

We are investigating the circumstellar material for a sample of B[e] stars using high spectral resolution data taken in the optical and near-infrared regions with ESO/FEROS and ESO/CRIRES spectrographs, respectively. B[e] stars are surrounded by dense disks of still unknown origin. While optical emission lines from [O I] and [Ca II] reflect the disk conditions close to the star (few stellar radii), the near-infrared data, especially the CO band emission, mirror the characteristics in the molecular part of the disk farther away from the star (several AU). Based on our high resolution spectroscopic data, we seek to derive the density and temperature structure of the disks, as well as their kinematics. This will allow us to obtain a better understanding of their structure, formation history and evolution. Here we present our preliminary results.

The appearance of the B[e] phenomenon in evolved massive stars such as B[e] supergiants is still a mystery. While these stars are generally found to have disks that are cool and dense enough for efficient molecule and dust condensation, the origin of the disk material is still unclear. We aim at studying the kinematics and origin of the disk in the eccentric binary system GG Car, whose primary component is proposed to be a B[e] supergiant. Based on medium- and high-resolution near-infrared spectra we analyzed the CO-band emission detected from GG Car. The complete CO-band structure delivers information on the density and temperature of the emitting region, and the detectable 13CO bands allow us to constrain the evolutionary phase. In addition, the kinematics of the CO gas can be extracted from the shape of the first 12CO band head. We find that the CO gas is located in a ring surrounding the eccentric binary system, and its kinematics agrees with Keplerian rotation with a velocity, projected to the line of sight, of (80\pm 1) km/s. The CO ring has a column density of (5\pm 3)x10^21 cm^-2 and a temperature of 3200\pm 500 K. In addition, the material is chemically enriched in 13CO, which agrees with the primary component being slightly evolved off the main sequence. We discuss two possible scenarios for the origin of the circumbinary disk: (i) non-conservative Roche lobe overflow, and (ii) the possibility that the progenitor of the primary component could have been a classical Be star. Neither can be firmly excluded, but for Roche lobe overflow to occur, a combination of stellar and orbital parameter extremawould be required.

Thanks to the high spatial resolution provided by long baseline interferometry, it is possible to understand the complex circumstellar geometry around stars with the B[e] phenomenon. These stars are composed by objects in different evolutionary stages, like high- and low-mass evolved stars, intermediate-mass pre-main sequence stars and symbiotic objects. However, up to now more than 50% of the confirmed B[e] stars are not well classified, being called unclassified B[e] stars. From instruments like VLTI/AMBER and VLTI/MIDI, we have identified the presence of gaseous and dusty circumstellar disks, which have provided us with some hints related to the nature of these objects. Here, we show our results for two galactic stars with the B[e] phenomenon, HD 50138 and CPD-529243, based on interferometric measurements.

LHA 115-S 18 is a very peculiar emission-line star exhibiting the B[e] phenomenon. Located in the Small Magellanic Cloud, its spectrum shows features of an extremely wide range of excitation and ionization stages, extending from highly ionized atomic lines (Si IV, C IV, He II) in the UV and optical regions to molecular emission bands of CO and TiO in the optical and IR regions. The most distinguishing spectral characteristic of LHA 115-S 18 is the high variability detected in the He II {\lambda}4686 emission line, which can be a very conspicuous or completely invisible feature. In this work, we report on another peculiarity of LHA 115-S 18. From high-resolution optical spectra taken between 2000 and 2008, we discovered the appearance and strengthening of two emission features at {\lambda}6825 \AA, and {\lambda}7082 \AA,, which we identified as Raman-scattered lines. This is the first time these lines have been detected in the spectrum of a massive luminous B[e] star. As the classification of LHA 115-S 18 is highly controversial, we discuss how the discovery of the appearance of Raman-scattered lines in this peculiar star might help us to solve this puzzle.

Molecular emission has been detected in several Magellanic Cloud B[e] supergiants. In this Letter, we report on the detection of CO band head emission in the B[e] supergiant LHA 115-S 65, and present a K-band near-infrared spectrum obtained with the Spectrograph for INtegral Field Observation in the Near-Infrared (SINFONI; R=4500) on the ESO VLT UT4 telescope. The observed molecular band head emission in S 65 is quite surprising in light of a previous non-detection by McGregor et al. 1989, as well as a high resolution (R=50000) Gemini/Phoenix spectrum of this star taken nine months earlier showing no emission. Based on analysis of the optical spectrum by Kraus et al. 2010, we suspect that the sudden appearance of molecular emission could be due to density build up in an outflowing viscous disk, as seen for Be stars. This new discovery, combined with variability in two other similar evolved massive stars, indicates an evolutionary link between B[e] supergiants and LBVs.

We calculated the influence of the limb-darkened finite disk correction factor in the theory of radiation-driven winds from massive stars. We solved the 1-D m-CAK hydrodynamical equation of rotating radiation-driven winds for all three known solutions, i.e., fast, \Omega-slow and \delta-slow. We found that for the fast solution, the mass loss rate is increased by a factor \sim 10%, while the terminal velocity is reduced about 10%, when compared with the solution using a finite disk correction factor from a uniformly bright star. For the other two slow solutions the changes are almost negligible. Although, we found that the limb darkening has no effects on the wind momentum luminosity relationship, it would affect the calculation of synthetic line profiles and the derivation of accurate wind parameters.

The theory of radiation-driven winds succeeded in describing terminal velocities and mass loss rates of massive stars. However, for A-type supergiants the standard m-CAK solution predicts values of mass loss and terminal velocity higher than the observed values. Based on the existence of a slow wind solution in fast rotating massive stars, we explore numerically the parameter space of radiation-driven flows to search for new wind solutions in slowly rotating stars, that could explain the origin of these discrepancies. We solve the 1-D hydrodynamical equation of rotating radiation-driven winds at different stellar latitudes and explore the influence of ionization's changes throughout the wind in the velocity profile. We have found that for particular sets of stellar and line-force parameters, a new slow solution exists over the entire star when the rotational speed is slow or even zero. In the case of slow rotating A-type supergiant stars the presence of this novel slow solution at all latitudes leads to mass losses and wind terminal velocities which are in agreement with the observed values. The theoretical Wind Momentum-Luminosity Relationship derived with these slow solutions shows very good agreement with the empirical relationship. In addition, the ratio between the terminal and escape velocities, which provides a simple way to predict stellar wind energy and momentum input into the interstellar medium, is also properly traced.

Since the external regions of the envelopes of rapidly rotating early-type stars are unstable to convection, a coupling may exist between the convection and the internal rotation. We explore what can be learned from spectroscopic and interferometric observations about the properties of the rotation law in the external layers of these objects. Using simple relations between the entropy and specific rotational quantities, some of which are found to be efficient at accounting for the solar differential rotation in the convective region, we derived analytical solutions that represent possible differential rotations in the envelope of early-type stars. A surface latitudinal differential rotation may not only be an external imprint of the inner rotation, but induces changes in the stellar geometry, the gravitational darkening, the aspect of spectral line profiles, and the emitted spectral energy distribution. By studying the equation of the surface of stars with non-conservative rotation laws, we conclude that objects undergo geometrical deformations that are a function of the latitudinal differential rotation able to be scrutinized both spectroscopically and by interferometry. The combination of Fourier analysis of spectral lines with model atmospheres provides independent estimates of the surface latitudinal differential rotation and the inclination angle. Models of stars at different evolutionary stages rotating with internal conservative rotation laws were calculated to show that the Roche approximation can be safely used to account for the gravitational potential. The surface temperature gradient in rapid rotators induce an acceleration to the surface angular velocity. A non-zero differential rotation parameter may indicate that the rotation is neither rigid nor shellular underneath the stellar surface.

The line driven- and rotation modulated-wind theory predicts an alternative slow solution, besides from the standard m-CAK solution, when the rotational velocity is close to the critical velocity. We study the behaviour of the winds of A-type supergiants (Asg) and show that under particular conditions, e.g., when the \delta line-force parameter is about 0.25, the slow solution could exist over the whole star, even for the cases when the rotational speed is slow or zero. We discuss density and velocity profiles as well as possible observational counterparts.

We have incorporated the oblate distortion of the shape of the star due to the stellar rotation, which modifies the finite disk correction factor (f_D) in the m-CAK hydrodynamical model. We implement a simplified version for the f_D allowing us to solve numerically the non-linear m- CAK momentum equation.We solve this model for a classical Be star in the polar and equatorial directions. The star's oblateness modifies the polar wind, which is now much faster than the spherical one, mainly because the wind receives radiation from a larger (than the spherical) stellar surface. In the equatorial direction we obtain slow solutions, which are even slower and denser than the spherical ones. For the case when the stellar rotational velocity is about the critical velocity, the most remarkable result of our calculations is that the density contrast between the equatorial density and the polar one, is about 100. This result could explain a long-standing problem on Be stars.

Effective temperatures of early-type supergiants are important to test stellar atmosphere- and internal structure-models of massive and intermediate mass objects at different evolutionary phases. However, these Teff values are more or less discrepant depending on the method used to determine them. We aim to obtain a new calibration of the Teff parameter for early-type supergiants as a function of observational quantities that are highly sensitive to the ionization balance in the photosphere and its gas pressure, independent of the interstellar extinction, and as much as possible model-independent. The observational quantities that best address our aims are the (lambda1, D) parameters of the BCD spectrophotometric system. They describe the energy distribution around the Balmer discontinuity, which is highly sensitive to Teff and logg. We perform a calibration of the (lambda1, D) parameters into Teff using effective temperatures derived with the bolometric-flux method for 217 program stars, whose individual absolute relative uncertainties are on average 0.05. The final comparison of calculated with obtained Teff values in the (lambda1,D) calibration further show that the latter have total relative uncertainties, which on average are of the order of 0.05 for all spectral types and luminosity classes. The effective temperatures of OB supergiants derived in this work agree on average within some 2000 K with other determinations found in the literature, except those issued from wind-free non-LTE plane-parallel models of stellar atmospheres, which produce effective temperatures that can be overestimated by up to more than 5000 K near Teff = 25000 K.

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.

We investigate the nature of OGLE05155332-6925581, one of the brightest members of the enigmatic group of Double Periodic Variables (DPVs) recently found in the Magellanic Clouds. The modeling of archival orbital light curves (LCs), along with the analysis of the radial velocities suggest that this object is a semi--detached binary with the less massive star transferring matter to the more massive and less evolved star, in an Algol--like configuration. We find evidence for additional orbital variability and H$\alpha$ emission, likely caused by an accretion disc around the primary star. As in the case of $\beta Lyr$ the circumprimary disc seems to be more luminous than the primary, but we do not detect orbital period changes. We find that the LC follows a loop in the color--magnitude diagram during the long cycle; the system is redder when brighter and the rising phase is bluer than during decline. Infrared excess is also present. The source of the long--term periodicity is not eclipsed, indicating its circumbinary origin. Strong asymmetries, discrete absorption components (DACs) and a $\gamma$ shift are new and essential observational properties in the infrared H I lines. The DACs strength and RV follow a saw--teeth pattern during the orbital cycle. We suggest that the system experiences supercycles of mass outflow feeding a circumbinary disc. Mass exchange and mass loss could produce comparable but opposite effects in the orbital period on a long time scale, resulting in a quasi--constancy of this parameter.

We carried out low resolution spectroscopic observations in the wavelength range 3400-4700 A of 20 He-weak and 8 He-strong stars to determine their fundamental parameters by means of the Divan-Chalonge-Barbier (BCD) spectrophotometric system. For a few He-weak stars we also estimate the effective temperatures and the angular diameters by integrating absolute fluxes observed over a wide spectral range. Non-LTE model calculations are carried out to study the influence of the He/H abundance ratio on the emergent radiation of He-strong stars and on their Teff determination. We find that the effective temperatures, surface gravities and bolometric absolute magnitudes of He-weak stars estimated with the BCD system and the integrated flux method are in good agreement between each other, and they also agree with previous determinations based on several different methods. The mean discrepancy between the visual absolute magnitudes derived using the Hipparcos parallaxes and the BCD values is on average 0.3 mag for He-weak stars, while it is 0.5 mag for He-strong stars. For He-strong stars, we note that the BCD calibration, based on stars in the solar environment, leads to overestimated values of Teff. By means of model atmosphere calculations with enhanced He/H abundance ratios we show that larger He/H ratios produce smaller BD which naturally explains the Teff overestimation. We take advantage of these calculations to introduce a method to estimate the He/H abundance ratio in He-strong stars. The BD of HD 37479 suggests that the Teff of this star remains fairly constant as the star spectrum undergoes changes in the intensity of H and He absorption lines. Data for the He-strong star HD 66765 are reported for the first time.

We show that Be stars belong to a high velocity tail of a single B-type star rotational velocity distribution in the MS. This implies that: 1) the number fraction N(Be)/N(Be+B) is independent of the mass; 2) Bn stars having ZAMS rotational velocities higher than a given limit might become Be stars.

The effects of rapid rotation and bi--stability upon the density contrast between the equatorial and polar directions of a B[e] supergiant are investigated. Based on a new slow solution for different high rotational radiation--driven winds and the fact that bi--stability allows a change in the line--force parameters ($\alpha$, $k$, and $\delta$), the equatorial densities are about $10^2$--$10^3$ times higher than the polar ones. These values are in qualitative agreement with the observations. This calculation also permits to obtain the aperture angle of the disk.

Galactic field Be stars were studied by taking into account the effects induced by the fast rotation on their fundamental parameters. Fractional ages $\tau/\tau\_{\rm MS}$ ($\tau\_{\rm MS}$ = time spent in the MS) against stellar mass reveal that: a) Be stars spread over the whole interval $0<\tau/\tau\_{\rm MS}<1$; b) the Be phenomenon in massive stars ($M>12M\_{\odot}$) is present at smaller age ratios than for less massive stars ($M<12M\_{\odot}$); c) there is a lack of Be stars with $M<7M\_{\odot}$ in the first half of the MS. Low mass fast rotators ($M<7M\_{\odot}$), called Bn stars, could be "becoming" Be stars.

We study FeII and Balmer hydrogen emission lines observed simultaneously of 18 southern Be stars. We use the self-absorption-curve method (SAC) to determine the optical depth regime of FeII emission lines and to derive first insights on the physical properties of their forming regions.

A sample of 97 galactic field Be stars were studied by taking into account the effects induced by the fast rotation on their fundamental parameters. All program stars were observed in the BCD spectrophotometric system in order to minimize the perturbations produced by the circumstellar environment on the spectral photospheric signatures. This is one of the first attempts at determining stellar masses and ages by simultaneously using model atmospheres and evolutionary tracks, both calculated for rotating objects. The stellar ages ($\tau$) normalized to the respective inferred time that each rotating star can spend in the main sequence phase ($\tau\_{\rm MS}$) reveal a mass-dependent trend. This trend shows that: a) there are Be stars spread over the whole interval $0 \la \tau/\tau\_{\rm MS} \la 1$ of the main sequence evolutionary phase; b) the distribution of points in the ($\tau/\tau\_{\rm MS},M/M\_{\odot}$) diagram indicates that in massive stars ($M \ga 12M\_{\odot}$) the Be phenomenon is present at smaller $\tau/\tau\_{\rm MS}$ age ratios than for less massive stars ($M \la 12M\_{\odot}$). This distribution can be due to: $i$) higher mass-loss rates in massive objets, which can act to reduce the surface fast rotation; $ii$) circulation time scales to transport angular momentum from the core to the surface, which are longer the lower the stellar mass.

The effects of rapid rotation and bi-stability upon the density contrast between the equatorial and polar directions of a B[e] supergiant are re-investigated. Based upon a new slow solution for different high rotational radiation driven winds (Cur\'e 2004) and the fact that bi--stability allows a change in the line--force parameters ($\alpha$, $k$, and $\delta$), the equatorial densities are about $10^2$--$10^4$ times higher than the polar ones. These values are in qualitative agreement with the observations.