New long Chandra grating observations of the O supergiant $\zeta$ Pup show not only a brightening in the x-ray emission line flux of a little more than 25 per cent in the 18 years since Chandra's first observing cycle, but also clear evidence of increased wind absorption signatures in its Doppler-broadened x-ray emission line profiles. Model fitting of ten lines in three day-long HETGS datasets taken during August 2018 (Chandra cycle 19) finds a mass-loss rate of $2.68 \pm 0.12 \times 10^{-6}$ Msun/yr which is a 50 per cent increase over the value obtained from fitting the cycle 1 data. The increase in the individual emission line fluxes is greater for short-wavelength lines than long-wavelength lines, as would be expected if a uniform increase in line emission is accompanied by an increase in the wavelength-dependent absorption by the cold wind in which the shock-heated plasma is embedded.

In this paper, we present results from the spectropolarimetric follow-up of photometrically selected candidate magnetic B stars from the MOBSTER project. Out of four observed targets, one (HD 38170) is found to host a detectable surface magnetic field, with a maximum longitudinal field measurement of 105$\pm$14 G. This star is chemically peculiar and classified as an $\alpha^2$ CVn variable. Its detection validates the use of TESS to perform a photometric selection of magnetic candidates. Furthermore, upper limits are derived for the remaining three stars, and we report the discovery of a previously unknown spectroscopic binary system, HD 25709. Finally, we use our non-detections as case studies to further inform the criteria to be used for the selection of a larger sample of stars to be followed up using high-resolution spectropolarimetry.

We review recent developments in determining the evolution of magnetic massive OB stars. One of the important scientific questions is the completeness and the detection limits of contemporaneous spectropolarimetric surveys across the HR diagram. We present the characteristics of the MiMeS survey of O-type stars, and how the limits of the current available observations warrant the design of new high-precision surveys that target the older O-type star population. Another important question is whether the presence of the magnetic fields changes stellar evolutionary tracks in a significant way, hence leading to a wrong determination of stellar parameters. We review new evolution models that include the effect of magnetic wind quenching, and suppression of convection in the iron opacity peak zone.

Observational and theoretical work has now established that the fossil fields of magnetic massive stars are surviving remnants from an earlier event, or an earlier evolutionary phase. However, many important questions remain regarding the effects of these fields on the late stages of stellar evolution, as well as their impact on the core-collapse mechanism and the formation of exotic compact objects such as magnetars and gravitational wave progenitors. There is currently a critical need to incorporate the impact of fossil fields in models of the structure and evolution of magnetic stars, and to determine the evolutionary history of magnetic massive stars. We present a preliminary population study of a cluster of co-evolving stars based on MESA evolutionary tracks that account for the effect of magnetic mass-loss quenching.

The time evolution of angular momentum and surface rotation of massive stars is strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields ($M_\star=5-60$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2-1.0$, $B_{\rm p} =1-20$ kG). We then employ a representative grid of B-type star models ($M_\star=5, 10, 15$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2 , 0.5, 0.8$, $B_{\rm p} = 1, 3 ,10, 30$ kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero age main sequence with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly-rotating stars.

Recent spectropolarimetric surveys (MiMeS, BOB) have revealed that approximately 7% of massive stars host stable, surface dipolar magnetic fields with strengths on the order of kG. These fields channel the dense radiatively driven stellar wind into a circumstellar magnetosphere. Wind-sensitive UV spectral lines can probe the density and velocity structure of massive star magnetospheres, providing insight into wind-field interactions. To date, large-scale magnetohydrodynamic modeling of this phenomenon has been limited by the associated computational cost. Our analysis, using the Analytic Dynamical Magnetosphere model, solves this problem by applying a simple analytic prescription to efficiently calculate synthetic UV spectral lines. It can therefore be applied in the context of a larger parameter study to derive the wind properties for the population of known magnetic O stars. We also present the latest UV spectra of the magnetic O star NGC 1624-2 obtained with HST/COS, which test the limits of our models and suggest a particularly complex magnetospheric structure for this archetypal object.

Magnetic massive and intermediate-mass stars constitute a separate population whose properties are still not fully understood. Increasing the sample of known objects of this type would help answer fundamental questions regarding the origins and characteristics of their magnetic fields. The MOBSTER Collaboration seeks to identify candidate magnetic A, B and O stars and explore the incidence and origins of photometric rotational modulation using high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. In this contribution, we present an overview of our methods and planned targeted spectropolarimetric follow-up surveys.

In this contribution, we present the MOBSTER Collaboration, a large community effort to leverage high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) in order to characterize the variability of magnetic massive and intermediate-mass stars. These data can be used to probe the varying column density of magnetospheric plasma along the line of sight for OB stars, thus improving our understanding of the interaction between surface magnetic fields and massive star winds. They can also be used to map out the brightness inhomogeneities present on the surfaces of Ap/Bp stars, informing present models of atomic diffusion in their atmospheres. Finally, we review our current and ongoing studies, which lead to new insights on this topic.

It was recently claimed that the magnetic O-type star HD 54879 exhibits important radial velocity variability indicative of its presence in a spectroscopic binary. More remarkably, it was furthermore reported that the star underwent a short, sudden variation in spectral type and magnetic field. In this Letter we examine new Narval and ESPaDOnS data of this star in addition to the previously-published FORS2 data and conclude that both the reported velocity variations and the sudden spectral and magnetic changes are spurious.

Our multi-view metric learning framework enables robust characterization of star categories by directly learning to discriminate in a multi-faceted feature space, thus, eliminating the need to combine feature representations prior to fitting the machine learning model. We also demonstrate how to extend standard multi-view learning, which employs multiple vectorized views, to the matrix-variate case which allows very novel variable star signature representations. The performance of our proposed methods is evaluated on the UCR Starlight and LINEAR datasets. Both the vector and matrix-variate versions of our multi-view learning framework perform favorably --- demonstrating the ability to discriminate variable star categories.

We present the construction of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern detection algorithm. We focus on the targeted identification of eclipsing binaries that demonstrate a feature known as the O'Connell effect. Our proposed methodology maps stellar variable observations to a new representation known as distribution fields (DFs). Given this novel representation, we develop a metric learning technique directly on the DF space that is capable of specifically identifying our stars of interest. The metric is tuned on a set of labeled eclipsing binary data from the Kepler survey, targeting particular systems exhibiting the O'Connell effect. The result is a conservative selection of 124 potential targets of interest out of the Villanova Eclipsing Binary Catalog. Our framework demonstrates favorable performance on Kepler eclipsing binary data, taking a crucial step in preparing the way for large-scale data volumes from next-generation telescopes such as LSST and SKA.

HD 62658 (B9p V) is a little-studied chemically peculiar star. Light curves obtained by the Kilodegree Extremely Little Telescope (KELT) and Transiting Exoplanet Survey Satellite (TESS) show clear eclipses with a period of about 4.75 d, as well as out-of-eclipse brightness modulation with the same 4.75 d period, consistent with synchronized rotational modulation of surface chemical spots. High-resolution ESPaDOnS circular spectropolarimetry shows a clear Zeeman signature in the line profile of the primary; there is no indication of a magnetic field in the secondary. PHOEBE modelling of the light curve and radial velocities indicates that the two components have almost identical masses of about 3 M$_\odot$. The primary's longitudinal magnetic field $\langle B_z \rangle$ varies between about $+100$ and $-250$ G, suggesting a surface magnetic dipole strength $B_{\rm d} = 850$~G. Bayesian analysis of the Stokes $V$ profiles indicates $B_{\rm d} = 650$~G for the primary and $B_{\rm d} < 110$ G for the secondary. The primary's line profiles are highly variable, consistent with the hypothesis that the out-of-eclipse brightness modulation is a consequence of rotational modulation of that star's chemical spots. We also detect a residual signal in the light curve after removal of the orbital and rotational modulations, which might be pulsational in origin; this could be consistent with the weak line profile variability of the secondary. This system represents an excellent opportunity to examine the consequences of magnetic fields for stellar structure via comparison of two stars that are essentially identical with the exception that one is magnetic. The existence of such a system furthermore suggests that purely environmental explanations for the origin of fossil magnetic fields are incomplete.

Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$\alpha$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed study of the magnetic and rotational evolution of this population has yet been performed. Using newly determined rotational periods, modern magnetic measurements, and atmospheric parameters determined via spectroscopic modelling, we have derived fundamental parameters, dipolar oblique rotator models, and magnetospheric parameters for 56 early B-type stars. Comparison to magnetic A- and O-type stars shows that the range of surface magnetic field strength is essentially constant with stellar mass, but that the unsigned surface magnetic flux increases with mass. Both the surface magnetic dipole strength and the total magnetic flux decrease with stellar age, with the rate of flux decay apparently increasing with stellar mass. We find tentative evidence that multipolar magnetic fields may decay more rapidly than dipoles. Rotational periods increase with stellar age, as expected for a magnetic braking scenario. Without exception, all stars with H$\alpha$ emission originating in a CM are 1) rapid rotators, 2) strongly magnetic, and 3) young, with the latter property consistent with the observation that magnetic fields and rotation both decrease over time.

