Cyg OB2 #8a is a massive O-type binary displaying strong non-thermal radio emission. Owing to the compactness of this binary, emission of non-thermal X-ray photons via inverse Compton scattering is expected. We first revised the orbital solution for Cyg OB2 #8a using new optical spectra. We then reduced and analysed X-ray spectra obtained with XMM-Newton, Swift, INTEGRAL, and NuSTAR. The analysis of the XMM-Newton and Swift data allows us to better characterise the X-ray emission from the stellar winds and colliding winds region at energies below 10 keV. We confirm the variation of the broad-band light curve of Cyg OB2 #8a along the orbit with, for the first time, the observation of the maximum emission around phase 0.8. The minimum ratio of the X-ray to bolometric flux of Cyg OB2 #8a remains well above the level expected for single O-type stars, indicating that the colliding wind region is not disrupted during the periastron passage. The analysis of the full set of publicly available INTEGRAL observations allows us to refine the upper limit on the non-thermal X-ray flux of the Cyg OB2 region between 20 and 200 keV. Two NuSTAR observations (phases 0.028 and 0.085) allow us to study the Cyg OB2 #8a spectrum up to 30 keV. These data do not provide evidence of the presence of non-thermal X-rays, but bring more stringent constraints on the flux of a putative non-thermal component. Finally, we computed, thanks to a new dedicated model, the anisotropic inverse Compton emission generated in the wind shock region. The theoretical non-thermal emission appears to be compatible with observational limits and the kinetic luminosity computed from these models is in good agreement with the unabsorbed flux observed below 10 keV.

We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions.

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.

In this paper, we investigate the photometric variability of magnetic O-type stars. Such stars possess oblique, predominantly dipolar magnetic fields that confine their winds roughly axisymmetrically about the magnetic equator, thus forming a magnetosphere. We interpret their photometric variability as phase-dependent magnetospheric occultations. For massive star winds dominated by electron scattering opacity in the optical and NIR, we can compute synthetic light curves from simply knowing the magnetosphere's mass density distribution. We exploit the newly-developed Analytical Dynamical Magnetosphere model (ADM) in order to obtain the predicted circumstellar density structures of magnetic O-type stars. The simplicity in our light curve synthesis model allows us to readily conduct a parameter space study. For validation purposes, we first apply our algorithm to HD 191612, the prototypical Of?p star. Next, we attempt to model the photometric variability of the Of?p-type stars identified in the Magellanic Clouds using OGLE photometry. We evaluate the compatibility of the ADM predictions with the observed photometric variations, and discuss the magnetic field properties that are implied by our modelling.

The Cyg OB2 #5 system is thought to consist of a short-period (6.6 d) eclipsing massive binary orbited by an OB-star orbiting with a period of ~6.7 yr; these stars in turn are orbited by a distant early B-star with a period of thousands of years. However, while the inner binary has been studied many times, information is missing on the other stars, in particular the third star whose presence was indirectly postulated from recurrent modulations in the radio domain. Besides, to this date, the X-ray light curve could not be fully interpreted, for example in the framework of colliding-wind emission linked to one of the systems. We obtained new optical and X-ray observations of Cyg OB2 #5, which we combined to archival data. We performed a thorough and homogeneous investigation of all available data, notably revisiting the times of primary minimum in photometry. In the X-ray domain, XMM-Newton provides scattered exposures over ~5000 d whilst Swift provides a nearly continuous monitoring for the last couple of years. Although the X-ray light curve reveals clear variability, no significant period can be found hence the high-energy emission cannot be explained solely in terms of colliding winds varying along either the short or intermediate orbits. The optical data reveal for the first time clear signs of reflex motion. The photometry indicates the presence of a 2366 d (i.e. 6.5 yr) period while the associated radial velocity changes are detected at the 3 sigma level in the systemic velocity of the He II 4686 emission line. With the revised period, the radio light curve is interpreted consistently in terms of a wind interaction between the inner binary and the tertiary star. From these optical and radio data, we derive constraints on the physical properties of the tertiary star and its orbit.

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

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.

