Massive, early type stars have been detected as radio sources for many decades. Their thermal winds radiate free-free continuum and in binary systems hosting a colliding-wind region, non-thermal emission has also been detected. To date, the most abundant data have been collected from frequencies higher than 1 GHz. We present here the results obtained from observations at 325 and 610 MHz, carried out with the Giant Metrewave Radio Telescope, of all known Wolf-Rayet and O-type stars encompassed in area of ~15 sq degrees centred on the Cygnus region. We report on the detection of 11 massive stars, including both Wolf-Rayet and O-type systems. The measured flux densities at decimeter wavelengths allowed us to study the radio spectrum of the binary systems and to propose a consistent interpretation in terms of physical processes affecting the wide-band radio emission from these objects. WR 140 was detected at 610 MHz, but not at 325 MHz, very likely because of the strong impact of free-free absorption. We also report - for the first time - on the detection of a colliding-wind binary system down to 150 MHz, pertaining to the system of WR 146, making use of complementary information extracted from the TIFR GMRT Sky Survey. Its spectral energy distribution clearly shows the turnover at a frequency of about 600 MHz, that we interpret to be due to free-free absorption. Finally, we report on the identification of two additional particle-accelerating colliding-wind binaries, namely Cyg OB2 12 and ALS 15108 AB.

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.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

We investigated the long-term behaviour in X-rays of the colliding wind binary WR 25, using archival data obtained with Suzaku, Swift, XMM-Newton, and NuSTAR spanning over ~16 years. Our analysis reveals phase-locked variations repeating consistently over many consecutive orbits, in agreement with an X-ray emission fully explained by thermal emission from the colliding winds in the 208-d orbit. We report on a significant deviation of the X-ray flux with respect to the 1/D trend (expected for adiabatic shocked winds) close to periastron passage. The absence of a drop in post-shock plasma temperature close to periastron suggests this break in trend cannot be explained in terms of reduced pre-shock velocities in this part of the orbit. Finally, NuSTAR data reveal a lack of hard X-ray emission (above 10.0 keV) above the background level. Upper limits on a putative non-thermal emission strongly suggest that the sensitivity of present hard X-ray observatories is not sufficient to detect non-thermal emission from massive binaries above 10 keV, unless the wind kinetic power is large enough to significantly feed particle acceleration in the wind-wind interaction.

Several tens of massive binary systems display indirect, or even strong evidence for non-thermal radio emission, hence their particle accelerator status. These objects are referred to as particle-accelerating colliding-wind binaries (PACWBs). WR 133 is one of the shortest period Wolf-Rayet + O systems in this category, and is therefore critical to characterize the boundaries of the parameter space adequate for particle acceleration in massive binaries. Our methodology consists in analyzing JVLA observations of WR 133 at different epochs to search for compelling evidence for a phase-locked variation attributable to synchrotron emission produced in the colliding-wind region. New data obtained during two orbits reveal a steady and thermal emission spectrum, in apparent contradiction with the previous detection of non-thermal emission. The thermal nature of the radio spectrum along the 112.4-d orbit is supported by the strong free-free absorption by the dense stellar winds, and shows that the simple binary scenario cannot explain the non-thermal emission reported previously. Alternatively, a triple system scenario with a wide, outer orbit would fit with the observational facts reported previously and in this paper, albeit no hint for the existence of a third component exists to date. The epoch-dependent nature of the identification of synchrotron radio emission in WR 133 emphasizes the issue of observational biases in the identification of PACWBs, that undoubtedly affect the present census of PACWB among colliding-wind binaries.

The massive binary system WR11 has been recently proposed as the counterpart of a Fermi source. If correct, it would be the second colliding wind binary detected in GeV gamma-rays. However, the reported flux measurements from 1.4 to 8.64GHz fail to establish the presence of non-thermal (synchrotron) emission from this source. Moreover, WR11 is not the only radio source within the Fermi detection box. Other possible counterparts have been identified in archival data, some of which present strong non-thermal radio emission. We conducted "-resolution observations towards WR11 at very low frequencies (150 to 1400~MHz) where the NT emission is expected to dominate, and present a catalog of more than 400 radio-emitters, among which a significant part is detected at more than one frequency, including limited spectral index information. A search for counterparts for this last group pointed at MOST0808-471, a source 2' away from WR11, as a promising candidate for high-energy emission, with resolved structure along 325 - 1390 MHz. For it, we reprocessed archive interferometric data up to 22.3 GHz and obtained a non-thermal radio spectral index of -0.97 +- 0.09. However, multiwavelength observations of this source are required to establish its nature and to assess whether it can produce (part of) the observed gamma-rays. WR11 spectrum follows a spectral index of 0.74 +- 0.03 from 150 to 230 GHz, consistent with thermal emission. We interpret that any putative synchrotron radiation from the colliding-wind region of this relatively short-period system is absorbed in the photospheres of the individual components. Notwithstanding, the new radio data allowed to derive a mass loss rate of 0.000025 Mo/yr, which, according to the latest models for gamma-ray emission in WR 11, would suffice to provide the required kinetic power to feed non-thermal radiation processes.

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

Massive runaway stars produce bow shocks through the interaction of their winds with the interstellar medium, with the prospect for particle acceleration by the shocks. These objects are consequently candidates for non-thermal emission. Our aim is to investigate the X-ray emission from these sources. We observed with XMM-Newton a sample of 5 bow shock runaways, which constitutes a significant improvement of the sample of bow shock runaways studied in X-rays so far. A careful analysis of the data did not reveal any X-ray emission related to the bow shocks. However, X-ray emission from the stars is detected, in agreement with the expected thermal emission from stellar winds. On the basis of background measurements we derive conservative upper limits between 0.3 and 10 keV on the bow shocks emission. Using a simple radiation model, these limits together with radio upper limits allow us to constrain some of the main physical quantities involved in the non-thermal emission processes, such as the magnetic field strength and the amount of incident infrared photons. The reasons likely responsible for the non-detection of non-thermal radiation are discussed. Finally, using energy budget arguments, we investigate the detectability of inverse Compton X-rays in a more extended sample of catalogued runaway star bow shocks. From our analysis we conclude that a clear identification of non-thermal X-rays from massive runaway bow shocks requires one order of magnitude (or higher) sensitivity improvement with respect to present observatories.

Particle-accelerating colliding-wind binaries (PACWBs) are multiple systems made of early-type stars able to accelerate particles up to relativistic velocities. The relativistic particles can interact with different fields (magnetic or radiation) in the colliding-wind region and produce non-thermal emission. In many cases, non-thermal synchrotron radiation might be observable and thus constitute an indicator of the existence of a relativistic particle population in these multiple systems. To date, the catalogue of PACWBs includes about 40 objects spread over many stellar types and evolutionary stages, with no clear trend pointing to privileged subclasses of objects likely to accelerate particles. This paper aims at discussing critically some criteria for selecting new candidates among massive binaries. The subsequent search for non-thermal radiation in these objects is expected to lead to new detections of particle accelerators. On the basis of this discussion, some broad ideas for observation strategies are formulated. At this stage of the investigation of PACWBs, there is no clear reason to consider particle acceleration in massive binaries as an anomaly or even as a rare phenomenon. We therefore consider that several PACWBs will be detected in the forthcoming years, essentially using sensitive radio interferometers which are capable of measuring synchrotron emission from colliding-wind binaries. Prospects for high-energy detections are also briefly addressed.

Our long term aim is to derive model-independent stellar masses and distances for long period massive binaries by combining apparent astrometric orbit with double-lined radial velocity amplitudes (SB2). We follow-up ten O+O binaries with AMBER, PIONIER and GRAVITY at the VLTI. Here, we report about 130 astrometric observations over the last 7 years. We combine this dataset with distance estimates to compute the total mass of the systems. We also compute preliminary individual component masses for the five systems with available SB2 radial velocities. Nine over the ten binaries have their three dimensional orbit well constrained. Four of them are known colliding wind, non-thermal radio emitters, and thus constitute valuable targets for future high angular resolution radio imaging. Two binaries break the correlation between period and eccentricity tentatively observed in previous studies. It suggests either that massive star formation produce a wide range of systems, or that several binary formation mechanisms are at play. Finally, we found that the use of existing SB2 radial velocity amplitudes can lead to unrealistic masses and distances. If not understood, the biases in radial velocity amplitudes will represent an intrinsic limitation for estimating dynamical masses from SB2+interferometry or SB2+Gaia. Nevertheless, our results can be combined with future Gaia astrometry to measure the dynamical masses and distances of the individual components with an accuracy of 5 to 15\%, completely independently of the radial velocities.