In this paper, we describe an analysis of the MiMeS Survey of O-type stars to explore the range of dipolar field strengths permitted by the polarisation spectra that do not yield a magnetic detection. We directly model the Stokes V profiles with a dipolar topology model using Bayesian inference. The noise statistics of the Stokes V profiles are in excellent agreement with those of the null profiles. Using a Monte-Carlo approach we conclude that a model in which all the stars in our sample were to host a 100\,G, dipolar magnetic field can be ruled out by the MiMeS data. Furthermore, if all the stars with no detection were to host a magnetic field just below their detection limit, the inferred distribution in strength of these undetected fields would be distinct from the known distribution in strength of the known magnetic O-type stars. This indicates that the Initial magnetic (B-)field Function (IBF) is likely bimodal -- young O-type stars are expected to either have weak/absent magnetic fields, or strong magnetic fields. We also find that better upper limits, by at least a factor of 10, would have been necessary to rule out a detection bias as an explanation for the apparent lack of evolved main-sequence magnetic O-type stars reported in the literature, and we conclude that the MiMeS survey cannot confirm or refute a magnetic flux decay in O-type stars.

It has been suggested that the absence of macroturbulence in the atmosphere of NGC 1624 - 2 is due its strong magnetic field (the strongest known for a massive O star) suppressing convection in its outer layers, removing the mechanism thought responsible for the observed macroturbulence in stars with lower field strengths. Here, we develop and apply a criterion for a uniform magnetic field to suppress convection in stellar envelopes in which radiation pressure is a significant contributor to hydrostatic balance. We find upper mass limits of ~55 Msun and ~30 Msun for magnetic suppression to be possible in Zero Age Main Sequence and Terminal Age Main Sequence stars, respectively. For evolved stars, magnetic suppression of convection can significantly alter the stars' evolution. For NGC 1624 - 2, we find a polar dipole strength of 16.5 +/- 5.9 kG is required to suppress convection, in good agreement with the value ~20 kG measured by spectropolarimetry.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

Surface magnetic fields have a strong impact on stellar mass loss and rotation and, as a consequence, on the evolution of massive stars. In this work we study the influence of an evolving dipolar surface fossil magnetic field with an initial field strength of 4 kG on the characteristics of 15 M$_{\odot}$ solar metallicity models using the Geneva stellar evolution code. Non-rotating and rotating models considering two different scenarios for internal angular momentum transport are computed, including magnetic field evolution, mass-loss quenching, and magnetic braking. Magnetic field evolution results in weakening the initially strong magnetic field, however, in our models an observable magnetic field is still maintained as the star evolves towards the red supergiant phase. At the given initial mass of the models, mass-loss quenching is modest. Magnetic braking greatly enhances chemical element mixing if radial differential rotation is allowed for, on the other hand, the inclusion of surface magnetic fields yields a lower surface enrichment in the case of near solid-body rotation. Models including surface magnetic fields show notably different trends on the Hunter diagram (plotting nitrogen abundance vs $v \sin i$) compared to those that do not. The magnetic models agree qualitatively with the anomalous `Group 2 stars', showing slow surface rotation and high surface nitrogen enhancement on the main sequence.

Atmospheric parameters determined via spectral modelling are unavailable for many of the known magnetic early B-type stars. We utilized high-resolution spectra together with NLTE models to measure effective temperatures $T_{\rm eff}$ and surface gravities $\log{g}$ of stars for which these measurements are not yet available. We find good agreement between our $T_{\rm eff}$ measurements and previous results obtained both photometrically and spectroscopically. For $\log{g}$, our results are compatible with previous spectroscopic measurements; however, surface gravities of stars previously determined photometrically have been substantially revised. We furthermore find that $\log{g}$ measurements obtained with HARPSpol are typically about 0.1 dex lower than those from comparable instruments. Luminosities were determined using Gaia Data Release 2 parallaxes. We find Gaia parallaxes to be unreliable for bright stars ($V<6$ mag) and for binaries; in these cases we reverted to Hipparcos parallaxes. In general we find luminosities systematically lower than those previously reported. Comparison of $\log{g}$ and $\log{L}$ to available rotational and magnetic measurements shows no correlation between either parameter with magnetic data, but a clear slow-down in rotation with both decreasing $\log{g}$ and increasing $\log{L}$, a result compatible with the expectation that magnetic braking should lead to rapid magnetic spindown that accelerates with increasing mass-loss.

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars.

NU Ori is a massive spectroscopic and visual binary in the Orion Nebula Cluster, with 4 components: Aa, Ab, B, and C. The B0.5 primary (Aa) is one of the most massive B-type stars reported to host a magnetic field. We report the detection of a spectroscopic contribution from the C component in high-resolution ESPaDOnS spectra, which is also detected in a Very Large Telescope Interferometer (VLTI) dataset. Radial velocity (RV) measurements of the inner binary (designated Aab) yield an orbital period of 14.3027(7) d. The orbit of the third component (designated C) was constrained using both RVs and interferometry. We find C to be on a mildly eccentric 476(1) d orbit. Thanks to spectral disentangling of mean line profiles obtained via least-squares deconvolution we show that the Zeeman Stokes $V$ signature is clearly associated with C, rather than Aa as previously assumed. The physical parameters of the stars were constrained using both orbital and evolutionary models, yielding $M_{\rm Aa} = 14.9 \pm 0.5 M_\odot$, $M_{\rm Ab} = 3.9 \pm 0.7 M_\odot$, and $M_{\rm C} = 7.8 \pm 0.7 M_\odot$. The rotational period obtained from longitudinal magnetic field $\langle B_z \rangle$ measurements is $P_{\rm rot} = 1.09468(7)$ d, consistent with previous results. Modeling of $\langle B_z \rangle$ indicates a surface dipole magnetic field strength of $\sim 8$ kG. NU Ori C has a magnetic field strength, rotational velocity, and luminosity similar to many other stars exhibiting magnetospheric H$\alpha$ emission, and we find marginal evidence of emission at the expected level ($\sim$1% of the continuum).

One of the defining processes which govern massive star evolution is their continuous mass loss via dense, supersonic line-driven winds. In the case of those OB stars which also host a surface magnetic field, the interaction between that field and the ionized outflow leads to complex circumstellar structures known as magnetospheres. In this contribution, we review recent developments in the field of massive star magnetospheres, including current efforts to characterize the largest magnetosphere surrounding an O star: that of NGC 1624-2. We also discuss the potential of the `analytic dynamical magnetosphere' (ADM) model to interpret multi-wavelength observations. Finally, we examine the possible effects of -- heretofore undetected -- small-scale magnetic fields on massive star winds and compare their hypothetical consequences to existing, unexplained observations.

HD 156324 is an SB3 (B2V/B5V/B5V) system in the Sco OB4 association. The He-strong primary possesses both a strong magnetic field, and H$\alpha$ emission believed to originate in its Centrifugal Magnetosphere (CM). We analyse a large spectroscopic and high-resolution spectropolarimetric dataset. The radial velocities (RVs) indicate that the system is composed of two sub-systems, which we designate A and B. Period analysis of the RVs of the three components yields orbital periods $P_{\rm orb} = 1.5806(1)$~d for the Aa and Ab components, and 6.67(2)~d for the B component, a PGa star. Period analysis of the longitudinal magnetic field \bz~and H$\alpha$ equivalent widths, which should both be sensitive to the rotational period $P_{\rm rot}$ of the magnetic Aa component, both yield $\sim$1.58~d. Since $P_{\rm orb} = P_{\rm rot}$ Aa and Ab must be tidally locked. Consistent with this, the orbit is circularized, and the rotational and orbital inclinations are identical within uncertainty, as are the semi-major axis and the Kepler corotation radius. The star's H$\alpha$ emission morphology differs markedly from both theoretical and observational expectations in that there is only one, rather than two, emission peaks. We propose that this unusual morphology may be a consequence of modification of the gravitocentrifugal potential by the presence of the close stellar companion. We also obtain upper limits on the magnetic dipole strength $B_{\rm d}$ for the Ab and B components, respectively finding $B_{\rm d} < 2.6$~kG and $<0.7$~kG.

In this contribution, we present BRITE observations of the early-B supergiants $\epsilon$ Ori and $\kappa$ Ori. We perform a preliminary analysis of the data acquired over the first two Orion observing runs. We evaluate whether they are compatible with co-rotating bright spots and discuss the challenges of such an approach.

