We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

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.

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.

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.

Although initially thought to be promising for production of the r-process nuclei, standard models of neutrino-heated winds from proto-neutron stars (PNSs) do not reach the requisite neutron-to-seed ratio for production of the lanthanides and actinides. However, the abundance distribution created by the r-, rp-, or $\nu p$-processes in PNS winds depends sensitively on the entropy and dynamical expansion timescale of the flow, which may be strongly affected by high magnetic fields. Here, we present results from magnetohydrodynamic simulations of non-rotating neutrino-heated PNS winds with strong dipole magnetic fields from $10^{14}-10^{16}$ G, and assess their role in altering the conditions for nucleosynthesis. The strong field forms a closed zone and helmet streamer configuration at the equator, with episodic dynamical mass ejections in toroidal plasmoids. We find dramatically enhanced entropy in these regions and conditions favorable for third-peak r-process nucleosynthesis if the wind is neutron-rich. If instead the wind is proton-rich, the conditions will affect the abundances from the $\nu p$-process. We quantify the distribution of ejected matter in entropy and dynamical expansion timescale, and the critical magnetic field strength required to affect the entropy. For $B\gtrsim10^{15}$ G, we find that $\gtrsim10^{-6}$ M$_{\odot}$ and up to $\sim10^{-5}$ M$_{\odot}$ of high entropy material is ejected per highly-magnetized neutron star birth in the wind phase, providing a mechanism for prompt heavy element enrichment of the universe. Former binary companions identified within (magnetar-hosting) supernova remnants, the remnants themselves, and runaway stars may exhibit overabundances. We provide a comparison with a semi-analytic model of plasmoid eruption and discuss implications and extensions.

The investigation of the emergence of life is a major endeavour of science. Astronomy is contributing to it in three fundamental manners: (1) by measuring the chemical enrichment of the Universe, (2) by investigating planet formation and searching for exoplanets with signatures of life and, (3) by determining the abundance of aminoacids and the chemical routes to aminoacid and protein growth in astronomical bodies. This proposal deals with the first two. In the Voyage to 2050, the world-wide scientific community is getting equipped with large facilities for the investigation of the emergence of life in the Universe (i.e. VLT, JWST, ELT, GMT, TMT, ALMA, FAST, VLA, ATHENA, SKA) including the ESA's CHEOPS, PLATO and ARIEL missions. This white paper is a community effort to call for the development of a large ultraviolet optical observatory to gather fundamental data for this investigation that will not be accessible through other ranges of the electromagnetic spectrum. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the pathway to life in the Universe.

Surface magnetic fields have been detected in 5 to 10% of isolated massive stars, hosting outer radiative envelopes. They are often thought to have a fossil origin, resulting from the stellar formation phase. Yet, magnetic massive stars are scarcer in (close) short-period binaries, as reported by the BinaMIcS (Binarity and Magnetic Interaction in various classes of Stars) collaboration. Different physical conditions in the molecular clouds giving birth to isolated stars and binaries are commonly invoked. In addition, we propose that the observed lower magnetic incidence in close binaries may be due to nonlinear tides. Indeed, close binaries are probably prone to tidal instability, a fluid instability growing upon the equilibrium tidal flow via nonlinear effects. Yet, stratified effects have hitherto been largely overlooked. We theoretically and numerically investigate tidal instability in rapidly rotating, stably stratified fluids permeated by magnetic fields. We use the short-wavelength stability method to propose a comprehensive (local) theory of tidal instability at the linear onset, discussing damping effects. Then, we propose a mixing-length theory for the mixing generated by tidal instability in the nonlinear regime. We successfully assess our theoretical predictions against proof-of-concept, direct numerical simulations. Finally, we compare our predictions with the observations of short-period, double-lined spectroscopic binary systems. Using new analytical results, cross-validated by a direct integration of the stability equations, we show that tidal instability can be generated by nonlinear couplings of inertia-gravity waves with the equilibrium tidal flow in short-period massive binaries, even against the Joule diffusion. In the nonlinear regime, a fossil magnetic field can be dissipated by the turbulent magnetic diffusion induced by the saturated tidal flows. We predict that the turbulent Joule diffusion of fossil fields would occur in a few million years for several short-period massive binaries. Therefore, turbulent tidal flows could explain the observed dearth of some short-period magnetic binaries.

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.

