abstract: The "new standard cosmology", based on the theory of inflation, has very impressive observational support. I review some outstanding problems of the new cosmology and the global view of the universe – the multiverse – that it suggests. I focus in particular on prospects for further observational tests of inflation and of the multiverse.
address:  Institute of Cosmology, Department of Physics and Astronomy,
  Tufts University, Medford, MA 02155, USA
author: Alexander Vilenkin
title: Perspectives in cosmology

# Introduction

We live at a very interesting juncture in cosmology, when a new paradigm appears to be shaping up. Its true significance will only be clear from the perspective of the future. In this talk I will try to outline where I think we stand and where we might be going.[^1] This is, of course, a tricky business. When speaking about perspectives, it's good to remember the quote from George Orwell, who wrote that even "to see what is in front of one's nose needs a constant struggle."

There are three basic questions in cosmology that we would like to answer. The first is: What is out there?, or What is the content of the universe? The second is: How did it get there?, or How did the universe evolve to its present state? And finally, What does it all mean? What is the big picture and how is it related to fundamental physics? There have been remarkable developments in all three of these areas over the last few decades.

The picture of the universe which has emerged from these developments can be called "the new standard cosmology." According to this picture, the universe has a nearly flat geometry, and its main three ingredients are dark energy (72%), dark matter (23%), and baryonic matter (less than 5%). The key role in the new cosmology is played by inflation – a period of super-fast, accelerated expansion in the early universe. Density perturbations that led to structure formation originated as quantum fluctuations during this inflationary period. Baryonic matter was formed in non-equilibrium, baryon-number-violating processes after inflation. After that, the evolution is along the lines of the hot big bang model. The new cosmology is still under construction, with details remaining largely uncertain. I will discuss some of these details in the rest of the talk.

# What is out there?

With some luck, we may be on the verge of discovering what dark matter is. Some of the dark matter candidates are well motivated from the particle physics perspective. For example, weakly interacting massive particles (WIMPs) could be produced as thermal relics in the early universe and are well suited for direct and indirect searches and for collider experiments. Axions could be produced by vacuum misalignment during inflation. There is no shortage of other candidates: supermassive particles, gravitinos, Q-balls, vortons, etc. I will stop at that, since dark matter will be discussed extensively at this meeting.

Turning now to dark energy, in Frank Wilczek's words, it is "the most mysterious fact in all of physical science." We have several dark energy candidates, but all of them have serious problems. First, the dark energy may simply be the vacuum energy, or the cosmological constant, $\rho_v$. The vacuum should then have pressure $P_v=-\rho_v$, with both $P_v$ and $\rho_v$ constant in time. A widely discussed alternative is a dynamical dark energy. Phenomenologically, it can be characterized by an equation of state, $$P=w\rho,$$ where $w\neq -1$ and may be time-dependent. In quintessence models, the dark energy is attributed to a scalar field with a slowly declining potential, e.g., $$V(\phi)\propto \phi^{-\alpha}.$$ As the field rolls down the potential, its energy density gradually gets smaller. The dark energy is so small today because the universe is so old, and the field $\phi$ had time to roll to very large values. Finally, there are modified gravity models, which explain the observed accelerated expansion by deviations from Einstein's gravity. Of all three alternatives, I think the front runner is the cosmological constant, for the reasons that I will now explain.

As I have mentioned, the cosmological constant has its problems. The first of these is the old cosmological constant problem, which has been with us for nearly half a century. All quantum fields, like the electromagnetic field, contribute to the vacuum energy. This contribution is given by the sum of zero-point oscillation energies over modes with all possible frequencies. This is a divergent sum, and we should cut it off at some mass scale $M$, which may be set, for example, by the supersymmetry breaking scale, $$\rho_v = {1\over{V}}\sum{1\over{2}} \omega \sim M^4.$$ (I use the system of units where $\hbar=c=1$.) With $M\sim 1$ TeV, this gives $\rho_v\sim 10^{60}\rho_m$, where $\rho_m$ is the average density of matter, while observations indicate that $\rho_v\sim 3\rho_m$. Fermions and bosons contribute to $\rho_v$ with opposite signs, but we would need cancellation up to the 60th decimal point in order to explain the observed value.

Another problem is related to the fact that the vacuum energy density remains constant in time, while the matter density is diluted by the expansion. Why then do we happen to live at the very special epoch when the two densities are comparable? This is the so-called coincidence problem.

As of now, the only viable resolution of these problems is provided by the anthropic approach. It assumes that $\rho_v$ takes a wide range of values in different parts of the universe. Of course, this can work only if the underlying particle physics admits many vacua with different values of $\rho_v$. This picture has some support from string theory, and I will get back to it toward the end of the talk. For now, let us assume that $\rho_v$ is indeed variable from one region to another.

The vacuum energy acts as a repulsive force, and it prevents galaxies from forming once it comes to dominate the universe. So, a large $\rho_v$ means no galaxies. We can calculate the probability distribution for an observer randomly picked in the universe to measure a given value of $\rho_v$. In regions with large values of $\rho_v$ there are no galaxies and therefore no observers, so the probability is strongly suppressed. Very small values, where $\rho_v$ is much smaller than it needs to be for galaxy formation, are also rather unlikely. One finds that the observed value of $\rho_v$ is well within the $2\sigma$ range of the resulting distribution. It is worth pointing out that the anthropic prediction of $\rho_v \sim
\rho_m$ was made well before the observation.[^2]

The main drawback of the alternative dynamical dark energy models is that they do not address the cosmological constant problems. For example, quintessence models assume that the dark energy density approaches zero as $\phi\to\infty$, while in general it should approach a nonzero vacuum energy density $\rho_v$. These models also offer no solution to the coincidence problem. The same criticisms apply to the modified gravity approach. Moreover, the equation of state parameter $w$ derived from the WMAP data combined with supernova observations is $$w=-0.97\pm 0.07,$$ which is very close to the cosmological constant value $w=-1$.

Apart from the three main components of the universe, what else may be out there? We know there are ultrahigh-energy cosmic rays. Until recently their origin appeared mysterious, but now it seems that they may have a rather mundane astrophysical origin.

There is much to be learned about cosmic magnetic fields. We have a rather sketchy information about the field distribution on the largest scales, and the origin of the magnetic fields remains a mystery.

We believe there should be gravitational waves, and we expect to detect them in the coming decade.

There are also some more exotic things which may or may not be there. For example, topological defects, such as cosmic strings.

Of course, there may also be some surprises – things we have not thought about which will change the direction of current research in new and unexpected ways.

I would like to discuss cosmic strings in a bit more detail. Strings are linear defects which could be formed at a phase transition in the early universe. They are predicted in a wide class of elementary particle models. Some superstring-inspired models suggest that fundamental strings may also have astronomical dimensions and play the role of cosmic strings.

Once they are formed, strings evolve in a scale-invariant manner. A typical horizon volume at any time has some long strings stretching across it and a large number of small closed loops. The average number of long strings per horizon is[^3] $$N\sim 10/p,$$ where $p$ is the probability for two intersecting strings to reconnect. Field-theory, solitonic strings always reconnect, so $p=1$, while for superstrings $p\ll 1$. This opens an exciting possibility of directly testing superstring theory. If we detect $N\gg 1$ long strings in the sky, we will know that we are looking at superstrings.

Strings can be detected through a variety of observational effects. Oscillating loops of string emit gravitational waves – both bursts and a stochastic background. They can also be sources of ultrahigh-energy cosmic rays. Long strings can act as gravitational lenses and can produce characteristic signatures in the CMB: discontinuous temperature jumps along the strings and a B-mode polarization pattern due to rotational perturbations induced by strings. If cosmic strings are discovered, this would open a window into physics of very high energies, which are not likely ever to be accessed in accelerators.

# Cosmic history

There is a great deal yet to be learned about the structure formation in the universe, from the dark ages, through star and galaxy formation, and the formation of supermassive black holes. We are now in the era of large surveys and massive computer simulations, and we can expect a lot of progress in this area in the coming decade. Here, I will focus on the theory of inflation and its upcoming tests.

A period of inflation in the early universe explains why the universe is large, expanding, homogeneous, isotropic, and flat. It is remarkable that no viable alternatives to inflation have been suggested, despite some heroic attempts. This is paralleled only by Darwin's theory in biology. And like Darwin's theory, inflation is not a specific model, but rather a paradigm encompassing a wide class of models.