Magnetic confinement of the winds of hot, massive stars has far-reaching consequences on timescales ranging from hours to Myr. Understanding the long-term effects of this interplay has already led to the identification of two new evolutionary pathways to form `heavy' stellar mass black holes and pair-instability supernova even at galactic metallicity. We are performing 1D stellar evolution model calculations that, for the first time, account for the surface effects and the time evolution of fossil magnetic fields. These models will be thoroughly confronted with observations and will potentially lead to a significant revision of the derived parameters of observed magnetic massive stars.

We present a new analysis of the rapidly oscillating Ap (roAp) star, 2MASS J$19400781-4420093$ (J1940; $V=13.1$). The star was discovered using SuperWASP broadband photometry to have a frequency of 176.39 d$^{-1}$ (2041.55 $\mu$Hz; $P = 8.2$ min; Holdsworth et al. 2014a) and is shown here to have a peak-to-peak amplitude of 34 mmag. J1940 has been observed during three seasons at the South African Astronomical Observatory, and has been the target of a Whole Earth Telescope campaign. The observations reveal that J1940 pulsates in a distorted quadrupole mode with unusual pulsational phase variations. A higher signal-to-noise ratio spectrum has been obtained since J1940's first announcement, which allows us to classify the star as A7 Vp Eu(Cr). The observing campaigns presented here reveal no pulsations other than the initially detected frequency. We model the pulsation in J1940 and conclude that the pulsation is distorted by a magnetic field of strength $1.5$ kG. A difference in the times of rotational maximum light and pulsation maximum suggests a significant offset between the spots and pulsation axis, as can be seen in roAp stars.

Magnetic massive stars comprise approximately 10% of the total OB star population. Modern spectropolarimetry shows these stars host strong, stable, large-scale, often nearly dipolar surface magnetic fields of 1 kG or more. These global magnetic fields trap and deflect outflowing stellar wind material, forming an anisotropic magnetosphere that can be probed with wind-sensitive UV resonance lines. Recent HST UV spectra of NGC 1624-2, the most magnetic O star observed to date, show atypically unsaturated P-Cygni profiles in the Civ resonant doublet, as well as a distinct variation with rotational phase. We examine the effect of non-radial, magnetically-channeled wind outflow on P-Cygni line formation, using a Sobolev Exact Integration (SEI) approach for direct comparison with HST UV spectra of NGC 1624-2. We demonstrate that the addition of a magnetic field desaturates the absorption trough of the P-Cygni profiles, but further efforts are needed to fully account for the observed line profile variation. Our study thus provides a first step toward a broader understanding of how strong magnetic fields affect mass loss diagnostics from UV lines.

Recent spectropolarimetric surveys of bright, hot stars have found that ~10% of OB-type stars contain strong (mostly dipolar) surface magnetic fields (~kG). The prominent paradigm describing the interaction between the stellar winds and the surface magnetic field is the magnetically confined wind shock (MCWS) model. In this model, the stellar wind plasma is forced to move along the closed field loops of the magnetic field, colliding at the magnetic equator, and creating a shock. As the shocked material cools radiatively it will emit X-rays. Therefore, X-ray spectroscopy is a key tool in detecting and characterizing the hot wind material confined by the magnetic fields of these stars. Some B-type stars are found to have very short rotational periods. The effects of the rapid rotation on the X-ray production within the magnetosphere have yet to be explored in detail. The added centrifugal force due to rapid rotation is predicted to cause faster wind outflows along the field lines, leading to higher shock temperatures and harder X-rays. However, this is not observed in all rapidly rotating magnetic B-type stars. In order to address this from a theoretical point of view, we use the X-ray Analytical Dynamical Magnetosphere (XADM) model, originally developed for slow rotators, with an implementation of new rapid rotational physics. Using X-ray spectroscopy from ESA's XMM-Newton space telescope, we observed 5 rapidly rotating B-type stars to add to the previous list of observations. Comparing the observed X-ray luminosity and hardness ratio to that predicted by the XADM allows us to determine the role the added centrifugal force plays in the magnetospheric X-ray emission of these stars.

Large-scale dipolar surface magnetic fields have been detected in a fraction of OB stars, however only few stellar evolution models of massive stars have considered the impact of these fossil fields. We are performing 1D hydrodynamical model calculations taking into account evolutionary consequences of the magnetospheric-wind interactions in a simplified parametric way. Two effects are considered: i) the global mass-loss rates are reduced due to mass-loss quenching, and ii) the surface angular momentum loss is enhanced due to magnetic braking. As a result of the magnetic mass-loss quenching, the mass of magnetic massive stars remains close to their initial masses. Thus magnetic massive stars - even at Galactic metallicity - have the potential to be progenitors of `heavy' stellar mass black holes. Similarly, at Galactic metallicity, the formation of pair instability supernovae is plausible with a magnetic progenitor.

Near-solar metallicity (and low-redshift) Pair-Instability Supernova (PISN) candidates challenge stellar evolution models. Indeed, at such a metallicity, even an initially very massive star generally loses so much mass by stellar winds that it will avoid the electron-positron pair-creation instability. We use recent results showing that a magnetic field at the surface of a massive star can significantly reduce its effective mass-loss rate to compute magnetic models of very massive stars (VMSs) at solar metallicity and explore the possibility that such stars end as PISNe. We implement the quenching of the mass loss produced by a surface dipolar magnetic field into the Geneva stellar evolution code and compute new stellar models with an initial mass of $200\,M_\odot$ at solar metallicity, with and without rotation. It considerably reduces the total amount of mass lost by the star during its life. For the non-rotating model, the total (CO-core) mass of the models is $72.8\,M_\odot$ ($70.1\,M_\odot$) at the onset of the electron-positron pair-creation instability. For the rotating model, we obtain $65.6\,M_\odot$ ($62.4\,M\odot$). In both cases, a significant fraction of the internal mass lies in the region where pair instability occurs in the $\log(T)-\log(\rho)$ plane. The interaction of the reduced mass loss with the magnetic field efficiently brakes the surface of the rotating model, producing a strong shear and hence a very efficient mixing that makes the star evolve nearly homogeneously. The core characteristics of our models indicate that solar metallicity models of magnetic VMSs may evolve to PISNe (and pulsation PISNe).

The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun.

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa.

We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four stars were detected at 3cm, whereas three were detected at 1cm and only one star was detected at 13cm. The detected B-type stars are significantly more radio luminous than the non-detected ones, which is not the case for O-type stars. The non-detections at 13cm are interpreted as due to thermal free-free absorption. Mass-loss rates were estimated using 3cm flux densities and were compared with theoretical mass-loss rates, which assume free-free emission. For HD 37742, the two values of the mass-loss rates were in good agreement, possibly suggesting that the radio emission for this star is mainly thermal. For the other three stars, the estimated mass-loss rates from radio observations were much higher than those expected from theory, suggesting either a possible contribution from non- thermal emission from the magnetic star or thermal or non-thermal emission due to interacting winds of the binary system, especially for HD 47129. All the detected stars are predicted to host centrifugal magnetospheres except HD 37742, which is likely to host a dynamical magnetosphere. This suggests that non-thermal radio emission is favoured in stars with centrifugal magnetospheres.

We present the analysis performed on spectropolarimetric data of 97 O-type targets included in the framework of the MiMeS (Magnetism in Massive Stars) Survey. Mean Least-Squares Deconvolved Stokes I and V line profiles were extracted for each observation, from which we measured the radial velocity, rotational and non-rotational broadening velocities, and longitudinal magnetic field. The investigation of the Stokes I profiles led to the discovery of 2 new multi-line spectroscopic systems (HD46106, HD204827) and confirmed the presence of a suspected companion in HD37041. We present a modified strategy of the Least-Squares Deconvolution technique aimed at optimising the detection of magnetic signatures while minimising the detection of spurious signatures in Stokes V. Using this analysis, we confirm the detection of a magnetic field in 6 targets previously reported as magnetic by the MiMeS collaboration (HD108, HD47129A2, HD57682, HD148937, CPD-28 2561, and NGC 1624-2), as well as report the presence of signal in Stokes V in 3 new magnetic candidates (HD36486, HD162978, HD199579). Overall, we find a magnetic incidence rate of 7+/-3%, for 108 individual O stars (including all O-type components part of multi-line systems), with a median uncertainty of the longitudinal field measurements of about 50\,G. An inspection of the data reveals no obvious biases affecting the incidence rate or the preference for detecting magnetic signatures in the magnetic stars. Similar to A- and B-type stars, we find no link between the stars' physical properties (e.g. Teff, mass, age) and the presence of a magnetic field. However, the Of?p stars represent a distinct class of magnetic O-type stars.

We present 16 new, and confirm 7 previously identified, DA white dwarfs in the Kepler field through ground-based spectroscopy with the Hale 200", Kitt Peak 4-meter, and Bok 2.3-meter telescopes. Using atmospheric models we determine their effective temperatures and surface gravities to constrain their position with respect to the ZZ Ceti (DA pulsator) instability strip, and look for the presence or absence of pulsation with Kepler's unprecedented photometry. Our results are as follows: i) From our measurements of temperature and surface gravity, 12 of the 23 DA white dwarfs from this work fall well outside of the instability strip. The Kepler photometry available for 11 of these WDs allows us to confirm that none are pulsating. One of these eleven happens to be a presumed binary, KIC 11604781, with a period of ~5 days. ii) The remaining 11 DA white dwarfs are instability strip candidates, potentially falling within the current, empirical instability strip, after accounting for uncertainties. These WDs will help constrain the strip's location further, as eight are near the blue edge and three are near the red edge of the instability strip. Four of these WDs do not have Kepler photometry, so ground-based photometry is needed to determine the pulsation nature of these white dwarfs. The remaining seven have Kepler photometry available, but do not show any periodicity on typical WD pulsation timescales.

We employ high resolution spectroscopy and spectropolarimetry to derive the physical properties and magnetic characteristics of the multiple system HD 164492C, located in the young open cluster M20. The spectrum reveals evidence of 3 components: a broad-lined early B star (HD 164492C1), a narrow-lined early B star (HD 164492C2), and a late B star (HD 164492C3). Components C2 and C3 exhibit significant ($>100$ km/s) bulk radial velocity variations with a period of $12.5351(7)$ d that we attribute to eccentric binary motion around a common centre-of-mass. Component C1 exhibits no detectable radial velocity variations. Using constraints derived from modeling the orbit of the C2+C3 binary and from synthesis of the combined spectrum, we determine the approximate physical characteristics of the components. We conclude that a coherent evolutionary solution consistent with the published age of M20 implies a distance to the system of $0.9\pm 0.2$ kpc, corresponding to the smallest published values. We confirm the detection of a strong magnetic field in the combined spectrum. The field is clearly associated with the broad-lined C1 component of the system. Repeated measurement of the longitudinal magnetic field allows the derivation of the rotation period of the magnetic star, $P_{\rm rot}=1.36986(6)$ d. We derive the star's magnetic geometry, finding $i=63\pm 6$ deg, $\beta=33\pm 6$ deg and a dipole polar strength $B_{\rm d}=7.9^{+1.2}_{-1.0}$ kG. Strong emission - varying according to the magnetic period - is detected in the H$\alpha$ profile. This is consistent with the presence of a centrifugal magnetosphere surrounding the rapidly rotating magnetic C1 component.

