Issue 
A&A
Volume 545, September 2012



Article Number  L1  
Number of page(s)  5  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201219472  
Published online  28 August 2012 
The VLTFLAMES Tarantula Survey^{⋆}
VI. Evidence for rotation of the young massive cluster R136
^{1} Scottish Universities Physics Alliance (SUPA), Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK
email: vhb@roe.ac.uk
^{2} Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
^{3} UK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK
^{4} Astronomical Institute “Anton Pannekoek”, University of Amsterdam, Postbus 94249, 1090 GE, Amsterdam, The Netherlands
^{5} Excellence Cluster Universe, Technische Universität München, Boltzmannstr. 2, 85748 Garching, Germany
^{6} Instituto de Astrofísica de AndalucíaCSIC, Glorieta de la Astronomía s/n, 18008 Granada, Spain
^{7} Institute of Astronomy with NAO, Bulgarian Academy of Sciences, PO Box 136, 4700 Smoljan, Bulgaria
^{8} School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, UK
^{9} European Southern Observatory, KarlSchwarzschildStrasse 2, 87548 Garching bei München, Germany
^{10} HarvardSmithsonian CfA, 60 Garden Street, Cambridge, MA 02138, USA
^{11} Astrophysics Group, LennardJones Laboratories, Keele University, Staffordshire ST5 5BG, UK
Received: 24 April 2012
Accepted: 31 July 2012
Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multiepoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud obtained as part of the VLTFLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single Otype stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ~3 km s^{1}. When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ~20% of the kinetic energy in rotation.
Key words: galaxies: star clusters: individual: R136 / Magellanic Clouds / stars: kinematics and dynamics / globular clusters: general
© ESO, 2012
1. Introduction
Despite their spherical shape, Milky Way globular clusters (GCs) rotate with amplitudes up to half the 1D velocity dispersion (0 ≲ V_{rot}sini /σ_{1D} ≲ 0.5, e.g. Meylan & Heggie 1997), so their amount of rotational energy is typically not dominant but also not negligible. Numerical studies have shown that this rotation has an important influence on star clusters by accelerating their dynamical evolution, for example by speeding up the collapse of the core through the gravogyro instability or by significantly increasing the escape rate for clusters in a strong tidal field (e.g. Einsel & Spurzem 1999; Kim et al. 2002; Ernst et al. 2007).
Most of the rotation signatures are found through radial velocity (RV) studies, but rotation has also been confirmed in the plane of the sky for ω Centauri and 47 Tucanae (van Leeuwen et al. 2000; Anderson & van der Marel 2010, respectively). RV studies are now able to measure rotational amplitudes in GCs below 1 km s^{1} and despite these precise measurements rotation has not been detected in some clusters (e.g. Lane et al. 2010), although this could also be an inclination effect.
It is unclear what the origin of the rotation is in some of these old clusters. It could be the result of the merging of two clusters (Baumgardt et al. 2003) or imprinted during the formation process. Observational input is now getting sufficiently abundant to look for correlations between rotational amplitude and other cluster properties. Bellazzini et al. (2012) report a correlation between horizontal branch (HB) morphology and V_{rot}sini /σ_{1D} in a sample of 20 GCs. Given that metallicity is the first parameter determining HB morphology, this in turn suggests a correlation between V_{rot}sini /σ_{1D} and metallicity, such that clusters with higher metallicity have greater fractions of energy in rotation. Since a higher metallicity in a gas implies a higher efficiency of energy dissipation by atomic transitions, this then hints at a significant role of dissipation in the process of cluster formation (e.g. Bekki 2010). Rotation may therefore be intimately linked to the formation of clusters.
Little is known about rotation in young clusters, partially because it is very challenging to measure accurately σ_{1D} given the high multiplicity fraction of massive stars (e.g. HénaultBrunet et al. 2012). Frenk & Fall (1982) found an ageellipticity relation for clusters in the Large Magellanic Cloud (LMC) with older clusters presenting less elongated shapes, which was interpreted as internal evolution erasing any asymmetry stemming from the violent relaxation process of the formation (although see Goodwin 1997). However, rotation and ellipticity are not necessarily equivalent. Ellipticity can be due to rotation (e.g. ω Cen; Meylan & Mayor 1986) but also to velocity anisotropy (Stephens et al. 2006; Hénon 1973), and rotating clusters can be spherical (LyndenBell 1960; Meza 2002).
Marginal evidence for rotation was found for the young (few 100 Myr) Galactic cluster GLIMPSEC01 with an amplitude of V_{rot}sini /σ_{1D} ≃ 0.2 (Davies et al. 2011). A rotational signal in the RVs was also detected in the ~100 Myr cluster NGC 1866 (Fischer et al. 1992) and in the ~50 Myr binary cluster NGC 1850 (Fischer et al. 1993), both in the LMC. To really confirm whether clusters form with a significant amount of angular momentum, we need to look for an even younger cluster. The young massive cluster (YMC) R136 in the 30 Doradus star forming region in the LMC is an ideal target to establish this. With an estimated mass of about 10^{5} M_{⊙} (Andersen et al. 2009) and its sub solar metallicity, it may at some stage resemble a typical metalrich GC as we find them in the Milky Way Bulge. With an age of less than 2 Myr (de Koter et al. 1998; Massey & Hunter 1998; Crowther et al. 2010), it is so young that any rotation needs to be attributed to the formation process, be it from merging of subclusters or directly from the angular momentum of the progenitor cloud. A rough estimate of the halfmass relaxation time (t_{rh}) of R136 can be obtained by assuming N = 10^{5} stars and a halfmass radius of 2.27 pc, which is found from multiplying the halflight radius of 1.7 pc (e.g. HénaultBrunet et al. 2012) by 4/3 (Spitzer 1987). Following the formula of Spitzer & Hart (1971), we obtain t_{rh} ≃ 366 Myr, so relaxation would not have had time to erase the original signature of rotation.
Here we report on evidence for rotation of R136 deduced from RV measurements of massive stars obtained as part of the VLTFLAMES Tarantula Survey (VFTS; Evans et al. 2011). We briefly present the data in Sect. 2 and describe our analysis of the rotational signature in Sect. 3. We discuss the implications of the rotation of R136 for cluster evolution in Sect. 4, and present our conclusions in Sect. 5.
Fig. 1 Illustration of the positions and RVs of the stars considered in this study. Symbol sizes denote the magnitude of the stellar velocities with respect to the average cluster velocity. The solid, dotted, and dashed lines correspond to the optimal rotation axis determined for models with a constant rotational velocity, a constant rotation rate, and a more realistic rotation curve (see Sect. 3), respectively. 
2. Data
The data used in this work consist of RV measurements and their uncertainties for 36 apparently single Otype stars within a projected radius of 10 pc from the centre of R136 (adopted here as the position of the star R136a1: α = 5^{h}38^{m}4239, δ = −69°06′0291, J2000). These are based on observations from the VFTS, for which at least five epochs were obtained for five different pointings of the FLAMESARGUS integralfield unit in the central arcminute of 30 Dor, in addition to nine configurations of the Medusa fibrefeed to the Giraffe spectrograph in the surrounding Tarantula Nebula over a 25 arcmin diameter fieldofview (Evans et al. 2011). The RV and variability analysis, based on Gaussian fitting of selected helium stellar absorption lines, is presented in HénaultBrunet et al. (2012) for the ARGUS data and in Sana et al. (2012) for the Medusa data. We retain here only the RVs of stars showing no significant variability throughout all epochs, and we apply a 3σ clipping centered on the mean velocity of the cluster, with σ = 6 km s^{1}, the observed lineofsight velocity dispersion of the cluster (i.e. before correction of the velocity dispersion for undetected binaries, see HénaultBrunet et al. 2012). This yields a total of 16 ARGUS sources, all within 5 pc from the centre, and 20 Medusa sources, all between 5 and 10 pc except two between 4 and 5 pc. The identification numbers, coordinates, spectral types, and RVs of these stars are listed in Table 3 of HénaultBrunet et al. (2012). In Fig. 1, we schematically present their positions and distribution of RVs with respect to the mean RV of the cluster.
3. Analysis
To look for the signature of rotation, we fit rotating models to our set of measured RVs by maximizing the logarithm of the likelihood function (Bevington 1969) (1)Given our relatively small data set, we only consider three simple models: constant rotational velocity amplitude, constant rotation rate, and finally a more realistic model with solid body rotation in the inner parts peaking near 2 halfmass radii followed by a decline. For these models, the probability density function P_{i} of a measurement RV_{i} with uncertainty σ_{i} given a position angle PA_{0} for the rotation axis and a lineofsight velocity dispersion σ_{1D} can be written as (2)where V_{rot}sini is a constant for a model with fixed rotational velocity. For a model with constant rotation rate, V_{rot}sini depends on the rotation rate (Ω) and the distance from the rotation axis (X_{j}) such that [V_{rot}sini] _{j} = Ω X_{j}. For the physically motivated model, we adopt a function of the form [V_{rot}sini] _{j} = A/2·X_{j}/(1 + (X_{j}/4)^{2}), where X_{j} is the distance to the rotation axis in pc and A is the maximum rotational velocity which we assume is at 4 pc for simplicity (the halfmass radius of R136 is about 2 pc). This rotation curve captures the general behaviour seen in simulations of rotating clusters (e.g. Kim et al. 2002). Note that because the velocity dispersion profile of R136 is relatively flat (HénaultBrunet et al. 2012), this function also describes V_{rot}/σ. We define the position angle (with respect to the centre of the cluster) as increasing anticlockwise in the plane of the sky from North (PA = 0°) towards East (PA = 90°). We adopt negative rotational velocities for position angles between PA_{0} and PA_{0}+180°. For simplicity, we assume that σ_{1D} is constant across the radius range considered. This σ_{1D} is largely due to cluster dynamics, but contains a small contribution from the orbital motions of undetected binaries (1–2 km s^{1}; HénaultBrunet et al. 2012), effectively adding noise to the rotation signature we are trying to detect.
Fig. 2 Top: RV in the system of the cluster as a function of distance X from the optimal rotation axis for models with constant rotation rate (left), constant rotational velocity (centre), and the more realistic rotation curve discussed above (right). The bestfit rotation curves are shown as solid lines, and ± σ_{1D} envelopes are represented by dashed lines. Bottom: histogram of amplitudes Ω (left), V_{rot} (centre) and A (right) from 10 000 Monte Carlo simulations of RV distributions with no rotation. The optimal values for the observed RV configuration are indicated by dashed vertical lines. Confidence levels of 94.6%, 93.6% and 95.6% are found respectively for the bestfit amplitude of each model. 
The values of the parameters that maximize the likelihood function are PA_{0} = 57 ± 15°, V_{rot}sini = 2.9 ± 1.0 km s^{1} and σ_{1D} = 5.9 ± 0.8 km s^{1} for the constant rotational velocity model, PA_{0} = 44 ± 18°, Ω = 0.75 ± 0.22 Myr^{1} and σ_{1D} = 6.0 ± 1.1 km s^{1} for the constant rotation rate model, and PA_{0} = 36 ± 14°, A = 3.7 ± 1.3 km s^{1} and σ_{1D} = 5.9 ± 0.8 km s^{1} for the final model. The bestfit rotation curves are shown in Fig. 2. The uncertainties were estimated using Monte Carlo calculations on simulated data sets comparable to our measured data (i.e. with the same rotational signal). We also obtained consistent uncertainties by considering the parameter change necessary to decrease M by from its value at the maximum (Bevington 1969). The maximum value of M is slightly higher for the constant rotational velocity model compared to the two other models, but likelihood ratio tests showed that the difference is not significant. Thus, we currently cannot favour one model over the others.
To establish the significance of the detected rotational signal, we performed Monte Carlo simulations and applied the maximum likelihood method above to 10 000 random distributions of velocities (i.e. nonrotating systems). The adopted size and spatial distribution of the simulated populations were taken to be the same as the observed sample in order to be sensitive to possible biases introduced by our nonuniform spatial sampling. The velocities were drawn from a Gaussian distribution with σ determined by the observed lineofsight velocity dispersion. For each simulated star, observational noise was added based on the RV uncertainty of the observed star at the same location, again to take into account the biases introduced by our nonuniform data set. The distribution of Ω, V_{rot}sini and A from 10 000 such simulations is shown in Fig. 2. The distributions do not peak at zero because the limited number of stars and the measurement noise generally result in a nonzero amplitude. Our bestfit values are located in the right tail of these distributions and the corresponding confidence level of the detection is around 95% for all three models.
The analysis outlined above was also performed on a subsample from which supergiant candidates and stars with possible composite spectra were excluded as their RVs might be inaccurate (see HénaultBrunet et al. 2012). Very similar results were obtained, with bestfit parameters all within 4% of those previously determined and a confidence level about 1% lower because of the smaller number of measurements. As an additional check, we ran the analysis on different subsamples of apparently single Otype stars from VFTS using the measurements of Sana et al. (2012), including stars farther away from R136 than in our main sample. The analysis of these subsamples suggests that the rotation signal is dominated by the stars in the inner regions and does not result from a velocity gradient across the field on a larger scale in the surrounding OB association. For example, when considering the sample of stars between 10 and 20 pc from the centre (37 stars), we find that V_{rot}sini goes down to 1.8 km s^{1} for a constant rotational velocity model and the confidence level of the rotational signature is only 51%. Note that stars that are part of the surrounding 30 Doradus region and not members of R136 could contaminate our sample. If we assume that the outer component of the doublecomponent EFF fit (Elson et al. 1987) to the light profile of R136 (Mackey & Gilmore 2003) is due to the surrounding association, then we may expect these stars to contribute to >50% of the sample beyond 5 pc from the centre and to dilute the rotational signal in the outer regions of the cluster.
We also looked for the signature of rotation with a method commonly used in studies of GCs. This method consists of dividing the sample with a line passing through the centre of the cluster and computing the difference in the mean RV between the two subsamples of stars as a function of the position angle of the dividing line (e.g. Côté et al. 1995; Bellazzini et al. 2012). Results very similar to those reported above were obtained for the position of the rotation axis and the rotational amplitude, but the maximum likelihood method has the advantage of directly comparing the data to models without binning.
4. Discussion
Given that σ_{1D} ≃ 5 ± 1 km s^{1} for R136 after correcting for undetected binaries (HénaultBrunet et al. 2012) and that the mean rotational velocity in the radius range considered is ≃ 3 ± 1 km s^{1} for the three bestfit rotation curves shown above, our analysis implies that V_{rot}sini /σ_{1D} ≃ 0.6 ± 0.3, which is somewhat larger than what is typically observed in GCs (e.g. Meylan & Heggie 1997). We have to keep in mind that V_{rot}sin i is itself a projected quantity, and so is σ_{1D}. We would like to know how these quantities relate to the 3D quantities ⟨ v_{φ} ⟩ and v_{rms}. Here ⟨ v_{φ} ⟩ is the mean tangential velocity component in the φ direction, used as the rotational component of the velocity vector in spherical coordinates for a cluster rotating about the z axis, and v_{rms} is the rootmean square of the velocities.
We can estimate that V_{rot}/σ_{1D} is typically about 10% larger than ⟨ v_{φ} ⟩ /v_{rms}. This can be understood as follows. For a circular orbit with unit velocity and the line of sight along the orbital plane (sini = 1) the mean RV component is 2/π. For an isotropic velocity dispersion the observed component is . This means that . If sini = 1 then we would have ⟨ v_{φ} ⟩ /v_{rms} ≃ 0.9V_{rot}/σ_{1D} and a lower limit on ⟨ v_{φ} ⟩ /v_{rms} of 0.5 ± 0.2. A ratio of ⟨ v_{φ} ⟩ /v_{rms} = 0.5 implies a ratio of kinetic energy in rotation over kinetic energy in random motions of 0.25, which in turn implies that at least 20% of the total kinetic energy is in rotation. If instead we assume that sini = 2/π (the average of sini assuming a random distribution of inclination angles), we have ⟨ v_{φ} ⟩ /v_{rms} = 0.9 and 45% of the total kinetic energy in rotation. The approximate criterion for stability against nonaxisymmetric perturbations provided by Ostriker & Peebles (1973) states that the ratio of the rotational kinetic energy over the potential energy should not exceed  T_{rot}/W  < 0.14. If we assume virial equilibrium, then T = −0.5W, so for T_{rot} = 0.20 T and T_{rot} = 0.45 T this implies  T_{rot}/W  = 0.1 and 0.225, respectively. Although not a rigorous test, this suggests that a low inclination is perhaps more reasonable.
If the cluster was flattened by rotation (although recall our previous words of caution about ellipticity and rotation), we would expect to see peaks in the azimuthal density profile, with density minima coinciding with the rotation axis (i.e. PA_{0} ~ 45° and 225°) and density maxima 90° away from the rotation axis (i.e. PA_{0} ~ 135 and 315°). However, only one minimum is seen at a position angle of ~100−120° in the K_{s} and Hband azimuthal density profiles of R136 by Campbell et al. (2010).
Given the very young age of R136, the cause of the rotation needs to be looked for in the details of the formation process of the cluster. This is a short and complicated phase in which various physical processes, with their respective timescales, operate simultaneously. A logical starting point is to see whether giant molecular clouds (GMCs), the birth sites of YMCs, rotate. Rosolowsky et al. (2003) found that GMCs in M 33 have nonzero angular momentum. From simple arguments based on differential rotation in a galactic potential and self gravity we expect that the rotation of GMCs should be prograde with the orbit in the galaxy. However, Rosolowsky et al. (2003) found that 40% of the GMCs have retrograde motions, which supports a scenario in which GMCs form through both agglomeration and self gravity and the angular momentum is the result of the clumpiness of the gas (Dobbs et al. 2011). As we pointed out above, a merger is another way to give rise to rotation. Interestingly, Sabbi et al. (2012) found a dual structure in the density of low mass stars in R136 that possibly hints at a relatively recent merger event of the main core of R136 and a second clump or cluster, but only about three of our targets are located in this second clump.
To establish whether the rotation in old GCs is a remnant of their formation, we need to know if the angular momentum can survive for a Hubble time of dynamical evolution. During the evolution, angular momentum is diffused outward (Fall & Frenk 1985; Einsel & Spurzem 1999) and ultimately lost through the escape of stars with high angular momentum (Agekian 1958; Shapiro & Marchant 1976). This process operates on a relaxation time and the angular momentum reduces after a fixed number of elapsed relaxation times. The relaxation time of expanding clusters grows roughly linearly in time, which makes the number of elapsed relaxation times grow slowly, namely as a logarithm of the age. Because the majority of GCs are in this expansion phase (Gieles et al. 2011), we may expect rotational signatures to still be present after a Hubble time because they have not evolved enough. For clusters that have entered the “massloss” phase, we do not expect the rotation to survive.
5. Conclusion
We presented evidence that the young massive cluster R136 is rotating with a rotational velocity amplitude of about 3 km s^{1}, which implies that at least ~20% of its total kinetic energy is in rotation. Obviously, RV measurements of more stars in this cluster would be desirable to better populate the rotation curve and confirm the rotational signal with a confidence level higher than the current 95%. Given the young age of R136, our results suggests that star clusters may form with a significant amount of angular momentum. This will place useful constraints on models of cluster formation. We finally argued that the rotation of GGs could originate from their formation, but this is clearly a topic where more detailed numerical investigations are welcome.
Acknowledgments
We would like to thank the referee for constructive feedback. We also wish to thank Sergey Koposov for useful discussions. VHB acknowledges support from the Scottish Universities Physics Alliance (SUPA) and from the Natural Science and Engineering Research Council of Canada (NSERC). M.G. acknowledges financial support from the Royal Society. N.B. was supported by the DFG cluster of excellence “Origin and Structure of the Universe” (www.universecluster.de). J.M.A. acknowledges support from [a] the Spanish Government Ministerio de Educación y Ciencia through grants AYA201015081 and AYA201017631 and [b] the Consejería de Educación of the Junta de Andalucía through grant P08TIC4075. N.M. was supported by the Bulgarian NSF (DO 0285).
References
 Agekian, T. A. 1958, Sov. Ast., 2, 22 [NASA ADS] [Google Scholar]
 Anderson, J., & van der Marel, R. P. 2010, ApJ, 710, 1032 [NASA ADS] [CrossRef] [Google Scholar]
 Andersen, M., Zinnecker, H., Moneti, A., et al. 2009, ApJ, 707, 1347 [NASA ADS] [CrossRef] [Google Scholar]
 Baumgardt, H., Makino, J., Hut, P., McMillan, S., & Portegies Zwart, S. 2003, ApJ, 589, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Bekki, K. 2010, ApJ, 724, L99 [NASA ADS] [CrossRef] [Google Scholar]
 Bellazzini, M., Bragaglia, A., Carretta, E., et al. 2012, A&A, 538, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bevington, P. R. 1969, Data reduction and error analysis for the physical sciences (New York: McGrawHill) [Google Scholar]
 Campbell, M. A., Evans, C. J., Mackey, A. D., et al. 2010, MNRAS, 405, 421 [NASA ADS] [Google Scholar]
 Côté, P., Welch, D. L., Fischer, P., & Gebhardt, K. 1995, ApJ, 454, 788 [NASA ADS] [CrossRef] [Google Scholar]
 Crowther, P. A., Schnurr, O., Hirschi, R., et al. 2010, MNRAS, 408, 731 [NASA ADS] [CrossRef] [Google Scholar]
 Davies, B., Bastian, N., Gieles, M., et al. 2011, MNRAS, 411, 1386 [NASA ADS] [CrossRef] [Google Scholar]
 de Koter, A., Heap, S. R., & Hubeny, I. 1998, ApJ, 509, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Dobbs, C. L., Burkert, A., & Pringle, J. E. 2011, MNRAS, 417, 1318 [NASA ADS] [CrossRef] [Google Scholar]
 Einsel, C., & Spurzem, R. 1999, MNRAS, 302, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Elson, R. A. W., Fall, S. M., & Freeman, K. C. 1987, ApJ, 323, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Ernst, A., Glaschke, P., Fiestas, J., Just, A., & Spurzem, R. 2007, MNRAS, 377, 465 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Evans, C. J., Taylor, W. D., HénaultBrunet, V., et al. 2011, A&A, 530, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fall, S. M., & Frenk, C. S. 1985, in Dynamics of Star Clusters, eds. J. Goodman, & P. Hut, IAU Symp., 113, 285 [Google Scholar]
 Fischer, P., Welch, D. L., Cote, P., Mateo, M., & Madore, B. F. 1992, AJ, 103, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Fischer, P., Welch, D. L., & Mateo, M. 1993, AJ, 105, 938 [NASA ADS] [CrossRef] [Google Scholar]
 Frenk, C. S., & Fall, S. M. 1982, MNRAS, 199, 565 [NASA ADS] [Google Scholar]
 Gieles, M., Heggie, D. C., & Zhao, H. 2011, MNRAS, 413, 2509 [NASA ADS] [CrossRef] [Google Scholar]
 Goodwin, S. P. 1997, MNRAS, 286, L39 [NASA ADS] [CrossRef] [Google Scholar]
 HénaultBrunet, V., Evans, C. J., Sana, H., et al. 2012, A&A, in press, DOI 10.1051/00046361/201219471 [Google Scholar]
 Hénon, M. 1973, A&A, 24, 229 [Google Scholar]
 Kim, E., Einsel, C., Lee, H. M., Spurzem, R., & Lee, M. G. 2002, MNRAS, 334, 310 [NASA ADS] [CrossRef] [Google Scholar]
 Lane, R. R., Kiss, L. L., Lewis, G. F., et al. 2010, MNRAS, 406, 2732 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D. 1960, MNRAS, 120, 204 [NASA ADS] [Google Scholar]
 Mackey, A. D., & Gilmore, G. F. 2003, MNRAS, 338, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Massey, P., & Hunter, D. A. 1998, ApJ, 493, 180 [NASA ADS] [CrossRef] [Google Scholar]
 Meylan, G., & Heggie, D. C. 1997, A&ARv, 8, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Meylan, G., & Mayor, M. 1986, A&A, 166, 122 [NASA ADS] [Google Scholar]
 Meza, A. 2002, A&A, 395, 25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ostriker, J. P., & Peebles, P. J. E. 1973, ApJ, 186, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Rosolowsky, E., Engargiola, G., Plambeck, R., & Blitz, L. 2003, ApJ, 599, 258 [NASA ADS] [CrossRef] [Google Scholar]
 Sabbi, E., Lennon, D. J., Gieles, M., et al. 2012, ApJ, 754, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Sana, H., de Koter, A., de Mink, S. E., et al. 2012, A&A, submitted [Google Scholar]
 Shapiro, S. L., & Marchant, A. B. 1976, ApJ, 210, 757 [NASA ADS] [CrossRef] [Google Scholar]
 Spitzer, L. 1987, Dynamical evolution of globular clusters (Princeton, NJ: Princeton University Press) [Google Scholar]
 Spitzer, Jr., L., & Hart, M. H. 1971, ApJ, 164, 399 [NASA ADS] [CrossRef] [Google Scholar]
 Stephens, A. W., Catelan, M., & Contreras, R. P. 2006, AJ, 131, 1426 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F., Le Poole, R. S., Reijns, R. A., Freeman, K. C., & de Zeeuw, P. T. 2000, A&A, 360, 472 [NASA ADS] [Google Scholar]
All Figures
Fig. 1 Illustration of the positions and RVs of the stars considered in this study. Symbol sizes denote the magnitude of the stellar velocities with respect to the average cluster velocity. The solid, dotted, and dashed lines correspond to the optimal rotation axis determined for models with a constant rotational velocity, a constant rotation rate, and a more realistic rotation curve (see Sect. 3), respectively. 

In the text 
Fig. 2 Top: RV in the system of the cluster as a function of distance X from the optimal rotation axis for models with constant rotation rate (left), constant rotational velocity (centre), and the more realistic rotation curve discussed above (right). The bestfit rotation curves are shown as solid lines, and ± σ_{1D} envelopes are represented by dashed lines. Bottom: histogram of amplitudes Ω (left), V_{rot} (centre) and A (right) from 10 000 Monte Carlo simulations of RV distributions with no rotation. The optimal values for the observed RV configuration are indicated by dashed vertical lines. Confidence levels of 94.6%, 93.6% and 95.6% are found respectively for the bestfit amplitude of each model. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.