Recently we have demonstrated that high-precision polarization observations can detect the polarization resulting from the rotational distortion of a rapidly rotating B-type star. Here we investigate the extension of this approach to an A-type star. Linear-polarization observations of $\alpha$ Oph (A5IV) have been obtained over wavelengths from 400 to 750 nm. They show the wavelength dependence expected for a rapidly-rotating star combined with a contribution from interstellar polarization. We model the observations by fitting rotating-star polarization models and adding additional constraints including a measured $v_e \sin{i}$. However, we cannot fully separate the effects of rotation rate and inclination, leaving a range of possible solutions. We determine a rotation rate $\omega = \Omega/\Omega_ c$ between 0.83 and 0.98 and an axial inclination i > 60 deg. The rotation-axis position angle is found to be 142 $\pm$ 4 deg, differing by 16 deg from a value obtained by interferometry. This might be due to precession of the rotation axis due to interaction with the binary companion. Other parameters resulting from the analysis include a polar temperature Tp = 8725 $\pm$ 175 K, polar gravity $\log{g_p} = 3.93 \pm 0.08$ (dex cgs), and polar radius $R_{\rm p} = 2.52 \pm 0.06$ Rsun. Comparison with rotating-star evolutionary models indicates that $\alpha$ Oph is in the later half of its main-sequence evolution and must have had an initial $\omega$ of 0.8 or greater. The interstellar polarization has a maximum value at a wavelength ($\lambda_{\rm max}$) of $440 \pm 110$ nm, consistent with values found for other nearby stars.

Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, organised (fossil) magnetic fields at their surfaces, little is known about the fields of very massive stars, nor the magnetic properties of stars outside our Galaxy. We aim to continue searching for strong magnetic fields in a diverse set of massive and very massive stars (VMS) in the Large and Small Magellanic Clouds (LMC/SMC), and we evaluate the overall capability of FORS2 to usefully search for and detect stellar magnetic fields in extra-galactic environments. We have obtained FORS2 spectropolarimetry of a sample of 41 stars, which principally consist of spectral types B, O, Of/WN, WNh, and classical WR stars in the LMC and SMC. Four of our targets are Of?p stars; one of them was just recently discovered. Each spectrum was analysed to infer the longitudinal magnetic field. No magnetic fields were formally detected in our study, although Bayesian statistical considerations suggest that the Of?p star SMC159-2 is magnetic with a dipolar field of the order of 2.4 to 4.4kG. In addition, our first constraints of magnetic fields in VMS provide interesting insights into the formation of the most massive stars in the Universe.

The ON stars are a rare subtype of O stars, of uncertain origin. We report two new, rapidly-rotating ON stars found in data acquired with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST. LS I +61 28 is an ON8.5 Vn dwarf with a projected equatorial rotational velocity of $v_{\textrm e}\sin{i} \simeq 298$ km$\;$s$^{-1}$, while HDE 236672 is an ON9 IVn subgiant with $v_{\textrm e}\sin{i} \simeq 253$ km$\;$s$^{-1}$ The former is the first rapidly-rotating ON dwarf to be found, and the latter is only the third ON subgiant. The luminosity classes of non-supergiant ON stars appear to be influenced by the axial inclination angle $i$: the rapidly-rotating giants are close to equator-on, while ON dwarfs with lower $v_{\textrm e}\sin{i}$ values are viewed more nearly pole-on. Combining parallaxes and proper motions from Gaia DR2 with radial-velocity measurements, we investigate the kinematics of non-supergiant ON stars, and infer that the dynamics, rapid rotation, and surface-nitrogen characteristics may all be consequences of binary interaction.

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

We scrutinize the Hipparcos parallax for the bright O supergiant zeta Pup, and confirm that the implied distance of 332+/-11pc appears to be reliable. We then review the implications for the star's physical parameters, and the consequences for the interpretation of P(phot), the 1.78-d photometric period. The equatorial rotation period is <3.7d (with 95% confidence), ruling out a proposed ~5.1d value. If the photometric period is the rotation period then i, the inclination of the rotation axis to the line of sight, is 33.2+/-1.8 degrees. The inferred mass, radius, and luminosity are securely established to be less than canonical values for the spectral type, and are not in agreement with single-star evolution models. The runaway status, rapid rotation, and anomalous physical properties are all indicative of an evolutionary history involving binary (or multiple-star) interaction. We perform simple starspot modelling to show that the low axial inclination required if P(rot) = 1.78d has testable spectroscopic consequences, which have not been identified in existing time series. If P(phot) is directly related to drivers of systematic, high-velocity stellar-wind variability (`discrete absorption components') in zeta Pup, antisolar differential rotation is required. Model line profiles calculated on that basis are at variance with observations.

We report new spectroscopic observations of the magnetic Of?p star HD 148937 obtained since 2015 that differ qualitatively from its extensive historical record of weak, periodic spectral variations. This remarkable behaviour represents clear evidence for an unprecedented change in the character of variability of the star. In this paper we describe the new spectral properties and compare them to the previous line profiles. Based on measurements of the radial velocities of the C iii/N iii} emission lines near 4640 \AA\ and the C iv absorption lines near 5800 \AA, we infer that HD 148937 is likely a high-mass, double-lined spectroscopic binary. Combining the spectroscopic orbit with an archival interferometric measurement of the apparent separation of the equal-brightness components, we tentatively conclude that HD 148937 consists of two O-type stars with masses of approximately 34 and $49~M_\odot$, orbiting in an eccentric ($e=0.75$), long-period ($P_{\rm orb}\sim 26$ y) orbit. We discuss the potential relationship of the binary system to the peculiar properties of HD 148937, and propose future observations to refine the orbital and stellar properties.

We examine claims of strong gravity-darkening effects in the O-star components of WR+O binaries. We generate synthetic spectra for a wide range of parameters, and show that the line-width results are consistent with extensive measurements of O stars that are either single or are members of `normal' binaries. By contrast, the WR+O results are at the extremes of, or outside, the distributions of both models and other observations. Remeasurement of the WR+O spectra shows that they can be reconciled with other results by judicious choice of pseudo-continuum normalization. With this interpretation, the supersynchronous rotation previously noted for the O-star components in the WR+O binaries with the longest orbital periods appears to be unexceptional. Our investigation is therefore consistent with the aphorism that if the title of a paper ends with a question mark, the answer is probably `no'.