(abridged version) Context : The detection of bright, hard, and variable X-ray emission in Tr16-22 prompted spectropolarimetric observations of this star, which in turn led to the discovery of a surface magnetic field. Aims : We want to further constrain the properties of this star, in particular to verify whether X-ray variations are correlated to changes in optical emission lines and magnetic field strength, as expected from the oblique rotator model that is widely accepted for magnetic O stars. Methods: We have obtained new low-resolution spectropolarimetric and long-term high-resolution spectroscopic monitoring of Tr16-22, and we also analyse new, serendipitous X-ray data. Results: The new X-ray observations are consistent with previous data, but their addition does not help to solve the ambiguity in the variation timescale because of numerous aliases. No obvious periodicity or any large variations are detected in the spectropolarimetric data of Tr16-22 obtained over three months. The derived field values appear to be in line with previous measurements, suggesting constancy of the field (though the possibility of small, short-term field variations cannot be excluded). Variations in the equivalent widths of Halpha are very small, and they do not appear to be related to the X-ray timescale; the overall lack of large variations in optical emission lines is consistent with the magnetic field constancy. In addition, variations of the radial velocities indicate that Tr16-22 is probably a SB1 binary with a very long period.

Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spatially complex and temporally variable, making it difficult to derive observational signatures and discern their overall scaling trends.Within a simplified, steady-state analysis based on overall conservation principles, we present here an "analytic dynamical magnetosphere" (ADM) model that provides explicit formulae for density, temperature and flow speed in each of these three components -- wind outflow, hot post-shock gas, and cooled inflow -- as a function of colatitude and radius within the closed (presumed dipole) field lines of the magnetosphere. We compare these scalings with time-averaged results from MHD simulations, and provide initial examples of application of this ADM model for deriving two key observational diagnostics, namely hydrogen H-alpha emission line profiles from the cooled infall, and X-ray emission from the hot post-shock gas. We conclude with a discussion of key issues and advantages in applying this ADM formalism toward derivation of a broader set of observational diagnostics and scaling trends for massive stars with such dynamical magnetospheres.

A significant fraction of massive main-sequence stars show strong, large-scale magnetic fields. The origin of these fields, their lifetimes, and their role in shaping the characteristics and evolution of massive stars are currently not well understood. We compile a catalogue of 389 massive main-sequence stars, 61 of which are magnetic, and derive their fundamental parameters and ages. The two samples contain stars brighter than magnitude 9 in the V band and range in mass between 5 and 100 Msun. We find that the fractional main-sequence age distribution of all considered stars follows what is expected for a magnitude limited sample, while that of magnetic stars shows a clear decrease towards the end of the main sequence. This dearth of old magnetic stars is independent of the choice of adopted stellar evolution tracks, and appears to become more prominent when considering only the most massive stars. We show that the decreasing trend in the distribution is significantly stronger than expected from magnetic flux conservation. We also find that binary rejuvenation and magnetic suppression of core convection are unlikely to be responsible for the observed lack of older magnetic massive stars, and conclude that its most probable cause is the decay of the magnetic field, over a time span longer than the stellar lifetime for the lowest considered masses, and shorter for the highest masses. We then investigate the spin-down ages of the slowly rotating magnetic massive stars and find them to exceed the stellar ages by far in many cases. The high fraction of very slowly rotating magnetic stars thus provides an independent argument for a decay of the magnetic fields.

Magnetically confined winds of early-type stars are expected to be sources of bright and hard X-rays. To clarify the systematics of the observed X-ray properties, we have analyzed a large series of Chandra and XMM observations, corresponding to all available exposures of known massive magnetic stars (over 100 exposures covering ~60% of stars compiled in the catalog of Petit et al. 2013). We show that the X-ray luminosity is strongly correlated with the stellar wind mass-loss-rate, with a power-law form that is slightly steeper than linear for the majority of the less luminous, lower-Mdot B stars and flattens for the more luminous, higher-Mdot O stars. As the winds are radiatively driven, these scalings can be equivalently written as relations with the bolometric luminosity. The observed X-ray luminosities, and their trend with mass-loss rates, are well reproduced by new MHD models, although a few overluminous stars (mostly rapidly rotating objects) exist. No relation is found between other X-ray properties (plasma temperature, absorption) and stellar or magnetic parameters, contrary to expectations (e.g. higher temperature for stronger mass-loss rate). This suggests that the main driver for the plasma properties is different from the main determinant of the X-ray luminosity. Finally, variations of the X-ray hardnesses and luminosities, in phase with the stellar rotation period, are detected for some objects and they suggest some temperature stratification to exist in massive stars' magnetospheres.

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 observed NGC 1624-2, the O-type star with the largest known magnetic field Bp~20 kG), in X-rays with the ACIS-S camera onboard the Chandra X-ray Observatory. Our two observations were obtained at the minimum and maximum of the periodic Halpha emission cycle, corresponding to the rotational phases where the magnetic field is the closest to equator-on and pole-on, respectively. With these observations, we aim to characterise the star's magnetosphere via the X-ray emission produced by magnetically confined wind shocks. Our main findings are: (i) The observed spectrum of NGC 1624-2 is hard, similar to the magnetic O-type star Theta 1 Ori C, with only a few photons detected below 0.8 keV. The emergent X-ray flux is 30% lower at the Halpha minimum phase. (ii) Our modelling indicated that this seemingly hard spectrum is in fact a consequence of relatively soft intrinsic emission, similar to other magnetic Of?p stars, combined with a large amount of local absorption (~1-3 x 10^22 cm^-2). This combination is necessary to reproduce both the prominent Mg and Si spectral features, and the lack of flux at low energies. NGC 1624-2 is intrinsically luminous in X-rays (log LX emission ~ 33.4) but 70-95% of the X-ray emission produced by magnetically confined wind shocks is absorbed before it escapes the magnetosphere (log LX ISM corrected ~ 32.5). (iii) The high X-ray luminosity, its variation with stellar rotation, and its large attenuation are all consistent with a large dynamical magnetosphere with magnetically confined wind shocks.

In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 12 obtained from high resolution Stokes $V$ spectra using the ESPaDOnS (CFHT) and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes $V$ spectra obtained with the DimaPol spectropolarimeter (DAO). HD 23478 was one of two rapidly rotating stars identified as potential "centrifugal magnetosphere" hosts based on IR observations from the Apache Point Observatory Galactic Evolution Experiment survey. We derive basic physical properties of this star including its mass ($M=6.1^{+0.8}_{-0.7}\,M_\odot$), effective temperature ($T_{\rm eff}=20\pm2\,$kK), radius ($R=2.7^{+1.6}_{-0.9}\,R_\odot$), and age ($\tau_{\rm age}=3^{+37}_{-1}\,$Myr). We repeatedly detect weakly-variable Zeeman signatures in metal, He and H lines in all our observations corresponding to a longitudinal magnetic field of $\langle B_z\rangle\approx-2.0\,$kG. The rotational period is inferred from Hipparcos photometry ($P_{\rm rot}=1.0498(4)\,$d). Under the assumption of the Oblique Rotator Model, our obsevations yield a surface dipole magnetic field of strength $B_d\geq9.5\,$kG that is approximately aligned with the stellar rotation axis. We confirm the presence of strong and broad H$\alpha$ emission and gauge the volume of this star's centrifugal magnetosphere to be consistent with those of other H$\alpha$ emitting centrifugal magnetosphere stars based on the large inferred Alfv\'en to Kepler radius ratio.

We have carried out a study of radio emission from a small sample of magnetic O- and B-type stars using the Giant Metrewave Radio Telescope, with the goal of investigating their magnetospheres at low frequencies. These are the lowest frequency radio measurements ever obtained of hot magnetic stars. The observations were taken at random rotational phases in the 1390 and the 610 MHz bands. Out of the 8 stars, we detect five B-type stars in both the 1390 and the 610 MHz bands. The O-type stars were observed only in the 1390 MHz band, and no detections were obtained. We explain this result as a consequence of free-free absorption by the free-flowing stellar wind exterior to the closed magnetosphere. We also study the variability of individual stars. One star - HD 133880 - exhibits remarkably strong and rapid variability of its low frequency flux density. We discuss the possibility of this emission being coherent emission as reported for CU Vir by Trigilio et al. (2000).