We have performed 3D isothermal MHD simulation of a magnetic rotating massive star with a non-zero dipole obliquity and predicted the radio/sub-mm observable lightcurves and continuum spectra for a frequency range compatible with ALMA. From these results we also compare the model input mass-loss to that calculated from the synthetic thermal emission. Spherical and cylindrical symmetry is broken due to the obliquity of the stellar magnetic dipole resulting in an inclination and phase dependence of both the spectral flux and inferred mass-loss rate, providing testable predictions of variability for oblique rotator. Both quantities vary by factors between 2 and 3 over a full rotational period of the star, demonstrating that the role of rotation as critical in understanding the emission. This illustrates the divergence from a symmetric wind, resulting in a two armed spiral structure indicative of a oblique magnetic rotator. We show that a constant spectral index, $\alpha$, model agrees well with our numerical prediction for a spherical wind for $\nu~<~10^{3} \ \mathrm{GHz}$, however it is unable to capture the behavior of emission at $\nu~>~10^{3} \ \mathrm{GHz}$. As such we caution the use of such constant $\alpha$ models for predicting emission from non-spherical winds such as those which form around magnetic massive stars.

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.

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.

Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dipole field. For characteristic stellar and disc parameters of a near-critically rotating B2e star, we find that a polar surface field strength of just 10 G can significantly disrupt the disc, while a field of 100 G, near the observational upper limit inferred for most Be stars, completely destroys the disc over just a few days. Our parameter study shows that the efficacy of this magnetic disruption of a disc scales with the characteristic plasma beta {(defined as the ratio between thermal and magnetic pressure)} in the disc, but is surprisingly insensitive to other variations, e.g. in stellar rotation speed, or the mass loss rate of the star's radiatively driven wind. The disc disruption seen here for even a modest field strength suggests that the presumed formation of such Be discs by decretion of material from the star would likely be strongly inhibited by such fields; this provides an attractive explanation for why no large-scale fields are detected from such Be stars.

Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $\Omega$-slow solution characterized by a dense and slow wind. Here, in this work we study the transition region of the solutions where the fast solution changes to the $\Omega$-slow. Using both time-steady and time-dependent numerical codes, we study this transition region for different equatorial models of B-type stars. In all the cases, at certain range of rotational speeds, we found a region where the fast and $\Omega$-slow solution can co-exist. We find that the type of solution obtained in this co-existence region depends heavily on the initial conditions of our models. We also test the stability of the solutions within the co-existence region by performing base density perturbations in the wind. We find that under certain conditions, the fast solution can switch to a $\Omega$-slow solution, or vice versa. Such solution switching may be a possible contributor of material injected into the circumstellar environment of Be stars, without requiring rotational speeds near critical values.

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.

Context. Red-giant stars may engulf planets. This may increase the rotation rate of their convective envelope, which could lead to strong dynamo-triggered magnetic fields. Aims. We explore the possibility of generating magnetic fields in red giants that have gone through the process of a planet engulfment. We compare them with similar models that evolve without any planets. We discuss the impact of stellar wind magnetic braking on the evolution of the surface velocity of the parent star. Methods. With rotating stellar models with and without planets and an empirical relation between the Rossby number and the surface magnetic field, we deduce the evolution of the surface magnetic field along the red-giant branch. The effects of wind magnetic braking is explored using a relation deduced from MHD simulations. Results. The stellar evolution model of a 1.7 M$_\odot$ without planet engulfment and that has a time-averaged rotation velocity during the Main-Sequence equal to 100 km s$^{-1}$, shows a surface magnetic field triggered by convection larger than 10 G only at the base of the red giant branch, that means for gravities log $g > 3$. When a planet engulfment occurs, such magnetic field can also appear at much lower gravities, i.e. at much higher luminosities along the red giant branch. Typically the engulfment of a 15 M$_J$ planet produces a dynamo triggered magnetic field larger than 10 G for gravities between 2.5 and 1.9. We show that for reasonable wind magnetic braking laws, the high surface velocity reached after a planet engulfment may be maintained sufficiently long for being observable. Conclusions. High surface magnetic fields for red giants in the upper part of the red giant branch is a strong indication of a planet engulfment or of an interaction with a companion. Our theory can be tested by observing fast rotating red giants and check whether they show magnetic fields.

High-resolution spectra of the magnetic star HD191612 were acquired using the Chandra X-ray observatory at both maximum and minimum emission phases. We confirm the flux and hardness variations previously reported with XMM-Newton, demonstrating the great repeatability of the behavior of HD191612 over a decade. The line profiles appear typical for magnetic massive stars: no significant line shift, relatively narrow lines for high-Z elements, and formation radius at about 2R*. Line ratios confirm the softening of the X-ray spectrum at the minimum emission phase. Shift or width variations appear of limited amplitude at most (slightly lower velocity and slightly increased broadening at minimum emission phase, but within 1--2 sigma of values at maximum). In addition, a fully self-consistent 3D magnetohydrodynamic (MHD) simulation of the confined wind in HD191612 was performed. The simulation results were directly fitted to the data leading to a remarkable agreement overall between them.

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.

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.

