The formation of the most massive stars in the Universe remains an unsolved problem. Are they able to form in relative isolation in a manner similar to the formation of solar-type stars, or do they necessarily require a clustered environment? In order to shed light on this important question, we study the origin of two very massive stars (VMS): the O2.5If*/WN6 star RFS7 ($\sim$100 $M_{\odot}$), and the O3.5If* star RFS8 ($\sim$70 $M_{\odot}$), found within $\approx$ 53 and 58 pc respectively from the Galactic massive young cluster NGC 3603, using Gaia data. RFS7 is found to exhibit motions resembling a runaway star from NGC 3603. This is now the most massive runaway star candidate known in the Milky Way. Although RFS8 also appears to move away from the cluster core, it has proper-motion values that appear inconsistent with being a runaway from NGC 3603 at the $3\sigma$ level (but with substantial uncertainties due to distance and age). Furthermore, no evidence for a bow-shock or a cluster was found surrounding RFS8 from available near-infrared photometry. In summary, whilst RFS7 is likely a runaway star from NGC 3603, making it the first VMS runaway in the Milky Way, RFS8 is an extremely young ($\sim$2 Myr) VMS, which might also be a runaway, but this would need to be established from future spectroscopic and astrometric observations, as well as precise distances. If RFS8 were still not meeting the criteria for being a runaway from NGC 3603 from such future data, this would have important ramifications for current theories of massive star formation, as well as the way the stellar initial mass function (IMF) is sampled.

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories.

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

This paper presents optical ($ugri$H$\alpha$)-infrared ($JHK$s,3.6--8.0$\mu$m) photometry, and $Gaia$ astrometry of 55 Classical T-Tauri stars (CTTS) in the star-forming region Sh 2-012, and it's central cluster NGC 6383. The sample was identified based on photometric H$\alpha$ emission line widths, and has a median age of 2.8$\pm$1.6 Myr, with a mass range between 0.3-1 $M_{\odot}$. 94% of CTTS with near-infrared cross-matches fall on the near-infrared T-Tauri locus, with all stars having mid-infrared photometry exhibiting evidence for accreting circumstellar discs. CTTS are found concentrated around the central cluster NGC 6383, and towards the bright rims located at the edges of Sh 2-012. Stars across the region have similar ages, suggestive of a single burst of star formation. Mass accretion rates ($\dot{M}_{\textrm{acc}}$) estimated via H$\alpha$ and $u$-band line intensities show a scatter (0.3 dex) similar to spectroscopic studies, indicating the suitability of H$\alpha$ photometry to estimate $\dot{M}_{\textrm{acc}}$. Examining the variation of $\dot{M}_{\textrm{acc}}$ with stellar mass ($M_{\ast}$), we find a smaller intercept in the $\dot{M}_{\textrm{acc}}$-$M_{\ast}$ relation than oft-quoted in the literature, providing evidence to discriminate between competing theories of protoplanetary disc evolution.

We present high-resolution (0.17\arcsec $\times$ 0.14\arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO\,(6-5) line, and 435\um\ dust continuum emission within a $\sim$9\arcsec $\times$ 9\arcsec\ area centered on the nucleus of the galaxy NGC\,5135. NGC\,5135 is a well-studied luminous infrared galaxy that also harbors a Compton-thick active galactic nucleus (AGN). At the achieved resolution of 48 $\times$ 40\,pc, the CO\,(6-5) and dust emissions are resolved into gas "clumps" along the symmetrical dust lanes associated with the inner stellar bar. The clumps have radii between $\sim$45-180\,pc and CO\,(6-5) line widths of $\sim$60-88\,\kms. The CO\,(6-5) to dust continuum flux ratios vary among the clumps and show an increasing trend with the \FeII/Br-$\gamma$ ratios, which we interpret as evidence for supernova-driven shocked gas providing a significant contribution to the \co65\ emission. The central AGN is undetected in continuum, nor in CO\,(6-5) if its line velocity width is no less than $\sim$\,40\,\kms. We estimate that the AGN contributes at most 1\% of the integrated CO\,(6-5) flux of 512 $\pm$ 24$\,$Jy\kms\ within the ALMA field of view, which in turn accounts for $\sim$32\% of the CO\,(6-5) flux of the whole galaxy.

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

Using Gaia DR2 data, combined with OmegaCAM ground-based optical photometry from the AD-HOC survey, and detailed Radial Velocity measurements from ESO-Gaia, we analyse in detail a 10x5 deg region around the Wolf-Rayet star $\gamma^2$ Vel, including the previously known clusters Gamma Vel and NGC2547. Using clustering analysis that considers positions, proper motions and parallax, we discover 6 clusters or associations -- 4 of which appear new. Analysis of the colour-magnitude diagram for these clusters show that 4 of them formed coevally from the same molecular clouds 10 Myr ago, while NGC 2547 formed together with a newly discovered cluster 30 Myr ago. This study shows the incredible wealth of data provided by Gaia for the study of young stellar clusters.

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged]