We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares Deconvolution yield a Zeeman signature detected in the mean Stokes V profile corresponding to phase 0.5 of the spectral ephemeris. Interpreting the 73.41 d period as the stellar rotational period, we have phased the equivalent widths and inferred longitudinal field measurements. The phased magnetic data exhibit a weak sinusoidal variation, with maximum of about 565 G at phase 0.5, and a minimum of about -335 G at phase 0.0, with extrema approximately in phase with the (double-wave) Halpha equivalent width variation. Modeling of the Halpha equivalent width variation assuming a quasi-3D magnetospheric model produces a unique solution for the ambiguous couplet of inclination and magnetic obliquity angles: (i, beta) or (beta, i)=(35 deg,90 deg). Adopting either geometry, the longitudinal field variation yields a dipole polar intensity Bd=2.6\pm 0.9~kG, consistent with that obtained from direct modelling of the Stokes V profiles. We derive a wind magnetic confinement parameter eta*\simeq 100, leading to an Alfv\'en radius RA\simeq 3-5~R*, and a Kepler radius RK\simeq 20~R*. This supports a physical scenario in which the Halpha emission and other line variability have their origin in an oblique, co-rotating 'dynamical magnetosphere' structure resulting from a magnetically channeled wind. Nevertheless, the details of the formation of spectral lines and their variability within this framework remain generally poorly understood.

In the context of the MiMeS survey of magnetism in massive stars, 85 classical Be stars were observed in circular polarization with the aim of detecting magnetic fields at their surfaces. No direct evidence of such fields is found, in contrast to the significant rate of detection (5-10\%) in non-Be B-type stars. In this paper we describe the sample properties, the methodology and the data quality. We describe a novel method, previously applied to Herbig Ae/Be stars, that allows us to infer upper limits on organized (dipolar) magnetic fields present in the photospheres of our targets. We review the characteristics and robustness of this null result, and discuss its implications.

The evolution of massive stars is still partly unconstrained. Mass, metallicity, mass loss and rotation are the main drivers of stellar evolution. Binarity and magnetic field may also significantly affect the fate of massive stars. Our goal is to investigate the evolution of single O stars in the Galaxy. For that, we use a sample of 74 objects comprising all luminosity classes and spectral types from O4 to O9.7. We rely on optical spectroscopy obtained in the context of the MiMeS survey of massive stars. We perform spectral modelling with the code CMFGEN. We determine the surface properties of the sample stars, with special emphasis on abundances of carbon, nitrogen and oxygen. Most of our sample stars have initial masses in the range 20 to 50 Msun. We show that nitrogen is more enriched and carbon/oxygen more depleted in supergiants than in dwarfs, with giants showing intermediate degrees of mixing. CNO abundances are observed in the range of values predicted by nucleosynthesis through the CNO cycle. More massive stars, within a given luminosity class, appear to be more chemically enriched than lower mass stars. We compare our results with predictions of three types of evolutionary models and show that, for two sets of models, 80% of our sample can be explained by stellar evolution including rotation. The effect of magnetism on surface abundances is unconstrained. Our study indicates that, in the 20-50 Msun mass range, the surface chemical abundances of most single O stars in the Galaxy are fairly well accounted for by stellar evolution of rotating stars.

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.

We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive cooling, the method presented here exploits the idea that the cooling post-shock plasma systematically passes through the temperature characteristic of distinct emission lines in the X-ray spectrum. In this way, the observed flux distribution among these X-ray lines can be used to construct the cumulative probability distribution of shock strengths that a typical wind parcel encounters as it advects through the wind. We apply this new method (Gayley 2014) to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. Correcting for wind absorption of the X-ray line emission is a crucial component of our analysis, and we use wind optical depth values derived from X-ray line-profile fitting (Cohen et al. 2014) in order to make that correction. The shock-heating rate results we derive for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least $10^6$ K as it advects through the wind, and the cumulative distribution of shock strengths is a strongly decreasing function of temperature, consistent with a negative power-law of index $n \approx 3$, implying a marginal distribution of shock strengths that scales as $T^{-4}$, and with hints of an even steeper decline or cut-off above $10^7$ K.

Using an extensive X-ray dataset, we analyze the X-ray emission of the massive O-star Tr16-22, which was recently found to be magnetic. Its bright X-ray emission is found to be modulated with a ~54d period. This timescale should represent the rotational timescale of the star, as for other magnetic massive stars. In parallel, new spectropolarimetric data confirm the published magnetic detection.

We have obtained high-quality magnetic field measurements of 19 sharp-lined B-type stars with precisely-measured N/C abundance ratios. Our primary goal is to test the idea that a magnetic field may explain extra drag (through the wind) on the surface rotation, thus producing more internal shear and mixing, and thus could provide an explanation for the appearance of slowly rotating N-rich main sequence B stars.

Magnetic B-type stars are often host to Centrifugal Magnetospheres (CMs). Here we describe the results of a population study encompassing the full sample of known magnetic early B-type stars, focusing on those with detectable CMs. We present revised rotational and magnetic parameters for some stars, clarifying their positions on the rotation-confinement diagram, and? find that plasma densities within their CMs are much lower than those predicted by centrifugal breakout.

OB stars exhibit various types of spectral variability associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These are proposed to be caused by either magnetic fields or non-radial pulsations (NRPs). In this paper, we evaluate the possible relation between large-scale, dipolar magnetic fields and the DAC phenomenon by investigating the magnetic properties of a sample of 13 OB stars exhibiting well-documented DAC behaviour. Using high-precision spectropolarimetric data acquired in part in the context of the Magnetism in Massive Stars (MiMeS) project, we find no evidence for surface dipolar magnetic fields in any of these stars. Using Bayesian inference, we compute upper limits on the strengths of the fields and use these limits to assess two potential mechanisms by which the field may influence wind outflow: magnetic wind confinement and local photospheric brightness enhancements. Within the limits we derive, both mechanisms fail to provide a systematic process capable of producing DACs in all of the stars of our sample. Therefore, this implies that dipolar fields are highly unlikely to be responsible for these structures in all massive stars, meaning that some other mechanism must come into play.

We obtained 128 high signal-to-noise ratio Stokes $V$ spectra of the B3V star $\iota$~Her on 5 consecutive nights in 2012 with the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, with the aim of searching for the presence of weak and/or complex magnetic fields. Least-Squares Deconvolution mean profiles were computed from individual spectra, averaged over individual nights and over the entire run. No Zeeman signatures are detected in any of the profiles. The longitudinal magnetic field in the grand average profile was measured to be $-0.24 \pm 0.32 $ G, as compared to $-0.22 \pm 0.32$ G in the null profile. Our observations therefore provide no evidence for the presence of Zeeman signatures analogous to those observed in the A0V star Vega by Lignieres et al. (2009). We interpret these observations in three ways. First, we compare the LSD profiles with synthetic Stokes $V$ profiles corresponding to organised (dipolar) magnetic fields, for which we find an upper limit of about 8~G on the polar strength of any surface dipole present. Secondly, we compare the grand average profile with calculations corresponding to the random magnetic spot topologies of Kochukhov & Sudnik (2013), inferring that spots, if present, of 2 degrees radius with strengths of 2-4 G and a filling factor of 50% should have been detected in our data. Finally, we compare the observations with synthetic $V$ profiles corresponding to the surface magnetic maps of Vega (Petit et al. 2010) computed for the spectral characteristics of $\iota$ Her. We conclude that while it is unlikely we would have detected a magnetic field identical to Vega's, we would have likely detected one with a peak strength of about 30 G, i.e. approximately four times as strong as that of Vega.

To understand the origin of the magnetic fields in massive stars as well as their impact on stellar internal structure, evolution, and circumstellar environment, within the MiMeS project, we searched for magnetic objects among a large sample of massive stars, and build a sub-sample for in-depth follow-up studies required to test the models and theories of fossil field origins, magnetic wind confinement and magnetospheric properties, and magnetic star evolution. We obtained high-resolution spectropolarimetric observations of a large number of OB stars thanks to three large programs that have been allocated on the high-resolution spectropolarimeters ESPaDOnS, Narval, and the polarimetric module HARPSpol of the HARPS spectrograph. We report here on the methods and first analysis of the HARPSpol magnetic detections. We identified the magnetic stars using a multi-line analysis technique. Then, when possible, we monitored the new discoveries to derive their rotation periods, which are critical for follow-up and magnetic mapping studies. We also performed a first-look analysis of their spectra and identified obvious spectral anomalies (e.g., abundance peculiarities, Halpha emission), which are also of interest for future studies. In this paper, we focus on eight of the 11 stars in which we discovered or confirmed a magnetic field from the HARPSpol LP sample (the remaining three were published in a previous paper). Seven of the stars were detected in early-type Bp stars, while the last star was detected in the Ap companion of a normal early B-type star. We report obvious spectral and multiplicity properties, as well as our measurements of their longitudinal field strengths, and their rotation periods when we are able to derive them. We also discuss the presence or absence of Halpha emission with respect to the theory of centrifugally-supported magnetospheres. (Abriged)

We use 2D MHD simulations to examine the effects of radiative cooling and inverse Compton (IC) cooling on X-ray emission from magnetically confined wind shocks (MCWS) in magnetic massive stars with radiatively driven stellar winds. For the standard dependence of mass loss rate on luminosity $\Mdot \sim L^{1.7} $, the scaling of IC cooling with $L$ and radiative cooling with $\Mdot$ means that IC cooling become formally more important for lower luminosity stars. However, because the sense of the trends is similar, we find the overall effect of including IC cooling is quite modest. More significantly, for stars with high enough mass loss to keep the shocks radiative, the MHD simulations indicate a linear scaling of X-ray luminosity with mass loss rate; but for lower luminosity stars with weak winds, X-ray emission is reduced and softened by a {\em shock retreat} resulting from the larger post-shock cooling length, which within the fixed length of a closed magnetic loop forces the shock back to lower pre-shock wind speeds. A semi-analytic scaling analysis that accounts both for the wind magnetic confinement and this shock retreat yields X-ray luminosities that have a similar scaling trend, but a factor few higher values, compared to time-averages computed from the MHD simulations. The simulation and scaling results here thus provide a good basis for interpreting available X-ray observations from the growing list of massive stars with confirmed large-scale magnetic fields.

