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.

Through the modelling of near-infrared photometry of star plus disk systems with the codes BEDISK/BERAY, we successfully describe the WISE photometric characteristics of Be stars in five young open clusters, NGC 663, NGC 869, NGC 884, NGC 3766 and NGC 4755, broadly studied in the literature. WISE photometry allows previously known Be stars to be detected and to find new Be candidates which could be confirmed spectroscopically. The location of Be stars in the WISE colour-magnitude diagram, separates them in two groups; active (Be stars hosting a developed circumstellar disk) and quiescent objects (Be stars in a diskless phase), and this way, we can explore how often stars are observed in these different stages. The variability observed in most active variable Be stars is compatible with a disk dissipation phase. We find that 50 percent of Be stars in the studied open clusters are in an active phase. We can interpret this as Be stars having a developed circumstellar disk half of the time. The location of Be stars with a developed disk in the CMD require mass loss rates in agreement with values recently reported in the literature. For these objects, we expect to have a tight relation between the equivalent width of the H{\alpha} line and the mass of the disk, if the inclination is known. Also, near-infrared photometry of Be stars in stellar clusters has the potential of being useful to test whether there is a preferential viewing angle.

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.

The post main-sequence evolution of massive stars is very sensitive to many parameters of the stellar models. Key parameters are the mixing processes, the metallicity, the mass-loss rate and the effect of a close companion. We study how the red supergiant lifetimes, the tracks in the Hertzsprung-Russel diagram (HRD), the positions in this diagram of the pre-supernova progenitor as well as the structure of the stars at that time change for various mass-loss rates during the red supergiant phase (RSG), and for two different initial rotation velocities. The surface abundances of RSGs are much more sensitive to rotation than to the mass-loss rates during that phase. A change of the RSG mass-loss rate has a strong impact on the RSG lifetimes and therefore on the luminosity function of RSGs. At solar metallicity, the enhanced mass-loss rate models do produce significant changes on the populations of blue, yellow and red supergiants. When extended blue loops or blue ward excursions are produced by enhanced mass-loss, the models predict that a majority of blue (yellow) supergiants are post RSG objects. These post RSG stars are predicted to show much smaller surface rotational velocities than similar blue supergiants on their first crossing of the HR gap. The position in the HRD of the end point of the evolution depends on the mass of the hydrogen envelope. More precisely, whenever, at the pre-supernova stage, the H-rich envelope contains more than about 5\% of the initial mass, the star is a red supergiant, and whenever the H-rich envelope contains less than 1\% of the total mass the star is a blue supergiant. For intermediate situations, intermediate colors/effective temperatures are obtained. Yellow progenitors for core collapse supernovae can be explained by the enhanced mass-loss rate models, while the red progenitors are better fitted by the standard mass-loss rate models.

The database of intermediate mass rotating stellar models presented in the past years by the Geneva Stellar Evolution Group can be used to build synthetic stellar populations that fully account for the effects of stellar rotation. However, up to now we still lacked stellar evolutionary tracks that rotate close to the critical limit during the whole MS phase. This occurs because the flat internal profile of rotation imposed at the Zero-Age MS (ZAMS) is modified by the action of meridional currents immediately after the ZAMS, causing the surface rotational velocity to decrease abruptly until it reaches a quasi-stationary state. We compute stellar models with non-solid rotation at the ZAMS to obtain stellar evolutionary tracks with a larger content of angular momentum, that attain rotational equatorial velocities close to the critical limit throughout their MS phase. These models have a longer MS lifetime and a higher surface chemical enrichment already at the end of the MS, particularly at Z=0.002. Stellar models with solid rotation at the ZAMS adequately represent the overall characteristics and evolution of differentially rotating models of identical angular momentum content, but with a lower initial surface rotational velocity. For these models we recommend to use as the initial rotational rate the values derived once the quasi-stationary state is reached, after the abrupt decrease in surface velocity. Interestingly, the initial equatorial rotational velocities are virtually metallicity independent for the stellar models we computed with the same mass and angular momentum content at the ZAMS. If, as some observational evidence indicates, B-type stars at Z=0.002 rotate with a higher equatorial velocity at the ZAMS than stars with Z=0.014, our finding would indicate that the angular momentum content of B-type stars in the SMC is higher than their Galactic counterparts.

