Context. After leaving the main sequence, massive stars undergo complex evolution, which is still poorly understood. With a population of hundreds of OB stars, the starburst cluster Westerlund 1 offers an unparallelled environment to study their evolutionary tracks. Aims. We characterise a large sample of evolved OB stars in the cluster, with the aim of determining cluster parameters and place stars in an evolutionary sequence. Methods. We used the FORS2 instrument on the VLT to obtain intermediate-resolution spectroscopy over the range 5800-9000Å of about a hundred stars selected as likely members of the cluster based on their photometry. We developed criteria for their spectral classification using only spectral features in the range observed. We discuss these criteria, useful for spectral classification of earlytype stars in the GAIA spectral region, in the appendix. Using these criteria, we obtain spectral classifications, probably accurate to one subtype, for 57 objects, most of which had no previous classification or a generic classification. Results. We identify more than 50 objects as OB supergiants. We find almost 30 luminous early-B supergiants and a number of less luminous late-O supergiants. In addition, we find a few mid B supergiants with very high luminosity, some of them displaying signs of heavy mass loss. All these stars form a sequence compatible with theoretical evolutionary tracks. In addition, two early B supergiants also show indication of heavy mass loss and may represent the evolutionary phase immediately prior to the Wolf-Rayet stage. We investigate cluster properties using the spectral types and existing photometry. We find that the reddening law to the cluster does not deviate strongly from standard, even though extinction is quite variable, with an average value A V = 10.8. Though evolutionary tracks for high-mass stars are subject to large uncertainties, our data support an age of > ∼ 5 Myr and a distance d ≈ 5 kpc for Westerlund 1. Conclusions. The spectral types observed are compatible with a single burst of star formation (the age range is very unlikely to be > 1 Myr). Westerlund 1 shows its potentiality as a laboratory for massive star evolution, which can be fulfilled by detailed study of the population presented here.

Aims. Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs) of spectral type B1-4 Ia + ; Cyg OB2 #12, ζ 1 Sco, HD 190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. Methods. Quantitative analysis of UV-radio photometric and spectroscopic datasets was undertaken with a non-LTE model atmosphere code in order to derive physical parameters for comparison with apparently closely related objects, such as B supergiants (BSGs) and luminous blue variables (LBVs), and theoretical evolutionary predictions. Results. The long-term photospheric and spectroscopic datasets compiled for the early-B HGs revealed that they are remarkably stable over long periods (≥40 yrs), with the possible exception of ζ 1 Sco prior to the 20th century; in contrast to the typical excursions that characterise LBVs. Quantitative analysis of ζ 1 Sco, HD 190603 and BP Cru yielded physical properties intermediate between BSGs and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars, for initial masses in the range ∼30−60 M . Comparison of the properties of ζ 1 Sco with the stellar population of its host cluster/association NGC 6231/Sco OB1 provides further support for such an evolutionary scenario. In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (∼110 M ) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, instead traversing an O4-6Ia→O4-6Ia + →WN7-9ha pathway, the properties of Cyg OB2 #12 are therefore difficult to understand under current evolutionary paradigms. Finally, we note that as with AG Car in its cool phase, despite exceeding the HD limit, the properties of Cyg OB2 #12 imply that it lies below the Eddington limit -thus we conclude that the HD limit does not define a region of the HR diagram inherently inimical to the presence of massive stars.

Smoothed particle hydrodynamics (SPH) method has been increasingly used for simulating fluid flows, however its ability to simulate evaporating flow requires significant improvements. This paper proposes an SPH method for evaporating multiphase flows. The present SPH method can simulate the heat and mass transfers across the liquid-gas interfaces. The conservation equations of mass, momentum and energy were reformulated based on SPH, then were used to govern the fluid flow and heat transfer in both the liquid and gas phases. The continuity equation of the vapor species was employed to simulate the vapor mass fraction in the gas phase. The vapor mass fraction at the interface was predicted by the Clausius-Clapeyron correlation. A new evaporation rate was derived to predict the mass transfer from the liquid phase to the gas phase at the interface. Because of the mass transfer across the liquid-gas interface, the mass of an SPH particle was allowed to change. New particle splitting and merging techniques were developed to avoid large mass difference between SPH particles of the same phase. The proposed method was tested by simulating three problems, including the Stefan problem, evaporation of a static drop, and evaporation of a drop impacting on a hot surface. For the Stefan problem, the SPH results of the evaporation rate at the interface agreed well with the analytical solution. For drop evaporation, the SPH result was compared with the result predicted by a level-set method from literature. In the case of drop impact on a hot surface, the evolution of the shape of the drop, temperature, and vapor mass fraction were predicted.

Aims. We aim to characterise the properties of a third massive, red supergiant dominated galactic cluster. Methods. To accomplish this we utilised a combination of near/mid-IR photometry and spectroscopy to identify and classify the properties of cluster members, and statistical arguments to determine the mass of the cluster. Results. We found a total of 16 strong candidates for cluster membership, for which formal classification of a subset yields spectral types from K3-M4 Ia and luminosities between log(L/L ) ∼ 4.5-4.8 for an adopted distance of 6 ± 1 kpc. For an age in the range of 16-20 Myr, the implied mass is 2-4×10 4 M , making it one of the most massive young clusters in the Galaxy. This discovery supports the hypothesis that a significant burst of star formation occurred at the base of Scutum-Crux arm between 10-20 Myr ago, yielding a stellar complex comprising at least ∼10 5 M of stars (noting that since the cluster identification criteria rely on the presence of RSGs, we suspect that the true stellar yield will be significantly higher). We highlight the apparent absence of X-ray binaries within the star formation complex and finally, given the physical association of at least two pulsars with this region, discuss the implications of this finding for stellar evolution and the production and properties of neutron stars.