Magnetic wind confinement has been proposed as one explanation for the complex wind structures of supergiant stars of spectral types B and A. Observational investigation of this hypothesis was undertaken using high-resolution ({\lambda}/{\Delta}{\lambda} {\sim} 65,000) circular polarization (Stokes V ) spectra of six late B and early A type supergiants ({\beta} Ori, B8Iae; 4 Lac, B9Iab; {\eta} Leo, A0Ib; HR1040, A0Ib; {\alpha} Cyg, A2Iae; {\nu} Cep, A2Iab), obtained with the instruments ESPaDOnS and Narval at the Canada-France-Hawaii Telescope and the Bernard Lyot Telescope. Least Squares Deconvolution (LSD) analysis of the Stokes V spectra of all stars yields no evidence of a magnetic field, with best longitudinal field 1{\sigma} error bars ranging from {\sim}0.5 to {\sim}4.5 G for most stars. Spectrum synthesis analysis of the LSD profiles using Bayesian inference yields an upper limit with 95.4% credibility on the polar strength of the (undetected) surface dipole fields of individual stars ranging from 3 to 30 G. These results strongly suggest that magnetic wind confinement due to organized dipolar magnetic fields is not the origin of the wind variability of BA supergiant stars. Upper limits for magnetic spots may also be inconsistent with magnetic wind confinement in the limit of large spot size and filling factor, depending on the adopted wind parameters. Therefore, if magnetic spots are responsible for the wind variability of BA supergiant stars, they likely occupy a small fraction of the photosphere.

The Magnetism in Massive Stars (MiMeS) project represents the largest systematic survey of stellar magnetism ever undertaken. Based on a sample of over 550 Galactic B and O-type stars, the MiMeS project has derived the basic characteristics of magnetism in hot, massive stars. Herein we report preliminary results.

OB stars are known to exhibit various types of wind variability, as detected in their ultraviolet spectra, amongst which are the ubiquitous discrete absorption components (DACs). These features have been associated with large-scale azimuthal structures extending from the base of the wind to its outer regions: corotating interaction regions (CIRs). There are several competing hypotheses as to which physical processes may perturb the star's surface and generate CIRs, including magnetic fields and non radial pulsations (NRPs), the subjects of this paper with a particular emphasis on the former. Although large-scale magnetic fields are ruled out, magnetic spots deserve further investigation, both on the observational and theoretical fronts.

Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfv\'en radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organisation of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff > 16 kK. Using these parameters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot of RK vs. RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK, vs. dynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Halpha emission, depend on a star's position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Halpha emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall timescale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown timescales. [abbreviated]

We present the first fully 3D MHD simulation for magnetic channeling and confinement of a radiatively driven, massive-star wind. The specific parameters are chosen to represent the prototypical slowly rotating magnetic O star \theta^1 Ori C, for which centrifugal and other dynamical effects of rotation are negligible. The computed global structure in latitude and radius resembles that found in previous 2D simulations, with unimpeded outflow along open field lines near the magnetic poles, and a complex equatorial belt of inner wind trapping by closed loops near the stellar surface, giving way to outflow above the Alfv\'{e}n radius. In contrast to this previous 2D work, the 3D simulation described here now also shows how this complex structure fragments in azimuth, forming distinct clumps of closed loop infall within the Alfv\'{e}n radius, transitioning in the outer wind to radial spokes of enhanced density with characteristic azimuthal separation of $15-20 \degr$. Applying these results in a 3D code for line radiative transfer, we show that emission from the associated 3D `dynamical magnetosphere' matches well the observed H\alpha emission seen from \theta^1 Ori C, fitting both its dynamic spectrum over rotational phase, as well as the observed level of cycle to cycle stochastic variation. Comparison with previously developed 2D models for Balmer emission from a dynamical magnetosphere generally confirms that time-averaging over 2D snapshots can be a good proxy for the spatial averaging over 3D azimuthal wind structure. Nevertheless, fully 3D simulations will still be needed to model the emission from magnetospheres with non-dipole field components, such as suggested by asymmetric features seen in the H\alpha equivalent-width curve of \theta^1 Ori C.

We report the detection of a strong, organized magnetic field in the secondary component of the massive O8III/I+O7.5V/III double-lined spectroscopic binary system HD 47129 (Plaskett's star), in the context of the Magnetism in Massive Stars (MiMeS) survey. Eight independent Stokes $V$ observations were acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope and the Narval spectropolarimeter at the T\'elescope Bernard Lyot. Using Least-Squares Deconvolution we obtain definite detections of signal in Stokes $V$ in 3 observations. No significant signal is detected in the diagnostic null ($N$) spectra. The Zeeman signatures are broad and track the radial velocity of the secondary component; we therefore conclude that the rapidly-rotating secondary component is the magnetized star. Correcting the polarized spectra for the line and continuum of the (sharp-lined) primary, we measured the longitudinal magnetic field from each observation. The longitudinal field of the secondary is variable and exhibits extreme values of $-810\pm 150$ G and $+680\pm 190$ G, implying a minimum surface dipole polar strength of $2850\pm 500$ G. In contrast, we derive an upper limit ($3\sigma$) to the primary's surface magnetic field of 230 G. The combination of a strong magnetic field and rapid rotation leads us to conclude that the secondary hosts a centrifugal magnetosphere fed through a magnetically confined wind. We revisit the properties of the optical line profiles and X-ray emission - previously interpreted as a consequence of colliding stellar winds - in this context. We conclude that HD 47129 represents a heretofore unique stellar system - a close, massive binary with a rapidly rotating, magnetized component - that will be a rich target for further study.

This paper presents a first observational investigation of the faint Of?cp star NGC 1624-2, yielding important new constraints on its spectral and physical characteristics, rotation, magnetic field strength, X-ray emission and magnetospheric properties. Modeling the spectrum and spectral energy distribution, we conclude that NGC 1624-2 is a main sequence star of mass M {\simeq} 30 M{\odot}, and infer an effective temperature of 35 {\pm} 2 kK and log g = 4.0 {\pm} 0.2. Based on an extensive time series of optical spectral observations we report significant variability of a large number of spectral lines, and infer a unique period of 157.99 {\pm} 0.94 d which we interpret as the rotational period of the star. We report the detection of a very strong - 5.35 {\pm} 0.5 kG - longitudinal magnetic field <Bz>, coupled with probable Zeeman splitting of Stokes I profiles of metal lines confirming a surface field modulus <B> of 14 {\pm} 1 kG, consistent with a surface dipole of polar strength >~ 20 kG. This is the largest magnetic field ever detected in an O-type star, and the first report of Zeeman splitting of Stokes I profiles in such an object. We also report the detection of reversed Stokes V profiles associated with weak, high-excitation emission lines of O iii, which we propose may form in the close magnetosphere of the star. We analyze archival Chandra ACIS-I X-ray data, inferring a very hard spectrum with an X-ray efficiency log Lx/Lbol = -6.4, a factor of 4 larger than the canonical value for O-type stars and comparable to that of the young magnetic O-type star {\theta}1 Ori C and other Of?p stars. Finally, we examine the probable magnetospheric properties of the star, reporting in particular very strong magnetic confinement of the stellar wind, with {\eta}* {\simeq} 1.5 {\times} 10^4, and a very large Alfven radius, RAlf = 11.4 R*.

The young O9.5 V spectroscopic binary Theta 2 Ori A shows moderately hard X-ray emission and relatively narrow X-ray lines, suggesting that it may be a Magnetically Confined Wind Shock (MCWS) source, similar to its more massive analogue Theta 1 Ori C. X-ray flares occurring near periastron led to the suggestion that the flares are produced via magnetic reconnection as magnetospheres on both components of the Theta 2 Ori A binary interact at closest approach. We use a series of high-resolution spectropolarimetric observations of Theta 2 Ori A to place an upper limit on the magnetic field strength of 135 G (95% credible region). Such a weak dipole field would not produce magnetic confinement, or a large magnetosphere. A sub-pixel analysis of the Chandra ACIS images of Theta 2 Ori A obtained during quiescence and flaring show that the hard, flaring X-rays are offset from the soft, quiescent emission by 0.4 arcsec. If the soft emission is associated with the A1/A2 spectroscopic binary, the offset and position angle of the hard, flaring source place it at the location of the intermediate-mass A3 companion, discovered via speckle interferometry. The spectropolarimetric and X-ray results taken together point to the A3 companion, not the massive A1/A2 binary, as the source of hard, flaring X-ray emission. We also discuss a similar analysis performed for the magnetic Bp star Sigma Ori E. We find a similar origin for its X-ray flaring.

X-ray surveys of the Carina nebula have revealed a few hard and luminous sources associated with early-type stars. Such unusual characteristics for the high-energy emission may be related to magnetically-confined winds. To search for the presence of magnetic fields in these objects, we performed a limited spectropolarimetric survey using the FORS instrument. The multi-object mode was used, so that a total of 21 OB stars could be investigated during a one-night-long run. A magnetic field was detected in two objects of the sample, with a 6 sigma significance; Tr16-22 and 13. Such a detection was expected for Tr16-22, as its X-ray emission is too bright, variable and hard, compared to other late-type O or O+OB systems. It is more surprising for Tr16-13, a poorly known star which so far has never shown any peculiar characteristics. Subsequent monitoring is now needed to ascertain the physical properties of these objects and enable a full modelling of their magnetic atmospheres and winds.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and complex magnetic field. Recently, Petit et al. 2011 presented the discovery of the first two tau Sco analogues -- HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS and Narval spectropolarimetric observations were obtained by the Magnetism in Massive Stars CFHT and TBL Large Programs, in order to characterize the stellar and magnetic properties of these stars. A magnetic field of similar surface strength was found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. We present additional phase-resolved observations secured by the MiMeS collaboration for HD 66665 in order to measure its magnetic geometry, and correlate that geometry with diagnostics of mass-loss.

Magnetic massive stars -- which are being discovered with increasing frequency -- represent a new category of wind-shaping mechanism for O and B stars. Magnetic channeling of these stars' radiation-driven winds, the Magnetically Confined Wind Shock paradigm, leads to the formation of a shock-heated magnetosphere, which can radiate X-rays, modify UV resonance lines, and create disks of Halpha emitting material. The dynamical properties of these magnetospheres are well understood from a theoretical point of view as an interplay between the magnetic wind confinement and rotation. However, the manifestations of magnetospheres across the spectrum may be more complex and various than first anticipated. On the other hand, recent advances in modeling these magnetospheres provide a key to better understand massive star winds in general. We will summarize the coordinated observational, theoretical, and modeling efforts from the Magnetism in Massive Star Project, addressing key outstanding questions regarding magnetosphere manifestations across the spectral domain.

In this paper we describe a Bayesian statistical method designed to infer the magnetic properties of stars observed using high-resolution circular spectropolarimetry in the context of large surveys. This approach is well suited for analysing stars for which the stellar rotation period is not known, and therefore the rotational phases of the observations are ambiguous. The model assumes that the magnetic observations correspond to a dipole oblique rotator, a situation commonly encountered in intermediate and high-mass stars. Using reasonable assumptions regarding the model parameter prior probability density distributions, the Bayesian algorithm determines the posterior probability densities corresponding to the surface magnetic field geometry and strength by performing a comparison between the observed and computed Stokes V profiles. Based on the results of numerical simulations, we conclude that this method yields a useful estimate of the surface dipole field strength based on a small number (i.e. 1 or 2) of observations. On the other hand, the method provides only weak constraints on the dipole geometry. The odds ratio, a parameter computed by the algorithm that quantifies the relative appropriateness of the magnetic dipole model versus the non-magnetic model, provides a more sensitive diagnostic of the presence of weak magnetic signals embedded in noise than traditional techniques. To illustrate the application of the technique to real data, we analyse seven ESPaDOnS and Narval observations of the early B-type magnetic star LP Ori. Insufficient information is available to determine the rotational period of the star and therefore the phase of the data; hence traditional modelling techniques fail to infer the dipole strength. In contrast, the Bayesian method allows a robust determination of the dipole polar strength, $B_d=911^{+138}_{-244}$ G.

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.)