This paper presents the results of the analysis of the very first dedicated X-ray observation with XMM-Newton of WR106. This carbon-rich WC9d Wolf-Rayet star belongs to the category of persistent dust makers (WCd stars). The issue of the multiplicity of these dust makers is pivotal to understand the dust formation process, and in this context X-ray observations may allow to reveal an X-ray emission attributable to colliding-winds in a binary system. The main result of this analysis is the lack of detection of X-rays coming from WR106. Upper limits on the X-ray flux are estimated, but the derived numbers are not sufficient to provide compelling constraints on the existence or not of a colliding-wind region. Detailed inspection of archive data bases reveals that persistent dust makers have been poorly investigated by the most sensitive X-ray observatories. Certainly, the combination of several approaches to indirectly constrain their multiplicity should be applied to lift a part of the veil on the nature of these persistent dust makers.

The long-term (over more than one decade) X-ray emission from two massive stellar systems known to be particle accelerators is investigated using XMM-Newton. Their X-ray properties are interpreted taking into account recent information about their multiplicity and orbital parameters. The two targets, HD168112 and HD167971 appear to be overluminous in X-rays, lending additional support to the idea that a significant contribution of the X-ray emission comes from colliding-wind regions. The variability of the X-ray flux from HD168112 is interpreted in terms of varying separation expected to follow the 1/D rule for adiabatic shocked winds. For HD167971, marginal decrease of the X-ray flux in September 2002 could tentatively be explained by a partial wind eclipse in the close pair. No long-term variability could be demonstrated despite the significant difference of separation between 2002 and 2014. This suggests the colliding-wind region in the wide orbit does not contribute a lot to the total X-ray emission, with a main contribution coming from the radiative shocked winds in the eclipsing pair. The later result provides evidence that shocks in a colliding-wind region may be efficient particle accelerators even in the absence of bright X-ray emission, suggesting particle acceleration may operate in a wide range of conditions. Finally, in hierarchical triple O-type systems, thermal X-rays do not necessarily constitute an efficient tracer to detect the wind-wind interaction in the long period orbit.

Aims. Having established the binary status of nineteen O-type stars located in four Cygnus OB associations, we now determine their fundamental parameters to constrain their properties and their evolutionary status. We also investigate their surface nitrogen abundances, which we compare with other results from the literature obtained for galactic O-type stars. Methods. Using optical spectra collected for each object in our sample and some UV data from the archives, we apply the CMFGEN atmosphere code to determine their main properties. For the binary systems, we have disentangled the components to obtain their individual spectra and investigate them as if they were single stars. Results. We find that the distances of several presumably single O-type stars seem poorly constrained because their luminosities are not in agreement with the "standard" luminosities of stars with similar spectral types. The ages of these O-type stars are all less than 7 Myrs. Therefore, the ages of these stars agree with those, quoted in the literature, of the four associations, except for CygOB8 for which the stars seem older than the association itself. However, we point out that the distance of certain stars is debatable relative to values found in the literature. The N content of these stars put in perspective with N contents of several other galactic O-type stars seems to draw the same five groups as found in the "Hunter" diagram for the O and B-type stars in the LMC even though their locations are obviously different. We determine mass-loss rates for several objects from the Halpha line and UV spectra. Finally, we confirm the "mass discrepancy" especially for O stars with masses smaller than 30 Msun. .

Radio observations are an effective tool to discover particle acceleration regions in colliding-wind binaries, through detection of synchrotron radiation; these regions are natural laboratories for the study of relativistic particles. Wind-collision region (WCR) models can reproduce the radio continuum spectra of massive binaries that contain both thermal and non-thermal radio emission; however, key constraints for models come from high-resolution imaging. Only five WCRs have been resolved to date at radio frequencies at milliarcsec (mas) angular scales. The source HD 93129A, prototype of the very few known O2 I stars, is a promising target for study: recently, a second massive, early-type star about 50 mas away was discovered, and a non-thermal radio source detected in the region. Preliminary long-baseline array data suggest that a significant fraction of the radio emission from the system comes from a putative WCR. We sought evidence that HD 93129A is a massive binary system with colliding stellar winds that produce non-thermal radiation, through spatially resolved images of the radio emitting regions. We completed observations with the Australian Long Baseline Array (LBA) to resolve the system at mas angular resolutions and reduced archival Australia Telescope Compact Array (ATCA) data to derive the total radio emission. We also compiled optical astrometric data of the system in a homogeneous way. We reduced historical Hubble Space Telescope data and obtained absolute and relative astrometry with milliarcsec accuracy. The astrometric analysis leads us to conclude that the two stars in HD 93129A form a gravitationally bound system. The LBA data reveal an extended arc-shaped non-thermal source between the two stars, indicative of a WCR. The wind momentum-rate ratio of the two stellar winds is estimated. The ATCA data show a point source with a change in flux level ...

Some colliding-wind massive binaries, called particle-accelerating colliding-wind binaries (PACWB), exhibit synchrotron radio emission, which is assumed to be generated by a stellar magnetic field. However, no measurement of magnetic fields in these stars has ever been performed. We aim at quantifying the possible stellar magnetic fields present in PACWB to provide constraints for models. We gathered 21 high-resolution spectropolarimetric observations of 9 PACWB available in the ESPaDOnS, Narval and HarpsPol archives. We analysed these observations with the Least Squares Deconvolution method. We separated the binary spectral components when possible. No magnetic signature is detected in any of the 9 PACWB stars and all longitudinal field measurements are compatible with 0 G. We derived the upper field strength of a possible field that could have remained hidden in the noise of the data. While the data are not very constraining for some stars, for several stars we could derive an upper limit of the polar field strength of the order of 200 G. We can therefore exclude the presence of strong or moderate stellar magnetic fields in PACWB, typical of the ones present in magnetic massive stars. Weak magnetic fields could however be present in these objects. These observational results provide the first quantitative constraints for future models of PACWB.

We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind.

Among evolved massive stars likely in transition to the Wolf-Rayet phase, IRC +10420 is probably one of the most enigmatic. It belongs to the category of yellow hypergiants and it is characterized by quite high mass loss episodes. Even though IRC +10420 benefited of many observations in several wavelength domains, it has never been a target for an X-ray observatory. We report here on the very first dedicated observation of IRC +10420 in X-rays, using the XMM-Newton satellite. Even though the target is not detected, we derive X-ray flux upper limits of the order of 1--3 10^-14 erg cm^-2 s^-1 (between 0.3 and 10.0 keV), and we discuss the case of IRC +10420 in the framework of emission models likely to be adequate for such an object. Using the Optical/UV Monitor on board XMM-Newton, we present the very first upper limits of the flux density of IRC +10420 in the UV domain (between 1800 and 2250 A, and between 2050 and 2450 A). Finally, we also report on the detection in this field of 10 X-ray and 7 UV point sources, and we briefly discuss their properties and potential counterparts at longer wavelengths.

Astrochemistry is a discipline that studies physico-chemical processes in astrophysical environments. Such environments are characterized by conditions that are substantially different from those existing in usual chemical laboratories. Models which aim to explain the formation of molecular species in interstellar environments must take into account various factors, including many that are directly, or indirectly related to the populations of massive stars in galaxies. The aim of this paper is to review the influence of massive stars, whatever their evolution stage, on the physico-chemical processes at work in interstellar environments. These influences include the ultraviolet radiation field, the production of high energy particles, the synthesis of radionuclides and the formation of shocks that permeate the interstellar medium.

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.

Massive systems made of two or more stars are known to be the site for interesting physical processes -- including at least in some cases -- particle acceleration. Over the past decade, this topic motivated a particular effort to unveil the properties of these systems and characterize the circumstances responsible for the acceleration of particles and the potential role of pre-supernova massive stars in the production of high energy particles in our Galaxy. Although previous studies on this topic were mostly devoted to processes in general, or to a few individual objects in particular, a unified target-oriented census of particle-accelerating colliding-wind binaries (hereafter PACWBs) does not exist yet. This paper aims at making a general and unified census of these systems, emphasizing their main properties. A general discussion includes energetic considerations along with wind properties in relation with non-thermal emission processes that are likely at work in colliding-wind binaries. Finally, some guidelines for future observational and theoretical studies are drawn.

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.