(abridged) We present a new tool for building synthetic colour-magnitude diagrams of coeval stellar populations. We study, from a theoretical point of view, the impact of axial rotation of stars on various observed properties of single-aged stellar populations: magnitude at the turnoff, photometric properties of evolved stars, surface velocities, surface abundances, and the impact of rotation on the age determination of clusters by an isochrone fitting. Stellar models for different initial masses, metallicities, and zero-age main sequence (ZAMS) rotational velocities are used for building interpolated stellar tracks, isochrones, and synthetic clusters for various ages and metallicities. The synthetic populations account for the effects of the initial distribution of the rotational velocities on the ZAMS, the impact of the inclination angle and the effects of gravity and limb darkening, unresolved binaries and photometric errors. Interpolated tracks, isochrones, and synthetic clusters can be computed through a public web interface. For clusters with a metallicity in the range [0.002, 0.014] and an age between 30 Myr and 1 Gyr, the fraction of fast rotators on the main sequence (MS) band is the largest just below the turnoff. This is a natural consequence of the increase in the MS lifetime due to rotation. The fraction of fast rotators one magnitude below the turnoff also increases with the age of the cluster between 30 Myr and 1 Gyr. There is an increase in the fraction of enriched stars when the metallicity decreases. We show that the use of isochrones computed from rotating stellar models with an initial rotation that is representative of the average initial rotation of the stars in clusters provides a reasonable estimate of the age, even though stars in a real cluster did not start their evolution with an identical initial rotation.

(shortened) We provide a grid of single star models covering a mass range from 0.8 to 120 Msun with an initial metallicity Z = 0.002 with and without rotation. We discuss the impact of a change in the metallicity by comparing the current tracks with models computed with exactly the same physical ingredients but with a metallicity Z = 0.014 (solar). We show that the width of the main-sequence (MS) band in the upper part of the Hertzsprung-Russell diagram (HRD), for luminosity above log(L/Lsun) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. We confirm, but here for the first time on the whole mass range, that surface enrichments are stronger at low metallicity provided that comparisons are made for equivalent initial mass, rotation and evolutionary stage. We show that the enhancement factor due to a lowering of the metallicity (all other factors kept constant) increases when the initial mass decreases. Present models predict an upper luminosity for the red supergiants (RSG) of log (L/Lsun) around 5.5 at Z = 0.002 in agreement with the observed upper limit of RSG in the Small Magellanic Cloud. We show that models using shear diffusion coefficient calibrated to reproduce the surface enrichments observed for MS B-type stars at Z = 0.014 can also reproduce the stronger enrichments observed at low metallicity. In the framework of the present models, we discuss the factors governing the timescale of the first crossing of the Hertzsprung gap after the MS phase. We show that any process favouring a deep localisation of the H-burning shell (steep gradient at the border of the H-burning convective core, low CNO content) and/or the low opacity of the H-rich envelope favour a blue position in the HRD for the whole or at least a significant fraction of the core He-burning phase.

We study how tides in a binary system induce some specific internal shear mixing, able to substantially modify the evolution of close binaries prior to mass transfer. We construct numerical models accounting for tidal interactions, meridional circulation, transport of angular momentum, shears and horizontal turbulence and consider a variety of orbital periods and initial rotation velocities. Depending on orbital periods and rotation velocities, tidal effects may spin down (spin down Case) or spin up (spin up Case) the axial rotation. In both cases, tides may induce a large internal differential rotation. The resulting tidally induced shear mixing (TISM) is so efficient that the internal distributions of angular velocity and chemical elements are greatly influenced. The evolutionary tracks are modified, and in both cases of spin down and spin up, large amounts of nitrogen can be transported to the stellar surfaces before any binary mass transfer. Meridional circulation, when properly treated as an advection, always tends to counteract the tidal interaction, tending to spin up the surface when it is braked down and vice versa. As a consequence, the times needed for the axial angular velocity to become equal to the orbital angular velocity may be larger than given by typical synchronization timescales. Also, due to meridional circulation some differential rotation remains in tidally locked binary systems.

B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars. In this paper, we provide stellar models covering the mass range between 1.7 to 15 Msun, which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence. We used the Geneva stellar-evolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 Msun present at the beginning of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. The values presently obtained are in agreement with the very few values obtained for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 Msun with a mean velocity on the MS of the order of 150 km/s show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.