Polarization in stars was first predicted by Chandrasekhar [1] who calculated a substantial linear polarization at the stellar limb for a pure electron-scattering atmosphere. This polarization will average to zero when integrated over a spherical star but could be detected if the symmetry is broken, for example by the eclipse of a binary companion. Nearly 50 years ago, Harrington and Collins [2] modeled another way of breaking the symmetry and producing net polarization - the distortion of a rapidly rotating hot star. Here we report the first detection of this effect. Observations of the linear polarization of Regulus, with two different high-precision polarimeters, range from +42 parts-per-million (ppm) at a wavelength of 741 nm to -22 ppm at 395 nm. The reversal from red to blue is a distinctive feature of rotation-induced polarization. Using a new set of models for the polarization of rapidly rotating stars we find that Regulus is rotating at 96.5$\substack{+0.6-0.8}$% of its critical angular velocity for breakup, and has an inclination greater than 76.5 degrees. The rotation axis of the star is at a position angle of 79.5$\pm$0.7 degrees. The conclusions are independent of, but in good agreement with, the results of previously published interferometric observations of Regulus [3]. The accurate measurement of rotation in early-type stars is important for understanding their stellar environments [4], and course of their evolution [5].

Previously, the X-ray emission of zeta Puppis was found to be variable with light curves harbouring "trends" with a typical timescale longer than the exposure length. The origin of these changes was proposed to be linked to large-scale structures in the wind, but further characterization of the variability at high energies was needed. Since then, a number of new X-ray observations have become available. Furthermore, a cyclic behaviour with a 1.78d period was identified in long optical photometric runs, which is thought to be associated with the launching mechanism of large-scale wind structures. We analysed these new X-ray data, revisited the old data, and compared X-ray with optical data, including when simultaneous. We found that the behaviour in X-rays cannot be explained in terms of a perfect clock because the amplitude and shape of its variations change with time. For example, zeta Puppis was much more strongly variable between 2007 and 2011 than before and after this interval. Comparing the X-ray spectra of the star at maximum and minimum brightness yields no compelling difference beyond the overall flux change: the temperatures, absorptions, and line shapes seem to remain constant, well within errors. The only common feature between X-ray datasets is that the variation amplitudes appear maximum in the medium (0.6-1.2keV) energy band. Finally, no clear and coherent correlation can be found between simultaneous X-ray and optical data. Only a subgroup of observations may be combined coherently with the optical period of 1.78d, although the simultaneous optical behaviour is unknown. The currently available data do not reveal any obvious, permanent, and direct correlation between X-ray and optical variations. The origin of the X-ray variability therefore still needs to be ascertained, highlighting the need for long-term monitoring in multiwavelengths, i.e. X-ray, UV, and optical.

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision can be obtained in modelling only if other astrophysical effects, including the stellar limb-darkening, are accounted for properly. In this paper, we explore the limits on precision due to the mathematical formulas currently adopted to approximate the stellar limb-darkening, and to the use of limb-darkening coefficients obtained either from stellar-atmosphere models or empirically. We propose a new limb-darkening law with two coefficients, `power-2', which outperforms other two-coefficient laws adopted in the literature in most cases, and particularly for cool stars. Empirical limb-darkening based on two-coefficient formulas can be significantly biased, even if the light-curve residuals are nearly photon-noise limited. We demonstrate an optimal strategy to fitting for the four-coefficients limb-darkening in the visible, using prior information on the exoplanet orbital parameters to break some of the degeneracies that otherwise would prevent the convergence of the fit. Infrared observations taken with the James Webb Space Telescope (JWST) will provide accurate measurements of the exoplanet orbital parameters with unprecedented precision, which can be used as priors to improve the stellar limb-darkening characterization, and therefore the inferred exoplanet parameters, from observations in the visible, such as those taken with Kepler/K2, JWST, other past and future instruments.

We analyse Kepler light-curves of the exoplanet KOI-13b transiting its moderately rapidly rotating (gravity-darkened) parent star. A physical model, with minimal ad hoc free parameters, reproduces the time-averaged light-curve at the ca. 10 parts per million level. We demonstrate that this Roche-model solution allows the absolute dimensions of the system to be determined from the star's projected equatorial rotation speed, v(e)sin(i), without any additional assumptions; we find a planetary radius 1.33+/-0.05 R(Jup), stellar polar radius 1.55+/-0.06 R(sun), combined mass M(*) + M(P) (\simeq M*) = 1.47 +/- 0.17 M(sun), and distance d \simeq 370+/-25 pc, where the errors are dominated by uncertainties in relative flux contribution of the visual-binary companion KOI-13B. The implied stellar rotation period is within ca. 5% of the non-orbital, 25.43-hr signal found in the Kepler photometry. We show that the model accurately reproduces independent tomographic observations, and yields an offset between orbital and stellar-rotation angular-momentum vectors of 60.25+/-0.05 degrees.

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.

Massive O-type stars play a dominant role in our Universe, but many of their properties remain poorly constrained. In the last decade magnetic fields have been detected in all Galactic members of the distinctive Of?p class, opening the door to a better knowledge of all O-type stars. With the aim of extending the study of magnetic massive stars to nearby galaxies, to better understand the role of metallicity in the formation of their magnetic fields and magnetospheres, and to broaden our knowledge of the role of magnetic fields in massive star evolution, we have carried out spectropolarimetry of five extra-Galactic Of?p stars, as well as a couple of dozen neighbouring stars. We have been able to measure magnetic fields with typical error bars from 0.2 to 1.0 kG, depending on the apparent magnitude and on weather conditions. No magnetic field has been firmly detected in any of our measurements, but we have been able to estimate upper limits to the field values of our target stars. One of our targets, 2dFS 936, exhibited an unexpected strengthening of emission lines. We confirm the unusual behaviour of BI 57, which exhibits a 787 d period with two photometric peaks and one spectroscopic maximum. The observed strengthening of the emission lines of 2dFS 936, and the lack of detection of a strong magnetic field in a star with such strong emission lines is at odd with expectations. Together with the unusual periodic behaviour of BI 57, it represents a challenge for the current models of Of?p stars. The limited precision that we obtained in our field measurements (in most cases as a consequence of poor weather) has led to field-strength upper limits that are substantially larger than those typically measured in Galactic magnetic O stars. Further higher precision observations and monitoring are clearly required.

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

Some O and B stars show unusually strong or weak lines of carbon and/or nitrogen. These objects are classified as OBN or OBC stars. It has recently been shown that nitrogen enrichment and carbon depletion are the most likely explanations for the existence of the ON class. We investigate OC stars (all being supergiants) to check that surface abundances are responsible for the observed anomalous line strengths. We perform a spectroscopic analysis of three OC supergiants using atmosphere models. A fourth star was previously studied by us. Our sample thus comprises all OC stars known to date in the Galaxy. We determine the stellar parameters and He, C, N, and O surface abundances. We show that all stars have effective temperatures and surface gravities fully consistent with morphologically normal O supergiants. However, OC stars show little, if any, nitrogen enrichment and carbon surface abundances consistent with the initial composition. OC supergiants are thus barely chemically evolved, unlike morphologically normal O supergiants.

Published photometry of fading events in the PTFO 8-8695 system is modelled using improved treatments of stellar geometry, surface intensities, and, particularly, gravity darkening, with a view to testing the planetary-transit hypothesis. Variability in the morphology of fading events can be reproduced by adopting convective-envelope gravity darkening, but near-critical stellar rotation is required. This leads to inconsistencies with spectroscopic observations; the model also predicts substantial photometric variability associated with stellar precession, contrary to observations. Furthermore, the empirical ratio of orbital to rotational angular momenta is at odds with physically plausible values. An exoplanet transiting a precessing, gravity-darkened star may not be the correct explanation of periodic fading events in this system.

The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the T\'elescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications.

We present grids of limb-darkening coefficients computed from non-LTE, line-blanketed TLUSTY model atmospheres, covering effective-temperature and surface-gravity ranges of 15--55kK and 4.75 dex (cgs) down to the effective Eddington limit, at 1x, 1x, 0.5x (LMC), 0.2x (SMC), and 0.1x solar. Results are given for the Bessell UBVRIJKHL, Sloan ugriz, Stromgren ubvy, WFCAM ZYJHK, Hipparcos, Kepler, and Tycho passbands, in each case characterized by several different limb-darkening `laws'. We examine the sensitivity of limb darkening to temperature, gravity, metallicity, microturbulent velocity, and wavelength, and make a comparison with LTE models. The dependence on metallicity is very weak, but limb darkening is a moderately strong function of log(g) in this temperature regime.

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

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.