Simple models of inflation make very definite predictions: (i) a flat universe, with the total density very close to the critical density, $$\Omega_{total}=1\pm 10^{-5};$$ (ii) a nearly flat spectrum of adiabatic, Gaussian density fluctuations, $$P_s(k)\propto k^{-n_s},$$ with the spectral index in the range $$0.95 < n_s < 0.98;$$ (iii) a stochastic background of gravitational waves with a flat spectrum. The amplitude of these gravitational waves is model-dependent, but in the simplest model with a quadratic inflaton potential, $$V(\phi)= {1\over{2}}m^2\phi^2,
\label{V}$$ the ratio of the power in tensor and scalar perturbations is $$r\equiv P_t/P_s \approx 0.14.$$ A gravitational wave background of this magnitude may be within the reach of Planck satellite. Current observations yield the values $$\Omega_{total}=1\pm 0.01, ~~~~~~ n_s=0.96\pm 0.01, ~~~~~~ r<0.22,$$ and are consistent with the simplest model.

If the simplest quadratic model () is supported by subsequent tests, this would be a spectacular triumph of inflationary cosmology. On the other hand, the CMB data show some tantalizing hints of non-Gaussianity. The parameter $f_{NL}$, which is related to the three-point correlation function and should vanish for Gaussian fluctuations, is $$f_{NL}=30\pm 15.$$ This is consistent with the Gaussian value $f_{NL}=0$ at the $2\sigma$ level, but barely. If indeed $f_{NL}\sim 30$, then Planck will detect it at $6\sigma$ level. That would be a very dramatic development, as it would rule out a wide class of models of inflation, including most single-field models.

Moreover, a number of anomalies have been reported in the CMB data: (i) the CMB quadrupole is much smaller than expected; (ii) a strange alignment between the low multipoles; (iii) an asymmetry between the correlation functions in the Northern and Southern hemispheres; (iv) a statistically unlikely cold spot in the WMAP microwave map of the sky. Are these anomalies due to chance? If you look hard enough at any stochastic pattern, you will start seeing things. For example, you can see my initials "AV" in the left upper part of the map.[^4] But it is hard to believe that all anomalies are due to chance, so some new physics may well be forthcoming from all this.

# The big picture

Inflationary cosmology has profound implications for the global structure of the universe, particularly for what the universe is like beyond the horizon, and for the beginning and the end of the universe.

The key point is that inflation is generically eternal. Even though it has ended in our local neighborhood, it still continues in remote parts of the universe. The reason is that the end of inflation involves quantum, probabilistic processes and does not occur everywhere at once. The details are somewhat model-dependent. For example, in models where inflation ends through bubble nucleation, we have an inflating high-energy false vacuum. Low-energy bubbles nucleate and expand in this inflating background. We live in one of the bubbles and can see only a small part of it. Even though the bubbles expand at speeds close to the speed of light, they rarely collide, because the false vacuum that separates them is expanding even faster. This makes room for more bubbles to form, and the process continues forever.

The physical properties of all low-energy bubbles do not have to be the same. In fact, string theory, which is at present our best candidate for the fundamental theory, appears to have a multitude of solutions describing vacua with different values of the low-energy constants of nature. These solutions are characterized by different compactifications of extra dimensions, by branes wrapped around extra dimensions is different ways, by different values of the fluxes, etc. The number of possibilities is combinatorial and can be as high as $10^{1000}$. Bubbles of all possible vacua will nucleate in the course of eternal inflation. The resulting multiverse, with bubbles within bubbles within bubbles, provides a natural arena for anthropic applications, such as the one I described for the cosmological constant.

The multiverse picture makes many people rather unhappy. There have been complaints that it cannot be tested, even in principle. If only a small part of one bubble is accessible to observation, how can we ever verify that other bubbles really exist?

I think such a test may in fact be possible. We can look for observational signatures of bubble collisions. Even though collisions are rare, our bubble will collide with an infinite number of other bubbles in the course of its expansion. A bubble collision could produce an imprint in the CMB; its detection would provide direct evidence for eternal inflation. There is, however, no guarantee that a bubble collision has occurred within our cosmic horizon.

Another approach is to use theoretical models of the multiverse to make statistical predictions for what we can expect to see in our observable part of the universe. One such prediction, for the cosmological constant, has already been confirmed. This may be our first evidence that there is indeed a huge multiverse out there. We can also derive predictions for other parameters, such as the density parameter $\Omega$, the ratio of baryonic and dark matter densities, the neutrino masses, etc. With more successful predictions, we may be able to prove the case for the multiverse beyond reasonable doubt.[^5]

If inflation has no end, could it also have no beginning? This would allow us to avoid many perplexing questions associated with the beginning of the universe. However, a theorem that I proved with Arvind Borde and Alan Guth states that, under very general conditions, inflationary spacetimes are incomplete to the past. This indicates that inflation must have had some sort of a beginning.

What could have happened before inflation? In my (biased) view, the most attractive possibility is the spontaneous nucleation of a small closed universe out of nothing. The nucleation process can be elegantly described in terms of a Euclidean instanton solution of Einstein's equations. Going beyond this semiclassical description will require further progress in quantum gravity.

I would like to conclude with some important news about the end of the universe. It is often said that if the dark energy is a cosmological constant, then the universe will expand forever. This is true for the universe as a whole, but not necessarily for our local region. In the multiverse picture, there must be a large number of negative-energy vacua, and bubbles of such vacua will inevitably form within our (nearly zero-energy) vacuum. At some point, probably in a very distant future, our neighborhood will be engulfed by a negative-energy bubble. The expansion will then locally turn into contraction, and our region will collapse to a big crunch. The good news is that the bubble will arrive at nearly the speed of light, so we will have no time to worry about it.

[^1]: Limitations of space do not allow me to provide adequate references to the literature, so I chose to give no references at all.

[^2]: I would like to emphasize that this prediction is not based on the trivial statement, often called "the anthropic principle", that observers can live only in regions where the constants of nature are consistent with their existence. Rather, it is based on "the principle of mediocrity", which asserts that we are typical representatives in a class of similar observers in the universe. Although not always explicitly stated, this principle is always assumed in statistical predictions, such as the prediction of the multipoles of the CMB.

[^3]: Both the normalization and the $p$-dependence in this equation are somewhat uncertain. The dynamics of a cosmic string network involves highly non-linear processes on a wide range of length scales; as a result, a quantitative description of string evolution is an extremely difficult problem. This problem is now successfully being addressed in new, high-resolution numerical simulations of cosmic strings.

[^4]: I am grateful to Andrei Linde for pointing this out.

[^5]: A long-standing problem with making predictions in the multiverse stems from the fact that any event having a nonzero probability will happen in the course of eternal inflation, and it will happen an infinite number of times. Statistical predictions are based on relative frequencies of events in the limit of $t\to\infty$. One finds however that the outcome sensitively depends on the limiting procedure. In the last few years there has been a significant progress towards resolution of this problem (which is known as "the measure problem").


author: P.L. Dufton, N. Langer, P.R. Dunstall, C.J. Evans, I. Brott, S.E. de Mink[^1], I.D. Howarth, M. Kennedy, C. McEvoy, A.T. Potter, O.H. Ramírez-Agudelo, H. Sana, S. Simón-Díaz, W. Taylor, J.S. Vink
bibliography: ./literature.bib
date: Received; accepted
title: The VLT-FLAMES Tarantula Survey X: Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Projected rotational velocities ($v_{\rm{e}} \sin i$) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Projected rotational velocities range up to approximately 450 km s$^{-1}$and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$$v_{\rm{e}}$$\leq$`<!-- -->`{=html}100 km s$^{-1}$and the high velocity component having $v_{\rm{e}}$$\sim 250$ km s$^{-1}$. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations.

# Introduction

Rotation is thought to have a major influence on the evolution of massive stars. It may affect both how the star evolves during hydrogen core burning and the nature of any supernova explosion . It can also provide insights into other phenomena such as star formation and how magnetic fields affect stellar evolution .