Astrochemistry aims at studying chemical processes in astronomical environments. This discipline -- located at the crossroad between astrophysics and chemistry -- is rapidly evolving and explores the issue of the formation of molecules of increasing complexity in particular physical conditions that deviate significantly from those frequently encountered in chemistry laboratories. The main goal of this paper is to provide an overview of this discipline. So far, about 170 molecules have been identified in the interstellar medium (ISM). The presence of this molecular diversity constitutes a firm evidence that efficient formation processes are at work in the interstellar medium. This paper aims at summarizing most of present ideas that are explored by astrochemists to investigate the chemistry taking place in various astronomical environments, with emphasis on the particular conditions which are met in space (including radiation fields, cosmic-rays, low densities...). The more ambitious question of the molecular complexity is addressed following two approaches presented to be converging. The first approach considers the growing complexity starting from the most simple chemical species in interstellar environments, and the second approach envisages successive precursors of the most complex species commonly found on Earth, and in particular in our biochemistry. The issue of molecular complexity constitutes one of the main modern scientific questions addressed by astrochemistry, and it is used as a guideline across this paper.

A rare group of high mass X-ray binaries (HMXBs) are known that also exhibit MeV, GeV, and/or TeV emission ("gamma-ray binaries"). Expanding the sample of gamma-ray binaries and identifying unknown Fermi sources are currently of great interest to the community. Based upon their positional coincidence with the unidentified Fermi sources 1FGL J1127.7-6244c and 1FGL J1808.5-1954c, the Be stars HD 99771 and HD 165783 have been proposed as gamma-ray binary candidates. During Fermi Cycle 4, we have performed multiwavelength observations of these sources using XMM-Newton and the CTIO 1.5m telescope. We do not confirm high energy emission from the Be stars. Here we examine other X-ray sources in the field of view that are potential counterparts to the Fermi sources.

The members of the scarce category of Of^+ supergiants present properties that are intermediate between regular O-stars and Wolf-Rayet (WR) stars. Significant similarities between these transitional stars and WN-type objects are now clearly established, at least in the visible and near-infrared domains, pointing to common stellar wind properties. In this study, we report on the first dedicated X-ray observations of HD16691 (O4If^+) and HD14947 (O5f^+), revealing a soft thermal spectrum in agreement with the expected X-ray emission from a single O-type star. However, the X-ray luminosity of our targets is slightly lower than expected for single O-type stars, suggesting that the particular properties of their stellar wind has also a significant impact on the X-ray emission of these objects on the way to the WN category. We argue that the X-ray under-luminosity of HD16691 and HD14947 may be interpreted as the signature in X-rays of the intermediate stage between O and WR stars, as a consequence of enhanced wind density.

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

Context. Establishing the multiplicity of O-type stars is the first step towards accurately determining their stellar parameters. Moreover, the distribution of the orbital parameters provides observational clues to the way that O-type stars form and to the interactions during their evolution. Aims. Our objective is to constrain the multiplicity of a sample of O-type stars belonging to poorly investigated OB associations in the Cygnus complex and for the first time to provide orbital parameters for binaries identified in our sample. Such information is relevant to addressing the issue of the binarity in the context of O-type star formation scenarios. Methods. We performed a long-term spectroscopic survey of nineteen O-type stars. We searched for radial velocity variations to unveil binaries on timescales from a few days up to a few years, on the basis of a large set of optical spectra. Results. We confirm the binarity for four objects: HD193443, HD228989, HD229234 and HD194649. We derive for the first time the orbital solutions of three systems, and we confirm the values of the fourth, showing that these four systems all have orbital periods shorter than 10 days. Besides these results, we also detect several objects that show non-periodic line profile variations in some of their spectral lines. These variations mainly occur in the spectral lines, that are generally affected by the stellar wind and are not likely to be related to binarity. Conclusions. The minimal binary fraction in our sample is estimated to be 21%, but it varies from one OB association to the next. Indeed, 33% of O stars belonging to CygOB1 are binary systems, 0% in CygOB3, 0% in CygOB8, and 33% in CygOB9. Our spectroscopic investigation also stresses the absence of long-period systems among the stars in our sample...

Cyg OB2 #9 is one of a small set of non-thermal radio emitting massive O-star binaries. The non-thermal radiation is due to synchrotron emission in the colliding-wind region. Cyg OB2 #9 was only recently discovered to be a binary system and a multi-wavelength campaign was organized to study its 2011 periastron passage. We report here on the results of the radio observations obtained in this monitoring campaign. We used the Expanded Very Large Array (EVLA) radio interferometer to obtain 6 and 20 cm continuum fluxes. The observed radio light curve shows a steep drop in flux sometime before periastron. The fluxes drop to a level that is comparable to the expected free-free emission from the stellar winds, suggesting that the non-thermal emitting region is completely hidden at that time. After periastron passage, the fluxes slowly increase. We introduce a simple model to solve the radiative transfer in the stellar winds and the colliding-wind region, and thus determine the expected behaviour of the radio light curve. From the asymmetry of the light curve, we show that the primary has the stronger wind. This is somewhat unexpected if we use the astrophysical parameters based on theoretical calibrations. But it becomes entirely feasible if we take into account that a given spectral type - luminosity class combination covers a range of astrophysical parameters. The colliding-wind region also contributes to the free-free emission, which can help to explain the high values of the spectral index seen after periastron passage. Combining our data with older Very Large Array (VLA) data allows us to derive a period P = 860.0 +- 3.7 days for this system. With this period, we update the orbital parameters that were derived in the first paper of this series.

HD 195592 is an O-type super-giant star, known as a well-established runaway. Recently, a Fermi gamma-ray source (2FGL J2030.7+4417) with a position compatible with that of HD 195592 has been reported. Our goal is to explore the scenario where HD 195592 is the counterpart of the Fermi gamma-ray source modeling the non-thermal emission produced in the bowshock of the runaway star. We calculate the spectral energy distribution of the radiation produced in the bowshock of HD 195592 and we compare it with Fermi observations of 2FGL J2030.7+4417. We present relativistic particle losses and the resulting radiation of the bowshock of HD 195592 and show that the latter is compatible with the detected gamma-ray emission. We conclude that the Fermi source 2FGL J2030.7+4417 might be produced, under some energetic assumptions, by inverse Compton up-scattering of photons from the heated dust in the bowshock of the runaway star. HD 195592 might therefore be the very first object detected belonging to the category of gamma-ray emitting runaway massive stars, whose existence has been recently predicted.

The O-type object 9 Sgr is a well-known synchrotron radio emitter. This feature is usually attributed to colliding-wind binary systems, but 9 Sgr was long considered a single star. We have conducted a long-term spectroscopic monitoring of this star to investigate its multiplicity and search for evidence for wind-wind interactions. Radial velocities are determined and analysed using various period search methods. Spectral disentangling is applied to separate the spectra of the components of the binary system. We derive the first ever orbital solution of 9 Sgr. The system is found to consist of an O3.5 V((f+)) primary and an O5- 5.5 V((f)) secondary moving around each other on a highly eccentric (e = 0.7), 8.6 year orbit. The spectra reveal no variable emission lines that could be formed in the wind interaction zone in agreement with the expected properties of the interaction in such a wide system. Our results provide further support to the paradigm of synchrotron radio emission from early-type stars being a manifestation of interacting winds in a binary system.

Using optical long baseline interferometry, we resolved for the first time the two wide components of HD167971, a candidate hierarchical triple system known to efficiently accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a gravitational link between the O8 supergiant and the close eclipsing O + O binary. The separation varies from 8 to 15 mas over the three-year baseline of our observations, suggesting that the components evolve on a wide and very eccentric orbit (most probably e>0.5). These results provide evidence that the wide orbit revealed by our study is not coplanar with the orbit of the inner eclipsing binary. From our measurements of the near-infrared luminosity ratio, we constrain the spectral classification of the components in the close binary to be O6-O7, and confirm that these stars are likely main-sequence objects. Our results are discussed in the context of the bright non-thermal radio emission already reported for this system, and we provide arguments in favour of a maximum radio emission coincident with periastron passage. HD167971 turns out to be an efficient O-type particle accelerator that constitutes a valuable target for future high angular resolution radio imaging using VLBI facilities.