Optical spectroscopic monitoring has been conducted of two O stars in the Small and one in the Large Magellanic Cloud, the spectral characteristics of which place them in the Of?p category, which has been established in the Galaxy to consist of oblique magnetic rotators. All of these Magellanic stars show systematic spectral variations typical of the Of?p class, further strengthening their magnetic candidacy to the point of virtual certainty. The spectral variations are related to photometric variations derived from OGLE data by Naze et al. (2015) in a parallel study, which yields rotational periods for two of them. Now circular spectropolarimetry is required to measure their fields, and ultraviolet spectroscopy to further characterize their low-metallicity, magnetically confined winds, in support of hydrodynamical analyses.

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.

Massive stars burn hydrogen through the CNO cycle during most of their evolution. When mixing is efficient, or when mass transfer in binary systems happens, chemically processed material is observed at the surface of O and B stars. ON stars show stronger lines of nitrogen than morphologically normal counterparts. Whether this corresponds to the presence of material processed through the CNO cycle or not is not known. Our goal is to answer this question. We perform a spectroscopic analysis of a sample of ON stars with atmosphere models. We determine the fundamental parameters as well as the He, C, N, and O surface abundances. We also measure the projected rotational velocities. We compare the properties of the ON stars to those of normal O stars. We show that ON stars are usually helium-rich. Their CNO surface abundances are fully consistent with predictions of nucleosynthesis. ON stars are more chemically evolved and rotate - on average - faster than normal O stars. Evolutionary models including rotation cannot account for the extreme enrichment observed among ON main sequence stars. Some ON stars are members of binary systems, but others are single stars as indicated by stable radial velocities. Hence, mass transfer is not a simple explanation for the observed chemical properties. We conclude that ON stars show extreme chemical enrichment at their surface, consistent with nucleosynthesis through the CNO cycle. Its origin is not clear at present.

We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 {\mu}m of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 {\mu}m transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1{\sigma}. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 {\sigma} significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) >30 s.

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.

We report a time-series analysis of the O4 I(n)fp star zeta Pup, based on optical photometry obtained with the SMEI instrument on the Coriolis satellite, 2003--2006. A single astrophysical signal is found, with P = (1.780938 \pm 0.000093) d and a mean semi-amplitude of (6.9 \pm 0.3) mmag. There is no evidence for persistent coherent signals with semi-amplitudes in excess of ca. 2~mmag on any of the timescales previously reported in the literature. In particular, there is no evidence for a signature of the proposed rotation period, ca. 5.1~days; zeta Pup is therefore probably not an oblique magnetic rotator. The 1.8-day signal varies in amplitude by a factor ca. 2 on timescales of 10--100d (and probably by more on longer timescales), and exhibits modest excursions in phase, but there is no evidence for systematic changes in period over the 1000-d span of our observations. Rotational modulation and stellar-wind variability appear to be unlikely candidates for the underlying mechanism; we suggest that the physical origin of the signal may be pulsation associated with low-l oscillatory convection modes.

Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$\mu$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$\mu$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year.

We report a time-series analysis of satellite photometry of the non-radially pulsating Oe star zeta Oph, principally using data from SMEI obtained 2003--2008, but augmented with MOST and WIRE results. Amplitudes of the strongest photometric signals, at 5.18, 2.96, and 2.67/d, each vary independently over the 6-year monitoring period (from ca. 30 to <2 mmag at 5.18/d), on timescales of hundreds of days. Signals at 7.19/d and 5.18/d have persisted (or recurred) for around two decades. Supplementary spectroscopic observations show an H-alpha emission episode in 2006; this coincided with small increases in amplitudes of the three strongest photometric signals.

[Aims] Accurate astrometry is required to reliably cross-match 20th-century photographic catalogues against 21st-century digital surveys. The present work provides modern-era identifications and astrometry for the 801 emission-line objects "of stellar appearance" in the Armagh survey (the largest of its nature to date). [Methods] Targets have been individually identified in digital images using the Armagh Atlas and, in most cases, unambiguously matched to entries in the UCAC astrometric catalogues. [Results] Astrometry with sub-arcsecond precision is now available for all the major photographic spectroscopic surveys of the LMC. The results are used to compile an annotated metacatalogue of 1675 individual, spectroscopically identified candidate H-alpha-emission stars, including detailed cross-matching between catalogues, and resolving many (though not all) identification ambiguities in individual primary sources.

