The Large Underground Xenon (LUX) dark matter search was a 250-kg active mass dual-phase time projection chamber that operated by detecting light and ionization signals from particles incident on a xenon target. In December 2015, LUX reported a minimum 90% upper C.L. of 6e-46 cm^2 on the spin-independent WIMP-nucleon elastic scattering cross section based on a 1.4e4 kg*day exposure in its first science run. Tension between experiments and the absence of a definitive positive detection suggest it would be prudent to search for WIMPs outside the standard spin-independent/spin-dependent paradigm. Recent theoretical work has identified a complete basis of 14 independent effective field theory (EFT) operators to describe WIMP-nucleon interactions. In addition to spin-independent and spin-dependent nuclear responses, these operators can produce novel responses such as angular-momentum-dependent and spin-orbit couplings. Here we report on a search for all 14 of these EFT couplings with data from LUX's first science run. Limits are placed on each coupling as a function of WIMP mass.

We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

The scattering of dark matter (DM) particles with sub-GeV masses off nuclei is difficult to detect using liquid xenon-based DM search instruments because the energy transfer during nuclear recoils is smaller than the typical detector threshold. However, the tree-level DM-nucleus scattering diagram can be accompanied by simultaneous emission of a Bremsstrahlung photon or a so-called "Migdal" electron. These provide an electron recoil component to the experimental signature at higher energies than the corresponding nuclear recoil. The presence of this signature allows liquid xenon detectors to use both the scintillation and the ionization signals in the analysis where the nuclear recoil signal would not be otherwise visible. We report constraints on spin-independent DM-nucleon scattering for DM particles with masses of 0.4-5 GeV/c$^2$ using 1.4$\times10^4$ kg$\cdot$day of search exposure from the 2013 data from the Large Underground Xenon (LUX) experiment for four different classes of mediators. This analysis extends the reach of liquid xenon-based DM search instruments to lower DM masses than has been achieved previously.

We report results from an extensive set of measurements of the \b{eta}-decay response in liquid xenon.These measurements are derived from high-statistics calibration data from injected sources of both $^{3}$H and $^{14}$C in the LUX detector. The mean light-to-charge ratio is reported for 13 electric field values ranging from 43 to 491 V/cm, and for energies ranging from 1.5 to 145 keV.

The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged]