There have been several large studies of rotation in massive main sequence stars both in our own Galaxy and in the Magellanic Clouds. For the former, one of the earliest comprehensive samples was provided by . As this sample was selected without any previous knowledge about the stellar spectra, it provides an unbiased sample. Estimates of the projected rotational velocities of 123 early-type stars were presented together with upper limits for 63 stars with relatively narrow absorption lines. investigated 1100 B-type stars selected from the bright star catalogue and this is one of the largest samples for which the biases are well understood. For late B-type stars they find evidence for a bi-modal distribution but early-B type stars do not seem to show this bi-modality, although the number of such stars in their sample is small. combined data from their observing campaigns to provide projected rotational velocities for approximately 660 cluster stars and 470 field stars. They find that the average projected rotational velocity is significantly lower for field stars, which they attribute to the fact that they are systematically more evolved and spun down than their cluster counterparts.

For the Magellanic Clouds, there have been two investigations using ESO large programmes, viz. and . These observed predominantly field stars, while as discussed above for the Galactic studies, the behaviour of cluster stars may be different. Also given the limited time cadence of the observations, it was not always possible to reliably distinguish between single and multiple systems.

Here we present projected rotational velocity estimates for a sample of 334 early B- and late O-type stars located in the 30 Doradus region . The rotational velocities of the more massive O-type stars will be discussed elsewhere (Ramirez-Agudelo et al. in preparation). Our sample has been selected as having no significant radial velocity variations, thereby reducing the contamination from multiple systems. We have also inferred the intrinsic rotational velocity distribution and discuss its bi-modal structure.

# Observations

The spectroscopic data were obtained as part of the European Southern Observatory (ESO) large programme using the Fibre Large Array Multi-Element Spectrograph, FLAMES on the Very Large Telescope (VLT), primarily with the Giraffe spectrograph, but also supplemented with data from the Ultraviolet and Visual Echelle Spectrograph, UVES . All observations have been based in and around the LMC open association 30 Doradus, with the core region R136 having a 15 arcsecond exclusion radius as it is too densely populated for use of the Medusa fibres. This region has been observed with five pointings using the ARGUS integral field unit (IFU). The current work will use the observations obtained using the two bluer region Giraffe spectra to obtain estimates of the projected rotational velocity and the red spectra to search for Be-type emission. The observational configurations are summarized in Table , with details of the target selection, observations and data reduction being available in Paper I.

The sample is that defined in Paper I as potentially B-type stars and comprises approximately 540 targets. Ongoing spectral classification indicates that some targets are late-O type stars with spectral types of O9.5 or O9.7. It is believed that our sample is effectively complete as regards the O9.7 and later spectral types but is incomplete for the O9.5 spectral type as approximately half these targets were included only in the dataset of potentially O-type stars. We return to this in Sect. . Additionally the magnitude cutoff employed in the survey implies that subject to spatial variations in extinction our sample should not include (near) zero-age main sequence stars with spectral types later than approximately B3.

have undertaken a radial velocity analysis for this dataset to identify binaries. They find 37% of their sample show unambiguous velocity variations with an additional 6% having an uncertain binary status. We excluded the former from our analysis but included the latter, although we consider these results separately. Additionally thirty stars without significant radial velocity variations have been identified as supergiants and a preliminary model atmosphere analysis, similar to that employed by implied that they all had logarithmic gravities between 2.5 and 3.2 dex. As these objects will normally have evolved from O-type stars, they will be discussed elsewhere (McEvoy et al, in preparation; Simón-Díaz et al., in preparation).

Hence our sample of 334 stars (of which 45 are of uncertain status) is a subset of the Tarantula B-type sample that lie close to the main sequence. They show no convincing evidence of any radial velocity variation typically at the 10 km s$^{-1}$ level over six epochs spread over one year. Assuming a spectral type range of B3 to O9.5, leads to an effective temperature range of approximately 20 000 to 34 000 K and a mass range of approximately 6 to 16  $\rm{M}_\odot$ for an LMC metallicity.

have discussed completeness and estimate that up to 50% of the apparent B-type main sequence sample could be undetected binaries. Indeed the failure to identify long period systems (typically $P>1000$ days) is an intrinsic limitation of the VLT-FLAMES Tarantula Survey (VFTS) given the one year baseline of the observations. Assuming that this analysis is correct, our sample could contain a significant fraction of binaries, albeit preferentially with long orbital period and/or low-mass companions. The light and spectrum from such binaries would be dominated by the main (B-type) component. Additionally the long orbital periods involved are unlikely to have led to any significant tidal interaction that would have modified the rotation of the objects.

As there were no significant radial velocity variations in our sample, spectra could be combined without including any wavelength shifts. Several methods were adopted with the earlier analysis being based on co-added spectra with exposures containing cosmic rays or having a low signal-to-noise (SNR) being excluded. Later, all the exposures were normalised using low order polynomials. Normally an extended wavelength region could be normalised with a single polynomial but for some exposures, there was significant variations in the continuum flux with wavelength (se Paper I for details). In these cases individual lines or groups of lines were normalised. Spectra were then combined with a simple median filter or a weighted mean with $\sigma$-clipping; in both cases low SNR exposures were again excluded. Given the relatively small range in SNR for a given star and wavelength region, these two approaches gave essentially identical spectra. Comparison of these data and of the earlier co-added spectra for five stars showed excellent agreement (at the 1-2% level) in the estimated mean projected rotational velocities indicating that the choice of reduction technique should not be a significant source of uncertainty.

|               |                 |       |                |
|:--------------|----------------:|------:|---------------:|
| Instrument    | $\lambda$-range |     R |      Exposures |
|               |             (Å) |       |      (seconds) |
| Giraffe LR02  |       3960-4560 |  7000 | 6 x (2 x 1815) |
| Giraffe LR03  |       4505-5050 |  8500 | 3 x (2 x 1815) |
| Giraffe HR15N |       6445-6815 | 16000 | 2 x (2 x 2265) |

Summary of the spectroscopy used in the rotational velocity analysis. For some stars additional LR02 observations were obtained when the initial observations were found to be of poor quality due to weather conditions.

# Projected rotational velocities

## Methodology

Several methods are available for estimating projected rotational velocities ($v_{\rm{e}} \sin i$). Those most commonly used for OB-type stars include a direct measure of the full-width-half-maximum of the spectral lines , profile fitting or cross correlation . have discussed the utility of the Fourier transform (FT) method for early type stars, as this should be able to separate the rotational broadening from other broadening mechanisms such as macroturbulence. Recently this approach has been widely used to study the different mechanisms contributing to the broadening of spectral lines in early-type supergiants .

We have adopted this Fourier transform method for all the stars in our sample. This approach relies on identifying the first minimum in the Fourier transform for a spectral line, which is assumed to be the first zero in the Fourier transform of the rotational broadening profile with the other broadening mechanisms exhibiting either no minima or only minima at higher frequencies. Fig. shows the line at 4026Å and its Fourier Transform for two targets, VFTS046 and VFTS699. Their spectra are of similar strength (implying a spectral type near B0) and hence the difference in their line profile reflects their different projected rotational velocities. Also shown are the rotational broadening profiles for a projected rotational velocity, corresponding to the first minimum in the Fourier Transform. As would be expected these agree well with the observations in the line centre but are narrower in the wings for VFTS046 due to the intrinsic Stark broadening of this diffuse line.

As discussed by , the first minimum in the Fourier Transform can be difficult to identify in spectra with, for example, low SNRs or significant nebular contamination. Therefore we have also profile fitted (PF) all the spectral features, assuming that rotation dominates the line broadening i.e. we have used a simple rotational profile and assumed that the intrinsic and instrumental profiles can be approximated by scaled $\delta$-functions. The neglect of other broadening mechanisms will lead to these estimates being upper limits but they provide a useful constraint on the actual projected rotational velocities. We discuss these results further in .

Two different sets of spectral lines were adopted for analysis depending on the degree of rotational broadening. For the narrower lined stars (with $v_{\rm{e}} \sin i$ $\la 150$km s$^{-1}$), a selection of metal lines and non-diffuse helium lines (Set 1 in Table ) was used. These had the advantage of being intrinsically narrow, whilst the former were not normally affected by nebula emission. For large rotational velocities, reliable results could only be obtained from stronger absorption lines and in these cases a selection of diffuse and non-diffuse helium lines (Set 2 in Table ) was considered. Some of the features are in fact close doublets or triplets, whilst additionally the helium diffuse lines have weak forbidden components (we note that blending of the 4026Å feature with a line should not be significant for our sample). For the former, the separations were less than the broadening due to rotation and indeed the instrumental profile, whilst the latter were only used for the more rapidly rotating stars where the profiles appear symmetric within the observational uncertainties.