We present for the first time phase-resolved UV spectroscopy of an Of?p star, namely, HD 191612. The observations were acquired with the Space Telescope Imaging Spectrograph (STIS) on-board the Hubble Space Telescope (HST). We report the variability observed in the main photospheric and wind features and compare the results with previous findings for the Of?p star HD 108. We show that UV line strengths, H(alpha), and longitudinal magnetic field, vary coherently according to the rotational period (P = 537.6d), providing additional support for the magnetic oblique rotator scenario. The stellar and wind parameters of HD 191612 are obtained based on NLTE expanding atmosphere models. The peculiar wind line profile variations revealed by the new STIS data - not reproduced by 1D atmosphere models - are addressed through non-spherical MHD simulations coupled with radiative transfer. The basic aspects of the UV variability observed are explained and the structure of the dynamical magnetosphere of HD 191612 is discussed.

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: Accurate astrometry is required to reliably cross-match 20th-century catalogues against 21st-century surveys. The present work aims to provide such astrometry for the 625 entries of the Bohannan & Epps (BE74) catalogue of H$\alpha$ emission-line stars. Methods: BE74 targets have been individually identified in digital images and, in most cases, unambiguously matched to entries in the UCAC4 astrometric catalogue. Results: Sub-arcsecond astrometry is now available for almost all BE74 stars. Several identification errors in the literature illustrate the perils of relying solely on positional coincidences using poorer-quality astrometry.

We report the discovery of two previously unknown WN3 stars in the Large Magellanic Cloud. Both are bright (15th magnitude), isolated, and located in regions covered in earlier surveys, although both are relatively weak-lined. We suggest that there may be $\mathcal{O}(10)$ remaining undiscovered WNE stars in the LMC.

The magnetic O-star HD191612 exhibits strongly variable, cyclic Balmer line emission on a 538-day period. We show here that its variable Halpha emission can be well reproduced by the rotational phase variation of synthetic spectra computed directly from full radiation magneto-hydrodynamical simulations of a magnetically confined wind. In slow rotators such as HD191612, wind material on closed magnetic field loops falls back to the star, but the transient suspension of material within the loops leads to a statistically overdense, low velocity region around the magnetic equator, causing the spectral variations. We contrast such "dynamical magnetospheres" (DMs) with the more steady-state "centrifugal magnetospheres" of stars with rapid rotation, and discuss the prospects of using this DM paradigm to explain periodic line emission from also other non-rapidly rotating magnetic massive stars.

We present the first ultraviolet spectrum of the peculiar, magnetic Of?p star HD 108 obtained in its spectroscopic low state. The new data, obtained with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, reveal significant changes compared to IUE spectra obtained in the high state: N V 1240, Si IV 1400, and C IV 1550 present weaker P-Cygni profiles (less absorption) in the new data, while N IV 1718 absorption is deeper, without the clear wind signature evident in the high state. Such changes contrast with those found in other magnetic massive stars, where more absorption is observed in the resonance doublets when the sightline is close to the plane of the magnetic equator. The new data show also that the photospheric Fe IV forest, at 1600--1700 angstroms, has strengthened compared to previous observations. The ultraviolet variability is large compared to that found in typical, non-magnetic O stars, but moderate when compared to the high-/low-state changes reported in the optical spectrum of HD 108 over several decades. We use non-LTE expanding-atmosphere models to analyze the new STIS observations. Overall, the results are in accord with a scenario in which the optical variability is mainly produced by magnetically constrained gas, close to the photosphere. The relatively modest changes found in the main ultraviolet wind lines suggest that the stellar wind is not substantially variable on a global scale. Nonetheless, multidimensional radiative-transfer models may be needed to understand some of the phenomena observed.

Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. The long-term datasets revealed that they are remarkably stable over long periods (>40yr), with the possible exception of zeta Sco prior to the 20th century, in contrast to the typical excursions that characterise luminous blue variables (LBVs). Zeta Sco, HD190603 and BP Cru possess physical properties intermediate between B supergiants and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars (for initial masses in the range ~30-60Msun). In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110Msun) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, the properties of Cyg OB2 #12 are difficult to understand under current evolutionary paradigms. [ABRIDGED]

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.

We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 within the context of the MiMeS LP at the CFHT. Thirty-two high signal-to-noise ratio circularly polarised (Stokes V) spectra and 13 unpolarised (Stokes I) spectra of HD 148937 were acquired in 2009 and 2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all observations (in both LSD mean profiles and individual spectral lines). The longitudinal magnetic field inferred from the Stokes V LSD profiles is consistently negative, in contrast to the essentially zero field strength measured from the diagnostic null profiles. A period search of equivalent width measurements confirms the previously-reported 7.03 d variability period. The variation of equivalent widths is not strictly periodic: we present evidence for evolution of the amount or distribution of circumstellar plasma. Interpreting the 7.03 d period as the stellar rotational period within the context of the ORM, we have phased the equivalent widths and longitudinal field measurements. The longitudinal field measurements show a weak sinusoidal variation of constant sign, with extrema out of phase with the H{\alpha} variation by about 0.25 cycles. The inferred magnetic configuration confirms the suggestion of Naz\'e et al (2010), who proposed that the weaker variability of HD 148937 as compared to other members of this class is a consequence of the stellar geometry. Based on the derived magnetic properties and published wind characteristics, we find a wind magnetic confinement parameter \eta\ast \simeq 20 and rotation parameter W = 0.12, supporting a picture in which the Halpha emission and other line variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. (Abridged.)