HD149404 is an evolved non-eclipsing O-star binary that has previously undergone a Roche lobe overflow interaction. Understanding some key properties of the system requires a determination of the orbital inclination and of the dimensions of the components. The BRITE-Heweliusz satellite was used to collect photometric data of HD149404. Additional photometry was retrieved from the SMEI archive. These data were analysed using a suite of period search tools. The orbital part of the lightcurve was modelled with the nightfall binary star code. The Gaia-DR2 parallax of HD149404 was used to provide additional constraints. The periodograms reveal a clear orbital modulation of the lightcurve with a peak-to-peak amplitude near 0.04 mag. The remaining non-orbital part of the variability is consistent with red noise. The lightcurve folded with the orbital period reveals ellipsoidal variations, but no eclipses. The minimum when the secondary star is in inferior conjunction is deeper than the other minimum due to mutual reflection effects between the stars. Combined with the Gaia-DR2 parallaxes, the photometric data indicate an orbital inclination in the range of 23{\deg} to 31{\deg} and a Roche lobe filling factor of the secondary larger than or equal to 0.96. The luminosity of the primary star is consistent with its present-day mass, whereas the more evolved secondary appears overluminous for its mass. We confirm that the primary's rotation period is about half the orbital period. Both features most probably stem from the past Roche lobe overflow episode.

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of the ESA Athena X-ray observatory. Over a field of view of 5' equivalent diameter, it will deliver X-ray spectra from 0.2 to 12 keV with a spectral resolution of 2.5 eV up to 7 keV on ~5 arcsecond pixels. The X-IFU is based on a large format array of super-conducting molybdenum-gold Transition Edge Sensors cooled at about 90 mK, each coupled with an absorber made of gold and bismuth with a pitch of 249 microns. A cryogenic anti-coincidence detector located underneath the prime TES array enables the non X-ray background to be reduced. A bath temperature of about 50 mK is obtained by a series of mechanical coolers combining 15K Pulse Tubes, 4K and 2K Joule-Thomson coolers which pre-cool a sub Kelvin cooler made of a 3He sorption cooler coupled with an Adiabatic Demagnetization Refrigerator. Frequency domain multiplexing enables to read out 40 pixels in one single channel. A photon interacting with an absorber leads to a current pulse, amplified by the readout electronics and whose shape is reconstructed on board to recover its energy with high accuracy. The defocusing capability offered by the Athena movable mirror assembly enables the X-IFU to observe the brightest X-ray sources of the sky (up to Crab-like intensities) by spreading the telescope point spread function over hundreds of pixels. Thus the X-IFU delivers low pile-up, high throughput (>50%), and typically 10 eV spectral resolution at 1 Crab intensities, i.e. a factor of 10 or more better than Silicon based X-ray detectors. In this paper, the current X-IFU baseline is presented, together with an assessment of its anticipated performance in terms of spectral resolution, background, and count rate capability. The X-IFU baseline configuration will be subject to a preliminary requirement review that is scheduled at the end of 2018.

A growing number of Be and Oe stars, named the gamma Cas stars, are known for their unusually hard and intense X-ray emission. This emission could either trace accretion by a compact companion or magnetic interaction between the star and its decretion disk. To test these scenarios, we carried out a detailed optical monitoring of HD45314, the hottest member of the class of gamma Cas stars, along with dedicated X-ray observations on specific dates. High-resolution optical spectra were taken to monitor the emission lines formed in the disk, while X-ray spectroscopy was obtained at epochs when the optical spectrum of the Oe star was displaying peculiar properties. Over the last four years, HD45314 has entered a phase of spectacular variations. The optical emission lines have undergone important morphology and intensity changes including transitions between single- and multiple-peaked emission lines as well as shell events, and phases of (partial) disk dissipation. Photometric variations are found to be anti-correlated with the equivalent width of the H-alpha emission. Whilst the star preserved its hard and bright X-ray emission during the shell phase, the X-ray spectrum during the phase of (partial) disk dissipation was significantly softer and weaker. The observed behaviour of HD45314 suggests a direct association between the level of X-ray emission and the amount of material simultaneously present in the Oe disk as expected in the magnetic star-disk interaction scenario.