It remains unclear whether massive star evolution is facilitated by mass loss through stellar winds only, or whether episodic mass loss during an eruptive luminous blue variable (LBV) phase is also significant. LBVs exhibit unique photometric and spectroscopic variability (termed S Dor variables). This may have tremendous implications for our understanding of the first stars, gravitational wave events, and supernovae. A key question here is whether all evolved massive stars passing through the blue supergiant phase are dormant S Dor variables transforming during a brief period, or whether LBVs are truly unique objects. By investigating the OGLE light-curves of 64 B supergiants (Bsgs) in the Small Magellanic Cloud (SMC) on a timescale of three years with a cadence of one night, the incidence of S Dor variables amongst the Bsgs population is investigated. From our sample, we find just one Bsg, AzV 261, that displays the photometric behaviour characteristic of S Doradus variables. We obtain and study a new VLT X-shooter spectrum of AzV 261 in order to investigate whether the object has changed its effective temperature over the last decade. We do not find any effective temperature variations indicating that the object is unlikely to be a LBV S Dor variable. As there is only one previous bona-fide S Dor variable known to be present in the SMC (R 40), we find the maximum duration of the LBV phase in the SMC to be at most $\sim$few 10$^3$yrs, or more likely that canonical Bsgs and S Dor LBVs are intrinsically different objects. We discuss the implications for massive star evolution in low metallicity environments, characteristic of the early Universe.

We present HK spectra of three sources located in the N66 region of the Small Magellanic Cloud. The sources display prominent stellar Br Gamma and extended H2 emission, and exhibit infrared excesses at lambda > 2 micron. Based on their spectral features, and photometric spectral energy distributions, we suggest that these sources are massive young stellar objects (mYSOs). The findings are interpreted as evidence of on-going high mass star formation in N66.

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.