A subset (~ 10%) of massive stars present strong, globally ordered (mostly dipolar) magnetic fields. The trapping and channeling of their stellar winds in closed magnetic loops leads to magnetically confined wind shocks (MCWS), with pre-shock flow speeds that are some fraction of the wind terminal speed. These shocks generate hot plasma, a source of X-rays. In the last decade, several developments took place, notably the determination of the hot plasma properties for a large sample of objects using XMM-Newton and Chandra, as well as fully self-consistent MHD modelling and the identification of shock retreat effects in weak winds. Despite a few exceptions, the combination of magnetic confinement, shock retreat and rotation effects seems to be able to account for X-ray emission in massive OB stars. Here we review these new observational and theoretical aspects of this X-ray emission and envisage some perspectives for the next generation of X-ray observatories.

The Of?p star CPD -28 2561 was monitored at high energies with XMM-Newton and HST. In X-rays, this magnetic oblique rotator displays bright and hard emission that varies by ~55% with rotational phase. These changes occur in phase with optical variations, as expected for magnetically confined winds; there are two maxima and two minima in X-rays during the 73d rotational period of CPD -28 2561. However, contrary to previously studied cases, no significant hardness variation is detected between minima and maxima, with the exception of the second minimum which is slightly distinct from the first one. In the UV domain, broad-band fluxes remain stable while line profiles display large variations. Stronger absorptions at low velocities are observed when the magnetic equator is seen edge-on, which can be reproduced by a detailed 3D model. However, a difference in absorption at high velocities in the CIV and NV lines is also detected for the two phases where the confined wind is seen nearly pole-on. This suggests the presence of strong asymmetries about the magnetic equator, mostly in the free-flowing wind (rather than in the confined dynamical magnetosphere).

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.

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.

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.

