A rare early-type star revealed in the Wing of the Small Magellanic Cloud

Sk 183 is the visually-brightest star in the N90 nebula, a young star-forming region in the Wing of the Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope which reveals Sk 183 to be one of the most massive O-type stars in the SMC. Classified as an O3-type dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption which suggests a later type. We propose that Sk 183 has a composite spectrum and that it is similar to another star in the SMC, MPG 324. This brings the number of rare O2- and O3-type stars known in the whole of the SMC to a mere four. We estimate physical parameters for Sk 183 from analysis of its spectrum. For a single-star model, we estimate an effective temperature of 46+/-2 kK, a low mass-loss rate of ~10^-7 Msun yr^-1, and a spectroscopic mass of 46^+9_-8 Msun (for an adopted distance modulus of 18.7 mag to the young population in the SMC Wing). An illustrative binary model requires a slightly hotter temperature (~47.5 kK) for the primary component. In either scenario, Sk 183 is the earliest-type star known in N90 and will therefore be the dominant source of hydrogen-ionising photons. This suggests Sk 183 is the primary influence on the star formation along the inner edge of the nebula.Comment: Accepted by ApJ, 10 pages, 7 figures, v2 after proof

Mass modelling globular clusters in the Gaia era: a method comparison using mock data from an N-body simulation of M 4

As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution function (DF)-based models, Jeans models, and a grid of N-body models) by applying them to mock observations from a star-by-star N-body simulation of the GC M 4 by Heggie. The mocks mimic existing and anticipated data for GCs: surface brightness or number density profiles, local stellar mass functions, line-of-sight velocities, and Hubble Space Telescope-and Gaia-like proper motions. We discuss the successes and limitations of the methods. We find that multimass DF-based models, Jeans, and N-body models provide more accurate mass profiles compared to single-mass DF-based models. We highlight complications in fitting the kinematics in the outskirts due to energetically unbound stars associated with the cluster ('potential escapers', captured neither by truncated DF models nor by N-body models of clusters in isolation), which can be avoided with DF-based models including potential escapers, or with Jeans models. We discuss ways to account for mass segregation. For example, three-component DF-based models with freedom in their mass function are a simple alternative to avoid the biases of single-mass models (which systematically underestimate the total mass, half-mass radius, and central density), while more realistic multimass DF-based models with freedom in the remnant content represent a promising avenue to infer the total mass and the mass function of remnants

Studying the kinematics of the giant star-forming region 30 Doradus. I. The data

We present high-quality VLT-FLAMES optical spectroscopy of the nebular gas in the giant star-forming region 30 Doradus. In this paper, the first of a series, we introduce our observations and discuss the main kinematic features of 30 Dor, as revealed by the spectroscopy of the ionized gas in the region. The primary data set consists of regular grid of nebular observations, which we used to produce a spectroscopic datacube of 30 Dor, centered on the massive star cluster R136 and covering a field-of-view of 10'x10'. The main emission lines present in the datacube are from Halpha and [NII]6548,6584. The Halpha emission-line profile varies across the region from simple single-peaked emission to complex, multiple-component profiles, suggesting that different physical mechanisms are acting on the excited gas. To analyse the gas kinematics we fit Gaussian profiles to the observed Halpha features. Unexpectedly, the narrowest Halpha profile in our sample lies close to the supernova remnant 30 Dor B. We present maps of the velocity field and velocity dispersion across 30 Dor, finding five previously unclassified expanding structures. These maps highlight the kinematic richness of 30 Dor (e.g. supersonic motions), which will be analysed in future papers.Comment: 24 pages, 12 figures. Accepted for publication in Astronomy and Astrophysic

The VLT-FLAMES Tarantula Survey. VII. A low velocity dispersion for the young massive cluster R136

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km/s on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4 and 5 km/s given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ~25 km/s if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk of interpreting velocity dispersion measurements of unresolved extragalactic young massive clusters.Comment: 18 pages, 7 figures, accepted by A&

Globular cluster number density profiles using Gaia DR2

Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data are supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters with a complete density profile covering the full tidal radius (and beyond) for each cluster. We model the density profiles using a set of single-mass models ranging from King and Wilson models to generalized lowered isothermal LIMEPY models and the recently introduced SPES models, which allow for the inclusion of potential escapers. We find that both King and Wilson models are too simple to fully reproduce the density profiles, with King (Wilson) models on average underestimating (overestimating) the radial extent of the clusters. The truncation radii derived from the LIMEPY models are similar to estimates for the Jacobi radii based on the cluster masses and their orbits. We show clear correlations between structural and environmental parameters, as a function of Galactocentric radius and integrated luminosity. Notably, the recovered fraction of potential escapers correlates with cluster pericentre radius, luminosity, and cluster concentration. The ratio of half mass over Jacobi radius also correlates with both truncation parameter and PE fraction, showing the effect of Roche lobe filling

Linking the rotation of a cluster to the spins of its stars: The kinematics of NGC6791 and NGC6819 in 3D

The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum cascading down to stellar scales. In this paper, we combine Gaia data and published line of sight velocities to explore if NGC6791 and NGC6819, two of the clusters for which an alignment of stellar spins has been reported, rotate in the same plane as their stars. We find evidence for rotation in NGC6791 using both proper motions and line of sight velocities. Our estimate of the inclination angle is broadly consistent with the mean inclination that has been determined for its stars, but the uncertainties are still substantial. Our results identify NGC6791 as a promising follow-up candidate to investigate the link between cluster and stellar rotation. We find no evidence for rotation in NGC6819

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are however found in multiple systems. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary sub-populations with one another as well as with that of presumed single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the massive stars spin rates. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (\vrot) for components of 114 spectroscopic binaries in 30 Doradus. The \vrot\ values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. \vrot\ calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. The overall \vrot\ distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at \vrot < 200 kms) and a shoulder at intermediate velocities (200 < \vrot < 300 kms). The distributions of binaries and single stars however differ in two ways. First, the main peak at \vrot \sim100 kms is broader and slightly shifted toward higher spin rates in the binary distribution compared to that of the presumed-single stars. Second, the \vrot distribution of primaries lacks a significant population of stars spinning faster than 300 kms while such a population is clearly present in the single star sample.Comment: 16 pages, 16 figures, paper accepted in Astronomy & Astrophysic

MUSE narrow field mode observations of the central kinematics of M15

We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providing the largest sample of radial velocities ever obtained for the innermost regions of this system. Combined with previous observations of M15 using MUSE in wide field mode and literature data, we find that the central kinematics of M15 are complex with the rotation axis of the core of M15 offset from the rotation axis of the bulk of the cluster. While this complexity has been suggested by previous work, we confirm it at higher significance and in more detail

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.Comment: Accepted for publication in A&A, 19 pages, 14 figures, 6 tables, (74 pages appendix, 68 figures, 4 tables