We report the discovery and analysis of a very strong magnetic field in the rapidly rotating early B-type star HR 5907, based on observations obtained as part of the Magnetism in Massive Stars (MiMeS) project. We infer a rotation period of 0.508276 +0.000015/-0.000012 d from photometric and H{\alpha} EW measurements, making this the shortest period, non-degenerate, magnetic massive star known to date. From the comparison of IUE UV and optical spectroscopy with LTE BRUCE/KYLIE models we find a solid-angle integrated, uniform black-body temperature of 17 000 \pm 1000 K, a projected rotational velocity of 290 \pm 10 km/s, an equatorial radius of 3.1 \pm 0.2 R_sun, a stellar mass of 5.5 \pm 0.5 M_sun, and an inclination angle of the rotation axis to our line-of-sight of 70 \pm 10\circ. Our measurements of the longitudinal magnetic field, which vary between -500 and -2000 G, phase coherently with the rotation period and imply a surface dipole field strength of \sim15.7 kG. On the other hand, from fits to mean Least-Squares Deconvolved Stokes V line profiles we infer a dipole field strength of \sim10.4 kG. This disagreement may result from a magnetic configuration more complex than our model, and/or from the non-uniform helium surface abundance distribution. In either case we obtain a magnetic obliquity nearly aligned with the rotation axis ({\beta} = 7+2/-1\circ). Our optical spectroscopy also shows weak variability in carbon, silicon and nitrogen lines. The emission variability in hydrogen Balmer and Paschen lines indicates the presence of a dense, highly structured magnetosphere, interpreted as a centrifugally supported, magnetically confined circumstellar disk.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive OB stars due to its high X-ray activity, peculiar wind diagnostics and highly complex magnetic field. This paper presents the discovery of the first two tau Sco analogues - HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to that of tau Sco. ESPaDOnS spectropolarimetric observations were secured by the Magnetism in Massive Stars CFHT Large Program, in order to characterize the stellar and magnetic properties of these stars. CMFGEN modelling of optical ESPaDOnS spectra and archived IUE UV spectra showed that these stars have stellar parameters similar to those of tau Sco. A magnetic field of similar surface strength is found on both stars, reinforcing the connection between the presence of a magnetic field and wind peculiarities. However, additional phase-resolved observations will be required in order to assess the potential complexity of the magnetic fields, and verify if the wind anomalies are linked to this property.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its high X-ray activity and remarkable wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and wind anomalies (Petit et al. 2010). We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars.

We will review our recent analysis of the magnetic properties of the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the Canada-France-Hawaii telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. We discuss our most recent determination of the rotational period from longitudinal magnetic field measurements and H_alpha variability - the latter obtained from over a decade's worth of professional and amateur spectroscopic observations. Lastly, we report on our investigation of the magnetic field geometry and the effects of the field on the circumstellar environment.

Spectropolarimetric observations of HD 155806 - the hottest Galactic Oe star - were obtained with CFHT/ESPaDOnS to test the hypothesis that disk signatures in its spectrum are due to magnetic channeling and confinement of its stellar wind. We did not detect a dipole field of sufficient strength to confine the wind, and could not confirm previous reports of a magnetic detection. It appears that stellar magnetism is not responsible for producing the disk of HD 155806.

In some massive stars, magnetic fields are thought to confine the outflowing radiatively-driven wind. Although theoretical models and MHD simulations are able to illustrate the dynamics of such a magnetized wind, the impact of this wind-field interaction on the observable properties of a magnetic star - X-ray emission, photometric and spectral variability - is still unclear. The aim of this study is to examine the relationship between magnetism, stellar winds and X-ray emission of OB stars, by providing empirical observations and confronting theory. In conjunction with the COUP survey of the Orion Nebula Cluster, we carried out spectropolarimatric ESPaDOnS observations to determine the magnetic properties of massive OB stars of this cluster.

The B0.2 V magnetic star tau Sco stands out from the larger population of massive magnetic OB stars due to its remarkable, superionized wind, apparently related to its peculiar magnetic field - a field which is far more complex than the mostly-dipolar fields usually observed in magnetic OB stars. tau Sco is therefore a puzzling outlier in the larger picture of stellar magnetism - a star that still defies interpretation in terms of a physically coherent model. Recently, two early B-type stars were discovered as tau Sco analogues, identified by the striking similarity of their UV spectra to that of tau Sco, which was - until now - unique among OB stars. We present the recent detection of their magnetic fields by the MiMeS collaboration, reinforcing the connection between the presence of a magnetic field and a superionized wind. We will also present ongoing observational efforts undertaken to establish the precise magnetic topology, in order to provide additional constrains for existing models attempting to reproduce the unique wind structure of tau Sco-like stars.

The presence of magnetic fields at the surfaces of many massive stars has been suspected for decades, to explain the observed properties and activity of OB stars. However, very few genuine high-mass stars had been identified as magnetic before the advent of a new generation of powerful spectropolarimeters that has resulted in a rapid burst of precise information about the magnetic properties of massive stars. During this talk, I will briefly review modern methods used to diagnose magnetic fields of higher-mass stars, and summarize our current understanding of the magnetic properties of OB stars.

The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual stars with spectral types from B5-O4, discovering new magnetic fields in a dozen hot, massive stars. The quality of this spectral and magnetic mat\'eriel is very high, and the Collaboration is keen to connect with colleagues capable of exploiting the data in new or unforeseen ways. In this paper we review the structure of the MiMeS observing programs and report the status of observations, data modeling and development of related theory.

I will review our recent analysis of the magnetic properties of the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the Canada-France-Hawaii telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. I discuss our most recent determination of the rotational period from longitudinal magnetic field measurements and Halpha variability - the latter obtained from over a decade's worth of professional and amateur spectroscopic observations. Lastly, I will report on our investigation of the magnetic field geometry and the effects of the field on the circumstellar environment.

Seventy-eight high-resolution Stokes V, Q and U spectra of the B8Iae supergiant Rigel were obtained with the ESPaDOnS spectropolarimeter at CFHT and its clone NARVAL at TBL in the context of the Magnetism in Massive Stars (MiMeS) Large Program, in order to scrutinize this core-collapse supernova progenitor for evidence of weak and/or complex magnetic fields. In this paper we describe the reduction and analysis of the data, the constraints obtained on any photospheric magnetic field, and the variability of photospheric and wind lines.