Context. HD150136 is one of the nearest systems harbouring an O3 star. Although this system was for a long time considered as binary, more recent investigations have suggested the possible existence of a third component. Aims. We present a detailed analysis of HD 150136 to confirm the triple nature of this system. In addition, we investigate the physical properties of the individual components of this system. Methods. We analysed high-resolution, high signal-to-noise data collected through multi-epoch runs spread over ten years. We applied a disentangling program to refine the radial velocities and to obtain the individual spectra of each star. With the radial velocities, we computed the orbital solution of the inner system, and we describe the main properties of the orbit of the outer star such as the preliminary mass ratio, the eccentricity, and the orbital-period range. With the individual spectra, we determined the stellar parameters of each star by means of the CMFGEN atmosphere code. Results. We offer clear evidence that HD 150136 is a triple system composed of an O3V((f\ast))-3.5V((f+)), an O5.5-6V((f)), and an O6.5-7V((f)) star. The three stars are between 0-3 Myr old. We derive dynamical masses of about 64, 40, and 35 Msun for the primary, the secondary and the third components by assuming an inclination of 49{\deg}. It currently corresponds to one of the most massive systems in our galaxy. The third star moves with a period in the range of 2950 to 5500 d on an outer orbit with an eccentricity of at least 0.3. This discovery makes HD 150136 the first confirmed triple system with an O3 primary star. However, because of the long orbital period, our dataset is not sufficient to constrain the orbital solution of the tertiary component with high accuracy.

Studies are under way to propose a new generation of post-VLTI interferometers. The Carlina concept studied at the Haute- Provence Observatory is one of the proposed solutions. It consists in an optical interferometer configured like a diluted version of the Arecibo radio telescope: above the diluted primary mirror made of fixed cospherical segments, a helium balloon (or cables suspended between two mountains), carries a gondola containing the focal optics. Since 2003, we have been building a technical demonstrator of this diluted telescope. First fringes were obtained in May 2004 with two closely-spaced primary segments and a CCD on the focal gondola. We have been testing the whole optical train with three primary mirrors. The main aim of this article is to describe the metrology that we have conceived, and tested under the helium balloon to align the primary mirrors separate by 5-10 m on the ground with an accuracy of a few microns. The servo loop stabilizes the mirror of metrology under the helium balloon with an accuracy better than 5 mm while it moves horizontally by 30 cm in open loop by 10-20 km/h of wind. We have obtained the white fringes of metrology; i.e., the three mirrors are aligned (cospherized) with an accuracy of {\approx} 1 micron. We show data proving the stability of fringes over 15 minutes, therefore providing evidence that the mechanical parts are stabilized within a few microns. This is an important step that demonstrates the feasibility of building a diluted telescope using cables strained between cliffs or under a balloon. Carlina, like the MMT or LBT, could be one of the first members of a new class of telescopes named diluted telescopes.

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.

As the brightest O-type X-ray source in the Carina nebula, HD 93250 (O4 III(fc)) is X-ray overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales of years suggesting either a single star, a binary system seen pole-one or a very long period and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 0.8+/-0.1 and a separation of 1.5+/-0.2 x 10E-03 arcsec. At the distance of Carina, this corresponds to a projected physical distance of 3.5 A.U. While a quantitative investigation would require a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively explain both its X-ray flux and hardness and its non-thermal radio emission in the framework of a colliding wind scenario. We also discuss various observational biases. We show that, due to line-blending of two similar spectral components, HD 93250 could have a period as short as one to several years despite the lack of measurable radial velocity variations.

Tracing the formation and evolution of all supermassive black holes, including the obscured ones, understanding how black holes influence their surroundings and how matter behaves under extreme conditions, are recognized as key science objectives to be addressed by the next generation of instruments. These are the main goals of the COSPIX proposal, made to ESA in December 2010 in the context of its call for selection of the M3 mission. In addition, COSPIX, will also provide key measurements on the non thermal Universe, particularly in relation to the question of the acceleration of particles, as well as on many other fundamental questions as for example the energetic particle content of clusters of galaxies. COSPIX is proposed as an observatory operating from 0.3 to more than 100 keV. The payload features a single long focal length focusing telescope offering an effective area close to ten times larger than any scheduled focusing mission at 30 keV, an angular resolution better than 20 arcseconds in hard X-rays, and polarimetric capabilities within the same focal plane instrumentation. In this paper, we describe the science objectives of the mission, its baseline design, and its performances, as proposed to ESA.

Some OB stars show variable non-thermal radio emission. The non-thermal emission is due to synchrotron radiation that is emitted by electrons accelerated to high energies. The electron acceleration occurs at strong shocks created by the collision of radiatively-driven stellar winds in binary systems. Here we present results of our modelling of two colliding wind systems: Cyg OB2 No. 8A and Cyg OB2 No. 9.

HD 259440 is a B0pe star that was proposed as the optical counterpart to the gamma-ray source HESS J0632+057. Here we present optical spectra of HD 259440 acquired to investigate the stellar parameters, the properties of the Be star disk, and evidence of binarity in this system. Emission from the H-alpha line shows evidence of a spiral density wave in the nearly edge-on disk. We find a best fit stellar effective temperature of 27500-30000 K and a log surface gravity of 3.75-4.0, although our fits are somewhat ambiguous due to scattered light from the circumstellar disk. We derive a mass of 13.2-19.0 M_sun and a radius of 6.0-9.6 R_sun. By fitting the spectral energy distribution, we find a distance between 1.1-1.7 kpc. We do not detect any significant radial velocity shifts in our data, ruling out orbital periods shorter than one month. If HD 259440 is a binary, it is likely a long period (> 100 d) system.

We study the non-thermal radio emission of the binary Cyg OB2 No. 8A, to see if it is variable and if that variability is locked to the orbital phase. We investigate if the synchrotron emission generated in the colliding-wind region of this binary can explain the observations and we verify that our proposed model is compatible with the X-ray data. We use both new and archive radio data from the Very Large Array (VLA) to construct a light curve as a function of orbital phase. We also present new X-ray data that allow us to improve the X-ray light curve. We develop a numerical model for the colliding-wind region and the synchrotron emission it generates. The model also includes free-free absorption and emission due to the stellar winds of both stars. In this way we construct artificial radio light curves and compare them with the observed one. The observed radio fluxes show phase-locked variability. Our model can explain this variability because the synchrotron emitting region is not completely hidden by the free-free absorption. In order to obtain a better agreement for the phases of minimum and maximum flux we need to use stellar wind parameters for the binary components which are somewhat different from typical values for single stars. We verify that the change in stellar parameters does not influence the interpretation of the X-ray light curve. Our model has trouble explaining the observed radio spectral index. This could indicate the presence of clumping or porosity in the stellar wind, which - through its influence on both the Razin effect and the free-free absorption - can considerably influence the spectral index. Non-thermal radio emitters could therefore open a valuable pathway to investigate the difficult issue of clumping in stellar winds.

After the first detection of its binary nature, the spectroscopic monitoring of the non-thermal radio emitter Cyg OB2 #9 (P=2.4yrs) has continued, doubling the number of available spectra of the star. Since the discovery paper of 2008, a second periastron passage has occurred in February 2009. Using a variety of techniques, the radial velocities could be estimated and a first, preliminary orbital solution was derived from the HeI5876 line. The mass ratio appears close to unity and the eccentricity is large, 0.7--0.75. X-ray data from 2004 and 2007 are also analyzed in quest of peculiarities linked to binarity. The observations reveal no large overluminosity nor strong hardness, but it must be noted that the high-energy data were taken after the periastron passage, at a time where colliding wind emission may be low. Some unusual X-ray variability is however detected, with a 10% flux decrease between 2004 and 2007. To clarify their origin and find a more obvious signature of the wind-wind collision, additional data, taken at periastron and close to it, are needed.

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

