Spectroscopic analyses of hydrogen-rich WN 5-6 stars within the young star clusters NGC 3603 and R136 are presented, using archival Hubble Space Telescope and Very Large Telescope spectroscopy, and high spatial resolution near-IR photometry, including Multi-Conjugate Adaptive Optics Demonstrator (MAD) imaging of R136. We derive high stellar temperatures for the WN stars in NGC 3603 (T * ∼ 42 ± 2 kK) and R136 (T * ∼ 53 ± 3 kK) plus clumping-corrected mass-loss rates of 2-5 × 10 −5 M yr −1 which closely agree with theoretical predictions from Vink et al. These stars make a disproportionate contribution to the global ionizing and mechanical wind power budget of their host clusters. Indeed, R136a1 alone supplies ∼7 per cent of the ionizing flux of the entire 30 Doradus region. Comparisons with stellar models calculated for the main-sequence evolution of 85-500 M accounting for rotation suggest ages of ∼1.5 Myr and initial masses in the range 105-170 M for three systems in NGC 3603, plus 165-320 M for four stars in R136. Our high stellar masses are supported by consistent spectroscopic and dynamical mass determinations for the components of NGC 3603A1. We consider the predicted X-ray luminosity of the R136 stars if they were close, colliding wind binaries. R136c is consistent with a colliding wind binary system. However, short period, colliding wind systems are excluded for R136a WN stars if mass ratios are of order unity. Widely separated systems would have been expected to harden owing to early dynamical encounters with other massive stars within such a high-density environment. From simulated star clusters, whose constituents are randomly sampled from the Kroupa initial mass function, both NGC 3603 and R136 are consistent with an tentative upper mass limit of ∼300 M . The Arches cluster is either too old to be used to diagnose the upper mass limit, exhibits a deficiency of very massive stars, or more likely stellar masses have been underestimated -initial masses for the most luminous stars in the Arches cluster approach 200 M according to contemporary stellar and photometric results. The potential for stars greatly exceeding 150 M within metal-poor galaxies suggests that such pair-instability supernovae could occur within the local universe, as has been claimed for SN 2007bi.

Abstract. We present new spectroscopic and photometric observations of the young Galactic open cluster Westerlund 1 (Wd 1) that reveal a unique population of massive evolved stars. We identify ∼200 cluster members and present spectroscopic classifications for ∼25% of these. We find that all stars so classified are unambiguously post-Main Sequence objects, consistent with an apparent lack of an identifiable Main Sequence in our photometric data to V ∼ 20. We are able to identify rich populations of Wolf Rayet stars, OB supergiants and short lived transitional objects. Of these, the latter group consists of both hot (Luminous Blue Variable and extreme B supergiant) and cool (Yellow Hypergiant and Red Supergiant) objects -we find that half the known Galactic population of YHGs resides within Wd 1. We obtain a mean V − M V ∼ 25 mag from the cluster Yellow Hypergiants, implying a Main Sequence turnoff at or below M V = −5 (O7 V or later). Based solely on the masses inferred for the 53 spectroscopically classified stars, we determine an absolute minimum mass of ∼1.5 × 10 3 M for Wd 1. However, considering the complete photometrically and spectroscopically selected cluster population and adopting a Kroupa IMF we infer a likely mass for Wd 1 of ∼10 5 M , noting that inevitable source confusion and incompleteness are likely to render this an underestimate. As such, Wd 1 is the most massive compact young cluster yet identified in the Local Group, with a mass exceeding that of Galactic Centre clusters such as the Arches and Quintuplet. Indeed, the luminosity, inferred mass and compact nature of Wd 1 are comparable with those of Super Star Clusters -previously identified only in external galaxies -and is consistent with expectations for a Globular Cluster progenitor.

We examine the luminosity and dynamical mass estimates for young massive stellar clusters. For many young (<50 Myr) clusters, the luminosity and dynamical mass estimates differ by a significant amount. We explain this as being due to many young clusters being out of virial equilibrium (which is assumed in dynamical mass estimates) because the clusters are undergoing violent relaxation after expelling gas not used in star formation. We show that, if we assume that luminous mass estimates are correct (for a standard IMF), at least 50 per cent of young clusters for which dynamical masses are known are likely to be destroyed within a few 10s Myr of their formation. Even clusters which will retain a bound core may lose a large fraction of their stellar mass. We also show that the core radius and other structural parameters change significantly during the violent relaxation that follows gas expulsion and that they should be considered instantaneous values only, not necessarily reflecting the final state of the cluster. In particular we note that the increasing core radii observed in young LMC/SMC clusters can be well explained as an effect of rapid gas loss.Comment: 8 pages, 5 figures. MNRAS, in pres

We discuss the observations and theory of star cluster formation to argue that clusters form dynamically cool (subvirial) and with substructure. We then perform an ensemble of simulations of cool, clumpy (fractal) clusters and show that they often dynamically mass segregate on timescales far shorter than expected from simple models. The mass segregation comes about through the production of a short-lived, but very dense core. This shows that in clusters like the Orion Nebula Cluster the stars ≥ 4M ⊙ can dynamically mass segregate within the current age of the cluster. Therefore, the observed mass segregation in apparently dynamically young clusters need not be primordial, but could be the result of rapid and violent early dynamical evolution.

Aims. We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. Methods. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particle's heating and cooling. The method captures -at minimal computational cost -the effects of (i) the rotational and vibrational degrees of freedom of H 2 ; (ii) H 2 dissociation and H o ionisation; (iii) opacity changes due to ice mantle melting, sublimation of dust, molecular lines, H − , bound-free and free-free processes and electron scattering; (iv) external irradiation; and (v) thermal inertia. Results. We use the new method to simulate the collapse of a 1 M cloud of initially uniform density and temperature. At first, the collapse proceeds almost isothermally (T ∝ ρ 0.08 ; cf. Larson 2005, MNRAS, 359, 211). The cloud starts heating fast when the optical depth to the cloud centre reaches unity (ρ C ∼ 7 × 10 −13 g cm −3 ). The first core forms at ρ C ∼ 4 × 10 −9 g cm −3 and steadily increases in mass. When the temperature at the centre reaches T C ∼ 2000 K, molecular hydrogen starts to dissociate and the second collapse begins, leading to the formation of the second (protostellar) core. The results mimic closely the detailed calculations of Masunaga & Inutsuka (2000, ApJ, 531, 350). We also simulate (i) the collapse of a 1.2 M cloud, which initially has uniform density and temperature, (ii) the collapse of a 1.2 M rotating cloud, with an m = 2 density perturbation and uniform initial temperature, and (iii) the smoothing of temperature fluctuations in a static, uniform density sphere. In all these tests the new algorithm reproduces the results of previous authors and/or known analytic solutions. The computational cost is comparable to a standard SPH simulation with a simple barotropic equation of state. The method is easy to implement, can be applied to both particle-and grid-based codes, and handles optical depths 0 < τ < ∼ 1011 .