We report the detection of a strong, organised magnetic field in the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. From the fitting of our spectra using NLTE model atmospheres we determined that HD 57682 is a $17^{+19}_{-9}$ M$_{\odot}$ star with a radius of $7.0^{+2.4}_{-1.8}$ R$_\odot$, and a relatively low mass-loss rate of $1.4^{+3.1}_{-0.95}\times10^{-9}$ M$_{\odot}$ yr$^{-1}$. The photospheric absorption lines are narrow, and we use the Fourier transform technique to infer $v\sin i=15\pm3$ km s$^{-1}$. This $v\sin i$ implies a maximum rotational period of 31.5 d, a value qualitatively consistent with the observed variability of the optical absorption and emission lines, as well as the Stokes $V$ profiles and longitudinal field. Using a Bayesian analysis of the velocity-resolved Stokes $V$ profiles to infer the magnetic field characteristics, we tentatively derive a dipole field strength of $1680^{+134}_{-356}$ G. The derived field strength and wind characteristics imply a wind that is strongly confined by the magnetic field.

We have obtained 40 high-resolution circular spectropolarimetric measurements of 12 slowly-pulsating B (SPB) stars, 8 Beta Cephei stars and two Be stars with the ESPaDOnS and NARVAL spectropolarimeters. The aim of these observations is to evaluate recent claims of a high incidence of magnetic field detections in stars of these types obtained using low-resolution spectropolarimetry by Hubrig (2006), Hubrig (2007) and Hubrig (2009). The precision achieved is generally comparable to or superior to that obtained by Hubrig et al., although our new observations are distinguished by their resolution of metallic and He line profiles, and their consequent sensitivity to magnetic fields of zero net longitudinal component. In the SPB stars we confirm the detection of magnetic field in one star (16 Peg), but find no evidence of the presence of fields in the remaining 11. In the Beta Cep stars, we detect a field in xi^1 CMa, but not in any of the remaining 7 stars. Finally, neither of the two B-type emission line stars shows any evidence of magnetic field. Based on our results, we conclude that fields are not common in SPB, Beta Cep and B-type emission line stars, consistent with the general rarity of fields in the broader population of main sequence B-type stars.

The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. The cornerstone of the project is the MiMeS Large Program at the Canada-France-Hawaii Telescope, which represents a dedication of 640 hours of telescope time from 2008-2012. The MiMeS Large Program will exploit the unique capabilities of the ESPaDOnS spectropolarimeter to obtain critical missing information about the poorly-studied magnetic properties of these important stars, to confront current models and to guide theory.

Massive stars are those stars with initial masses above about 8 times that of the sun, eventually leading to catastrophic explosions in the form of supernovae. These represent the most massive and luminous stellar component of the Universe, and are the crucibles in which the lion's share of the chemical elements are forged. These rapidly-evolving stars drive the chemistry, structure and evolution of galaxies, dominating the ecology of the Universe - not only as supernovae, but also during their entire lifetimes - with far-reaching consequences. Although the existence of magnetic fields in massive stars is no longer in question, our knowledge of the basic statistical properties of massive star magnetic fields is seriously incomplete. The Magnetism in Massive Stars (MiMeS) Project represents a comprehensive, multidisciplinary strategy by an international team of recognized researchers to address the "big questions" related to the complex and puzzling magnetism of massive stars. This paper present the first results of the MiMeS Large Program at the Canada-France-Hawaii Telescope.

In massive stars, magnetic fields are thought to confine the outflowing radiatively-driven wind, resulting in X-ray emission that is harder, more variable and more efficient than that produced by instability-generated shocks in non-magnetic winds. Although magnetic confinement of stellar winds has been shown to strongly modify the mass-loss and X-ray characteristics of massive OB stars, we lack a detailed understanding of the complex processes responsible. The aim of this study is to examine the relationship between magnetism, stellar winds and X-ray emission of OB stars. In conjunction with a Chandra survey of the Orion Nebula Cluster, we carried out spectropolarimatric ESPaDOnS observations to determine the magnetic properties of massive OB stars of this cluster.

Oe stars are a subset of the O-type stars that exhibit emission lines from a circumstellar disk. The recent detection of magnetic fields in some O-type stars suggests a possible explanation for the stability of disk-like structures around Oe stars. According to this hypothesis, the wind of the star is channeled by a dipolar magnetic field producing a disc in the magnetic equatorial plane. As a test of this model, we have obtained spectropolarimetric observations of the hottest Galactic Oe star HD 155806. Here we discuss the results and implications of those observations.

In massive stars, magnetic fields are thought to confine the outflowing radiatively-driven wind, resulting in X-ray emission that is harder, more variable and more efficient than that produced by instability-generated shocks in non-magnetic winds. Although magnetic confinement of stellar winds has been shown to strongly modify the mass-loss and X-ray characteristics of massive OB stars, we lack a detailed understanding of the complex processes responsible. The aim of this study is to examine the relationship between magnetism, stellar winds and X-ray emission of OB stars. In conjunction with a Chandra survey of the Orion Nebula Cluster, we carried out spectropolarimatric ESPaDOnS observations to determine the magnetic properties of massive OB stars of this cluster. We found of two new massive magnetic stars in the Orion Nebula Cluster: HD 36982 and HD 37061, for which the estimated dipole polar strengths are 1150 (+320 -200) G and 620 (+220 -170) G, respectively. However, the apparent lack of clear correlation between X-ray indicator and the presence of a magnetic fields brings forth new challenges for understanding the processes leading to X-ray emission in massive stars.

The young O-type star theta1 OriC, the brightest star of the Trapezium cluster in Orion, is one of only two known magnetic rotators among the O stars. However, not all spectroscopic variations of this star can be explained by the magnetic rotator model. We present results from a long-term monitoring to study these unexplained variations and to improve the stellar rotational period. We want to study long-term trends of the radial velocity of theta1 OriC, to search for unusual changes, to improve the established rotational period and to check for possible period changes. We combine a large set of published spectroscopic data with new observations and analyze the spectra in a homogeneous way. We study the radial velocity from selected photo-spheric lines and determine the equivalent width of the Halpha and HeII4686 lines. We find evidence for a secular change of the radial velocity of theta1 OriC that is consistent with the published interferometric orbit. We refine the rotational period of theta1 OriC and discuss the possibility of detecting period changes in the near future.

The origin of the magnetic fields in neutron stars, and the physical differences between magnetars and strongly magnetised radio pulsars are still under vigorous debate. It has been suggested that the properties of the progenitors of neutron stars (the massive OB stars), such as rotation, magnetic fields and mass, may play an important role in the outcome of core collapse leading to type II SNe. Therefore, knowing the magnetic properties of the progenitor OB stars would be an important asset for constraining models of stellar evolution leading to the birth of a neutron star. We present here the beginning of a broad study with the goal of characterising the magnetic properties of main sequence massive OB stars. We report the detection of two new massive magnetic stars in the Orion Nebula Cluster: Par 1772 (HD 36982) and NU Ori (HD 37061), for which the estimated dipole polar strengths, with 1 sigma error bars, are 1150 (+320,-200) G and 650 (+220,-170) G respectively.

Ferrario & Wickramasinghe (2006) explored the hypothesis that the magnetic fields of neutron stars are of fossil origin. In this context, they predicted the field distribution of the progenitor OB stars, finding that 5 per cent of main sequence massive stars should have fields in excess of 1kG. We have carried out sensitive ESPaDOnS spectropolarimetric observations to search for direct evidence of such fields in all massive B- and O-type stars in the Orion Nebula Cluster star-forming region. We have detected unambiguous Stokes V Zeeman signatures in spectra of three out of the eight stars observed (38%). Using a new state-of-the-art Bayesian analysis, we infer the presence of strong (kG), organised magnetic fields in their photospheres. For the remaining five stars, we constrain any dipolar fields in the photosphere to be weaker than about 200G. Statistically, the chance of finding three ~kG fields in a sample of eight OB stars is quite low (less than 1%) if the predictions of Ferrario & Wickramasinghe are correct. This implies that either the magnetic fields of neutron stars are not of fossil origin, that the flux-evolution model of Ferrario & Wickramasinghe is incomplete, or that the ONC has unusual magnetic properties. We are undertaking a study of other young star clusters, in order to better explore these possibilities.

Magnetic fields have been frequently invoked as a likely source of variability and confinement of the winds of massive stars. To date, the only magnetic field detected in O-type stars are those of theta 1 Ori C (HD 37022; Donati et al. 2002), the brightest and most massive member of the Orion Nebula Cluster (ONC), and HD 191612 (Donati et al. 2006). Notably, theta 1 Ori C is an intense X-ray emitter, and the source of these X-rays is thought to be strong shocks occurring in its magnetically-confined wind (Babel & Montmerle 1997a, Donati et al. 2002). Recently, Stelzer et al. (2005) found significant X-ray emission from all massive stars in the ONC. Periodic rotational modulation in X-rays and other indicators suggested that theta 1 Ori C may be but one of many magnetic B- and O-type stars in this star-forming region. In 2005B we carried out sensitive ESPaDOnS observations to search for direct evidence of such fields, detecting unambiguous Zeeman signatures in two objects.