All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citep{Naze} and spectroscopic \citep{Walborn} variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing for the synthesis of photometric observables based on the Analytic Dynamical Magnetosphere (ADM) model of \citet{Owocki}. We apply our model to OGLE photometry in order to constrain their magnetic geometries and surface dipole strengths. We predict that the field strengths for some of these candidate extra-Galactic magnetic stars may be within the detection limits of the FORS2 instrument.

Aims: Past observations of fast-rotating massive stars exhibiting normal nitrogen abundances at their surface have raised questions about the rotational mixing paradigm. We revisit this question thanks to a spectroscopic analysis of a sample of bright fast-rotating OB stars, with the goal of quantifying the efficiency of rotational mixing at high rotation rates. Method: Our sample consists of 40 fast rotators on the main sequence, with spectral types comprised between B0.5 and O4. We compare the abundances of some key element indicators of mixing (He, CNO) with the predictions of evolutionary models for single objects and for stars in interacting binary systems. Results: The properties of half of the sample stars can be reproduced by single evolutionary models, even in the case of probable or confirmed binaries that can therefore be true single stars in a pre-interaction configuration. The main problem for the rest of the sample is a mismatch for the [N/O] abundance ratio (we confirm the existence of fast rotators with a lack of nitrogen enrichment) and/or a high helium abundance that cannot be accounted for by models. Modifying the diffusion coefficient implemented in single- star models does not solve the problem as it cannot simultaneously reproduce the helium abundances and [N/O] abundance ratios of our targets. Since part of them actually are binaries, we also compared their chemical properties with predictions for post-mass transfer systems. We found that these models can explain the abundances measured for a majority of our targets, including some of the most helium-enriched, but fail to reproduce them in other cases. Our study thus reveals that some physical ingredients are still missing in current models.

Aims. The magnetic star HD148937 is the only Galactic Of?p star surrounded by a nebula. The structure of this nebula is particularly complex and is composed, from the center out outwards, of a close bipolar ejecta nebula (NGC6164/5), an ellipsoidal wind-blown shell, and a spherically symmetric Stromgren sphere. The exact formation process of this nebula and its precise relation to the star's evolution remain unknown. Methods. We analyzed infrared Spitzer IRS and far-infrared Herschel/PACS observations of the NGC6164/5 nebula. The Herschel imaging allowed us to constrain the global morphology of the nebula. We also combined the infrared spectra with optical spectra of the central star to constrain its evolutionary status. We used these data to derive the abundances in the ejected material. To relate this information to the evolutionary status of the star, we also determined the fundamental parameters of HD148937 using the CMFGEN atmosphere code. Results. The H$\alpha$ image displays a bipolar or "8"-shaped ionized nebula, whilst the infrared images show dust to be more concentrated around the central object. We determine nebular abundance ratios of N/O = 1.06 close to the star, and N/O = 1.54 in the bright lobe constituting NGC6164. Interestingly, the parts of the nebula located further from HD148937 appear more enriched in stellar material than the part located closer to the star. Evolutionary tracks suggest that these ejecta have occured $\sim$1.2-1.3 and $\sim$0.6 Myrs ago, respectively. In addition, we derive abundances of argon for the nebula compatible with the solar values and we find a depletion of neon and sulfur. The combined analyses of the known kinematics and of the new abundances of the nebula suggest either a helical morphology for the nebula, possibly linked to the magnetic geometry, or the occurrence of a binary merger.

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.