The very young open cluster NGC 6383 centered on the O-star binary HD 159176 is an interesting place for studying the impact of early-type stars with strong radiation fields and powerful winds on the formation processes of low-mass stars. To investigate this process, it is necessary to determine the characteristics (age, presence, or absence of circumstellar material) of the population of low-mass pre-main-sequence (PMS) stars in the cluster. We obtained deep U B V (R I)_c H-alpha photometric data of the entire cluster as well as medium-resolution optical spectroscopy of a subsample of X-ray selected objects. Our spectroscopic data reveal only very weak H-alpha emission lines in a few X-ray selected PMS candidates. We photometrically identify a number of H-alpha emission candidates but their cluster membership is uncertain. We find that the fainter objects in the field of view have a wide range of extinction (up to A_V = 20), one X-ray selected OB star having A_V ~ 8. Our investigation uncovers a population of PMS stars in NGC 6383 that are probably coeval with HD 159176. In addition, we detect a population of reddened objects that are probably located at different depths within the natal molecular cloud of the cluster. Finally, we identify a rather complex spatial distribution of H-alpha emitters, which is probably indicative of a severe contamination by foreground and background stars.

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to explore the mechanism by which it was ejected from the open cluster of its birth. Here we present recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic binary and physical parameters of the O star primary. We also present the first XMM-Newton observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place limits on the emission from a proposed compact companion in the system, and we rule out a quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved secondary companion is also not detected in our optical spectra of the binary, and it is difficult to conclude that a gravitational interaction could have ejected this runaway binary with a low mass optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent black hole with conditions unfavorable for accretion at the time of our observations.

With an orbital period of about 3.078 days, the double-lined spectroscopic binary HD 48099 is, still now, the only short-period O+O system known in the Mon OB2 association. Even though an orbital solution has already been derived for this system, few information are available about the individual stars. We present, in this paper, the results of a long-term spectroscopic campaign. We derive a new orbital solution and apply a disentangling method to recover the mean spectrum of each star. To improve our knowledge concerning both components, we determine their spectral classifications and their projected rotational velocities. We also constrain the main stellar parameters of both stars by using the CMFGEN atmosphere code and provide the wind properties for the primary star through the study of IUE spectra. This investigation reveals that HD 48099 is an O5.5 V((f))+O9 V binary with M_1 sin^3 i = 0.70 M_{\sun} and M_2 sin^3 i = 0.39 M_{\sun}, implying a rather low orbital inclination. This result, combined with both a large effective temperature and log g, suggests that the primary star (v sini ~ 91 km s^-1) is actually a fast rotator with a strongly clumped wind and a nitrogen abundance of about 8 times the solar value.

The transition from early Of stars to WN type objects is poorly understood. O-type supergiants with emission lines (OIf+) are considered to be intermediate between these two classes. The scope of this paper is to investigate the spectral variability of three Of+ supergiants. We constituted spectral time series of unprecedented quality for our targets (~ 200 spectra in total), essentially in the blue domain, covering time-scales from a few hours up to a few years. Time Variance Spectrum (TVS) and Fourier analyses were performed in order to characterize their spectral variability. We report on a correlated significant line profile variability in the prominent He II \lambda 4686 and H\beta lines most likely related to the strong stellar winds. The variability pattern is similar for the three stars investigated (HD14947, HD15570 and HD16691), and the main differences are more quantitative than qualitative. However, the reported time-scales are somewhat different, and the most striking variability pattern is reported for HD16691. We did not find any clear evidence for binarity, and we focus mainly on an interpretation based on a single star scenario. We show that the behaviour of the three stars investigated in this study present strong similarities, pointing to a putative common scenario, even though a few differences should be noted. Our preferred interpretation scheme is that of Large Scale Corotating Structures modulating the profile of the lines that are produced in the strong stellar wind.

The scope of this paper is to perform a detailed spectroscopic study of the northern O-type supergiant HD 195592. We use a large sample of high quality spectra in order to investigate its multiplicity, and to probe the line profile variability. Our analysis reveals a clear spectroscopic binary signature in the profile of the He {\sc i} $\lambda$ 6678 line, pointing to a probable O + B system. We report on low amplitude radial velocity variations in every strong absorption line in the blue spectrum of HD 195592. These variations are ruled by two time-scales respectively of 5.063 and about 20 days. The former is firmly established, whilst the latter is poorly constrained. We report also on a very significant line profile variability of the H $\beta$ line, with time scales strongly related to those of the radial velocities. Our results provide significant evidence that HD 195592 is a binary system, with a period that might be the variability time-scale of about 5 days. The second time scale may be the signature of an additional star moving along a wider orbit provided its mass is low enough, even though direct evidence for the presence of a third star is still lacking. Alternatively, the second time-scale may be the signature of a variability intrinsic to the stellar wind of the primary, potentially related to the stellar rotation.

Aims. We present the results obtained from a long-term spectroscopic campaign devoted to the multiplicity of O-type stars in the young open cluster NGC2244 and in the Mon OB2 association. Methods. Our spectroscopic monitoring was performed over several years, allowing us to probe different time-scales. For each star, several spectral diagnostic tools are applied, in order to search for line shifts and profile variations. We also measure the projected rotational velocity and revisit the spectral classification. Results. In our sample, several stars were previously considered as spectroscopic binaries, though only a few scattered observations were available. Our results now reveal a more complex situation. Our study identifies two new spectroscopic binaries (HD46149 in NGC2244 and HD46573 in MonOB2). The first object is a long-period double-lined spectroscopic binary, though the exact value of its period remains uncertain and the second object is classified as an SB1 system with a period of about 10.67 days but the time series of our observations do not enable us to derive a unique orbital solution for this system. We also classify another star as variable in radial velocity (HD46150) and we detect line profile variations in two rapid rotators (HD46056 and HD46485). Conclusions. This spectroscopic investigation places a firm lower limit (17%) on the binary fraction of O-stars in NGC2244 and reveals the lack of short-period O+OB systems in this cluster. In addition, a comparison of these new results with two other well-studied clusters (NGC6231 and IC1805) puts forward possible hints of a relation between stellar density and binarity, which could provide constraints on the theories about the formation and early evolution of hot stars.

In response to ESO's call for proposals for second generation instruments for the Very Large Telescope Interferometer (VLTI), a consortium is currently developing the VLTI-Spectro-Imager (VSI). In the context of the Phase A study, a science group has prepared a science case taking advantage of the expected performances of VSI. Among several science topics, the case of dust making Wolf-Rayet binaries producing the so-called pinwheel nebulae has been considered. Here, we review the main specifications of VSI, and we provide preliminary results expected to illustrate the imaging capabilities of VSI, and the interest for the study of pinwheel nebulae similar to those formed close to well-known systems such as WR98a and WR104.

Massive binaries are crucial laboratories that allow us to investigate processes occurring in quite extreme conditions, such as particle acceleration, high-energy emission, or even dust formation. All these processes are intimately dependent on binarity. Our understanding of the underlying physics strongly requires preliminary multiplicity studies likely to uncover still undetected binaries, and determine their orbital parameters. However, classical spectroscopic approaches sometimes fail to provide a solution to this issue. Long baseline interferometry turns out to be a promising complementary technique to address the question of the multiplicity of massive stars. A campaign has been initiated with the VLTI to take benefit of this technique.

Several colliding-wind massive binaries are known to be non-thermal emitters in the radio domain. This constitutes strong evidence for the fact that an efficient particle acceleration process is at work in these objects. The acceleration mechanism is most probably the Diffusive Shock Acceleration (DSA) process in the presence of strong hydrodynamic shocks due to the colliding-winds. In order to investigate the physics of this particle acceleration, we initiated a multiwavelength campaign covering a large part of the electromagnetic spectrum. In this context, the detailed study of the hard X-ray emission from these sources in the SIMBOL-X bandpass constitutes a crucial element in order to probe this still poorly known topic of astrophysics. It should be noted that colliding-wind massive binaries should be considered as very valuable targets for the investigation of particle acceleration in a similar way as supernova remnants, but in a different region of the parameter space.