This paper examines how to compare stellar limb-darkening coefficients evaluated from model atmospheres with those derived from photometry. Different characterizations of a given model atmosphere can give quite different numerical results (even for a given limb-darkening `law'), while light-curve analyses yield limb-darkening coefficients that are dependent on system geometry, and that are not directly comparable to any model-atmosphere representation. These issues are examined in the context of exoplanetary transits, which offer significant advantages over traditional binary-star eclipsing systems in the study of stellar limb darkening. `Like for like' comparisons between light-curve analyses and new model-atmosphere results, mediated by synthetic photometry, are conducted for a small sample of stars. Agreement between the resulting synthetic-photometry/atmosphere-model (SPAM) limb-darkening coefficients and empirical values ranges from very good to quite poor, even though the targets investigated show only a small dispersion in fundamental stellar parameters.

This paper reports high-precision Stokes V spectra of HD 191612 acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, in the context of the Magnetism in Massive stars (MiMeS) Project. Using measurements of the equivalent width of the Halpha line and radial velocities of various metallic lines, we have updated both the spectroscopic and orbital ephemerides of this star. We confirm the presence of a strong magnetic field in the photosphere of HD 191612, and detect its variability. We establish that the longitudinal field varies in a manner consistent with the spectroscopic period of 537.6 d, in an approximately sinusoidal fashion. This demonstrates a firm connection between the magnetic field and the processes responsible for the line and continuum variability. Interpreting the variation of the longitudinal magnetic field within the context of the dipole oblique rotator model we obtain a best-fit surface magnetic field model with obliquity beta=67\pm 5 deg and polar strength Bd=2450\pm 400 G . The inferred magnetic field strength implies an equatorial wind magnetic confinement parameter eta*~50, supporting a picture in which the Halpha emission and photometric variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. This interpretation is supported by our successful Monte Carlo radiative transfer modeling of the photometric variation, which assumes the enhanced plasma densities in the magnetic equatorial plane above the star implied by such a picture. Predictions of the continuum linear polarisation resulting from Thompson scattering from the magnetospheric material indicate that the Stokes Q and U variations are highly sensitive to the magnetospheric geometry, and that expected amplitudes are in the range of current instrumentation. (abridged)

We present high precision transit observations of the exoplanet WASP-21b, obtained with the RISE instrument mounted on 2.0m Liverpool Telescope. A transit model is fitted, coupled with an MCMC routine to derive accurate system parameters. The two new high precision transits allow to estimate the stellar density directly from the light curve. Our analysis suggests that WASP-21 is evolving off the main sequence which led to a previous overestimation of the stellar density. Using isochrone interpolation, we find a stellar mass of 0.86 \pm 0.04 Msun which is significantly lower than previously reported (1.01 \pm 0.03 Msun). Consequently, we find a lower planetary mass of $0.27 \pm 0.01 Mjup$. A lower inclination (87.4 \pm 0.3 degrees) is also found for the system than previously reported, resulting in a slightly larger stellar (R_* =1.10 \pm 0.03 Rsun) and planetary radius (R_p = 1.14 \pm 0.04 Rjup). The planet radius suggests a hydrogen/helium composition with no core which strengthens the correlation between planetary density and host star metallicity. A new ephemeris is determined for the system, i.e., t0 =2455084.51974 \pm 0.00020 (HJD) and P=4.3225060 \pm 0.0000031 days. We found no transit timing variations in WASP-21b.

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.

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.

New grids of Atlas9 models have been calculated using revised convection parameters and updated opacity-distribution functions, for chemical compositions intended to be representative of solar, [M/H] = +0.3, +0.5, Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) abundances. The grids cover T(eff) = 3.5-50kK, from log(g) = 5.0 to the effective Eddington limit. Limb-darkening coefficients and synthetic photometry are presented in the UBVRIJHKLM, uvby, ugriz, WFCAM, Hipparcos/Tycho, and Kepler passbands for these models, and for Castelli's comparable `new-ODF' grids. Flux distributions are given for the new models. The sensitivity of limb-darkening coefficients to the adopted physics is illustrated.

Only 5 Of?p stars have been identified in the Galaxy. Of these, 3 have been studied in detail, and within the past 5 years magnetic fields have been detected in each of them. The observed magnetic and spectral characteristics are indicative of organised magnetic fields, likely of fossil origin, confining their supersonic stellar winds into dense, structured magnetospheres. The systematic detection of magnetic fields in these stars strongly suggests that the Of?p stars represent a general class of magnetic O-type stars.

We report the detection of a magnetic field on the Of?p star HD108. Spectropolarimetric observations conducted in 2007, 2008 and 2009 respectively with NARVAL@TBL and ESPaDOnS@CFHT reveal a clear Zeeman signature in the average Stokes V profile, stable on timescales of days to months and slowly increasing in amplitude on timescales of years. We speculate that this timescale is the same as that on which Ha emission is varying and is equal to the rotation period of the star. The corresponding longitudinal magnetic field, measured during each of the three seasons, increases slowly from 100 to 150G, implying that the polar strength of the putatively-dipolar large-scale magnetic field of HD108 is at least 0.5kG and most likely of the order of 1-2 kG. The stellar and wind properties are derived through a quantitative spectroscopic analysis with the code CMFGEN. The effective temperature is difficult to constrain because of the unusually strong HeI4471 and HeI5876 lines. Values in the range 33000-37000 K are preferred. A mass loss rate of about 1e-7 Msun/yr (with a clumping factor f=0.01) and a wind terminal velocity of 2000 km/s are derived. The wind confinement parameter eta_star is larger than 100, implying that the wind of HD108 is magnetically confined. Stochastic short-term variability is observed in the wind-sensitive lines but not in the photospheric lines, excluding the presence of pulsations. Material infall in the confined wind is the most likely origin for lines formed in the inner wind. Wind-clumping also probably causes part of the Ha variability. The projected rotational velocity of HD108 is lower than 50 km/s, consistent with the spectroscopic and photometric variation timescales of a few decades. Overall, HD108 is very similar to the magnetic O star HD191612 except for an even slower rotation.

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.