Mass and momentum exchanges in close massive binaries play an important role in their evolution, and produce several observational signatures such as asynchronous rotation and altered chemical compositions, that remain after the stars detach again. We investigated these effects for the detached massive O-star binary HD 149404 (O7.5If + ON9.7I, P = 9.81 days), which is thought to have experienced a past episode of case A Roche-lobe overflow (RLOF). Using phase-resolved spectroscopy, we performed the disentangling of the optical spectra of the two stars. The reconstructed primary and secondary spectra were then analysed with the CMFGEN model atmosphere code to determine stellar parameters, such as the effective temperatures and surface gravities, and to constrain the chemical composition of the components. We complemented the optical study with the study of IUE spectra, which we compare to the synthetic binary spectra. The properties of the stars were compared to evolutionary models. We confirmed a strong overabundance in nitrogen ($\left[\rm N/C\right]\sim150\left[\rm N/C\right]_{\odot}$) for the secondary and a slight nitrogen overabundance ($\left[\rm N/C\right]\sim5\left[\rm N/C\right]_{\odot}$) for the primary star. Comparing the two stars, we found evidence for asynchronous rotation, with a rotational period ratio of $0.50 \pm 0.11$. The hypothesis of a past case A RLOF interaction in HD 149404 is most plausible to explain its chemical abundances and rotational asynchronicity. Some of the observed properties, such as the abundance pattern, are clearly a challenge for current case A binary evolution models, however.

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.

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phase=0.0 when the secondary $\delta$ Orionis Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

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.

Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system $\delta$ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the ${\rm \it Hipparcos}$ parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if $\delta$ Ori lies at about twice the ${\rm \it Hipparcos}$ distance, in the vicinity of the $\sigma$-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be $\Delta V \approx 2.\!\!^{\rm m}8$. The inferred parameters suggest the secondary is an early B-type dwarf ($\approx$ B1 V), while the tertiary is an early B-type subgiant ($\approx$ B0 IV). We find evidence for rapid turbulent velocities ($\sim 200$ km ${\rm s}^{-1}$) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind ($\log L_{\text{X}} / L_{\text{Bol}} \approx -6.85$), initiating close to the stellar surface at $R_0 \sim 1.1\,R_*$. Accounting for clumping, the mass-loss rate of the primary is found to be $\log \dot{M} \approx -6.4\,[M_\odot\,{\rm yr}^{-1}]$, which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio spectral domains.

X-ray and radio data recently acquired as part of a project to study Cyg OB2#9 are used to constrain physical models of the binary system, providing in-depth knowledge about the wind-wind collision and the thermal, and non-thermal, emission arising from the shocks. We use a three-dimensional, adaptive mesh refinement simulation (including wind acceleration, radiative cooling, and the orbital motion of the stars) to model the gas dynamics of the wind-wind collision. The simulation output is used as the basis for radiative transfer calculations considering the thermal X-ray emission and the thermal/non-thermal radio emission. To obtain good agreement with the X-ray observations, our initial mass-loss rate estimates require a down-shift by a factor of roughly 7.7 to $6.5\times10^{-7}$ and $7.5\times10^{-7}$ solar mass per year for the primary and secondary star, respectively. Furthermore, the low gas densities and high shock velocities in Cyg OB2#9 are suggestive of unequal electron and ion temperatures, and the X-ray analysis indicates that an (immediately post-shock) electron-ion temperature ratio of $\simeq 0.1$ is also required. The radio emission is dominated by (non-thermal) synchrotron emission. A parameter space exploration provides evidence against models assuming equipartition between magnetic and relativistic energy densities. However, fits of comparable quality can be attained with models having stark contrasts in the ratio of magnetic-to-relativistic energy densities. The radio models also reveal a subtle effect whereby inverse Compton cooling leads to an increase in emissivity as a result of the synchrotron characteristic frequency being significantly reduced. Finally, using the results of the radio analysis, we estimate the surface magnetic field strengths to be $\approx 0.3-52\;$G. (Abridged)

We obtained deep $UBVRI$ H$\alpha$ photometric data of the field situated to the west of the main Carina nebula and centered on WR~22. Medium-resolution optical spectroscopy of a subsample of X-ray selected objects along with archival data sets from $Chandra$, $XMM-Newton$ and 2MASS surveys were used for the present study. Different sets of color-color and color-magnitude diagrams are used to determine reddening for the region and to identify young stellar objects (YSOs) and estimate their age and mass. Our spectroscopic results indicate that the majority of the X-ray sources are late spectral type stars. The region shows a large amount of differential reddening with minimum and maximum values of $E(B-V)$ as 0.25 and 1.1 mag, respectively. Our analysis reveals that the total-to-selective absorption ratio $R \rm_{V}$ is $\sim$3.7 $\pm$ 0.1, suggesting an abnormal grain size in the observed region. We identified 467 YSOs and studied their characteristics. The ages and masses of the 241 optically identified YSOs range from $\sim$0.1 to 10 Myr and $\sim$0.3 to 4.8M$_\odot$, respectively. However, the majority of them are younger than 1 Myr and have masses below 2 M$_\odot$. The high mass star WR 22 does not seem to have contributed to the formation of YSOs in the CrW region. The initial mass function slope, $\Gamma$, in this region is found to be -1.13 $\pm$ 0.20 in the mass range of 0.5 < M/M$_\odot$ < 4.8. The $K$-band luminosity function slope ($\alpha$) is also estimated as 0.31 $\pm$ 0.01. We also performed minimum spanning tree analysis of the YSOs in this region, which reveals that there are at least ten YSO cores associated with the molecular cloud, and that leads to an average core radius of 0.43 pc and a median branch length of 0.28 pc.