The properties of the early-type binary Cyg OB2 #5 have been debated for many years and spectroscopic and photometric investigations yielded conflicting results. We have attempted to constrain the physical properties of the binary by collecting new optical and X-ray observations. We find that the orbital period of the system slowly changes though we are unable to discriminate between several possible explanations of this trend. The best fit solution of the continuum light curve reveals a contact configuration with the secondary star being significantly brighter and hotter on its leading side facing the primary. The mean temperature of the secondary star turns out to be only slightly lower than that of the primary, whilst the bolometric luminosity ratio is found to be 3.1. The solution of the light curve yields a distance of 925 +/- 25 pc much lower than the usually assumed distance of the Cyg OB2 association. Whilst we confirm the existence of episodes of higher X-ray fluxes, the data reveal no phase-locked modulation with the 6.6 day period of the eclipsing binary nor any clear relation between the X-ray flux and the 6.7 yr radio cycle. The bright region of the secondary star is probably heated by energy transfer in a common envelope in this contact binary system as well as by the collision with the primary's wind. The existence of a common photosphere probably also explains the odd mass-luminosity relation of the stars in this system. Most of the X-ray, non-thermal radio, and possibly gamma-ray emission of Cyg OB2 #5 is likely to arise from the interaction of the combined wind of the eclipsing binary with at least one additional star of this multiple system.

Classically, optical and near-infrared interferometry have relied on closure phase techniques to produce images. Such techniques allow us to achieve modest dynamic ranges. In order to test the feasibility of next generation optical interferometers in the context of the VLTI-spectro-imager (VSI), we have embarked on a study of image reconstruction and analysis. Our main aim was to test the influence of the number of telescopes, observing nights and distribution of the visibility points on the quality of the reconstructed images. Our results show that observations using six Auxiliary Telescopes (ATs) during one complete night yield the best results in general and is critical in most science cases; the number of telescopes is the determining factor in the image reconstruction outcome. In terms of imaging capabilities, an optical, six telescope VLTI-type configuration and ~200 meter baseline will achieve 4 mas spatial resolution, which is comparable to ALMA and almost 50 times better than JWST will achieve at 2.2 microns. Our results show that such an instrument will be capable of imaging, with unprecedented detail, a plethora of sources, ranging from complex stellar surfaces to microlensing events.

One of the aims of next generation optical interferometric instrumentation is to be able to make use of information contained in the visibility phase to construct high dynamic range images. Radio and optical interferometry are at the two extremes of phase corruption by the atmosphere. While in radio it is possible to obtain calibrated phases for the science objects, in the optical this is currently not possible. Instead, optical interferometry has relied on closure phase techniques to produce images. Such techniques allow only to achieve modest dynamic ranges. However, with high contrast objects, for faint targets or when structure detail is needed, phase referencing techniques as used in radio interferometry, should theoretically achieve higher dynamic ranges for the same number of telescopes. Our approach is not to provide evidence either for or against the hypothesis that phase referenced imaging gives better dynamic range than closure phase imaging. Instead we wish to explore the potential of this technique for future optical interferometry and also because image reconstruction in the optical using phase referencing techniques has only been performed with limited success. We have generated simulated, noisy, complex visibility data, analogous to the signal produced in radio interferometers, using the VLTI as a template. We proceeded with image reconstruction using the radio image reconstruction algorithms contained in AIPS IMAGR (CLEAN algorithm). Our results show that image reconstruction is successful in most of our science cases, yielding images with a 4 milliarcsecond resolution in K band. (abridged)

We review the properties of the very young (2 Myr) open cluster NGC6383. The cluster is dominated by the massive binary HD159176 (O7V + O7V). The distance to NGC6383 is consistently found to be 1.3 +- 0.1 kpc and the average reddening is determined to be E(B-V) = 0.32 +- 0.02. Several pre-main sequence candidates have been identified using different criteria relying on the detection of emission lines, infrared excesses, photometric variability and X-ray emission.

We present here the analysis of an extensive set of high resolution optical spectra of Plaskett's star (HD 47129). We use a disentangling method to separate the individual spectra of each star. We derive a new orbital solution and discuss the spectral classification of both components. A Doppler tomography technique applied to the emission lines H alpha and He II 4686 yields a Doppler map that illustrates the wind interactions in the system. Finally, an atmosphere code is used to determine the different chemical abundances of the system components and the wind parameters. HD 47129 appears to be an O8 III/I + O7.5 III binary system in a post RLOF evolutionary stage, where matter has been transferred from the primary to the secondary star. The He overabundance of the secondary supports this scenario. In addition, the N overabundance and C underabundance of the primary component confirm previous results based on X-ray spectroscopy and indicate that the primary is an evolved massive star. Furthermore, the secondary star has a large rotational velocity that deforms its surface, leading to a non-uniform distribution in effective temperature. This could explain the variations in the equivalent widths of the secondary lines with phase. We suggest that the wind of the secondary star is confined near the equatorial plane because of its high rotational velocity, affecting the ram pressure equilibrium in the wind interaction zone.

The VLTI Spectro Imager (VSI) was proposed as a second-generation instrument of the Very Large Telescope Interferometer providing the ESO community with spectrally-resolved, near-infrared images at angular resolutions down to 1.1 milliarcsecond and spectral resolutions up to R=12000. Targets as faint as K=13 will be imaged without requiring a brighter nearby reference object. The unique combination of high-dynamic-range imaging at high angular resolution and high spectral resolution enables a scientific program which serves a broad user community and at the same time provides the opportunity for breakthroughs in many areas of astrophysic including: probing the initial conditions for planet formation in the AU-scale environments of young stars; imaging convective cells and other phenomena on the surfaces of stars; mapping the chemical and physical environments of evolved stars, stellar remnants, and stellar winds; and disentangling the central regions of active galactic nuclei and supermassive black holes. VSI will provide these new capabilities using technologies which have been extensively tested in the past and VSI requires little in terms of new infrastructure on the VLTI. At the same time, VSI will be able to make maximum use of new infrastructure as it becomes available; for example, by combining 4, 6 and eventually 8 telescopes, enabling rapid imaging through the measurement of up to 28 visibilities in every wavelength channel within a few minutes. The current studies are focused on a 4-telescope version with an upgrade to a 6-telescope one. The instrument contains its own fringe tracker and tip-tilt control in order to reduce the constraints on the VLTI infrastructure and maximize the scientific return.

The line profile variability and photometric variability of the O9.5 Vp star HD93521 are examined in order to establish the properties of the non-radial pulsations in this star. Fourier techniques are used to characterize the modulations of the He I 5876, 6678 and H-alpha lines in several spectroscopic time series and to search for variations in a photometric time series. Our spectroscopic data confirm the existence of two periods of 1.75 and 2.89 hr. The line profiles, especially those affected by emission wings, exhibit also modulations on longer time scales, but these are epoch-dependent and change from line to line. Unlike previous claims, we find no unambiguous signature of the rotational period in our data, nor of a third pulsation period (corresponding to a frequency of 2.66 day$^{-1}$). HD 93521 very likely exhibits non-radial pulsations with periods of 1.75 and 2.89 hr with $l \simeq 8 \pm 1$ and $l \simeq 4 \pm 1$ respectively. No significant signal is found in the first harmonics of these two periods. The 2.89 hr mode is seen at all epochs and in all lines investigated, while the visibility of the 1.75 hr mode is clearly epoch dependent. Whilst light variations are detected, their connection to these periodicities is not straightforward.

Aims: Non-thermal radio emission associated with massive stars is believed to arise from a wind-wind collision in a binary system. However, the evidence of binarity is still lacking in some cases, notably Cyg OB2 #9 Methods: For several years, we have been monitoring this heavily-reddened star from various observatories. This campaign allowed us to probe variations both on short and long timescales and constitutes the first in-depth study of the visible spectrum of this object. Results: Our observations provide the very first direct evidence of a companion in Cyg OB2 #9, confirming the theoretical wind-wind collision scenario. These data suggest a highly eccentric orbit with a period of a few years, compatible with the 2yr-timescale measured in the radio range. In addition, the signature of the wind-wind collision is very likely reflected in the behaviour of some emission lines.

Nowadays, compact sources like surfaces of nearby stars, circumstellar environments of stars from early stages to the most evolved ones and surroundings of active galactic nuclei can be investigated at milli-arcsecond scales only with the VLT in its interferometric mode. We propose a spectro-imager, named VSI (VLTI spectro-imager), which is capable to probe these sources both over spatial and spectral scales in the near-infrared domain. This instrument will provide information complementary to what is obtained at the same time with ALMA at different wavelengths and the extreme large telescopes.