Even though it is broadly accepted that single Be stars are rapidly rotating stars surrounded by a flat rotating circumstellar disk, there is still a debate about how fast these stars rotate and also about the mechanisms involved in the angular-momentum and mass input in the disk. We study the properties of stars that rotate near their critical-rotation rate and investigate the properties of the disks formed by equatorial mass ejections. We used the most recent Geneva stellar evolutionary tracks for rapidly rotating stars that reach the critical limit and used a simple model for the disk structure. We obtain that for a 9 Msun star at solar metallicity, the minimum average velocity during the Main Sequence phase to reach the critical velocity is around 330 km/s, whereas it would be 390 km/s at the metallicity of the Small Magellanic Cloud (SMC). Red giants or supergiants originating from very rapid rotators rotate six times faster and show N/C ratios three times higher than those originating from slowly rotating stars. This difference becomes stronger at lower metallicity. It might therefore be very interesting to study the red giants in clusters that show a large number of Be stars on the MS band. On the basis of our single-star models, we show that the observed Be-star fraction with cluster age is compatible with the existence of a temperature-dependent lower limit in the velocity rate required for a star to become a Be star. The mass, extension, and diffusion time of the disks produced when the star is losing mass at the critical velocity, obtained from simple parametrized expressions, are not too far from those estimated for disks around Be-type stars. At a given metallicity, the mass and the extension of the disk increase with the initial mass and with age on the MS phase. Denser disks are expected in low-metallicity regions.

We review the significant role played by red supergiants (RSGs) in stellar populations, and some challenges and questions they raise for theoretical stellar evolution. We present how metallicity and rotation modify the way stars go to the red part of the Hertzsprung- Russell diagram or come back from it, and how RSGs might keep a trace of their main-sequence evolution. We compare theoretical popu- lation ratios with observed ones.

Mass-loss rates are one of the most relevant parameters determining the evolution of massive stars. In particular, the rates at which the star loses mass during the red-supergiant (RSG) phase is the least constrained by the observations or theory. In this paper, we show how the mass loss during the RSG phase affects the later evolution of the star, as well as the final type of supernova towards which it leads. We also discuss some possibilities to discriminate between blue stars that went through a RSG phase and those which remained in the blue part of the Hertzsprung-Russell diagram.

[abridged] Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z=0.014, spanning a wide mass range from 0.8 to 120 Msun. For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the ZAMS with a rotation rate Vini/Vcrit=0.4. The evolution is computed until the end of the central carbon-burning phase, the early AGB phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models. [...]

We argue that NSV 11749, an eruption observed in the early twentieth century, was a rare event known as "very late thermal pulse" (VLTP). To support our argument we compare the lightcurve of NSV 11749 with those of the two bonafide VLTP objects known to date, V4334 Sgr and V605 Aql, and with those predicted by state of the art stellar evolution models. Next, we explore the IPHAS and 2MASS catalogues for possible counterparts of the eruption. Our analysis shows that the VLTP scenario outperforms all other proposed scenarios as an explanation of NSV 11749. We identify an IPHAS/2MASS source at the eruption location of NSV 11749. The derived colors suggest that the object is not enshrouded in a thick dust shell as V605 Aql and V4334 Sgr. Also the absence of an apparent planetary nebula (PN) at the eruption location suggests differences with known VLTP objects which might be linked to the intensity of the eruption and the mass of the object. Further exploration of this source and scenario seems desirable. If NSV 11749 was a born again star, it would be the third event of its kind to have been observed and will strongly help us to increase our understanding on the later stages of stellar evolution and violent reactive convective burning.

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.

The Be star phenomenon is related to fast rotation, although the cause of this fast rotation is not yet clearly established. The basic effects of fast rotation on the stellar structure are reviewed: oblateness, mixing, anisotropic winds. The processes governing the evolution of the equatorial velocity of a single star (transport mechanisms and mass loss) are presented, as well as their metallicity dependence. The theoretical results are compared to observations of B and Be stars in the Galaxy and the Magellanic Clouds.

Be stars are thought to be fast rotating stars surrounded by an equatorial disc. The formation, structure and evolution of the disc are still not well understood. In the frame of single star models, it is expected that the surface of an initially fast rotating star can reach its keplerian velocity (critical velocity). The Geneva stellar evolution code has been recently improved, in order to obtain some estimates of the total mass loss and of the mechanical mass loss rates in the equatorial disc during the whole critical rotation phase. We present here the first results of the computation of a grid of fast rotating B stars evolving towards the Be phase, and discuss the first estimates we obtained.

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.

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.