NGC 796 is a massive young cluster located 59 kpc from us in the diffuse intergalactic medium of the 1/5-1/10 $Z_{\odot}$ Magellanic Bridge, allowing to probe variations in star formation and stellar evolution processes as a function of metallicity in a resolved fashion, providing a link between resolved studies of nearby solar-metallicity and unresolved distant metal-poor clusters located in high-redshift galaxies. In this paper, we present adaptive optics $gri$H$\alpha$ imaging of NGC 796 (at 0.5", which is ~0.14 pc at the cluster distance) along with optical spectroscopy of two bright members to quantify the cluster properties. Our aim is to explore if star formation and stellar evolution varies as a function of metallicity by comparing the properties of NGC 796 to higher metallicity clusters. We find from isochronal fitting of the cluster main sequence in the colour-magnitude diagram an age of 20$^{+12}_{-5}$ Myr. Based on the cluster luminosity function, we derive a top-heavy stellar initial mass function (IMF) with a slope $\alpha$ = 1.99$\pm$0.2, hinting at an metallicity and/or environmental dependence of the IMF which may lead to a top-heavy IMF in the early Universe. Study of the H$\alpha$ emission line stars reveals that Classical Be stars constitute a higher fraction of the total B-type stars when compared with similar clusters at greater metallicity, providing some support to the chemically homogeneous theory of stellar evolution. Overall, NGC 796 has a total estimated mass of 990$\pm200$ $M_{\odot}$, and a core radius of 1.4$\pm$0.3 pc which classifies it as a massive young open cluster, unique in the diffuse interstellar medium of the Magellanic Bridge.

New sensitive CO(2-1) observations of the 30 Doradus region in the Large Magellanic Cloud are presented. We identify a chain of three newly discovered molecular clouds we name KN1, KN2 and KN3 lying within 2--14 pc in projection from the young massive cluster R136 in 30 Doradus. Excited H$_2$ 2.12$\mu$m emission is spatially coincident with the molecular clouds, but ionized Br$\gamma$ emission is not. We interpret these observations as the tails of pillar-like structures whose ionized heads are pointing towards R136. Based on infrared photometry, we identify a new generation of stars forming within this structure.

As part of the Accretion Discs in H$\alpha$ with OmegaCAM (ADHOC) survey, we imaged in r, i and H-alpha a region of 12x8 square degrees around the Orion Nebula Cluster. Thanks to the high-quality photometry obtained, we discovered three well-separated pre-main sequences in the color-magnitude diagram. The populations are all concentrated towards the cluster's center. Although several explanations can be invoked to explain these sequences we are left with two competitive, but intriguing, scenarios: a population of unresolved binaries with an exotic mass ratio distribution or three populations with different ages. Independent high-resolution spectroscopy supports the presence of discrete episodes of star formation, each separated by about a million years. The stars from the two putative youngest populations rotate faster than the older ones, in agreement with the evolution of stellar rotation observed in pre-main sequence stars younger than 4 Myr in several star forming regions. Whatever the final explanation, our results prompt for a revised look at the formation mode and early evolution of stars in clusters.

We present a spectroscopic study of the dynamics of the ionized and neutral gas throughout the Lagoon nebula (M8), using VLT/FLAMES data from the Gaia-ESO Survey. We explore the connections between the nebular gas and the stellar population of the associated star cluster NGC6530. We characterize through spectral fitting emission lines of H-alpha, [N II] and [S II] doublets, [O III], and absorption lines of sodium D doublet, using data from the FLAMES/Giraffe and UVES spectrographs, on more than 1000 sightlines towards the entire face of the Lagoon nebula. Gas temperatures are derived from line-width comparisons, densities from the [S II] doublet ratio, and ionization parameter from H-alpha/[N II] ratio. Although doubly-peaked emission profiles are rarely found, line asymmetries often imply multiple velocity components along the line of sight. This is especially true for the sodium absorption, and for the [O III] lines. Spatial maps for density and ionization are derived, and compared to other known properties of the nebula and of its massive stars 9 Sgr, Herschel 36 and HD 165052 which are confirmed to provide most of the ionizing flux. The detailed velocity fields across the nebula show several expanding shells, related to the cluster NGC6530, the O stars 9 Sgr and Herschel 36, and the massive protostar M8East-IR. The origins of kinematical expansion and ionization of the NGC6530 shell appear to be different. We are able to put constrains on the line-of-sight (relative or absolute) distances between some of these objects and the molecular cloud. The large obscuring band running through the middle of the nebula is being compressed by both sides, which might explain its enhanced density. We also find an unexplained large-scale velocity gradient across the entire nebula. At larger distances, the transition from ionized to neutral gas is studied using the sodium lines.

We present the first extensive spectroscopic study of the global population in star clusters Trumpler~16, Trumpler~14 and Collinder~232 in the Carina Nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. In addition to the standard homogeneous Survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. We find about four hundred good candidate members ranging from OB types down to slightly sub-solar masses. About one-hundred heavily-reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically-derived temperatures for nearly 300 low-mass members allows the inference of individual extinction values, and the study of the relative placement of stars along the line of sight. We find a complex spatial structure, with definite clustering of low-mass members around the most massive stars, and spatially-variable extinction. By combining the new data with existing X-ray data we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters, and of their connection to bright and dark nebulosity, and UV sources. The identification of tens of background giants enables us also to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star-formation history of the region, with Trumpler~14 stars found to be systematically younger than stars in other sub-clusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range.

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$\alpha$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted.

We present observations from the Gaia-ESO Survey in the lines of H$\alpha$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$\alpha$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $\eta$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $\eta$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$\alpha$ are found to arise from expanding dust reflecting the $\eta$ Car spectrum.

The Gaia ESO Public Spectroscopic Survey (GES) is providing the astronomical community with high-precision measurements of many stellar parameters including radial velocities (RVs) of stars belonging to several young clusters and star-forming regions. One of the main goals of the young cluster observations is to study of their dynamical evolution and provide insight into their future, revealing if they will eventually disperse to populate the field, rather than evolve into bound open clusters. In this paper we report the analysis of the dynamical state of L1688 in the $\rho$~Ophiuchi molecular cloud using the dataset provided by the GES consortium. We performed the membership selection of the more than 300 objects observed. Using the presence of the lithium absorption and the location in the Hertzspung-Russell diagram, we identify 45 already known members and two new association members. We provide accurate RVs for all 47 confirmed members.A dynamical analysis, after accounting for unresolved binaries and errors, shows that the stellar surface population of L1688 has a velocity dispersion $\sigma \sim$1.14$\pm$0.35 km s$^{-1}$ that is consistent with being in virial equilibrium and is bound with a $\sim$80% probability. We also find a velocity gradient in the stellar surface population of $\sim$1.0 km s$^{-1}$pc$^{-1}$ in the northwest/southeast direction, which is consistent with that found for the pre-stellar dense cores, and we discuss the possibility of sequential and triggered star formation in L1688.

We estimate the accretion rates of 235 Classical T Tauri star (CTTS) candidates in the Lagoon Nebula using $ugri$H$\alpha$ photometry from the VPHAS+ survey. Our sample consists of stars displaying H$\alpha$-excess, the intensity of which is used to derive accretion rates. For a subset of 87 stars, the intensity of the $u$-band excess is also used to estimate accretion rates. We find the mean variation in accretion rates measured using H$\alpha$ and $u$-band intensities to be $\sim$ 0.17 dex, agreeing with previous estimates (0.04-0.4 dex) but for a much larger sample. The spatial distribution of CTTS align with the location of protostars and molecular gas suggesting that they retain an imprint of the natal gas fragmentation process. Strong accretors are concentrated spatially, while weak accretors are more distributed. Our results do not support the sequential star forming processes suggested in the literature.

O and early B stars are at the apex of galactic ecology, but in the Milky Way, only a minority of them may yet have been identified. We present the results of a pilot study to select and parametrise OB star candidates in the Southern Galactic plane, down to a limiting magnitude of $g=20$. A 2 square-degree field capturing the Carina Arm around the young massive star cluster, Westerlund 2, is examined. The confirmed OB stars in this cluster are used to validate our identification method, based on selection from the $(u-g, g-r)$ diagram for the region. Our Markov Chain Monte Carlo fitting method combines VPHAS+ $u, g, r, i$ with published $J, H, K$ photometry in order to derive posterior probability distributions of the stellar parameters $\log(\rm T_{\rm eff})$ and distance modulus, together with the reddening parameters $A_0$ and $R_V$. The stellar parameters are sufficient to confirm OB status while the reddening parameters are determined to a precision of $\sigma(A_0)\sim0.09$ and $\sigma(R_V)\sim0.08$. There are 489 objects that fit well as new OB candidates, earlier than $\sim$B2. This total includes 74 probable massive O stars, 5 likely blue supergiants and 32 reddened subdwarfs. This increases the number of previously known and candidate OB stars in the region by nearly a factor of 10. Most of the new objects are likely to be at distances between 3 and 6 kpc. We have confirmed the results of previous studies that, at these longer distances, these sight lines require non-standard reddening laws with $3.5<R_V<4$.

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

We present optical spectra of pre-main sequence (PMS) candidates around the H$\alpha$ region taken with the Southern African Large Telescope, SALT, in the low metallicity ($Z$) Galactic region Sh 2-284, which includes the open cluster Dolidze 25 with an atypical low metallicity of $Z$ $\sim$ 1/5 $Z_{\odot}$. It has been suggested on the basis of both theory and observations that PMS mass-accretion rates, $\dot M_{\rm{acc}}$, are a function of $Z$. We present the first sample of spectroscopic estimates of mass-accretion rates for PMS stars in any low-$Z$ star-forming region. Our data-set was enlarged with literature data of H$\alpha$ emission in intermediate-resolution R-band spectroscopy. Our total sample includes 24 objects spanning a mass range between 1 - 2 $M_{\odot}$ and with a median age of approximately 3.5 Myr. The vast majority (21 out of 24) show evidence for a circumstellar disk on the basis of 2MASS and Spitzer infrared photometry. We find $\dot M_{\rm{acc}}$ in the 1 - 2 $M_{\odot}$ interval to depend quasi-quadratically on stellar mass, with $\dot M_{\rm{acc}}$ $\propto$ $M_{\ast}^{2.4\,\pm\,0.35}$, and inversely with stellar age $\dot M_{\rm{acc}}$ $\propto$ $t_{\ast}^{-\,0.7\,\pm\,0.4}$. Furthermore, we compare our spectroscopic $\dot M_{\rm{acc}}$ measurements with solar $Z$ Galactic PMS stars in the same mass range, but, surprisingly find no evidence for a systematic change in $\dot M_{\rm{acc}}$ with $Z$. We show that literature accretion-rate studies are influenced by detection limits, and we suggest that $\dot M_{\rm{acc}}$ may be controlled by factors other than $Z_{\ast}$, $M_{\ast}$, and age.

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (log $L$ = 4.04 $\pm$ 0.25 $L_{\odot}$) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to measure mass accretion rates at low metallicity spectroscopically, to understand the effect of metallicity on accretion rates.