In this paper, I present a general discussion of several astrophysical processes likely to play a role in the production of non-thermal emission in massive stars, with emphasis on massive binaries. Even though the discussion will start in the radio domain where the non-thermal emission was first detected, the census of physical processes involved in the non-thermal emission from massive stars shows that many spectral domains are concerned, from the radio to the very high energies. First, the theoretical aspects of the non-thermal emission from early-type stars will be addressed. The main topics that will be discussed are respectively the physics of individual stellar winds and their interaction in binary systems, the acceleration of relativistic electrons, the magnetic field of massive stars, and finally the non-thermal emission processes relevant to the case of massive stars. Second, this general qualitative discussion will be followed by a more quantitative one, devoted to the most probable scenario where non-thermal radio emitters are massive binaries. I will show how several stellar, wind and orbital parameters can be combined in order to make some semi-quantitative predictions on the high-energy counterpart to the non-thermal emission detected in the radio domain. These theoretical considerations will be followed by a census of results obtained so far, and related to this topic... (see paper for full abstract)

An increasing number of early-type (O and Wolf-Rayet) colliding wind binaries (CWBs) is known to accelerate particles up to relativistic energies. In this context, non-thermal emission processes such as inverse Compton (IC) scattering are expected to produce a high energy spectrum, in addition to the strong thermal emission from the shock-heated plasma. SIMBOL-X will be the ideal observatory to investigate the hard X-ray spectrum (above 10 keV) of these systems, i.e. where it is no longer dominated by the thermal emission. Such observations are strongly needed to constrain the models aimed at understanding the physics of particle acceleration in CWB. Such systems are important laboratories for investigating the underlying physics of particle acceleration at high Mach number shocks, and probe a different region of parameter space than studies of supernova remnants.

We present a spectroscopic analysis of four massive binary systems that are known or are good candidates to display the Struve-Sahade effect (defined as the apparent strengthening of the secondary spectrum of the binary when the star is approaching, and the corresponding weakening of the lines when it is receding). We use high resolution optical spectra to determine new orbital solutions and spectral types of HD 165052, HD 100213, HD 159176 and DH Cep. As good knowledge of the fundamental parameters of the considered systems is necessary to examine the Struve-Sahade effect. We then study equivalent width variations in the lines of both components of these binaries during their orbital cycle. In the case of these four systems, variations appear in the equivalent widths of some lines during the orbital cycle, but the definition given above can any longer be valid, since it is now clear that the effect modifies the primary spectrum as much as the secondary spectrum. Furthermore, the lines affected, and the way in which they are affected, depend on the considered system. For at least two of them (HD 100213 and HD 159176) these variations probably reflect the ellipsoidal variable nature of the system.

Aims: We analyze INTEGRAL-ISGRI data in order to probe the hard X-ray emission (above 20 keV) from point sources in the Cyg OB2 region and to investigate the putative non-thermal high-energy emission from early-type stars (Wolf-Rayet and O-type stars). Among the targets located in the field of view, we focus on the still unidentified EGRET source 3EG 2033+4118 that may be related to massive stars known to produce non-thermal emission in the radio domain, and on the wide colliding-wind binary WR 140. Methods: Using a large set of data obtained with the IBIS-ISGRI imager onboard INTEGRAL, we run the OSA software package in order to find point sources in the fully coded field of view of the instrument. Results: Our data do not allow the detection of a lower-energy counterpart of 3EG J2033+4118 nor of any other new point sources in the field of view, and we derive upper limits on the high-energy flux for a few targets: 3EG J2033+4118, TeV J2032+4130, WR140, WR146 and WR147. The results are discussed in the context of the multiwavelength investigation of these objects. Conclusions: The upper limits derived are valuable constraints for models aimed at understanding the acceleration of particles in non-thermal emitting massive stars, and of the still unidentified very-high gamma-ray source TeV J2032+4130.

Aims. The properties of the early-type stars in the core of the Westerlund2 cluster are examined in order to establish a link between the cluster and the very massive Wolf-Rayet binary WR20a as well as the HII complex RCW49. Methods. Photometric monitoring as well as spectroscopic observations of Westerlund2 are used to search for light variability and to establish the spectral types of the early-type stars in the cluster core. Results. The first light curves of the eclipsing binary WR20a in B and V filters are analysed and a distance of 8kpc is inferred. Three additional eclipsing binaries, which are probable late O or early B-type cluster members, are discovered, but none of the known early O-type stars in the cluster displays significant photometric variability above 1% at the 1-sigma level. The twelve brightest O-type stars are found to have spectral types between O3 and O6.5, significantly earlier than previously thought. Conclusions. The distance of the early-type stars in Westerlund2 is established to be in excellent agreement with the distance of WR20a, indicating that WR20a actually belongs to the cluster. Our best estimate of the cluster distance thus amounts to 8.0pm1.4kpc. Despite the earlier spectral types, the currently known population of early-type stars in Westerlund2 does not provide enough ionizing photons to account for the radio emission of the RCW49 complex. This suggests that there might still exist a number of embedded early O-stars in RCW49.

HD 167971 is a triple system consisting of a 3.3-day eclipsing binary (O5-8 V + O5-8 V) and an O8 supergiant. It is also a well known non-thermal radio emitter. We observed the radio emission of HD 167971 with the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). By combining these data with VLA archive observations we constructed a radio lightcurve covering a 20-yr time-range. We searched for, but failed to find, the 3.3-day spectroscopic period of the binary in the radio data. This could be due to the absence of intrinsic synchrotron radiation at the colliding-wind region between the two components of the eclipsing binary, or due to the large amount of free-free absorption that blocks the synchrotron radiation. We are able to explain many of the observed characteristics of the radio data if the non-thermal emission is produced in a colliding-wind region between the supergiant and the combined winds of the binary. Furthermore, if the system is gravitationally bound, the orbital motion occurs over a period of ~ 20 years, or longer, as suggested by the long-term variability in the radio data. We argue that the variability is due to the free-free absorption that changes with orbital phase or may also in part be due to changes in separation, should the orbit be eccentric.

Aims: We address the issue of the multiplicity of the three brightest early-type stars of the young open cluster IC1805, namely HD15570, HD15629 and HD15558. Methods: For the three stars, we measured the radial velocity by fitting Gaussian curves to line profiles in the optical domain. In the case of the massive binary HD15558, we also used a spectral disentangling method to separate the spectra of the primary and of the secondary in order to derive the radial velocities of the two components. These measurements were used to compute orbital solutions for HD15558. Results: For HD15570 and HD15629, the radial velocities do not present any significant trend attributable to a binary motion on time scales of a few days, nor from one year to the next. In the case of HD15558 we obtained an improved SB1 orbital solution with a period of about 442 days, and we report for the first time on the detection of the spectral signature of its secondary star. We derive spectral types O5.5III(f) and O7V for the primary and the secondary of HD15558. We tentatively compute a first SB2 orbital solution although the radial velocities from the secondary star should be considered with caution. The mass ratio is rather high, i.e. about 3, and leads to very extreme minimum masses, in particular for the primary object. Minimum masses of the order of 150 \pm 50 and 50 \pm 15 M_\odot are found respectively for the primary and the secondary. Conclusions: We propose that HD15558 could be a triple system. This scenario could help to reconcile the very large minimum mass derived for the primary object with its spectral type. In addition, considering new and previously published results, we find that the binary frequency among O-stars in IC1805 has a lower limit of 20%, and that previously published values (80%) are probably overestimated.

We discuss the results of two XMM-Newton observations of the open cluster NGC6604 obtained in April and September 2002. We concentrate mainly on the multiple system HD167971 (O5-8V + O5-8V + (O8I)). The soft part of the EPIC spectrum of this system is thermal with typical temperatures of about 2 10^6 to 9 10^6 K. The nature (thermal vs non-thermal) of the hard part of the spectrum is not unambiguously revealed by our data. If the emission is thermal, the high temperature of the plasma (~ 2.3 10^7 to 4.6 10^7 K) would be typical of what should be expected from a wind-wind interaction zone within a long period binary system. This emission could arise from an interaction between the combined winds of the O5-8V + O5-8V close binary system and that of the more distant O8I companion. Assuming instead that the hard part of the spectrum is non-thermal, the photon index would be rather steep (~3). Moreover, a marginal variability between our two XMM-Newton pointings could be attributed to an eclipse of the O5-8V + O5-8V system. The overall X-ray luminosity points to a significant X-ray luminosity excess of about a factor 4 possibly due to colliding winds. Considering HD167971 along with several recent X-ray and radio observations, we propose that the simultaneous observation of non-thermal radiation in the X-ray (below 10.0 keV) and radio domains appears rather unlikely. Our investigation of our XMM-Newton data of NGC6604 reveals a rather sparse distribution of X-ray emitters. Including the two bright non-thermal radio emitters HD168112 and HD167971, we present a list of 31 X-ray sources... (see paper for full abstract)