The Onfp class of rotationally broadened, hot spectra was defined some time ago in the Galaxy, where its membership to date numbers only eight. The principal defining characteristic is a broad, centrally reversed He II $\lambda$4686 emission profile; other emission and absorption lines are also rotationally broadened. Recent surveys in the Magellanic Clouds (MCs) have brought the class membership there, including some related spectra, to 28. We present a survey of the spectral morphology and rotational velocities, as a first step toward elucidating the nature of this class. Evolved, rapidly rotating hot stars are not expected theoretically, because the stellar winds should brake the rotation. Luminosity classification of these spectra is not possible, because the principal criterion (He II $\lambda$4686) is peculiar; however, the MCs provide reliable absolute magnitudes, which show that they span the entire range from dwarfs to supergiants. The Onfp line-broadening distribution is distinct and shifted toward larger values from those of normal O dwarfs and supergiants with >99.99% confidence. All cases with multiple observations show line-profile variations, which even remove some objects from the class temporarily. Some of them are spectroscopic binaries; it is possible that the peculiar profiles may have multiple causes among different objects. The origin and future of these stars are intriguing; for instance, they could be stellar mergers and/or gamma-ray-burst progenitors.

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 follow the bright, highly energetic afterglow of Swift-discovered GRB 080721 at z=2.591 out to 36 days or 3e6 s since the trigger in the optical and X-ray bands. We do not detect a break in the late-time light curve inferring a limit on the opening angle of theta_j >= 7.3 deg and setting tight constraints on the total energy budget of the burst of E_gamma >= 9.9e51 erg within the fireball model. To obey the fireball model closure relations the GRB jet must be expanding into a homogeneous surrounding medium and likely lies behind a significant column of dust. The energy constraint we derive can be used as observational input for models of the progenitors of long gamma-ray bursts: we discuss how such high collimation-corrected energies could be accommodated with certain parameters of the standard massive star core-collapse models. We can, however, most probably rule out a magnetar progenitor for this GRB which would require 100% efficiency to reach the observed total energy.

We have previously analysed the spectra of 135 early B-type stars in the LMC and found several groups of stars that have chemical compositions that conflict with the theory of rotational mixing. Here we extend this study to Galactic and SMC metallicities with the analysis of ~50 Galactic and ~100 SMC early B-type stars with rotational velocities up to ~300km/s. The surface nitrogen abundances are utilised as a probe of the mixing process. In the SMC, we find a population of slowly rotating nitrogen-rich stars amongst the early B type core-hydrogen burning stars, similar to the LMC. In the Galactic sample we find no significant enrichment amongst the core hydrogen-burning stars, which appears to be in contrast with the expectation from both rotating single-star and close binary evolution models. However, only a small number of the rapidly rotating stars have evolved enough to produce a significant nitrogen enrichment, and these may be analogous to the non-enriched rapid rotators previously found in the LMC sample. Finally, in each metallicity regime, a population of highly enriched supergiants is observed, which cannot be the immediate descendants of core-hydrogen burning stars. Their abundances are, however, compatible with them having gone through a previous red supergiant phase. Together, these observations paint a complex picture of the nitrogen enrichment in massive main sequence and supergiant stellar atmospheres, where age and binarity cause crucial effects. Whether rotational mixing is required to understand our results remains an open question at this time, but could be answered by identifying the true binary fraction in those groups of stars that do not agree with single-star evolutionary models (abridged).

The paradigmatic Luminous Blue Variable R127 in the Large Magellanic Cloud has been found in the intermediate, peculiar early-B state, and substantially fainter in visual light, signaling the final decline from its major outburst that began between 1978 and 1980. This transformation was detected in 2008 January, but archival data show that it began between early 2005 and early 2007. In fact, significant changes from the maximum, peculiar A-type spectrum, which was maintained from 1986 through 1998, had already begun the following year, coinciding with a steep drop in visual light. We show detailed correspondences between the spectrum and light, in which the decline mimics the rise. Moreover, these trends are not monotonic but are characterized by multiple spikes and dips, which may provide constraints on the unknown outburst mechanism. Intensive photometric and spectroscopic monitoring of R127 should now resume, to follow the decline presumably back to the quiescent Ofpe/WN9 state, in order to fully document the remainder of this unique observational opportunity.

We report here the detection of a weak magnetic field of 50 - 100 G on the O9.7 supergiant zeta Ori A, using spectropolarimetric observations obtained with NARVAL at the 2m Telescope Bernard Lyot atop Pic du Midi (France). zeta Ori A is the third O star known to host a magnetic field (along with theta^1 Ori C and HD 191612), and the first detection on a 'normal' rapidly-rotating O star. The magnetic field of zeta Ori A is the weakest magnetic field ever detected on a massive star. The measured field is lower than the thermal equipartition limit (about 100 G). By fitting NLTE model atmospheres to our spectra, we determined that zeta Ori A is a 40 Msun star with a radius of 25 Rsun and an age of about 5 - 6 Myr, showing no surface nitrogen enhancement and losing mass at a rate of about 2x10^(-6) Msol/yr. The magnetic topology of zeta Ori A is apparently more complex than a dipole and involves two main magnetic polarities located on both sides of the same hemisphere; our data also suggest that zeta Ori A rotates in about 7.0 d and is about 40 degrees away from pole-on to an Earth-based observer. Despite its weakness, the detected magnetic field significantly affects the wind structure; the corresponding Alfven radius is however very close to the surface, thus generating a different rotational modulation in wind lines than that reported on the two other known magnetic O stars. The rapid rotation of zeta Ori A with respect to theta^1 Ori C appears as a surprise, both stars having similar unsigned magnetic fluxes (once rescaled to the same radius); it may suggest that the sub-equipartition field detected on zeta Ori A is not a fossil remnant (as opposed to that of theta^1 Ori C and HD 191612), but the result of an exotic dynamo action produced through MHD instabilities.