Important insights into the formation, structure, evolution and environment of all types of stars can be obtained through the measurement of their winds and possible magnetospheres. However, this has hardly been done up to now mainly because of the lack of UV instrumentation available for long periods of time. To reach this aim, we have designed UVMag, an M-size space mission equipped with a high-resolution spectropolarimeter working in the UV and visible spectral range. The UV domain is crucial in stellar physics as it is very rich in atomic and molecular lines and contains most of the flux of hot stars. Moreover, covering the UV and visible spectral domains at the same time will allow us to study the star and its environment simultaneously. Adding polarimetric power to the spectrograph will multiply tenfold the capabilities of extracting information on stellar magnetospheres, winds, disks, and magnetic fields. Examples of science objectives that can be reached with UVMag are presented for pre-main sequence, main sequence and evolved stars. They will cast new light onto stellar physics by addressing many exciting and important questions. UVMag is currently undergoing a Research and Technology study and will be proposed at the forthcoming ESA call for M-size missions. This spectropolarimeter could also be installed on a large UV and visible observatory (e.g. NASA's LUVOIR project) within a suite of instruments.

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.

We examine X-rays from radiatively cooled shocks, focusing on how their thin-shell instability reduces X-ray emission. For 2D simulations of collision between equal expanding winds, we carry out a parameter study of such instability as a function of the ratio of radiative vs. adiabatic-expansion cooling lengths. In the adiabatic regime, the extended cooling layer suppresses instability, leading to planar shock compression with X-ray luminosity that follows closely the expected (Lx ~ M^2) quadratic scaling with mass-loss rate M . In the strongly radiative limit, the X-ray emission now follows an expected linear scaling with mass loss (Lx ~ M), but the instability deforms the shock compression into extended shear layers with oblique shocks along fingers of cooled, dense material. The spatial dispersion of shock thermalization limits strong X-ray emission to the tips and troughs of the fingers, and so reduces the X-ray emission (here by about a factor 1/50) below what is expected from analytic radiative-shock models without unstable structure. Between these two limits, X-ray emission can switch between a high-state associated with extended shock compression, and a low-state characterized by extensive shear. Further study is needed to clarify the origin of this "shear mixing reduction factor" in X-ray emission, and its dependence on parameters like the shock Mach number.

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.

Oe stars possibly form an extension to higher temperatures of the Be phenomenon, but it is still unclear whether these stars have disks. X-ray spectra could provide hints for interactions of the star with a putative surrounding disk. We obtained XMM-Newton observations of two Oe stars, HD45314 and HD60848. Spectra and light curves were extracted and analysed. Optical spectra were also obtained to support the X-ray observations. We find that both stars display very different X-ray properties. Whilst HD60848 has an X-ray spectrum and emission level typical for its spectral type, HD45314 displays a very hard X-ray emission, dominated by a thermal plasma with kT ~ 21 keV. Furthermore, HD45314 displays count rate variations by a factor 2 on timescales of ~ 1000s and a high log(L_X/L_bol) = -6.10 +/- 0.03. The X-ray properties of HD45314 indicate that this star is a new member of the class of gamma Cas analogs, the first one among the original category of Oe stars.

The growth of luminous structures and the building blocks of life in the Universe began as primordial gas was processed in stars and mixed at galactic scales. The mechanisms responsible for this development are not well understood and have changed over the intervening 13 billion years. To follow the evolution of matter over cosmic time, it is necessary to study the strongest (resonance) transitions of the most abundant species in the Universe. Most of them are in the ultraviolet (UV; 950A-3000A) spectral range that is unobservable from the ground. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the Universe. Habitable planets grow in protostellar discs under ultraviolet irradiation, a by-product of the star-disk interaction that drives the physical and chemical evolution of discs and young planetary systems. The electronic transitions of the most abundant molecules are pumped by the UV field, providing unique diagnostics of the planet-forming environment that cannot be accessed from the ground. Earth's atmosphere is in constant interaction with the interplanetary medium and the solar UV radiation field. A 50-100 times improvement in sensitivity would enable the observation of the key atmospheric ingredients of Earth-like exoplanets (carbon, oxygen, ozone), provide crucial input for models of biologically active worlds outside the solar system, and provide the phenomenological baseline to understand the Earth atmosphere in context. In this white paper, we outline the key science that such a facility would make possible and outline the instrumentation to be implemented.

Stars over a wide range of masses and evolutionary stages are nowadays known to emit X-rays. This X-ray emission is a unique probe of the most energetic phenomena occurring in the circumstellar environment of these stars, and provides precious insight on magnetic phenomena or hydrodynamic shocks. Owing to its large collecting area, Athena+ will open up an entirely new window on these phenomena. Indeed, Athena+ will not only allow us to study many more objects with an unprecedented spectral resolution, but will also pioneer the study of the dynamics of these objects via time-resolved high-resolution spectroscopy. In this way, Athena+ will be a unique tool to study accretion processes in TTauri stars, flaring activity in young stars, dynamos in ultra-cool dwarfs, small and large-scale structures in the winds of single massive stars, wind interactions in massive binary systems, hot bubbles in planetary nebula... All these studies will lead to a deeper understanding of yet poorly understood processes which have profound impact in star and planetary system formation as well as in feedback processes on Galactic scale.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

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

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.

The O9IV star HD 57682, discovered to be magnetic within the context of the MiMeS survey in 2009, is one of only eight convincingly detected magnetic O-type stars. Among this select group, it stands out due to its sharp-lined photospheric spectrum. Since its discovery, the MiMeS Collaboration has continued to obtain spectroscopic and magnetic observations in order to refine our knowledge of its magnetic field strength and geometry, rotational period, and spectral properties and variability. In this paper we report new ESPaDOnS spectropolarimetric observations of HD 57682, which are combined with previously published ESPaDOnS data and archival H{\alpha} spectroscopy. This dataset is used to determine the rotational period (63.5708 \pm 0.0057 d), refine the longitudinal magnetic field variation and magnetic geometry (dipole surface field strength of 880\pm50 G and magnetic obliquity of 79\pm4\circ as measured from the magnetic longitudinal field variations, assuming an inclination of 60\circ), and examine the phase variation of various lines. In particular, we demonstrate that the H{\alpha} equivalent width undergoes a double-wave variation during a single rotation of the star, consistent with the derived magnetic geometry. We group the variable lines into two classes: those that, like H{\alpha}, exhibit non-sinusoidal variability, often with multiple maxima during the rotation cycle, and those that vary essentially sinusoidally. Based on our modelling of the H{\alpha} emission, we show that the variability is consistent with emission being generated from an optically thick, flattened distribution of magnetically-confined plasma that is roughly distributed about the magnetic equator. Finally, we discuss our findings in the magnetospheric framework proposed in our earlier study.

Due to computational requirements and numerical difficulties associated with coordinate singularity in spherical geometry, fully dynamic 3D magnetohydrodynamic (MHD) simulations of massive star winds are not readily available. Here we report preliminary results of the first such a 3D simulation using $\theta^1$ Ori C (O5.5 V) as a model. The oblique magnetic rotator $\theta^1$ Ori C is a source of hard X-ray emitting plasma in its circumstellar environment. Our numerical model can explain both the hardness and the location of the X-ray emission from this star confirming that magnetically confined wind shock (MCWS) is the dominating mechanism for hard Xrays in some massive stars.

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

The key empirical property of the X-ray emission from O stars is a strong correlation between the bolometric and X-ray luminosities. In the framework of the Chandra Carina Complex Project, 129 O and B stars have been detected as X-ray sources; 78 of those, all with spectral type earlier than B3, have enough counts for at least a rough X-ray spectral characterization. This leads to an estimate of the Lx/Lbol ratio for an exceptional number of 60 O stars belonging to the same region and triples the number of Carina massive stars studied spectroscopically in X-rays. The derived log(Lx/Lbol) is -7.26 for single objects, with a dispersion of only 0.21dex. Using the properties of hot massive stars listed in the literature, we compare the X-ray luminosities of different types of objects. In the case of O stars, the Lx/Lbol ratios are similar for bright and faint objects, as well as for stars of different luminosity classes or spectral types. Binaries appear only slightly harder and slightly more luminous in X-rays than single objects; the differences are not formally significant (at the 1% level), except for the Lx/Lbol ratio in the medium (1.0--2.5keV) energy band. Weak-wind objects have similar X-ray luminosities but they display slightly softer spectra compared to "normal" O stars with the same bolometric luminosity. Discarding three overluminous objects, we find a very shallow trend of harder emission in brighter objects. The properties of the few B stars bright enough to yield some spectral information appear to be different overall (constant X-ray luminosities, harder spectra), hinting that another mechanism for producing X-rays, besides wind shocks, might be at work. However, it must be stressed that the earliest and X-ray brightest amongst these few detected objects are similar to the latest O stars, suggesting a possibly smooth transition between the two processes.

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

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.

We obtained relatively high signal-to-noise X-ray spectral data of the early massive star tau Sco (B0.2V) with the Suzaku XIS instrument. This source displays several unusual features that motivated our study: (a) redshifted absorption in UV P Cygni lines to approximately +250 km/s suggestive of infalling gas, (b) unusually hard X-ray emission requiring hot plasma at temperatures in excess of 10 MK whereas most massive stars show relatively soft X-rays at a few MK, and (c) a complex photospheric magnetic field of open and closed field lines. In an attempt to understand the hard component better, X-ray data were obtained at six roughly equally spaced phases within the same epoch of tau Sco's 41 day rotation period. The XIS instrument has three operable detectors: XIS1 is back illuminated with sensitivity down to 0.2 keV; XIS0 and XIS2 are front illuminated with sensivitity only down to 0.4 keV and have overall less effective area than XIS1. The XIS0 and XIS3 detectors show relatively little variability. In contrast, there is a 4sigma detection of a 4% drop in the count rate of the XIS1 detector at one rotational phase. In addition, all three detectors show a 3% increase in count rate at the same phase. The most optimistic prediction of X-ray variability allows for a 40% change in the count rate, particularly near phases where we have pointings. Observed modulations in the X-ray light curve on the rotation cycle is an order of magnitude smaller than this, which places new stringent constraints on future modeling of this interesting magnetic massive star.

In recent years several in-depth investigations of the three Galactic Of?p stars were undertaken. These multiwavelength studies revealed the peculiar properties of these objects (in the X-rays as well as in the optical): magnetic fields, periodic line profile variations, recurrent photometric changes. However, many questions remain unsolved. To clarify some of the properties of the Of?p stars, we have continued their monitoring. A new XMM observation and two new optical datasets were obtained. Additional information for the prototypical Of?p trio has been found. HD108 has now reached its quiescent, minimum-emission state, for the first time in 50--60yrs. The echelle spectra of HD148937 confirm the presence of the 7d variations in the Balmer lines and reveal similar periodic variations (though of lower amplitudes) in the HeI5876 and HeII4686 lines, underlining its similarities with the other two prototypical Of?p stars. The new XMM observation of HD191612 was taken at the same phase in the line modulation cycle but at a different orbital phase as previous data. It clearly shows that the X-ray emission of HD191612 is modulated by the 538d period and not the orbital period of 1542d - it is thus not of colliding-wind origin and the phenomenon responsible for the optical changes appears also at work in the high-energy domain. There are however problems: our MHD simulations of the wind magnetic confinement predict both a harder X-ray flux of a much larger strength than what is observed (the modeled DEM peaks at 30-40MK, whereas the observed one peaks at 2MK) and narrow lines (hot gas moving with velocities of 100--200km/s, whereras the observed FWHM is ~2000km/s).

We present new U-band photometry of the magnetic Helium-strong star Sigma Ori E, obtained over 2004-2009 using the SMARTS 0.9-m telescope at Cerro Tololo Inter-American Observatory. When combined with historical measurements, these data constrain the evolution of the star's 1.19 d rotation period over the past three decades. We are able to rule out a constant period at the p_null = 0.05% level, and instead find that the data are well described (p_null = 99.3%) by a period increasing linearly at a rate of 77 ms per year. This corresponds to a characteristic spin-down time of 1.34 Myr, in good agreement with theoretical predictions based on magnetohydrodynamical simulations of angular momentum loss from magnetic massive stars. We therefore conclude that the observations are consistent with Sigma Ori E undergoing rotational braking due to its magnetized line-driven wind.

High-energy emissions are good indicators of peculiar behaviours in stars. We have therefore obtained an XMM-Newton observation of HD155806 and 1RXSJ171502.4-333344, and derived their spectral properties for the first time. The X-ray spectrum of HD155806 appears soft, even slightly softer than usual for O-type stars (as shown by a comparison with the O9 star HD155889 in the same XMM field). It is well-fitted with a two-component thermal model with low temperatures (0.2 and 0.6 keV), and it shows no overluminosity (log[LX/Lbol]=-6.75). The high-resolution spectrum, though noisy, reveals a few broad, symmetric X-ray lines (FWHM ~ 2500 km/s). The X-ray emission is compatible with the wind-shock model and therefore appears unaffected by the putative dense equatorial regions at the origin of the Oe classification. 1RXSJ171502.4-333344 is a nearby flaring source of moderate X-ray luminosity (log[LX/Lbol]=-3), with a soft thermal spectrum composed of narrow lines and presenting a larger abundance of elements (e.g. Ne) with a high first ionization potential (FIP) compared to lower-FIP elements. All the evidence indicates a coronal origin for the X-ray emission, in agreement with the dMe classification of this source.

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.

We examine the angular momentum loss and associated rotational spin-down for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, W(=Vrot/Vorb) and 'wind magnetic confinement', $\eta_\ast$ defined below. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss averaged over both solid angle and time follows closely the general WD scaling $\dot {J} \sim \dot {M} \Omega R_A^2$. The key distinction is that for a dipole field Alfv\`en radius $R_A$ is significantly smaller than for the monopole field WD used in their analyses. This leads to a slower stellar spin-down for the dipole field with typical spin-down times of order 1 Myr for several known magnetic massive stars.

We examine the angular momentum loss and associated rotational spindown for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, one characterizing the wind magnetic confinement ($\eta_{\ast} \equiv B_{eq}^{2} R_{\ast}^{2}/{\dot M} v_{\infty}$), and the other the ratio ($W \equiv V_{rot}/V_{orb}$) of stellar rotation to critical (orbital) speed. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss ${\dot J}$ averaged over both solid angle and time follows closely the general WD scaling ${\dot J} = (2/3) {\dot M} \Omega R_{A}^{2}$, where ${\dot M}$ is the mass loss rate, $\Omega$ is the stellar angular velocity, and $R_{A}$ is a characteristic Alfv\'{e}n radius. However, a key distinction here is that for a dipole field, this Alfv\'{e}n radius has a strong-field scaling $R_{A}/R_{\ast} \approx \eta_{\ast}^{1/4}$, instead of the scaling $R_{A}/R_{\ast} \sim \sqrt{\eta_{\ast}}$ for a monopole field. This leads to a slower stellar spindown time that in the dipole case scales as $\tau_{spin} = \tau_{mass} 1.5k/\sqrt{\eta_{\ast}}$, where $\tau_{mass} \equiv M/{\dot M}$ is the characteristic mass loss time, and $k$ is the dimensionless factor for stellar moment of inertia. The full numerical scaling relation we cite gives typical spindown times of order 1 Myr for several known magnetic massive stars.

Building upon our previous MHD simulation study of magnetic channeling in radiatively driven stellar winds, we examine here the additional dynamical effects of stellar {\em rotation} in the (still) 2-D axisymmetric case of an aligned dipole surface field. In addition to the magnetic confinement parameter $\eta_{\ast}$ introduced in Paper I, we characterize the stellar rotation in terms of a parameter $W \equiv V_{\rm{rot}}/V_{\rm{orb}}$ (the ratio of the equatorial surface rotation speed to orbital speed), examining specifically models with moderately strong rotation $W =$ 0.25 and 0.5, and comparing these to analogous non-rotating cases. Defining the associated Alfv\'{e}n radius $R_{\rm{A}} \approx \eta_{\ast}^{1/4} \Rstar$ and Kepler corotation radius $R_{\rm{K}} \approx W^{-2/3} \Rstar$, we find rotation effects are weak for models with $R_{\rm{A}} < R_{\rm{K}}$, but can be substantial and even dominant for models with $R_{\rm{A}} \gtwig R_{\rm{K}}$. In particular, by extending our simulations to magnetic confinement parameters (up to $\eta_{\ast} = 1000$) that are well above those ($\eta_{\ast} = 10$) considered in Paper I, we are able to study cases with $R_{\rm{A}} \gg R_{\rm{K}}$; we find that these do indeed show clear formation of the {\em rigid-body} disk predicted in previous analytic models, with however a rather complex, dynamic behavior characterized by both episodes of downward infall and outward breakout that limit the buildup of disk mass. Overall, the results provide an intriguing glimpse into the complex interplay between rotation and magnetic confinement, and form the basis for a full MHD description of the rigid-body disks expected in strongly magnetic Bp stars like $\sigma$ Ori E.

We introduce a new Rigid-Field Hydrodynamics approach to modeling the magnetospheres of massive stars in the limit of very-strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the 1-dimensional flow along each, subject to pressure, radiative, gravitational, and centrifugal forces. We solve these equations numerically for a large ensemble of field lines, to build up a 3-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star sigma Ori E. Since the flow along each field line can be solved for independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disk defined by the locus of minima of the effective (gravitational plus centrifugal) potential. But a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ~3,000 km/s) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disk can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

We present results of numerical simulations of wind variability in the magnetic B1 IVe star $\beta$ Cephei. 2D-MHD simulations are used to determine the structure of the wind. From these wind models we calculate line profiles for different aspect angles to simulate rotation. The results are compared with the observed UV wind line profiles.

We present 2D MHD simulations of the radiatively driven outflow from a rotating hot star with a dipole magnetic field aligned with the star's rotation axis. We focus primarily on a model with moderately rapid rotation (half the critical value), and also a large magnetic confinement parameter, $\eta_{\ast} \equiv B_{\ast}^2 R_{\ast}^{2} / \dot{M} V_{\infty} = 600$. The magnetic field channels and torques the wind outflow into an equatorial, rigidly rotating disk extending from near the Kepler corotation radius outwards. Even with fine-tuning at lower magnetic confinement, none of the MHD models produce a stable Keplerian disk. Instead, material below the Kepler radius falls back on to the stellar surface, while the strong centrifugal force on material beyond the corotation escape radius stretches the magnetic loops outwards, leading to episodic breakout of mass when the field reconnects. The associated dissipation of magnetic energy heats material to temperatures of nearly $10^{8}$K, high enough to emit hard (several keV) X-rays. Such \emph{centrifugal mass ejection} represents a novel mechanism for driving magnetic reconnection, and seems a very promising basis for modeling X-ray flares recently observed in rotating magnetic Bp stars like $\sigma$ Ori E.

We report on four Chandra grating observations of the oblique magnetic rotator theta^1 Ori C (O5.5 V) covering a wide range of viewing angles with respect to the star's 1060 G dipole magnetic field. We employ line-width and centroid analyses to study the dynamics of the X-ray emitting plasma in the circumstellar environment, as well as line-ratio diagnostics to constrain the spatial location, and global spectral modeling to constrain the temperature distribution and abundances of the very hot plasma. We investigate these diagnostics as a function of viewing angle and analyze them in conjunction with new MHD simulations of the magnetically channeled wind shock mechanism on theta^1 Ori C. This model fits all the data surprisingly well, predicting the temperature, luminosity, and occultation of the X-ray emitting plasma with rotation phase.

This talk summarizes results from recent MHD simulations of the role of a dipole magnetic field in inducing large-scale structure in the line-driven stellar winds of hot, luminous stars. Unlike previous fixed-field analyses, the MHD simulations here take full account of the dynamical competition between the field and the flow. A key result is that the overall degree to which the wind is influenced by the field depends largely on a single, dimensionless `wind magnetic confinement parameter', $\eta_\ast (= B_{eq}^2 R_{\ast}^2/\dot{M} v_\infty$), which characterizes the ratio between magnetic field energy density and kinetic energy density of the wind. For weak confinement, $\eta_\ast \le 1$, the field is fully opened by wind outflow, but nonetheless, for confinement as small as $\eta_\ast=1/10$ it can have significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement, $\eta_\ast > 1$, the magnetic field remains closed over limited range of latitude and height above the equatorial surface, but eventually is opened into nearly radial configuration at large radii. Within the closed loops, the flow is channeled toward loop tops into shock collisions that are strong enough to produce hard X-rays. Within the open field region, the equatorial channeling leads to oblique shocks that are again strong enough to produce X-rays and also lead to a thin, dense, slowly outflowing ``disk'' at the magnetic equator.

We summarize recent 2D MHD simulations of line-driven stellar winds from rotating hot-stars with a dipole magnetic field aligned to the star's rotation axis. For moderate to strong fields, much wind outflow is initially along closed magnetic loops that nearly corotate as a solid body with the underlying star, thus providing a torque that results in an effective angular momentum spin-up of the outflowing material. But instead of forming the ``magnetically torqued disk'' (MTD) postulated in previous phenemenological analyses, the dynamical simulations here show that material trapped near the tops of such closed loops tends either to fall back or break out, depending on whether it is below or above the Keplerian corotation radius. Overall the results raise serious questions about whether magnetic torquing of a wind outflow could naturally result in a Keplerian circumstellar disk. However, for very strong fields, it does still seem possible to form a %``magnetically confined, centrifugally supported, rigid-body disk'', centrifugally supported, ``magnetically rigid disk'' (MRD), in which the field not only forces material to maintain a rigid-body rotation, but for some extended period also holds it down against the outward centrifugal force at the loop tops. We argue that such rigid-body disks seem ill-suited to explain the disk emission from Be stars, but could provide a quite attractive paradigm for circumstellar emission from the magnetically strong Bp and Ap stars.

We carry out an extended analytic study of how the tilt and faster-than-radial expansion from a magnetic field affect the mass flux and flow speed of a line-driven stellar wind. A key motivation is to reconcile results of numerical MHD simulations with previous analyses that had predicted non-spherical expansion would lead to a strong speed enhancement. By including finite-disk correction effects, a dynamically more consistent form for the non-spherical expansion, and a moderate value of the line-driving power index $\alpha$, we infer more modest speed enhancements that are in good quantitative agreement with MHD simulations, and also are more consistent with observational results. Our analysis also explains simulation results that show the latitudinal variation of the surface mass flux scales with the square of the cosine of the local tilt angle between the magnetic field and the radial direction. Finally, we present a perturbation analysis of the effects of a finite gas pressure on the wind mass loss rate and flow speed in both spherical and magnetic wind models, showing that these scale with the ratio of the sound speed to surface escape speed, $a/v_{esc}$, and are typically 10-20% compared to an idealized, zero-gas-pressure model.

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.

We present numerical magnetohydrodynamic (MHD) simulations of the effect of stellar dipole magnetic fields on line-driven wind outflows from hot, luminous stars. Unlike previous fixed-field analyses, the simulations here take full account of the dynamical competition between field and flow, and thus apply to a full range of magnetic field strength, and within both closed and open magnetic topologies. A key result is that the overall degree to which the wind is influenced by the field depends largely on a single, dimensionless, `wind magnetic confinement parameter', $\eta_{\ast}$ ($ = B_{eq}^2 R_\ast^2/{\dot M} v_\infty$), which characterizes the ratio between magnetic field energy density and kinetic energy density of the wind. For weak confinement $\eta_{\ast} \le 1$, the field is fully opened by the wind outflow, but nonetheless for confinements as small as $\eta_{\ast}=1/10$ can have a significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement $\eta_{\ast} > 1$, the magnetic field remains closed over a limited range of latitude and height about the equatorial surface, but eventually is opened into a nearly radial configuration at large radii.

High-resolution X-ray spectra of high-mass stars and low-mass T-Tauri stars obtained during the first year of the Chandra mission are providing important clues about the mechanisms which produce X-rays on very young stars. For zeta Puppis (O4 If) and zeta Ori (O9.5 I), the broad, blue-shifted line profiles, line ratios, and derived temperature distribution suggest that the X-rays are produced throughout the wind via instability-driven wind shocks. For some less luminous OB stars, like theta^1 Ori C (O7 V) and tau Sco (B0 V), the line profiles are symmetric and narrower. The presence of time-variable emission and very high-temperature lines in theta^1 Ori C and tau Sco suggest that magnetically confined wind shocks may be at work. The grating spectrum of the classical T-Tauri star TW Hya is remarkable because the forbidden-line emission of He-like Ne IX and O VII is very weak, implying that the X-ray emitting region is very dense, n = 6E+12 cgs, or that the X-rays are produced very close to the ultraviolet hotspot at the base of an accretion funnel. ACIS light curves and spectra of flares and low-mass and high-mass young stellar objects in Orion and rho Ophiuchus further suggest that extreme magnetic activity is a general property of many very young stars.