We present a radio lightcurve of the O5.5 III(f+) star HD 168112, based on archive data from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). The fluxes show considerable variability and a negative spectral index, thereby confirming that HD 168112 is a non-thermal radio emitter. The non-thermal radio emission is believed to be due to synchrotron radiation from relativistic electrons that have been Fermi accelerated in shocks. For HD 168112, it is not known whether these shocks are due to a wind-wind collision in a binary system or to the intrinsic instability of the stellar wind driving mechanism. Assuming HD 168112 to be a single star, our synchrotron model shows that the velocity jump of the shocks should be very high, or there should be a very large number of shocks in the wind. Neither of these is compatible with time-dependent hydrodynamical calculations of O star winds. If, on the other hand, we assume that HD 168112 is a binary, the high velocity jump is easily explained by ascribing it to the wind-wind collision. By further assuming the star to be an eccentric binary, we can explain the observed radio variability by the colliding-wind region moving in and out of the region where free-free absorption is important. The radio data presented here show that the binary has a period between one and two years. By combining the radio data with X-ray data, we find that the most likely period is about 1.4 yr.

We present results from a radial velocity study of two runaway O-type stars, HD 14633 (ON8.5 V) and HD 15137 (O9.5 III(n)). We find that HD 14633 is a single-lined spectroscopic binary with an orbital period of 15.4083 days. The second target HD 15137 is a radial velocity variable and a possible single-lined spectroscopic binary with a period close to 1 month. Both binaries have large eccentricity, small semiamplitude, and a small mass function. We show the trajectories of the stars in the sky based upon an integration of motion in the Galactic potential, and we suggest that both stars were ejected from the vicinity of the open cluster NGC 654 in the Perseus spiral arm. The binary orbital parameters and runaway velocities are consistent with the idea that both these stars were ejected by supernova explosions in binaries and that they host neutron star companions. We find that the time-of-flight since ejection is longer than the predicted evolutionary timescales for the stars, which may indicate that they have a lower mass than normally associated with their spectral classifications and/or that their lives have been extended through rapid rotation.

We present the very first analysis of the spectroscopic variability of the three rapidly rotating Oef stars HD192281 (O5(ef)), HD14442 (O5.5ef) and HD14434 (O6.5(ef)). Radial velocities of the HeII lambda 4541 line reveal no evidence of binarity on time scales of a few days, or from one year to the next, for any of the targets. The HeII lambda 4686 double-peaked emission and, to some extent, the Hbeta absorption line display significant profile variability in the spectra of all three stars. Data gathered during different observing runs spread over six years reveal a rather stable time scale for HD192281 and HD14442, whereas the variability pattern changes significantly from one year to the other. The case of HD14434 is less clear as no obvious time scale emerges from our analysis. In a temptative way to interpret this variability, stellar rotation remains a possible clock for HD192281 and HD14442. However, currently available models addressing stellar rotation fail to explain some crucial aspects of the observed variability behaviour, which appear to be even more complex in the case of HD14434.

We present the results of a spectroscopic campaign revealing that the non-thermal radio emitter Cyg OB2 #8A is an O6 + O5.5 binary system. We propose the very first orbital solution indicating a period of about 21.9 days. The system appears to be eccentric (0.24 +/- 0.04) and is likely seen under a rather low inclination angle. The mass ratio of the components is close to unity. The impact of the binarity of this star in the framework of our understanding of non-thermal radio emission from early-type stars is briefly discussed.

We present the results of an intensive spectroscopic campaign in the optical waveband revealing that Cyg OB2 #8A is an O6 + O5.5 binary system with a period of about 21.9 d. Cyg OB2 #8A is a bright X-ray source, as well as a non-thermal radio emitter. We discuss the binarity of this star in the framework of a campaign devoted to the study of non-thermal emitters, from the radio waveband to gamma-rays. In this context, we attribute the non-thermal radio emission from this star to a population of relativistic electrons, accelerated by the shock of the wind-wind collision. These relativistic electrons could also be responsible for a putative gamma-ray emission through inverse Compton scattering of photospheric UV photons, thus contributing to the yet unidentified EGRET source 3EG J2033+4118.

We report the results of a multiwavelength study of the non-thermal radio emitter HD168112 (O5.5III(f^+)). The detailed analysis of two quasi-simultaneous XMM-Newton and VLA observations reveals strong variability of this star both in the X-ray and radio ranges. The X-ray observations separated by five months reveal a decrease of the X-ray flux of ~30%. The radio emission on the other hand increases by a factor 5-7 between the two observations obtained roughly simultaneously with the XMM-Newton pointings. The X-ray data reveal a hard emission that is most likely produced by a thermal plasma at kT ~2-3 keV while the VLA data confirm the non-thermal status of this star in the radio waveband. Comparison with archive X-ray and radio data confirms the variability of this source in both wavelength ranges over a yet ill defined time scale. The properties of HD168112 in the X-ray and radio domain point towards a binary system with a significant eccentricity and an orbital period of a few years. However, our optical spectra reveal no significant changes of the star's radial velocity suggesting that if HD168112 is indeed a binary, it must be seen under a fairly low inclination.

We analyse spectroscopic observations of WR20a revealing that this star is a massive early-type binary system with a most probable orbital period of \sim 3.675 days. Our spectra indicate that both components are most likely of WN6ha or O3If^*/WN6ha spectral type. The orbital solution for a period of 3.675 days yields extremely large minimum masses of 70.7 \pm 4.0 and 68.8 \pm 3.8 M_{\odot} for the two stars. These properties make WR20a a cornerstone system for the study of massive star evolution.

We report the analysis of an XMM-Newton observation of the close binary HD159176 (O7V + O7V). The observed L_X/L_bol ratio reveals an X-ray luminosity exceeding by a factor ~ 7 the expected value for X-ray emission from single O-stars, therefore suggesting a wind-wind interaction scenario. EPIC and RGS spectra are fitted consistently with a two temperature mekal optically thin thermal plasma model, with temperatures ranging from ~ 2 to 6 10^6 K. At first sight, these rather low temperatures are consistent with the expectations for a close binary system where the winds collide well before reaching their terminal velocities. We also investigate the variability of the X-ray light curve of HD159176 on various short time scales. No significant variability is found and we conclude that if hydrodynamical instabilities exist in the wind interaction region of HD159176, they are not sufficient to produce an observable signature in the X-ray emission. Hydrodynamic simulations using wind parameters from the literature reveal some puzzling discrepancies. The most striking one concerns the predicted X-ray luminosity which is one or more orders of magnitude larger than the observed one. A significant reduction of the mass loss rate of the components compared to the values quoted in the literature alleviates the discrepancy but is not sufficient to fully account for the observed luminosity. Because hydrodynamical models are best for the adiabatic case whereas the colliding winds in HD159176 are most likely highly radiative, a totally new approach has been envisaged... (see paper for complete abstract)

In 2001, using a large spectroscopic dataset from an extensive monitoring campaign, we discovered that the peculiar Of star HD108 displayed spectacular line variations. This strange behaviour could be attributed to a variety of models, and an investigation of the high energy properties of HD108 was crucially needed to test the predictions from these models. Our dedicated XMM-Newton observation of HD108 shows that its spectrum is well represented by a two temperature thermal plasma model with kT_1\sim0.2 keV and kT_2sim1.4 keV. In addition, we find that the star does not display any significant short-term changes during the XMM-Newton exposure. Compared to previous Einstein and ROSAT detections, it also appears that HD108 does not present long-term flux variations either. While the dramatic line variations continue to modify HD108's spectrum in the optical domain, the X-ray emission of the star appears thus surprisingly stable: no simple model is for the moment able to explain such an unexpected behaviour. Thanks to its high sensitivity, the XMM-Newton observatory has also enabled the serendipitous discovery of 57 new X-ray sources in the field of HD108. Their properties are also discussed in this paper.