We present radial velocities for 2045 stars in the Small Magellanic Cloud (SMC), obtained from the 2dF survey by Evans et al. (2004). The great majority of these stars are of OBA type, tracing the dynamics of the young stellar population. Dividing the sample into ad hoc `bar' and `wing' samples (north and south, respectively, of the line: $\delta$ = -77$^{\circ}$50' + [4$\alpha$]', where $\alpha$ is in minutes of time) we find that the velocities in the SMC bar show a gradient of 26.3 +/- 1.6 km/s/deg. at a position angle of 126 +/- 4 deg. The derived gradient in the bar is robust to the adopted line of demarcation between the two samples. The largest redshifts are found in the SMC wing, in which the velocity distribution appears distinct from that in the bar, most likely a consequence of the interaction between the Magellanic Clouds that is predicted to have occurred 0.2 Gyr ago. The mean velocity for all stars in the sample is +172.0 +/- 0.2 km/s (redshifted by ~20 km/s when compared to published results for older populations), with a velocity dispersion of 30 km/s.

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe-O8fp). The Balmer and HeI lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. Metal lines and HeII absorptions (including many selective emission lines but excluding He II 4686A emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with P(orb) = 1542d, e=0.45. We conduct a model-atmosphere analysis of the primary, and find that the system is consistent with a O8: giant with a B1: main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying `clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.

We present the analysis of phase-resolved X-ray and optical observations of the peculiar hot star HD191612 (Of?p). This star is known to display line-profile variations that are recurrent with a period of 538d and its spectrum was found to present the signature of a magnetic field. In the X-rays, it is slightly overluminous compared to the canonical L_X/L_BOL relation and appears brighter when the optical lines are strongest. Our XMM-Newton observations further reveal that the X-ray spectrum of \hd exhibits rather broad lines and is dominated by a `cool' (0.2-0.6 keV) thermal component, two characteristics at odds with the proposed magnetic rotator model. We also report for the first time the low-level variability of the metallic (abs/em) lines and HeII absorptions that appear to be associated with radial-velocity shifts. Finally, we compare our results with observations of the early-type stars and discuss several possible scenarios.

We report the discovery of a medium-strength (~0.5kG) magnetic field on the young, massive star tauSco (B0.2V), which becomes the third-hottest magnetic star known. Circularly polarized Zeeman signatures are clearly detected in observations collected mostly with the ESPaDOnS spectropolarimeter, recently installed on the 3.6-m Canada-France-Hawaii Telescope; temporal variability is also clearly established in the polarimetry, and can be unambiguously attributed to rotational modulation with a period close to 41d. Archival UV spectra confirm that this modulation repeats over timescales of decades. By reconstructing the large-scale structure of its magnetic topology, we find that the magnetic structure is unusually complex for a hot star. The surface topology is dominated by a potential field, although a moderate toroidal component is probably present. We fail to detect intrinsic temporal variability of the magnetic structure over the 1.5-yr period of our spectropolarimetric observations (in agreement with the stable temporal variations of the UV spectra), and infer that any differential surface rotation must be very small. The topology of the extended magnetic field that we derive from the photospheric magnetic maps is also more complex than a global dipole, and features in particular a significantly warped torus of closed magnetic loops encircling the star (tilted at about 90deg to the rotation axis), with additional, smaller, networks of closed field lines. This topology appears to be consistent with the exceptional Xray properties of tauSco and also provides a natural explanation of the variability observed in wind-formed UV lines. We conclude that its magnetic field is most probably a fossil remnant from the star-formation stage.

Recent results strongly challenge the canonical picture of massive star winds: various evidence indicates that currently accepted mass-loss rates, Mdot, may need to be revised downwards significantly. This is because the most commonly used mass-loss diagnostics are affected by ``clumping'' (small-scale density inhomogeneities), influencing our interpretation of observed spectra and fluxes. Such downward revisions would have dramatic consequences for the evolution of, and feedback from, massive stars, and thus robust determinations of the clumping properties and mass-loss rates are urgently needed. Here, we present a first attempt to constrain the radial stratification of the so-called clumping factor. To this end, we have analyzed a sample of 19 Galactic O-type supergiants/giants, by combining data for Halpha, IR, mm and radio fluxes, and using appropriate analysis methods. Clumping has been included into our analysis in the ``conventional'' way, by assuming the inter-clump matter to be void. Because (almost) all our diagnostics depends on the square of density, we cannot derive absolute clumping factors, but only factors normalized to a certain minimum. This minimum was usually found to be located in the outermost, radio-emitting region, i.e., the radio mass-loss rates are the lowest ones, compared to Mdot derived from Halpha and the IR. The radio rates agree well with those predicted by theory, but are only upper limits, due to unknown clumping in the outer wind. Our most important result concerns a (physical) difference between denser and thinner winds: for denser winds, the innermost region is more strongly clumped than the outermost one (with a normalized clumping factor of 4.1+/-1.4), whereas thinner winds have similar clumping properties in the inner and outer regions.

From observations made with the ESPaDOnS spectropolarimeter, recently installed on the 3.6-m Canada--France--Hawaii Telescope, we report the discovery of a strong magnetic field in the Of?p spectrum variable HD 191612 -- only the second known magnetic O star (following theta1 Ori C). The stability of the observed Zeeman signature over four nights of observation, together with the non-rotational shape of line profiles, argue that the rotation period of HD 191612 is significantly longer than the 9-d value previously proposed. We suggest that the recently identified 538-d spectral-variability period is the rotation period, in which case the observed line-of-sight magnetic field of -220+-38 G implies a large-scale field (assumed dipolar) with a polar strength of about -1.5 kG. If confirmed, this scenario suggests that HD 191612 is, essentially, an evolved version of the near-ZAMS magnetic O star theta1 Ori C, but with an even stronger field (about 15 kG at an age similar to that of theta1Ori C). We suggest that the rotation rate of HD 191612, which is exceptionally slow by accepted O-star standards, could be due to angular-momentum dissipation through a magnetically confined wind.

We have conducted a programme to determine the fundamental parameters of a substantial number of eclipsing binaries of spectral types O and B in the Small Magellanic Cloud. New spectroscopic data, obtained with the two-degree-field multi-object spectrograph on the 3.9-m Anglo-Australian Telescope, have been used in conjunction with photometry from the Optical Gravitational Lens Experiment (OGLE-II) database of SMC eclipsing binaries. Previously we reported results for 10 systems; in this second and concluding paper we present spectral types, masses, radii, temperatures, surface gravities and luminosities for the components of a further 40 binaries. The full sample of 50 OB-type eclipsing systems is the largest single set of fundamental parameters determined for high-mass binaries in any galaxy. We find that 21 of the systems studied are in detached configurations, 28 are in semi-detached post-mass-transfer states, and one is a contact binary. Each system provides a primary distance indicator. We find a mean distance modulus to the SMC of 18.91+/-0.03+/-0.1 (internal and external uncertainties; D=60.6+/-1.0 kpc). This value represents one of the most precise available determinations of the distance to the SMC.

The observational picture of the enigmatic O-type spectrum variable HD191612 has been sharpened substantially. A symmetrical, low-amplitude light curve with a period near 540 d has recently been reported from Hipparcos photometry. This period satisfies all of the spectroscopy since at least 1982, including extensive new observations during 2003 and 2004, and it has predicted the next transition during September--October 2004. Measurements of the H alpha equivalent width reveal a sharp emission peak in the phase diagram, in contrast to the apparently sinusoidal light curve. The He II absorption-line strength is essentially constant, while He I varies strongly, possibly filled in by emission in the O6 state, thus producing the apparent spectral-type variations. The O8 state appears to be the "normal" one. Two intermediate O7 observations have been obtained, which fall at the expected phases, but these are the only modern observations of the transitions so far. The period is too long for rotation or pulsation; although there is no direct evidence as yet for a companion, a model in which tidally induced oscillations drive an enhanced wind near periastron of an eccentric orbit appears promising. Further observations during the now predictable transitions may provide a critical test. Ultraviolet and X-ray observations during both states will likely also prove illuminating.

New constraints on the large-scale wind structures responsible for discrete absorption components are obtained from far-ultraviolet time-series observations of O supergiants in the Large Magellanic Cloud.

We present a catalogue of new spectral types for hot, luminous stars in the Small Magellanic Cloud. The catalogue contains 4161 objects, giving an order of magnitude increase in the number of SMC stars with published spectroscopic classifications. The targets are primarily B- and A-type stars (2862 and 853 objects respectively), with 1 Wolf-Rayet, 139 O-type, and 306 FG stars, sampling the main sequence to ~mid-B. The selection and classification criteria are described, and objects of particular interest are discussed, including UV-selected targets from the Ultraviolet Imaging Telescope (UIT) experiment, Be and B[e] stars, `anomalous A supergiants', and composite-spectrum systems. We examine the incidence of Balmer-line emission, and the relationship between H-gamma equivalent width and absolute magnitude for BA stars.

From a survey of the 3400 A region in the earliest O-type spectra, we have found that two of the four O2 giants observed in the LMC have the OIV lines there stronger than the NIV, while the other two have the opposite. A SMC counterpart also has the NIV stronger than OIV. Inspection of the blue spectra of these stars shows that the former pair have weaker N lines in all ionization states (III, IV, V) present as well as lines of CIV lambda 4658, while the latter three have stronger N lines and greater He/H. Space ultraviolet observations of two of the N-strong stars show NV wind profiles substantially stronger than those of CIV, while in the N-weak stars the CIV features are equal to or stronger than the NV. The N-strong stars are now reclassified as ON2 III(f*), newly defining that category. These characteristics strongly suggest a larger fraction of processed material in the atmospheres of the ON2 stars, which we confirm by modeling the optical spectra. In the context of current models, it is in turn implied that the ON2 stars are in a more advanced evolutionary state than the others, and/or that they had higher initial rotational velocities. The recent formulation of the effects of rotation on massive stellar evolution introduces an additional fundamental parameter, which the CNO abundances are in principle able to constrain. We present some illustrative comparisons with current Geneva evolutionary models for rotating massive stars. It is possible that these very hot, nitrogen-rich objects are products of homogeneous evolution. Our results will provide motivation for further physical modeling of the atmospheres and evolutionary histories of the most massive hot stars.

We argue that, in general, observational studies of Be-star rotation have paid insufficient attention to the effects of equatorial gravity darkening. We present new line-profile calculations that emphasize the insensitivity of line width to rotation for fast rotators. Coupled with a critical review of observational procedures, these calculations suggest that the observational parameter v sin(i) may systematically underestimate the true projected equatorial rotation velocity, ve sin(i), by some tens of per cent for rapid rotators. The crucial implication of this work is that Be stars may be rotating much closer to their critical velocities than is generally supposed, bringing a range of new processes into contention for the elusive physical mechanism responsible for the circumstellar disk thought to be central to the Be phenomenon.

We address the relationship between spectral type and physical properties for A-type supergiants in the SMC. We first construct a self-consistent classification scheme for A supergiants, employing the calcium K to H epsilon line ratio as a temperature-sequence discriminant. Following the precepts of the `MK process', the same morphological criteria are applied to Galactic and SMC spectra with the understanding there may not be a correspondence in physical properties between spectral counterparts in different environments. We then discuss the temperature scale, concluding that A supergiants in the SMC are systematically cooler than their Galactic counterparts at the same spectral type, by up to ~10%. Considering the relative line strengths of H gamma and the CH G-band we extend our study to F and early G-type supergiants, for which similar effects are found. We note the implications for analyses of extragalactic luminous supergiants, for the flux-weighted gravity-luminosity relationship and for population synthesis studies in unresolved stellar systems.

We present the first results of an observational programme to measure the fundamental parameters of over 100 eclipsing binaries in the Small Magellanic Cloud (SMC). The spectroscopic data were obtained by using the 2dF multi-object spectrograph on the 3.9-m Anglo-Australian Telescope, and are used in conjunction with photometry from the OGLE database of SMC eclipsing binaries. Ten systems are discussed in this first paper. Three are detached early-B binaries, six are in a semi-detached configuration, and one is in a marginal contact state. We conclude that the semi-detached systems are undergoing the slow mass-transfer phase of case-A binary evolution, in which the mass donor has reached its Roche lobe while still on the main sequence. Each system provides a primary distance indicator. By constructing a new calibration between spectral type and temperature for O and early B stars, we find a mean distance modulus to the SMC of 18.89 +- 0.04 (statistical) +- 0.10 (systematic). This value represents one of the most precise determinations to date of the distance to the SMC.

A stellar wind module has been developed for the PHOENIX stellar atmosphere code for the purpose of computing non-LTE, line-blanketed, expanding atmospheric structures and detailed synthetic spectra of hot luminous stars with winds. We apply the code to observations of Deneb, for which we report the first positive detections of mm and cm emission (obtained using the SCUBA and the VLA), as well a strong upper limit on the 850 micron flux (using the HHT). The slope of the radio spectrum shows that the stellar wind is partially ionized. We report a uniform-disk angular diameter measurement, 2.40 +/- 0.06 mas, from the Navy Prototype Optical Interferometer (NPOI). The measured bolometric flux and corrected NPOI angular diameter yield an effective temperature of 8600 +/- 500 K. Least-squares comparisons of synthetic spectral energy distributions from 1220 A to 3.6 cm with the observations provide estimates for the effective temperature and the mass-loss rate of 8400 +/- 100 K and 8 +/- 3 E-7 M_sun/yr, respectively. This range of mass-loss rates is consistent with that derived from high dispersion UV spectra when non-LTE metal-line blanketing is considered. We are unable achieve a reasonable fit to a typical Halpha P-Cygni profile with any model parameters over a reasonable range. This is troubling because the \ha profile is the observational basis for Wind Momentum-Luminosity Relationship.