We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest (($B-V$) $>$ 1 mag) and brightest ($V$ $<$ 16 mag) objects in the VFTS. We use optical and near-IR photometry to estimate their temperatures and luminosities, and introduce the luminosity-age diagram to estimate their ages. We derive physical parameters for our targets, including temperatures from a new calibration of $(J-K_{\rm s})_{0}$ colour for luminous cool stars in the LMC, luminosities from their $J$-band magnitudes (thence radii), and ages from comparisons with state-of-the-art evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12-24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. We argue that the effects of mass-transfer in binaries may lead to more massive and luminous RSGs (which we call `red stragglers') than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24$^{+5}_{-3}$ Myr for Hodge 301, 22$^{+6}_{-5}$ Myr for SL 639, and 23$^{+4}_{-2}$ Myr for NGC 2100.

Various dark matter models predict annual and diurnal modulations of dark matter interaction rates in Earth-based experiments as a result of the Earth's motion in the halo. Observation of such features can provide generic evidence for detection of dark matter interactions. This paper reports a search for both annual and diurnal rate modulations in the LUX dark matter experiment using over 20 calendar months of data acquired between 2013 and 2016. This search focuses on electron recoil events at low energies, where leptophilic dark matter interactions are expected to occur and where the DAMA experiment has observed a strong rate modulation for over two decades. By using the innermost volume of the LUX detector and developing robust cuts and corrections, we obtained a stable event rate of 2.3$\pm$0.2~cpd/keV$_{\text{ee}}$/tonne, which is among the lowest in all dark matter experiments. No statistically significant annual modulation was observed in energy windows up to 26~keV$_{\text{ee}}$. Between 2 and 6~keV$_{\text{ee}}$, this analysis demonstrates the most sensitive annual modulation search up to date, with 9.2$\sigma$ tension with the DAMA/LIBRA result. We also report no observation of diurnal modulations above 0.2~cpd/keV$_{\text{ee}}$/tonne amplitude between 2 and 6~keV$_{\text{ee}}$.

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

ERIS is an instrument that will both extend and enhance the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It will replace two instruments that are now being maintained beyond their operational lifetimes, combine their functionality on a single focus, provide a new wavefront sensing module that makes use of the facility Adaptive Optics System, and considerably improve their performance. The instrument will be competitive with respect to JWST in several regimes, and has outstanding potential for studies of the Galactic Center, exoplanets, and high redshift galaxies. ERIS had its final design review in 2017, and is expected to be on sky in 2020. This contribution describes the instrument concept, outlines its expected performance, and highlights where it will most excel.

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities <80 \kms and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars from a previous survey have been re-evaluated. Approximately 75-80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities <40 \kms and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. Adopting an enhancement threshold of 0.6 dex, we observed five stars, yet stellar evolution models with rotation predict only 1.25$\pm$1.11 for random stellar viewing inclinations. The excess of such objects is estimated to be 20-30% of all stars with current rotational velocities of less than 40 \kms and 2-4% of the total single B-type sample. These estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. 70% of the stars with rotational velocities less than 40 \kms and 6-8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking.

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]

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

We report an absolute calibration of the ionization yields($\textit{Q$_y$})$ and fluctuations for electronic recoil events in liquid xenon at discrete energies between 186 eV and 33.2 keV. The average electric field applied across the liquid xenon target is 180 V/cm. The data are obtained using low energy $^{127}$Xe electron capture decay events from the 95.0-day first run from LUX (WS2013) in search of Weakly Interacting Massive Particles (WIMPs). The sequence of gamma-ray and X-ray cascades associated with $^{127}$I de-excitations produces clearly identified 2-vertex events in the LUX detector. We observe the K- (binding energy, 33.2 keV), L- (5.2 keV), M- (1.1 keV), and N- (186 eV) shell cascade events and verify that the relative ratio of observed events for each shell agrees with calculations. The N-shell cascade analysis includes single extracted electron (SE) events and represents the lowest-energy electronic recoil $\textit{in situ}$ measurements that have been explored in liquid xenon.

The first searches for axions and axion-like particles with the Large Underground Xenon (LUX) experiment are presented. Under the assumption of an axio-electric interaction in xenon, the coupling constant between axions and electrons, gAe is tested, using data collected in 2013 with an exposure totalling 95 live-days $\times$ 118 kg. A double-sided, profile likelihood ratio statistic test excludes gAe larger than 3.5 $\times$ 10$^{-12}$ (90% C.L.) for solar axions. Assuming the DFSZ theoretical description, the upper limit in coupling corresponds to an upper limit on axion mass of 0.12 eV/c$^{2}$, while for the KSVZ description masses above 36.6 eV/c$^{2}$ are excluded. For galactic axion-like particles, values of gAe larger than 4.2 $\times$ 10$^{-13}$ are excluded for particle masses in the range 1-16 keV/c$^{2}$. These are the most stringent constraints to date for these interactions.

We present experimental constraints on the spin-dependent WIMP-nucleon elastic cross sections from the total 129.5 kg-year exposure acquired by the Large Underground Xenon experiment (LUX), operating at the Sanford Underground Research Facility in Lead, South Dakota (USA). A profile likelihood ratio analysis allows 90% CL upper limits to be set on the WIMP-neutron (WIMP-proton) cross section of $\sigma_n$ = 1.6$\times 10^{-41}$ cm$^{2}$ ($\sigma_p$ = 5$\times 10^{-40}$ cm$^{2}$) at 35 GeV$c^{-2}$, almost a sixfold improvement over the previous LUX spin-dependent results. The spin-dependent WIMP-neutron limit is the most sensitive constraint to date.

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.

In this Technical Design Report (TDR) we describe the LZ detector to be built at the Sanford Underground Research Facility (SURF). The LZ dark matter experiment is designed to achieve sensitivity to a WIMP-nucleon spin-independent cross section of three times ten to the negative forty-eighth square centimeters.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs.

We report constraints on spin-independent weakly interacting massive particle (WIMP)-nucleon scattering using a 3.35e4 kg-day exposure of the Large Underground Xenon (LUX) experiment. A dual-phase xenon time projection chamber with 250 kg of active mass is operated at the Sanford Underground Research Facility under Lead, South Dakota (USA). With roughly fourfold improvement in sensitivity for high WIMP masses relative to our previous results, this search yields no evidence of WIMP nuclear recoils. At a WIMP mass of 50 GeV/c^2, WIMP-nucleon spin-independent cross sections above 2.2e-46 cm^2 are excluded at the 90% confidence level. When combined with the previously reported LUX exposure, this exclusion strengthens to 1.1e-46 cm^2 at 50 GeV/c^2.

The Large Underground Xenon (LUX) experiment is a dual-phase liquid xenon time projection chamber (TPC) operating at the Sanford Underground Research Facility in Lead, South Dakota. A calibration of nuclear recoils in liquid xenon was performed $\textit{in situ}$ in the LUX detector using a collimated beam of mono-energetic 2.45 MeV neutrons produced by a deuterium-deuterium (D-D) fusion source. The nuclear recoil energy from the first neutron scatter in the TPC was reconstructed using the measured scattering angle defined by double-scatter neutron events within the active xenon volume. We measured the absolute charge ($Q_{y}$) and light ($L_{y}$) yields at an average electric field of 180 V/cm for nuclear recoil energies spanning 0.7 to 74 keV and 1.1 to 74 keV, respectively. This calibration of the nuclear recoil signal yields will permit the further refinement of liquid xenon nuclear recoil signal models and, importantly for dark matter searches, clearly demonstrates measured ionization and scintillation signals in this medium at recoil energies down to $\mathcal{O}$(1 keV).

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

We propose a new technique for the calibration of nuclear recoils in large noble element dual-phase time projection chambers used to search for WIMP dark matter in the local galactic halo. This technique provides an $\textit{in situ}$ measurement of the low-energy nuclear recoil response of the target media using the measured scattering angle between multiple neutron interactions within the detector volume. The low-energy reach and reduced systematics of this calibration have particular significance for the low-mass WIMP sensitivity of several leading dark matter experiments. Multiple strategies for improving this calibration technique are discussed, including the creation of a new type of quasi-monoenergetic 272 keV neutron source. We report results from a time-of-flight based measurement of the neutron energy spectrum produced by an Adelphi Technology, Inc. DD108 neutron generator, confirming its suitability for the proposed nuclear recoil calibration.

We present constraints on weakly interacting massive particles (WIMP)-nucleus scattering from the 2013 data of the Large Underground Xenon dark matter experiment, including $1.4\times10^{4}\;\mathrm{kg\; day}$ of search exposure. This new analysis incorporates several advances: single-photon calibration at the scintillation wavelength, improved event-reconstruction algorithms, a revised background model including events originating on the detector walls in an enlarged fiducial volume, and new calibrations from decays of an injected tritium $\beta$ source and from kinematically constrained nuclear recoils down to 1.1 keV. Sensitivity, especially to low-mass WIMPs, is enhanced compared to our previous results which modeled the signal only above a 3 keV minimum energy. Under standard dark matter halo assumptions and in the mass range above 4 $\mathrm{GeV}\,c^{-2}$, these new results give the most stringent direct limits on the spin-independent WIMP-nucleon cross section. The 90% C.L. upper limit has a minimum of 0.6 zb at 33 $\mathrm{GeV}\,c^{-2}$ WIMP mass.

We present measurements of the electron-recoil (ER) response of the LUX dark matter detector based upon 170,000 highly pure and spatially-uniform tritium decays. We reconstruct the tritium energy spectrum using the combined energy model and find good agreement with expectations. We report the average charge and light yields of ER events in liquid xenon at 180 V/cm and 105 V/cm and compare the results to the NEST model. We also measure the mean charge recombination fraction and its fluctuations, and we investigate the location and width of the LUX ER band. These results provide input to a re-analysis of the LUX Run3 WIMP search.

The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

Extremely broad emission wings at H$\beta$ and H$\alpha$ have been found in VFTS data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of H$\beta$ and the shortward wing of H$\alpha$. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables. No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.

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

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

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

The advent of extremely large telescopes will bring unprecedented light-collecting power and spatial resolution, but it will also lead to a significant increase in the size and complexity of focal-plane instruments. The use of freeform mirrors could drastically reduce the number of components in optical systems. Currently, manufacturing issues limit the common use of freeform mirrors at short wavelengths. This article outlines the use of freeform mirrors in astronomical instruments with a description of two efficient freeform optical systems. A new manufacturing method is presented which seeks to overcome the manufacturing issues through hydroforming of thin polished substrates. A specific design of an active array is detailed, which will compensate for residual manufacturing errors, thermoelastic deformation, and gravity-induced errors during observations. The combined hydroformed mirror and the active array comprise the Freeform Active Mirror Experiment, which will produce an accurate, compact, and stable freeform optics dedicated to visible and near-infrared observations.

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30\,Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{km s^{-1}}$. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300\,$\rm{km s^{-1}}$, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part -- and could potentially be completely due -- to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

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

The blue supergiant Sher25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula, which shows similarities to the triple-ring structure seen around SN1987A. From optical spectroscopy over six consecutive nights, we detect periodic radial velocity variations in the stellar spectrum of Sher25 with a peak-to-peak amplitude of ~12 km/s on a timescale of about 6 days, confirming the tentative detec-tion of similar variations by Hendry et al. From consideration of the amplitude and timescale of the signal, coupled with observed line profile variations, we propose that the physical origin of these variations is related to pulsations in the stellar atmosphere, rejecting the previous hypothesis of a massive, short-period binary companion. The radial velocities of two other blue supergiants with similar bipolar nebulae, SBW1 and HD 168625, were also monitored over the course of six nights, but these did not display any significant radial velocity variations.

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

Context: The commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+NIR photometry. Aims: We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical+NIR photometry to determine the effective temperature of OB stars. Methods: We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions such as the family of extinction laws. Results: We derive a new family of optical+NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405-5495) ~ 1.5 mag.

Detailed spectral classifications are presented for 352 O-B0 stars in the VLT-FLAMES Tarantula Survey, of which 213 O-type are of sufficient quality for further morphological analysis. Among them, six subcategories of special interest are distinguished. (1) Several new examples of the earliest spectral types O2-O3 have been found. (2) A group of extremely rapidly rotating main-sequence objects has been isolated, including the largest $v\sin i$ values known, the spatial and radial-velocity distributions of which suggest ejection from the two principal ionizing clusters. (3) Several new examples of the evolved, rapidly rotating Onfp class show similar evidence. (4) No fewer than 48 members of the Vz category, hypothesized to be on or near the ZAMS, are found in this sample; in contrast to the rapid rotators, they are strongly concentrated to the ionizing clusters, supporting their interpretation as very young objects, as do their relatively faint absolute magnitudes. (5) A surprisingly large fraction of the main-sequence spectra belong to the recently recognized V((fc)) class, with C III emission lines of similar strength to the usual N III in V((f)) spectra; there are also six objects with very high-quality data but no trace of either mission feature, presenting new challenges to physical interpretations. (6) Five spectra with morphologically enhanced nitrogen lines have been detected. Absolute visual magnitudes have been derived for each star with individual extinction laws, and composite HRDs provide evidence of the multiple generations present in this field. Associations with X-ray sources are noted. Further analyses of this unique dataset underway will provide new insights into the evolution of massive stars and starburst clusters.

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V stars, identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about the O Vz phenomenon. RESULTS. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2-4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. CONCLUSIONS. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{\Delta RV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun also showed high contributions to the global feedback, ~25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

Aims: Projected rotational velocities (\vsini) 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. Methods: 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. Results: Projected rotational velocities range up to approximately 450 \kms and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$\ve$\leq$100\,\kms and the high velocity component having \ve$\sim 250$\,\kms. 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. Conclusions: 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.

The Tarantula Nebula (30 Dor) is a spectacular star-forming region in the Large Magellanic Cloud, seen through gas in the Galactic Disc and Halo. Diffuse Interstellar Bands offer a unique probe of the diffuse, cool-warm gas in these regions. The aim is to use DIBs as diagnostics of the local interstellar conditions, whilst at the same time deriving properties of the yet-unknown carriers. Spectra of over 800 early-type stars from the VLT Flames Tarantula Survey (VFTS) were analysed. Maps were created, separately, for the Galactic and LMC absorption in the DIBs at 4428 and 6614 Ang and - in a smaller region near the central cluster R136 - neutral sodium (Na I D); we also measured the DIBs at 5780 and 5797 Ang. The maps show strong 4428 and 6614 Ang DIBs in the quiescent cloud complex to the south of 30 Dor but weak absorption in the harsher environments to the north (bubbles) and near the OB associations. The Na maps show at least five kinematic components in the LMC and a shell-like structure surrounding R136, and small-scale structure in the Milky Way. The strengths of the 4428, 5780, 5797 and 6614 Ang DIBs are correlated, also with Na absorption and visual extinction. The strong 4428 Ang DIB is present already at low Na column density but the 6614, 5780 and 5797 Ang DIBs start to be detectable at subsequently larger Na column densities. The relative strength of the 5780 and 5797 Ang DIBs clearly confirm the Tarantula Nebula and Galactic high-latitude gas to represent a harsh radiation environment. The resilience of the 4428 Ang DIB suggests its carrier is large, compact and neutral. Structure is detected in the distribution of cool-warm gas on scales between one and >100 pc in the LMC and as little as 0.01 pc in the Sun's vicinity. Stellar winds from the central cluster R136 have created an expanding shell; some infalling gas is also detected, reminiscent of a galactic "fountain".

Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^\pi\ (with 0.15 < log P < 3.5) and f_q ~ q^\kappa\ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are \pi = -0.45 +/- 0.30 and \kappa = -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

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

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 multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type 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. 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.

We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared photometry and X-ray imaging to characterise the system. We find radial velocity variations of around 400 km/s in the high excitation Si IV, N III and He II spectra, and photometric variability of ~0.6 mag with a period of 12.7 days. In addition, we detect long-term photometric variations of ~0.25 mag, which may be due to a longer-term variability with a period of ~400 days. We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbiting a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate from an optically-thick accretion disc directly surrounding the primary.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

The spectral and spatial characteristics of two special categories of O stars found in the VFTS dataset are presented. One of them comprises very rapid rotators, including several more extreme than any previously known. These objects are distributed around the peripheries of the main 30 Doradus clusters, suggesting a runaway nature for which their radial velocities already provide preliminary supporting evidence. The other category consists of a large number of Vz stars, previously hypothesized on spectroscopic grounds to be on or very near the ZAMS. Their distribution is the inverse of that of the rapid rotators: the Vz are strongly concentrated to the ionizing clusters, plus a newly recognized band of recent and current star formation to the north, which provides strong circumstantial evidence for their extreme youth.

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fastest stellar winds known. From analysis of the CIV 1548-51 doublet we find a terminal velocity, v_infty=3450 +/- 50km/s. Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, from which we rule out a massive companion (with 2d<P<1yr) to a confidence of 98%. The radial velocity of #016 is offset from the systemic value by -85km/s, suggesting that the star has traveled the 120pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions.

EAGLE is the multi-object, spatially-resolved, near-IR spectrograph instrument concept for the E-ELT, relying on a distributed Adaptive Optics, so-called Multi Object Adaptive Optics. This paper presents the results of a phase A study. Using 84x84 actuator deformable mirrors, the performed analysis demonstrates that 6 laser guide stars and up to 5 natural guide stars of magnitude R<17, picked-up in a 7.3' diameter patrol field of view, allow us to obtain an overall performance in terms of Ensquared Energy of 35% in a 75x75 mas^2 spaxel at H band, whatever the target direction in the centred 5' science field for median seeing conditions. The computed sky coverage at galactic latitudes |b|~60 is close to 90%.

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

We prove that inflation is forbidden in the most well understood class of semi-realistic type IIA string compactifications: Calabi-Yau compactifications with only standard NS-NS 3-form flux, R-R fluxes, D6-branes and O6-planes at large volume and small string coupling. With these ingredients, the first slow-roll parameter satisfies epsilon >= 27/13 whenever V > 0, ruling out both inflation (including brane/anti-brane inflation) and de Sitter vacua in this limit. Our proof is based on the dependence of the 4-dimensional potential on the volume and dilaton moduli in the presence of fluxes and branes. We also describe broader classes of IIA models which may include cosmologies with inflation and/or de Sitter vacua. The inclusion of extra ingredients, such as NS 5-branes and geometric or non-geometric NS-NS fluxes, evades the assumptions used in deriving the no-go theorem. We focus on NS 5-branes and outline how such ingredients may prove fruitful for cosmology, but we do not provide an explicit model. We contrast the results of our IIA analysis with the rather different situation in IIB.

We present an overview of the EAGLE science case, which spans spatially-resolved spectroscopy of targets from five key science areas - ranging from studies of heavily-obscured Galactic star clusters, right out to the first galaxies at the highest redshifts. Here we summarise the requirements adopted for study and also evaluate the availability of natural guide stars in example fields, which will impact on the adaptive optics performance and architecture.