We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse spectroscopic data. HR5171A exhibits a complex appearance. Our AMBER data reveal a surprisingly large star for a YHG R*=1315+/-260Rsun\ (~6.1AU) at the distance of 3.6+/-0.5kpc. The source is surrounded by an extended nebulosity, and these data also show a large level of asymmetry in the brightness distribution of the system, which we attribute to a newly discovered companion star located in front of the primary star. The companion's signature is also detected in the visual photometry, which indicates an orbital period of Porb=1304+/-6d. Modeling the light curve with the NIGHTFALL program provides clear evidence that the system is a contact or possibly over-contact eclipsing binary. A total current system mass of 39^{+40}_{-22} solar mass and a high mass ratio q>10 is inferred for the system. The low-mass companion of HR5171 A is very close to the primary star that is embedded within its dense wind. Tight constraints on the inclination and vsini of the primary are lacking, which prevents us from determining its influence precisely on the mass-loss phenomenon, but the system is probably experiencing a wind Roche-Lobe overflow. Depending on the amount of angular momentum that can be transferred to the stellar envelope, HR5171 A may become a fast-rotating B[e]/Luminous Blue Variable (LBV)/Wolf-Rayet star. In any case, HR5171 A highlights the possible importance of binaries for interpreting the unstable YHGs and for massive star evolution in general.

Context. Massive-binary evolution models predict that some systems will go through an evolutionary phase where the original primary has become a supernova and left a compact object behind that then orbits a Wolf-Rayet (hereafter, WR) star. WR 138 is an X-ray bright WR star that has been described as a triple system, including a compact companion in a short-period orbit. Aims. Our goal is to search for spectroscopic evidence of a compact companion around WR 138. Methods. We used optical and X-ray spectra to search for signatures of a compact companion, which can be revealed by systematic variations in WR optical spectral lines induced by orbital motion of the compact companion or by hard, luminous X-rays from accretion onto this companion. Results. The optical spectra display emission-line profile variations that are most probably caused by clumps inside the stellar winds. The radial velocities do not vary on a short time-scale compatible with the suggested orbital period of a putative compact companion. The X-ray spectra are found to be normal for a WN5-6+OB system with no indication of accretion by a compact companion. Conclusions. There is no evidence for the presence of a compact companion, and we therefore conclude that WR 138 is a normal long-period (P~1521 d) eccentric WR+OB system.

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.

Context. Zeta Pup is the X-ray brightest O-type star of the sky. This object was regularly observed with the RGS instrument aboard XMM-Newton for calibration purposes, leading to an unprecedented set of high-quality spectra. Aims. We have previously reduced and extracted this data set and combined it into the most detailed high-resolution X-ray spectrum of any early-type star so far. Here we present the analysis of this spectrum accounting for the presence of structures in the stellar wind. Methods. For this purpose, we use our new modeling tool that allows fitting the entire spectrum with a multi-temperature plasma. We illustrate the impact of a proper treatment of the radial dependence of the X-ray opacity of the cool wind on the best-fit radial distribution of the temperature of the X-ray plasma. Results. The best fit of the RGS spectrum of Zeta Pup is obtained assuming no porosity. Four plasma components at temperatures between 0.10 and 0.69 keV are needed to adequately represent the observed spectrum. Whilst the hardest emission is concentrated between ~3 and 4 R*, the softer emission starts already at 1.5 R* and extends to the outer regions of the wind. Conclusions. The inferred radial distribution of the plasma temperatures agrees rather well with theoretical expectations. The mass- loss rate and CNO abundances corresponding to our best-fit model also agree quite well with the results of recent studies of Zeta Pup in the UV and optical domain.