Normally it was not possible to obtain estimates from all the spectral lines in the adopted set for three main reasons. Firstly, some lines were too weak to provide reliable measurements. For example, the and features could not be used for the hotter stars, whilst the and features were weak in the cooler objects. Secondly, blending could affect some lines and in particular the line at 4121Å (blended with ) and the line at 4481Å (blended with ). However the most serious constraint was the presence of strong nebular emission in the spectrum of some stars, which could not be reliably removed for our multi-fibre observations. This emission could affect the feature but more seriously the spectrum. The following approach was taken for increasing amounts of contamination.

1.  The emission was snipped out of the spectra (see Fig. ), which was then regridded and analysed. The quality of the Fourier transform, the consistency between the estimates from the FT and PF approaches and the quality of the fit of the rotationally broadened profile to the observations was then used to assess the reliability of the estimate.

2.  The lines with the greatest nebular contamination (normally the diffuse triplet transitions) were excluded. As can be seen from Table 4, there are no systematic differences between the estimates obtained from different lines. Hence the use of a subset of lines is unlikely to introduce significant biases.

3.  When a reliable estimate could not be obtained, an attempt was made to assign the star to a velocity range covering 40 km s$^{-1}$ (0-40, 40-80, 80-120 km s$^{-1}$ etc.). The observed profiles were then compared with rotationally broadened profiles appropriate to this range to validate this decision. Hence although no reliable projected rotational velocity could be estimated, it would be possible to include such stars in any statistical study.

There were a number of narrow lined stars whose line width appeared to be dominated by the instrumental profile which corresponded to a FWHM of approximately 40 km s$^{-1}$. As would be expected, their Fourier Transforms did not show convincing minima with individual measurements showing a significant degree of scatter. Other stars showed a greater degree of broadening and yielded mean projected rotational velocities in the range of 30 to 40 km s$^{-1}$. Here we have taken a conservative approach and have estimated rotational velocities only for stars with individual measurements from the Fourier Transform methodology that were consistently greater than 40 km s$^{-1}$ and where the profile fitting implied a significant rotational broadening component.

For each star, a projected rotational velocity was estimated only if there was at least three independent measurements. For the narrow lined stars, the average number of measurements was 6.9, with approximately 90% having five or more measurements, with the corresponding values for the broader lined stars being 6.3 and 93%. For the latter, an additional 16 stars were assigned projected rotational velocity ranges whilst one star, VFTS729, appeared to have asymmetric line profiles, suggesting a binary or composite nature, and was excluded.

Tables 3 and 4 (only available online) lists all the individual measurements using both the FT and the PF methodologies, together with the percentage depth of the feature as estimated from the profile fitting. The former is for the metal line set and the latter for the helium line set. In both tables stars that are of an uncertain status i.e. could possibly be binaries (see Sect. ) are identified with an asterisk. The tables also give means and standard deviations with the FT estimates being weighted by the central depths. Using unweighted means would have lead to no systematic shift in the estimates with typical changes of 2% for line set 1, and a 1% decrease for line set 2 with typical changes again being 2%. Hence the choice of weighting is unlikely to be a significant source of error. We also list the mean values from the profile fitting; in this case the individual estimates were not weighted as these are better considered as providing upper bounds. Stars with narrow lined spectra that were instrumentally dominated have been assigned a projected rotational velocity of $\la 40$ km s$^{-1}$.

and have investigated the effects of equatorial gravity darkening due to angular rotation ($\Omega$) on observed line profiles. They find that for angular velocities near the critical limit ($\Omega_c$) the effects can be significant and lead to a systematic underestimation of the real projected rotational velocity. For example, considered a representative metal and helium diffuse line in stars with a spectral range B0-B3 (similar to that observed here) and for $\Omega$/$\Omega_c$=0.95 and $\sin i = 1$. They found the projected rotational velocity would be underestimated by between 9 and 22%. Unfortunately we do not know either the $\sin i$ values or atmospheric parameters of our targets and are hence unable to systematically correct for this effect.

However we can use the results of to obtain an insight into the magnitude of any systematic effects. They used the methodology of to correct observed projected rotational velocities for 47 Be-type stars towards the LMC cluster, NGC 2004. For $\Omega$/$\Omega_c=0.85$ (the lowest value that they considered), the correction were 5.5% for stars with $v_{\rm{e}} \sin i$$<200$km s$^{-1}$, 3.6% ($200\leq$$v_{\rm{e}} \sin i$$<300$km s$^{-1}$), 3.3% ($300\leq$$v_{\rm{e}} \sin i$$<400$km s$^{-1}$) and 3.0% ($400\leq$$v_{\rm{e}} \sin i$$<500$km s$^{-1}$). The systematic errors in our estimates should be considerably smaller as it is unlikely that the majority of our targets are rotating with such high values of $\Omega$/$\Omega_c$. For example, taking the median rotational velocity estimated from the de-convolution discussed in Sect. and the critical velocities from the grid of for early B-type stars leads to a value of $\Omega$/$\Omega_c\simeq 0.3-0.4$. Hence we do not expect that there will be a significant effect on our sample as a whole, although a small number of individual estimates (where the target is rotating with near critical velocity) could be affected.

|  Set 1  |               |  Set 2  |               |     |     |     |     |
|:-------:|:-------------:|:-------:|:-------------:|:---:|:---:|:---:|:---:|
| Species | $\lambda$ (Å) | Species | $\lambda$ (Å) |     |     |     |     |
|         |     4120      |         |     4009      |     |     |     |     |
|         |     4169      |         |     4026      |     |     |     |     |
|         |     4267      |         |     4120      |     |     |     |     |
|         |     4437      |         |     4143      |     |     |     |     |
|         |     4481      |         |     4387      |     |     |     |     |
|         |     4552      |         |     4471      |     |     |     |     |
|         |     4567      |         |     4713      |     |     |     |     |
|         |     4574      |         |     4921      |     |     |     |     |
|         |     4661      |   \-    |      \-       |     |     |     |     |
|         |     4713      |   \-    |      \-       |     |     |     |     |

Sets of absorption lines used for projected rotational velocity estimates

## Results and estimated uncertainties