We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonstrate that a cluster wind or wind-blown bubble is unlikely to provide a significant contribution to the X-ray emission detected from the central area of the cluster. We therefore suggest that X-ray emission at the cluster core originates from an ensemble of low- and solar-mass pre-main-sequence (PMS) stars, each of which would be too weak in X-rays to be detected individually. We attribute the X-ray emission from the dusty ridge to the embedded tight cluster of the new-born stars known in this area from infrared studies. Assuming that the levels of X-ray activity in young stars in the low-metallicity environment of NGC 602a are comparable to their Galactic counterparts, then the detected spatial distribution, spectral properties, and level of X-ray emission are largely consistent with those expected from low- and solar-mass PMS stars and young stellar objects (YSOs). This is the first discovery of X-ray emission attributable to PMS stars and YSOs in the SMC, which suggests that the accretion and dynamo processes in young, low-mass objects in the SMC resemble those in the Galaxy.

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

High-resolution X-ray spectra of O-type stars revealed less wind absorption than expected from smooth winds with conventional mass-loss rates. Various solutions have been proposed, including porous winds, optically thick clumps or an overall reduction of the mass-loss rates. The latter has a strong impact on the evolution of the star. Our final goal is to analyse high resolution X-ray spectra of O-type stars with a multi temperature plasma model in order to determine crucial wind parameters such as the mass loss rate, the CNO abundances and the X-ray temperature plasma distribution in the wind. In this context we are developing a modelling tool to calculate synthetic X-ray spectra. We present, here, the main ingredients and physics necessary for a such work. Our code uses the AtomDB emissivities to compute the intrinsic emissivity of the hot plasma as well as the CMFGEN model atmosphere code to evaluate the opacity of the cool wind. Following the comparison between two formalisms of stellar wind fragmentation, we introduce, for the first time in X-rays, the effects of a tenuous inter-clump medium. We then explore the quantitative impact of different model parameters on the X-ray spectra such as the position in the wind of the X-ray emitting plasma. For the first time, we show that the two formalisms of stellar wind fragmentation yield different results, although the differences for individual lines are small and can probably not be tested with the current generation of X-ray telescopes. As an illustration of our method, we compare various synthetic line profiles to the observed O VIII {\lambda} 18.97{\AA} line in the spectrum of {\zeta} Puppis. We illustrate how different combinations of parameters can actually lead to the same morphology of a single line, underlining the need to analyse the whole spectrum in a consistent way when attempting to constrain the parameters of the wind.

During recent years some non-thermal radio emitting OB stars have been discovered to be binary, or multiple systems. The non-thermal emission is due to synchrotron radiation that is emitted by electrons accelerated up to high energies. The electron acceleration occurs at the strong shocks created by the collision of radiatively-driven winds. Here we summarize the available radio data and more recent observations for the binary Cyg OB2 No. 9. We also show a new emission model which is being developed to compare the theoretical total radio flux and the spectral index with the observed radio light curves. This comparison will be useful in order to solve fundamental questions, such as the determination of the stellar mass loss rates, which are perturbed by clumping.

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

The distance of the very young open cluster Westerlund 2, which contains the very massive binary system WR 20a and is likely associated with a TeV source, has been the subject of much debate. We attempt a joint analysis of spectroscopic and photometric data of eclipsing binaries in the cluster to constrain its distance. A sample of 15 stars, including three eclipsing binaries (MSP 44, MSP 96, and MSP 223) was monitored with the FLAMES multi-object spectrograph. The spectroscopic data are analysed together with existing BV photometry. The analysis of the three eclipsing binaries clearly supports the larger values of the distance, around 8 kpc, and rules out values of about 2.4 - 2.8 kpc that have been suggested in the literature. Furthermore, our spectroscopic monitoring reveals no clear signature of binarity with periods shorter than 50 days in either the WN6ha star WR 20b, the early O-type stars MSP 18, MSP 171, MSP 182, MSP 183, MSP 199, and MSP 203, or three previously unknown mid O-type stars. The only newly identified candidate binary system is MSP 167. The absence of a binary signature is especially surprising for WR 20b and MSP 18, which were previously found to be bright X-ray sources. The distance of Westerlund 2 is confirmed to be around 8 kpc as previously suggested based on the spectrophotometry of its population of O-type stars and the analysis of the light curve of WR 20a. Our results suggest that short-period binary systems are not likely to be common, at least not among the population of O-type stars in the cluster.

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