Before discussing the significance of the results presented in Tables 3 and 4, it is important to validate their reliability. Several tests were undertaken as discussed below. Firstly one star (VFTS331) was inadvertently analysed twice with excellent agreement; the mean of the FT estimates agreed to within 2 km s$^{-1}$ and that of the PF estimates to within 3 km s$^{-1}$. Secondly four targets in the overlap region between the narrow and broad lines stars ($100\la$$v_{\rm{e}} \sin i$$\la 200$ km s$^{-1}$) were measured using both sets of lines and in Table the mean values and the estimated errors are summarized. The estimates are in reasonable agreement with no evidence of systematic differences. The t-statistic and degrees of freedom for unequal sample sizes and unequal variance (Welch's t-test) are also listed. None of these are significant at the 10% levels, but those for VFTS234 and VFTS452 are significant at at the 20% level.

| Star | Line set 1 | Line set 2 | $|t|$ | df |  |
|:---|:--:|:--:|:--:|:--:|:--:|
| 046 | 169$\pm$`<!-- -->`{=html}15 (7) | 171$\pm$`<!-- -->`{=html}15 (7) | 0.2 | 6 |  |
| 060 | 129$\pm$`<!-- -->`{=html}16 (9) | 129$\pm$`<!-- -->`{=html}8 (7) | 0.0 | 6 |  |
| 234 | 157$\pm$`<!-- -->`{=html}12 (7) | 147$\pm$`<!-- -->`{=html}6 (7) | 1.9 | 6 |  |
| 452 | 125$\pm$`<!-- -->`{=html}15 (7) | 134$\pm$`<!-- -->`{=html}5 (8) | 1.5 | 7 |  |

Comparison of mean projected rotational velocities (km s$^{-1}$) deduced using the Fourier Transform methodology from the two line sets. Also listed are the standard deviations and number (in brackets) of the estimates. The modulus of the t-statistic and number of degrees of freedom are also listed.

A comparison has also been made of the mean projected rotational velocities obtained from the FT and PF methodologies. These are illustrated in Fig. for the the narrow lined (line set 1) and broader lined (line set 2) subsamples. In general the results are in excellent agreement for both line sets. At lower projected rotational velocities (e.g. $\la 80$ km s$^{-1}$), the profile fitting methodology yields systematically higher estimates, which would be expected given that only rotational broadening was considered.

To test this hypothesis, we took the spectrum from a TLUSTY model , with an effective temperature of 25 000 K, a logarithmic gravity of 4.0 dex and an LMC composition . This was convolved with a Gaussian profile to represent instrumental broadening and with different amounts of projected rotational broadening. The profiles were then analysed using the PF method to estimate the projected rotational velocity. Three representative absorption features were considered, viz. the line at 4552Å, the non-diffuse helium line at 4713Å and the diffuse helium line at 4387Å. These results are also shown in Fig. and reproduce the trend in the differences between the observational estimates discussed above. Indeed although they only represent a subset of lines for a single model, they provide strong evidence for the reliability of our Fourier Transform results.

The sample standard deviations are typically 10% of the mean projected rotational velocities for both line sets. However the error in the mean values will be smaller and if the estimates are normally distributed, this would be the sample standard deviation divided by the square root of the number of observations. This would then lead to a typical uncertainty in the mean values of 4%.

An independent error estimate can be found at higher projected rotational velocities, as here the profile fitting methodology should provide reliable estimates. We have therefore compared the two methodologies in these regimes. For line set 1, we considered all stars that had estimated projected rotational velocities of greater than 80 km s$^{-1}$. For each star the ratio of the PF to the FT estimate was found, with the mean of these ratios being $1.01\pm 0.06$. For line set 2 and limiting our sample to those stars with estimates greater than 300 km s$^{-1}$ leads to a mean ratio of $1.02\pm 0.04$. Hence in these high projected rotational velocity regimes, the methodologies give excellent agreement, confirming the trends seen in Fig. . Additionally as the standard deviations for the ratios will include a contribution from uncertainties in the PF approach, they are consistent with our estimated error of 4% for the means of the FT estimates.

Rotational velocities have also been estimated for the apparently single O-type stars in the Tarantula Survey (Ramirez-Agudelo et al., in preparation). They also used the FT technique but a different sets of lines reflecting the different temperature range. Due to uncertainties in the preliminary spectral classification, late O- and early B-type targets were included in both datasets (see Sect. ). Excluding targets that have only been assigned here to a velocity bin, there are 42 stars with independent estimates. There is no significant systematic difference between the two sets of measurements, with the mean value of the absolute differences being 3$\pm$`<!-- -->`{=html}15 km s$^{-1}$and of the relative differences being 0$\pm$`<!-- -->`{=html}8%. For thirty seven of the targets, the difference is less than 10%, with the remaining five ranging from 10-30%. We discuss these discrepancy cases below:

**VFTS 226, 627:** These targets have narrow lined spectra and are hence close to the resolution limit for getting reliable projected rotational velocities. Additionally macroturbulence may be significant given their relatively large effective temperatures. Hence we do not consider the relatively small absolute discrepancies (7-11km s$^{-1}$) in the estimates to be significant.

**VFTS 373, 543:** Strong emission is present in both these spectra, making reliable estimates difficult. Indeed our sample standard deviations are approximately 30 km s$^{-1}$(see Table 4), implying that the estimates are consistent within the uncertainties.

**VFTS 141:** The spectrum of this star is well observed, with the nebular emission being weak. Although the two sets of measurements differ by 19%, when only lines that are common to both analyses are included this decreases to 7%. Hence we believe that the measurements are consistent.

## Distribution of projected radial velocities

A histogram of our projected rotational velocities, excluding the targets with an uncertain status, is shown in Fig. . A binsize of 40 km s$^{-1}$ has been adopted as this is consistent with both our spectrally unresolved targets and those where only a velocity range was assigned. The distribution appears to be double peaked with approximately 25% of the sample, having $v_{\rm{e}} \sin i$$\leq 80$km s$^{-1}$. The sample should be relatively unbiased as it was selected on the basis of magnitude (with no colour cut-off) and then on the ability to position the spectrograph fibres. The former could bias the sample towards targets with low extinction, whilst the latter could bias the sample away from regions with high star densities but neither of these would appear to be sufficient to explain the double peak. Another explanation of the bi-modal structure could be that many of the narrow-lined 'single' stars are in fact part of binary systems with synchronised orbits. We would expect such systems to have relatively small orbital periods with large radial velocity variations, whilst their radial velocity estimates will have a high accuracy. As such the radial velocity analysis of should have identified such targets and hence this bias should not exist.

As discussed in Sect. , we believe that our sample contains effectively all the targets with spectral types O9.7 or later but only partially samples those with spectral type O9.5. Hence in Fig. , we also show a histogram with the O9.5 stars excluded. We note that the two histograms are very similar with any differences being less than the random Poisson noise expected from finite sample sizes. Hence in the subsequent discussion of these data, we have considered the complete sample but excluded stars with an uncertain binary status. Our principle conclusions would remain unchanged if we had excluded the targets with a O9.5 spectral type and/or included those of uncertain status.

## Be-type stars

The classical Be-type star is a B-type star that shows or has previously shown prominent emission features in its Balmer line spectrum, indicating the presence of a geometrically flattened, circumstellar disc . Although the mechanism for the Be phenomena remains unclear, such stars are believed to rotate with velocities greater than 40% up to approximately 90% of the critical limit, where the centrifugal force balances that of gravity . Reliable identification of Be-type stars in this sample was difficult because of the strong nebular emission. Additionally sky subtraction did not work well for these fibre-fed spectra, because of the spatially highly variable nature of this emission. Hence to search for evidence of Be-type emission, we examined red-region (HR15N) spectra processed without sky subtraction. Although there was no prospect of identifying weak or narrow stellar emission, broad Be-type H$\alpha$ emission wings were identifiable in 52 out of the 289 constant-velocity stars , and 8 out of 45 stars of uncertain status – i.e., about 18% of the sample. This fraction is comparable with canonical figures for LMC *clusters* (e.g., ), even though we are only able to identify relatively strong emission. The Be-type candidates are identified in Tables 3 and 4 by crosses (+).

In the bottom panel of Fig. , the projected rotational velocity distribution is shown for the sample excluding the Be-type candidates. As might be expected, this leads to a decrease in the number of rapidly rotating stars, whilst the bimodality is still present. There remains a small number of rapidly rotating stars ($v_{\rm{e}} \sin i$$\ge$ 400 km s$^{-1}$), although we cannot exclude them having weak Be-type emission.

The median value of $v_{\rm{e}} \sin i$ for the Be-star sample (excluding VFTS 102) is 254 km s$^{-1}$ (238 km s$^{-1}$ excluding the stars of uncertain status), and the mean $249\pm11$ km s$^{-1}$ ($244\pm12$ km s$^{-1}$). Although any bias introduced by the undetectability of weak stellar emission would presumably have skewed the results in favour of faster rotators, these values are somewhat smaller than the mean of $\sim$`<!-- -->`{=html}270 km s$^{-1}$ for LMC field and cluster Be stars reported by .

## Comparison with previous studies

Previous observational investigations of rotation in B-type stars have been discussed in the introduction. Here we limit the discussion to surveys in the Large Magellanic Cloud that are directly comparable with the current results. The FLAMES-I survey observed targets towards two LMC clusters, N11 and NGC2004 and provided estimates of projected rotational velocities . These estimates were not corrected for the systematic effect at near critical rotation discussed in Sect. . The rotational velocity distribution inferred by had a modal value of 100km s$^{-1}$corresponding to $\Omega$/$\Omega_c\simeq 0.2$ and hence we would expect that most of their sample will not contain significant systematic effects. In Figs. and , the projected rotational velocity distribution and cumulative distribution function are shown for the 92 targets that had no evidence for binarity, were not supergiants and had a spectral type of O9.5-B3 (consistent with our current sample).

provided estimates for 106 B-type and 47 Be-type stars towards NGC2004. Supergiants and spectroscopic binaries were also excluded in this study, and most investigated stars had masses in the range 4 $\rm{M}_\odot$to 20 $\rm{M}_\odot$. Again limiting the sample to the 146 O9.5-B3 leads to the distributions shown in Fig. and .

Neither of these two studies identified a bi-modal distribution. For example, fitted a similar histogram (but with a different bin size and sample; their Fig. 9a) by a uni-modal distribution. However tailoring their sample to match our range of spectral types leads to evidence for a bi-modal distribution, although this must be considered as marginal given their small sample size. By contrast, the larger sample of shows no clear evidence for bi-modality, although this cannot be excluded.

To investigate whether these samples originated from the same parent population as our current sample, we have used the two sample Kolmogorov-Smirnov and Kuiper tests. These implied that for the sample of , this was unlikely (P $\sim 1\times 10^{-3}$ and $6\times 10^{-3}$ respectively) but for that of , the probability was relatively high (P $\sim 0.4$ and 0.2 ). This is consistent with cumulative distribution functions plotted in Fig. .

These comparisons must be treated with caution as the samples remain poorly matched. Firstly the selection criteria adopted for the three surveys will have introduced different observational biases. For example both the samples of and contained evolved main sequence stars (see their Figs. 7 and 5, respectively). As the 30 Doradus region is thought to be young, our sample should contain a higher fraction of relatively unevolved B-type stars, although we do not currently have spectral types or luminosity estimates to confirm this. Secondly the previous studies observed a mixture of cluster and field stars, with the latter dominating. For example, estimated that $\sim$`<!-- -->`{=html}75% of their sample were field stars; given the complex nature of the 30 Dor region, the fraction of field stars (in the traditional sense) in the current sample will be relatively low. Thirdly the sample of was obtained at two epochs and they only identified approximately 10% of their sample as binaries. By contrast for the sample of , had identified 23% of the NGC2004 targets and 36% of the N11 targets as binaries and these have been excluded from the sample considered here. Hence given that the Tarantula survey has a lower limit of 37% on the binary fraction , we conclude that the sample of probably contains a significant fraction of binaries.

# Rotational velocity distribution

The relatively large size of our sample has been exploited to investigate its intrinsic rotational velocity distribution. We first describe our methodology, and then outline our results and investigate their robustness.

## Deconvolution of the $v_{\rm{e}} \sin i$ distribution

If the stellar rotation axes are randomly distributed, it is possible to infer the probability density for the rotational velocity distribution (P($v_{\rm{e}}$)) from that for $v_{\rm{e}} \sin i$. Previously several authors have investigated this by assuming an analytical function for P($v_{\rm{e}}$) and then optimising the fit with observation after convolution to allow for random inclinations. For example, for the analysis of the FLAMES-I datasets, rectangular, Maxwellian and Gaussian distributions were adopted. Given the structure in the projected rotational velocity distribution discussed above, it is not clear what adopted function would be appropriate.

We have therefore attempted to deconvolve the projected rotational velocity distribution to directly infer that of the rotational velocities. We have adopted the iterative procedure of which if it converges will converge to the point of maximum likelihood. For the continuous function, representing the observed projected rotational velocities ($\tilde{\phi}$ in the notation of Lucy), we followed his recommended procedure. The initial estimate for the rotational velocity distribution ($\psi^{0}$ in the same notation) was taken to be a Gaussian with a positive abscissa; tests showed that this choice did not affect the final iterative solution, consistent with the discussion of .

The effectiveness of this procedure has been tested by taking simulated distributions for P($v_{\rm{e}}$) and then using a Monte-Carlo technique to generate P($v_{\rm{e}} \sin i$). These were then de-convoluted to obtain an estimate for P($v_{\rm{e}}$). Simulated distributions considered included Gaussian, triangular and step functions. In general the deconvolution retrieved the large scale structure, while the estimated *projected* rotational velocity distribution, ($\tilde{\phi}$ in Lucy's notation), agreed well with the simulated P($v_{\rm{e}} \sin i$). Additional small scale structure was present especially where the adopted rotational velocity distribution was discontinuous. Also structure varied from simulation to simulation and appeared to be connected with that due to the Poisson noise. The only significant concern to arise in these simulations was that on occasions a small excess of slowly rotating stars were found even though there was excellent agreement between the simulated projected rotational velocity distribution, P($v_{\rm{e}} \sin i$), and that inferred from the deconvolution procedure ($\tilde{\phi}$). In turn this implied that relatively small errors in the observed projected rotational velocity distribution could propagate through to that for the rotational velocities, P($v_{\rm{e}}$).

Our adopted procedure is similar to that used by to investigate the rotational velocity distributions for Galactic A-type stars. However we differ from in not preprocessing our observational projected rotational distributions using a Gaussian filter. Our approach will not remove spurious structure in our inferred rotational velocity distribution due to the Poisson noise in the observational distributions. However it will also illustrate the limitations of this approach and its sensitivity to the assumptions made.

The choices of the number of iterations and the bin size for the observed projected rotational velocity are important. For the former, the iteration initially converges rapidly and then tends to accentuate small scale structure that may not be real. in his two rotational velocity examples only considered two iterations. In Fig. , we show (for an observation bin size of 40 km s$^{-1}$), the originally assumed distribution of the rotational velocities and the solutions after two, four, six and eight iterations. The solution after two iterations shows the major features of the distribution with further iterations introducing additional structure at approximately at 70, 270 and 420km s$^{-1}$. Although the signal will not be band-limited, this structure would appear to contain frequencies higher than the Nyquist frequency (or equivalently in the rotational velocity domain smaller than half the bin size). Additionally the errors due to Poisson statistics map onto a typical uncertainty of approximately 5$\times 10^{-4}$ in P($v_{\rm{e}}$) and hence this additional structure is probably not significant.

In the case of the adopted observational bin size it is necessary to balance the resolution in the velocity domain against the (Poisson) noise in the star counts. We have taken a pragmatic approach by considering two different possible bin sizes of 20 and 40 km s$^{-1}$, with for the former the spectrally unresolved targets being equally distributed between the first two bins. The results are illustrated in Figs. and for different number of iterations. The large scale structure of the estimates of the rotational velocity attribution are similar but the smaller bin size leads to more small scale structure. Although this will to some extent reflect the higher sampling rate, we believe that mainly that it arises from random errors in the $v_{\rm{e}} \sin i$ statistics.

On the basis of the above tests, we have adopted the larger bin size together with 4 iterations in our analysis and discussion given below. We note that the major conclusions would not have been significantly affected by other choices.

## Results

In order to illustrate our results, Fig. shows the observed $v_{\rm{e}} \sin i$ histogram from Fig. , normalised to unity when integrated over velocity. Also shown are the estimates after four iterations of the probability density of the projected rotational velocity distribution ($\tilde{\phi^{4}}$ in Lucy's notation) and of the corresponding probability density for the rotational velocity, P($v_{\rm{e}}$) ($\psi^{4}$ in Lucy's notation). There is good agreement between the observed $v_{\rm{e}} \sin i$ distribution and that inferred ($\tilde{\phi^{4}}$) from the iterative process. Additionally the estimate of the rotational velocity distribution appears reasonable with the large scale structure of the $v_{\rm{e}} \sin i$ distributions being conserved but moved to larger velocities due to the effect of inclination.

In Table our best estimates of the probability distribution function are tabulated. For smaller rotational velocities, the estimates have been rounded to 5$\times 10^{-5}$ but for larger velocities a small interval was chosen. We emphasis that these values are subjective as they require a judgement about the validity of the smaller scale structure. Also listed in this table is the cumulative distribution function (cdf) corresponding to this best estimate and because this is truncated at $v_{\rm{e}}$$\leq 600$ km s$^{-1}$, it does not reach unity. We recommend that:

1.  either the range up to the critical velocity is populated

2.  or that the cdf and estimate of P($v_{\rm{e}}$) are renormalised

depending on the nature of the application.

| $v_{\rm{e}}$ | P($v_{\rm{e}}$)$\times 10^3$ |  cdf  |     |     |     |     |
|-------------:|:----------------------------:|:-----:|:---:|:---:|:----|:----|
|            0 |             2.20             |   0   |     |     |     |     |
|           20 |             2.40             | 0.046 |     |     |     |     |
|           40 |             2.45             | 0.094 |     |     |     |     |
|           60 |             2.55             | 0.144 |     |     |     |     |
|           80 |             2.25             | 0.192 |     |     |     |     |
|          100 |             1.40             | 0.239 |     |     |     |     |
|          120 |             1.00             | 0.253 |     |     |     |     |
|          140 |             0.75             | 0.270 |     |     |     |     |
|          160 |             1.05             | 0.288 |     |     |     |     |
|          180 |             2.40             | 0.322 |     |     |     |     |
|          200 |             3.00             | 0.377 |     |     |     |     |
|          220 |             2.85             | 0.435 |     |     |     |     |
|          240 |             2.75             | 0.492 |     |     |     |     |
|          260 |             2.90             | 0.548 |     |     |     |     |
|          280 |             3.20             | 0.609 |     |     |     |     |
|          300 |             3.30             | 0.674 |     |     |     |     |
|          320 |             3.20             | 0.739 |     |     |     |     |
|          340 |             2.50             | 0.796 |     |     |     |     |
|          360 |             2.00             | 0.841 |     |     |     |     |
|          380 |             1.80             | 0.879 |     |     |     |     |
|          400 |             1.50             | 0.912 |     |     |     |     |
|          420 |             1.10             | 0.938 |     |     |     |     |
|          440 |             0.75             | 0.956 |     |     |     |     |
|          460 |             0.70             | 0.971 |     |     |     |     |
|          480 |             0.50             | 0.983 |     |     |     |     |
|          500 |             0.21             | 0.990 |     |     |     |     |
|          520 |             0.09             | 0.993 |     |     |     |     |
|          540 |             0.08             | 0.995 |     |     |     |     |
|          560 |             0.06             | 0.996 |     |     |     |     |
|          580 |             0.04             | 0.997 |     |     |     |     |

Estimates of the probability density of the rotational velocity, P($v_{\rm{e}}$), and its cumulative distribution function, cdf.

All the de-convolutions show a double peak in the equatorial velocity distribution consistent with that observed in the projected rotational velocity distribution. Approximately one quarter of the sample have a rotational velocity of less than 100 km s$^{-1}$ and there appears to be another peak 250-300 km s$^{-1}$ although this is complicated by the varying degree of small scale structure found in the different de-convolutions.

We have attempted to fit the de-convolved distribution using analytical functions. Initially two Gaussian profiles were adopted but these were found to give a relatively poor fit to our estimated distribution. Two Maxwellian functions as discussed by were also considered but would not have reproduced the significant value of P($v_{\rm{e}}$) as $v_{\rm{e}}$$\rightarrow 0$. We therefore recommend that either the values listed in Table or analytical fits appropriate to the specific applications are adopted.

The most distinctive feature of the $v_{\rm{e}} \sin i$ distribution is its bi-modal nature, which is also present in the estimation of the rotational velocity distribution, illustrated in Fig. . Although our sample contains approximately 300 targets, it will still be subject to significant random sampling errors due to its finite size. In Fig. , the population of the two 40 km s$^{-1}$ bins with the lowest projected rotational velocities have values of approximately 40 corresponding to an estimated standard error of approximately 6. Hence even decreasing the populations of both these bins by twice this estimated error (note that if these bins had been overpopulated our estimate of the standard error would also be too high) would not remove this bi-modal behaviour. Additionally Kolmogorov-Smirnov tests using uni-modal distributions (for example a single Gaussian fit to our estimated rotational velocity distribution) lead to very small probabilities that they were the parent populations. Hence we conclude that the projected rotational velocity distribution (and hence the underlying rotational velocity distribution) is bi-modal, although we accept that our sample size limits the information on the structure of the two components.

In Sect. , we discussed the projected rotational velocities in two other LMC samples. Those of showed some evidence for a bimodal distribution but that of appeared unimodal. However given the differences in the samples in terms of age, fraction of field stars and undetected binaries, we do not consider this discrepancy to be significant.

# Origin of bi-modal distribution

The bi-modal distribution of rotational velocities implies that our stellar sample is composed of, at least, two different components. Which physical processes lead to the existence of these two components? Why do most B-type main sequence stars rotate fast, while about a quarter or so rotate slowly? Whatever the answers to these questions, they are not contained in the standard evolutionary theories of single stars, which predicts only small changes of the surface rotational velocity during core hydrogen burning .

The two components, which we find, could either correspond to different ages, different star formation conditions, or emerge as a consequence of differences in the evolution of stars (e.g., in close binaries). We emphasize again that age differences can only lead to the difference in the mean rotation rate of the two components unless the standard stellar evolution picture is wrong or incomplete. Differences in age and star formation conditions could produce stellar components with different spatial or kinematic properties. We inspect our data in this respect in the next two subsections. Thereafter, we discuss the possible impact of non-standard stellar evolution.

## Spatial variations

The 30 Doradus region has had a complex star-formation history over the past few tens of Myrs, as demonstrated by the five distinct stellar groups identified by . For instance, the VFTS observations include targets in the cluster Hodge 301, which is $\sim$`<!-- -->`{=html}20 Myr older than the central parts of 30 Doradus . Indeed, even in the central 30 Doradus cluster (NGC 2070), identified three distinct temporal components, with recently identifying two spatially-distinct components which appear to have an age difference of $\sim$`<!-- -->`{=html}1 Myr. An obvious test is thus to check whether the slow rotators concentrate in particular regions of the VFTS field. In Fig , the spatial positions of the sample of our targets are illustrated, together with their estimated projected rotational velocity. The concentrations of stars at specific positions follows directly from the morphology of the region.

From this figure, there appears to be no evidence for any spatial segregation with projected rotational velocity. For example, the percentage of targets with values of $v_{\rm{e}} \sin i$ less than 100 km s$^{-1}$ within 0.1 degrees of the central concentration of stars associated with Hodge 301 is 33%, in excellent agreement with the 33% found for the whole sample.

To further assess any spatial variation, the sample was split into five velocity bins, each of 100km s$^{-1}$, except the top bin which encompassed all objects with $v_{\rm{e}} \sin i$$\geq 400$km s$^{-1}$. Cumulative radial-distribution functions were calculated for each bin. Two-sample Kolmogorov-Smirnov (KS) tests were then performed on all possible pairs of the radial-distribution functions. No statistically significant variation was found between the pairings, suggesting the spatial distribution does not influence the $v_{\rm{e}} \sin i$ distribution.

In this context it is interesting to consider the results from for projected rotational velocities of 24 OB-type stars in R136. Their sample stars were located within $\sim$`<!-- -->`{=html}25$''$ of the centre of R136 (equivalent to a projected distance of $\sim$`<!-- -->`{=html}6 pc), i.e. the majority of their stars lie within the dense central part of the FLAMES field that was not observed with the Medusa fibres. From comparisons with results for other samples in the LMC (including those from Hunter et al. 2008b), Wolff et al. concluded that the mean rotational velocity of the massive stars in R136 is larger than that found from other samples in the LMC, both for stars in (less dense) clusters and in the field population.

noted that R136 lacks the population of very slowly-rotating stars seen elsewhere (and which also features in the results presented here). They divided their sample into high- and low-mass bins, with inferred evolutionary masses in the range of 15 $\le$ $M_\ast$ $\le$ 30 $M_\odot$ and 6 $\le$ $M_\ast$ $\le$ 12 $M_\odot$, respectively. The majority of the stars in the lower-mass bin lie below the faint magnitude cut-off for the VFTS, but the lack of slow rotators is also evident in their high-mass sub-sample (see their Fig. 11) – they argued that these results might reflect the initial conditions within which the stars formed, i.e. the relationship of rotational velocity with stellar density. As noted above, we see no compelling spatial trends in the results from our significantly larger VFTS sample, with the caveat that our results do not include stars in the central few pc of R136.

## Radial velocity variations

A bi-modality found in Galactic targets by has been interpreted as due to a field star and a cluster population. Although we see no evidence for this in the spatial distribution of targets, we have also investigated the estimated stellar radial velocities of our sample. These have been measured Kennedy et al. (in prep.) using profile fitting to selected spectral features with a methodology similar to that discussed by . We have divided our sample into two groups with projected rotational velocities less than or greater than 100 km s$^{-1}$. The former will be a mixture of slowly rotating stars and those viewed at low angles of inclination; however the contamination should be relatively small given that the mean value of $\sin i$ for random inclinations is approximately 0.78. The latter will be stars with rotational velocities in excess of 100 km s$^{-1}$.

|  |  |  |
|:---|:--:|:--:|
|  |  |  |
| Sample | Number | $\langle v_{\rm rad} \rangle$ km s$^{-1}$ |
|  |  |  |
| $v_{\rm{e}} \sin i$$< 100$ km s$^{-1}$ | 89 | 269$\pm$`<!-- -->`{=html}17 |
| $v_{\rm{e}} \sin i$$\geq100$ km s$^{-1}$ | 172 | 267$\pm$`<!-- -->`{=html}18 |
|  |  |  |
| All | 261 | 267$\pm$`<!-- -->`{=html}17 |

Average radial velocities for sub-samples with $vsini$ smaller and larger than 100 km/s (cf. Fig ), compared to the full sample.

In Table , we summarize the mean and standard deviations for the different samples, together with the sample sizes. Although reliable values were not available for approximately 10% of the targets, we do not expect this to introduce a significant bias. For the two samples, both the means and standard deviations are effectively identical. Coupled with the results presented in Sect. , we conclude that there is no evidence that the bi-modality arises from different populations, due to either sequential star formation or from contamination with field stars.

## Binary evolution

The fraction of Galactic massive close binaries is high, i.e., at least 50% . Similar high fractions have been found for the O- and B-type VFTS samples . Since the time evolution of the rotation of massive main sequence stars can be strongly affected by the interaction with a companion , this could in principle affect the distribution of rotation rates of our sample. The reason is that, while radial velocity variables have been eliminated from our sample, found that the vast majority of post-interaction objects are either single stars or appear as such.

De Mink et al. (2012b) showed that for constant star formation, the fraction of apparently single post-binary interaction objects can be as high as 30%. This fraction would be considerably larger in a sample in which radial velocity variables are removed, as except for a small fraction of Algol-type systems, those are almost exclusively pre-interaction binaries. For the current sample, for which the age distribution is not yet known, on one hand the fraction of apparently single close binary products is expected to be smaller than 30% due to the expected presence of early B-type stars which are coeval with the core region of 30 Doradus. On the other hand, this fraction is enhanced due to the removal of the radial velocity variables from our sample. In conclusion, at present a substantial contribution of stars that are the product of binary evolution, can not be excluded.

The binary population synthesis models of imply that binary interactions, i.e., mass and angular momentum transfer and stellar mergers, produce a population of fast rotators, with rotational velocities up to about 600 km/s. In principle, binary interaction can also produce slow rotators, for example, as a result of spin down by tides. However, find that the effect of this on a stellar population is very small. The low velocity peak of our distribution (see Fig. 8) contains about 25%, which is too many to be explained solely as the products of binary evolution. Indeed, to produce the low velocity component of our distribution by close binary effects could require the hypothesis that mass transfer and/or mergers lead to significant magnetic fields in the binary remnant (cf., Sect. 5.4) as supported by the recent discovery by of a strong magnetic field in the companion of Plaskett's star.

## Magnetic fields

have undertaken an extensive study of the rotational velocities in Galactic late B-type and A-type stars, using a similar methodology to estimate rotational velocity distributions. Their large dataset allowed them to study variations as a function of age and mass. They find that stars with masses $\la 2.5$  $\rm{M}_\odot$ have a uni-modal distribution of rotational velocities, whilst the more massive stars (2.5-3.9 $\rm{M}_\odot$) show a bi-modal structure. fit the latter with two Maxwellian functions and find that the percentage of targets in the low velocity Maxwellian varies from 8-20% depending on the mass range considered. Additionally they find that the peak of the higher velocity Maxwellian lies in the range 200$\la$$v_{\rm{e}}$$\la$`<!-- -->`{=html}240 km s$^{-1}$ for the highest masses they consider.

These results are not directly comparable with those presented here since sampled different field and cluster environments. Additionally our sample is not large enough to be subdivided into smaller mass ranges, while we do not have age estimates for our targets. However, both the fraction of slow rotators (25% for our fit; $8-20$% in the higher mass samples of ) and the peak value of the high velocity component (280 km s$^{-1}$ and 200-240 km s$^{-1}$  respectively) are similar.

discuss two possible explanations for the bi-modality in their distributions. The first is that a fraction of stars loses angular momentum via tidal breaking within a binary system. For our sample, we can exclude this scenario, as show that the number of affected stars is much too small. The second possibility discussed by is that all stars were born rotating rapidly, but that a fraction of them had or acquired magnetic fields. Indeed, the fraction of magnetic main sequence stars in the mass range 2.5-3.9 $\rm{M}_\odot$is believed to be about 15% , whilst their rotation may be slowed down by magnetic breaking.

It has become evident in recent years that magnetic fields could potentially play a significant role in the evolution of some massive stars. In the LMC sample of early B-type stars discussed in Sect. , found a sub-sample of 15% of their stars with slow rotation and high nitrogen enrichment. While direct magnetic field measurements in these stars are not yet feasible, and have identified Galactic analogues for which the incidence of magnetic fields is indeed high. find magnetic fields in $\sim$`<!-- -->`{=html}10% of their early B stars.

and have explored stellar evolution models for magnetic massive main sequence stars. Based on a model for magnetic braking of , they both find that magnetic massive main sequence models can in fact spin-down for reasonable field strength. performed a population synthesis matched to the B-type stellar sample of . They produced a bi-modal velocity distribution, with a significant component of slow rotators ($v_e < 60$ km/s). However their distribution of fast rotators peaks at $v_e \simeq 100$ km/s (possibly due to their choice of the initial rotational velocity distribution), compared to $v_e \simeq 250$ km/s for our sample (Fig. ).

Of course, if magnetic spin down was responsible for the low-velocity component in our rotational velocity distribution, one may wonder about the origin of the magnetic fields in these main sequence stars. Three possibilities are being discussed in the literature. The fields could be fossil, i.e. inherited from the star formation process . Alternatively, the fields could be generated by strong binary interaction, i.e. binary merger, mass transfer or common evelope evolution . Finally, the fields could be generated by a dynamo process inside the the main sequence star . We emphasize that these three alternatives predict different time-dependances of the stellar rotation and magnetic field incidence, which might be testable if the age distribution of our sample stars can be established.

An additional important braking mechanism could be mass loss via stellar winds. Winds are thought to be ubiquitous amongst massive stars, and despite the spectroscopic evidence for high-speed outflows being most prevalant amongst the most luminous O-type objects (in the form of P Cygni ultraviolet lines), they are also observed in the lower mass B-type stars. Although mass-loss rates are higher in O stars, wind braking could be relevant for B-type stars given their longer evolutionary timescales. One interesting possibility is that of bi-stability braking , which might lead to a population of both slow and fast rotators, and hence play a role in the observed bi-modality. This could be tested when the stellar parameters of our sample stars have been established.

# Conclusions

From the VFTS B-type stellar sample (approximately 540 stars), we have constructed a sample 289 early B-type stars (plus 45 stars of uncertain binary status) by removing radial velocity variable stars and supergiants. These will have masses in the range 6 to 16  $\rm{M}_\odot$and are centered on the 30 Doradus region in the LMC. We have estimated their projected rotational velocities, and find that their distribution shows a distinct bi-modality, which neither correlates with position in the field nor with radial velocity.

We have deconvolved the $v_{\rm{e}} \sin i$ distribution to obtain the most likely underlying distribution of intrinsic rotational velocities. This shows a low velocity component containing about one quarter of the sample stars, with $v_{\rm e} <$ 100 km/s, and a second component of stars rotating predominantly with velocities in the range 200 km/s$< v_{\rm e} <$ 350 km/s.

This distribution is not inconsistent with previous LMC surveys, despite a bi-modality having not been clearly demonstrated before. A conclusive interpretation of our finding is not yet possible, in particular because the age distribution and chemical composition of the sample stars are still unknown. An origin of the bi-modality due to evolutionary effects relating to binarity and/or magnetic fields or that the velocity distribution was inherited in the star formation process remain possibilities.

[^1]: Hubble Fellow