We present a 1.42 square degree mosaic of diffuse X-ray emission in the Great Nebula in Carina from the Chandra X-ray Observatory Advanced CCD Imaging Spectrometer camera. After removing >14,000 X-ray point sources from the field, we smooth the remaining unresolved emission, tessellate it into segments of similar apparent surface brightness, and perform X-ray spectral fitting on those tessellates to infer the intrinsic properties of the X-ray-emitting plasma. By modeling faint resolved point sources, we estimate the contribution to the extended X-ray emission from unresolved point sources and show that the vast majority of Carina's unresolved X-ray emission is truly diffuse. Line-like correlated residuals in the X-ray spectral fits suggest that substantial X-ray emission is generated by charge exchange at the interfaces between Carina's hot, rarefied plasma and its many cold neutral pillars, ridges, and clumps.

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.

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

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

We present the results of a spectroscopic monitoring campaign devoted to three O-type stars in the Carina nebula. We derive the full SB2 orbital solution of the binary system Tr 16-112, an exceptional dissymmetrical system consisting of an O5.5-6V((f$^+$?p)) primary and a B2V-III secondary. We also report on low-amplitude brightness variations in Tr 16-112 that are likely due to the ellipsoidal shape of the O5.5-6 primary revolving in an eccentric orbit around the system's centre of mass. We detect for the first time a clear SB2 binary signature in the spectrum of HD 93343 (O8 + O8), although our data are not sufficient to establish an orbital solution. This system also displays low amplitude photometric modulations. On the other hand, no indication of multiplicity is found in the optical spectra of HD 93250. Finally, we discuss the general properties of multiple massive stars in the Carina OB1 association.

The Gaia satellite will be launched at the end of 2011. It will observe at least 1 billion stars, and among them several million emission line stars and hot stars. Gaia will provide parallaxes for each star and spectra for stars till V magnitude equal to 17. After a general description of Gaia, we present the codes and methods, which are currently developed by our team. They will provide automatically the astrophysical parameters and spectral classification for the hot and emission line stars in the Milky Way and other close Local Group galaxies such as the Magellanic Clouds.

We present the results of an optical and X-ray monitoring campaign on the short-period massive SB2 binary HD 152218. Combining our HiRes spectroscopic data with previous observations, we unveil the contradictions between the published orbital solutions. In particular, we solve the aliasing on the period and derive a value close to 5.604 d. Our eccentricity e = 0.259 +/- 0.006 is slightly lower than previously admitted. We show that HD 152218 is probably undergoing a relatively rapid apsidal motion of about 3deg/yr and we confirm the O9IV + O9.7V classification. We derive minimal masses of 15.82 +/- 0.26 Msol operator and 12.00 +/- 0.19 Msol operator and constrain the radius of the components to R1 = 10.3 +/- 1.3 Rsol and R2 = 7.8 +/- 1.7 Rsol. We also report the results of an XMM-Newton monitoring of the HD 152218 X-ray emission throughout its orbital motion. The averaged X-ray spectrum is relatively soft and it is well reproduced by a 2-T optically thin thermal plasma model with component temperatures about 0.3 and 0.7 keV. The system presents an increase of its X-ray flux by about 30% near apastron compared to periastron, which is interpreted as the signature of an ongoing wind-wind interaction process occurring within the wind acceleration region.

We present new XMM-Newton results on the field around the NGC346 star cluster in the SMC. This continues and extends previously published work on Chandra observations of the same field. The two XMM-Newton observations were obtained, respectively, six months before and six months after the previously published Chandra data. Of the 51 X-ray sources detected with XMM-Newton, 29 were already detected with Chandra. Comparing the properties of these X-ray sources in each of our three datasets has enabled us to investigate their variability on times scales of a year. Changes in the flux levels and/or spectral properties were observed for 21 of these sources. In addition, we discovered long-term variations in the X-ray properties of the peculiar system HD5980, a luminous blue variable star, that is likely to be a colliding wind binary system, which displayed the largest luminosity during the first XMM-Newton observation.