The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

Recent and on-going large ground-based multi-object spectroscopic surveys allow to significantly increase the sample of spectroscopic binaries to get insight into their statistical properties. We investigate the repeated spectral observations of the Gaia-ESO Survey (GES) internal data release 5 to identify and characterize spectroscopic binaries with one visible component (SB1) in fields covering the discs, the bulge, the CoRot fields, and stellar clusters and associations. A statistical chi2-test is performed on spectra of the iDR5 sub-sample of approximately 43500 stars characterized by at least 2 observations and a S/N > 3. Our sample of RV variables is cleaned from contamination by pulsation/convection-induced variables using Gaia DR2 parallaxes and photometry. Monte-Carlo simulations using the SB9 catalogue of spectroscopic orbits allow to estimate our detection efficiency and to correct the SB1 rate to evaluate the GES SB1 binary fraction and its dependence with effective temperature and metallicity. We find 641 (resp., 803) FGK SB1 candidates at the 5 sigma (resp., 3 sigma) level. The orbital-period distribution is estimated from the RV standard-deviation distribution and reveals that the detected SB1 probe binaries with log(P[d]) < 4. We estimate the global GES SB1 fraction to be in the range 7-14% with a typical uncertainty of 4%. The GES SB1 frequency decreases with metallicity at a rate of -9+/-3%/dex in the metallicity range -2.7<FeH<+0.6. This anticorrelation is obtained with a confidence level higher than 93% on a homogeneous sample covering spectral types FGK and a large range of metallicities. When the present-day mass function is accounted for, this rate turns to 4+/-2%/dex with a confidence level higher than 88%. In addition we provide the variation of the SB1 fraction with metallicity separately for F, G, and K spectral types, as well as for dwarf and giant primaries.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

This paper presents the specification, design, and development of the Radial Velocity Spectrometer (RVS) on the European Space Agency's Gaia mission. Starting with the rationale for the full six dimensions of phase space in the dynamical modelling of the Galaxy, the scientific goals and derived top-level instrument requirements are discussed, leading to a brief description of the initial concepts for the instrument. The main part of the paper is a description of the flight RVS, considering the optical design, the focal plane, the detection and acquisition chain, and the as-built performance drivers and critical technical areas. After presenting the pre-launch performance predictions, the paper concludes with the post-launch developments and mitigation strategies, together with a summary of the in-flight performance at the end of commissioning.

The Gaia Data Release 2 contains the 1st release of radial velocities complementing the kinematic data of a sample of about 7 million relatively bright, late-type stars. Aims: This paper provides a detailed description of the Gaia spectroscopic data processing pipeline, and of the approach adopted to derive the radial velocities presented in DR2. Methods: The pipeline must perform four main tasks: (i) clean and reduce the spectra observed with the Radial Velocity Spectrometer (RVS); (ii) calibrate the RVS instrument, including wavelength, straylight, line-spread function, bias non-uniformity, and photometric zeropoint; (iii) extract the radial velocities; and (iv) verify the accuracy and precision of the results. The radial velocity of a star is obtained through a fit of the RVS spectrum relative to an appropriate synthetic template spectrum. An additional task of the spectroscopic pipeline was to provide 1st-order estimates of the stellar atmospheric parameters required to select such template spectra. We describe the pipeline features and present the detailed calibration algorithms and software solutions we used to produce the radial velocities published in DR2. Results: The spectroscopic processing pipeline produced median radial velocities for Gaia stars with narrow-band near-IR magnitude Grvs < 12 (i.e. brighter than V~13). Stars identified as double-lined spectroscopic binaries were removed from the pipeline, while variable stars, single-lined, and non-detected double-lined spectroscopic binaries were treated as single stars. The scatter in radial velocity among different observations of a same star, also published in DR2, provides information about radial velocity variability. For the hottest (Teff > 7000 K) and coolest (Teff < 3500 K) stars, the accuracy and precision of the stellar parameter estimates are not sufficient to allow selection of appropriate templates. [Abridged]

For Gaia DR2 (GDR2), 280 million spectra, collected by the RVS instrument on-board Gaia, were processed and median radial velocities were derived for 9.8 million sources brighter than Grvs = 12 mag. This paper describes the validation and properties of the median radial velocities published in GDR2. Quality tests and filters are applied to select, from the 9.8 million radial velocities, those with the quality to be published in GDR2. The accuracy of the selected sample is assessed with respect to ground-based catalogues. Its precision is estimated using both ground-based catalogues and the distribution of the Gaia radial velocity uncertainties. GDR2 contains median radial velocities for 7 224 631 stars, with Teff in the range [3550, 6900] K, which passed succesfully the quality tests. The published median radial velocities provide a full sky-coverage and have a completness with respect to the astrometric data of 77.2\% (for $G \leq 12.5$ mag). The median radial velocity residuals with respect to the ground-based surveys vary from one catalogue to another, but do not exceed a few 100s m/s. In addition, the Gaia radial velocities show a positive trend as a function of magnitude, which starts around Grvs $\sim 9$ mag and reaches about $+500$ m/s at Grvs $= 11.75$ mag. The overall precision, estimated from the median of the Gaia radial velocity uncertainties, is 1.05 km/s. The radial velocity precision is function of many parameters, in particular the magnitude and effective temperature. For bright stars, Grvs in [4, 8] mag, the precision is in the range 200-350 m/s, which is about 3 to 5 times more precise than the pre-launch specification of 1 km/s. At the faint end, Grvs = 11.75 mag, the precisions for Teff = 5000 K and 6500 K are respectively 1.4 km/s and 3.7 km/s.

Aims. The Radial Velocity Spectrometer (RVS) on board the ESA satellite mission Gaia has no calibration device. Therefore, the radial velocity zero point needs to be calibrated with stars that are proved to be stable at a level of 300 m/s during the Gaia observations. Methods. We compiled a dataset of ~71000 radial velocity measurements from five high-resolution spectrographs. A catalogue of 4813 stars was built by combining these individual measurements. The zero point was established using asteroids. Results. The resulting catalogue has seven observations per star on average on a typical time baseline of six years, with a median standard deviation of 15 m/s. A subset of the most stable stars fulfilling the RVS requirements was used to establish the radial velocity zero point provided in Gaia Data Release 2. The stars that were not used for calibration are used to validate the RVS data.

In the colliding-wind region of massive binaries, non-thermal radio emission occurs. This non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimetrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 +- 0.033 mJy. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer.

The Gaia-ESO Survey (GES) is a large spectroscopic survey that provides a unique opportunity to study the distribution of spectroscopic multiple systems among different populations of the Galaxy. We aim at detecting binarity/multiplicity for stars targeted by the GES from the analysis of the cross-correlation functions (CCFs) of the GES spectra with spectral templates. We develop a method based on the computation of the CCF successive derivatives to detect multiple peaks and determine their radial velocities, even when the peaks are strongly blended. The parameters of the detection of extrema (DOE) code have been optimized for each GES GIRAFFE and UVES setup to maximize detection. This code therefore allows to automatically detect multiple line spectroscopic binaries (SBn, n>1). We apply this method on the fourth GES internal data release and detect 354 SBn candidates (342 SB2, 11 SB3 and even one SB4), including only 9 SB2 known in the literature. This implies that about 98% of these SBn candidates are new (because of their faint visual magnitude that can reach V=19). Visual inspection of the SBn candidate spectra reveals that the most probable candidates have indeed a composite spectrum. Among SB2 candidates, an orbital solution could be computed for two previously unknown binaries: 06404608+0949173 (known as V642 Mon) in NGC 2264 and 19013257-0027338 in Berkeley 81. A detailed analysis of the unique SB4 (four peaks in the CCF) reveals that HD 74438 in the open cluster IC 2391 is a physically bound stellar quadruple system. The SB candidates belonging to stellar clusters are reviewed in detail to discard false detections. We warn against the use of atmospheric parameters for these system components rather than by SB-specific pipelines. Our implementation of an automatic detection of spectroscopic binaries within the GES has allowed an efficient discovery of many new multiple systems.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4).

Gaia is a space mission currently measuring the five astrometric parameters as well as spectrophotometry of at least 1 billion stars to G = 20.7 mag with unprecedented precision. The sixth parameter in phase space (radial velocity) is also measured thanks to medium-resolution spectroscopy being obtained for the 150 million brightest stars. During the commissioning phase, two fields, one around each ecliptic pole, have been repeatedly observed to assess and to improve the overall satellite performances as well as the associated reduction and analysis software. A ground-based photometric and spectroscopic survey was therefore initiated in 2007, and is still running in order to gather as much information as possible about the stars in these fields. This work is of particular interest to the validation of the Radial Velocity Spectrometer (RVS) outputs. The paper presents the radial velocity measurements performed for the Southern targets in the 12 - 17 R magnitude range on high- to mid-resolution spectra obtained with the GIRAFFE and UVES spectrographs.

Mass-loss in massive stars plays a critical role in their evolution, although the precise mechanism(s) responsible - radiatively driven winds, impulsive ejection and/or binary interaction -remain uncertain. In this paper we present ALMA line and continuum observations of the supergiant B[e] star Wd1-9, a massive post-Main Sequence object located within the starburst cluster Westerlund 1. We find it to be one of the brightest stellar point sources in the sky at millimetre wavelengths, with (serendipitously identified) emission in the H41alpha radio recombination line. We attribute these properties to a low velocity (~100 km/s) ionised wind, with an extreme mass-loss rate 6.4x10^-5(d/5kpc)^1.5 Msol/yr. External to this is an extended aspherical ejection nebula indicative of a prior phase of significant mass-loss. Taken together, the millimetre properties of Wd1-9 show a remarkable similarity to those of the highly luminous stellar source MWC349A.We conclude that these objects are interacting binaries evolving away from the main sequence and undergoing rapid case-A mass transfer. As such they - and by extension the wider class of supergiant B[e] stars - may provide a unique window into the physics of a process that shapes the life-cycle of ~70% of massive stars found in binary systems.

We present e-MERLIN 21 cm (L-band) observations of single luminous OB stars in the Cygnus OB2 association, from the COBRaS Legacy programme. The radio observations potentially offer the most straightforward, least model-dependent, determinations of mass-loss rates, and can be used to help resolve current discrepancies in mass-loss rates via clumped and structured hot star winds. We report here that the 21 cm flux densities of O3 to O6 supergiant and giant stars are less than ~ 70 microJy. These fluxes may be translated to `smooth' wind mass-loss upper limits of ~ 4.4 - 4.8 x 10^(-6) M_sol/yr for O3 supergiants and < 2.9 x 10^(-6) M_sol/yr for B0 to B1 supergiants. The first ever resolved 21 cm detections of the hypergiant (and LBV candidate) Cyg OB2 #12 are discussed; for multiple observations separated by 14 days, we detect a ~ 69% increase in its flux density. Our constraints on the upper limits for the mass-loss rates of evolved OB stars in Cyg OB2 support the model that the inner wind region close to the stellar surface (where Halpha forms) is more clumped than the very extended geometric region sampled by our radio observations.

The long-period, highly eccentric O-star binary 9 Sgr, known for its non-thermal radio emission and its relatively bright X-ray emission, went through its periastron in 2013. Such an event can be used to observationally test the predictions of the theory of colliding stellar winds over a broad range of wavelengths. We have conducted a multi-wavelength monitoring campaign of 9 Sgr around the 2013 periastron. In this paper, we focus on X-ray observations and optical spectroscopy. The optical spectra allow us to revisit the orbital solution of 9 Sgr and to refine its orbital period to 9.1 years. The X-ray flux is maximum at periastron over all energy bands, but with clear differences as a function of energy. The largest variations are observed at energies above 2 keV, whilst the spectrum in the soft band (0.5 - 1.0 keV) remains mostly unchanged indicating that it arises far from the collision region, in the inner winds of the individual components. The level of the hard emission at periastron clearly deviates from the 1/r relation expected for an adiabatic wind interaction zone, whilst this relation seems to hold at the other phases covered by our observations. The spectra taken at phase 0.946 reveal a clear Fe xxv line at 6.7 keV, but no such line is detected at periastron (phi = 0.000) although a simple model predicts a strong line that should be easily visible in the data. The peculiarities of the X-ray spectrum of 9 Sgr could reflect the impact of radiative inhibition as well as a phase-dependent efficiency of particle acceleration on the shock properties.

Context. Luminous blue variables (LBVs) are rare massive stars with very high luminosity. They are characterized by strong photo-metric and spectroscopic variability related to transient eruptions. The mechanisms at the origin of these eruptions is not well known. In addition, their formation is still problematic and the presence of a companion could help to explain how they form. Aims. This article presents a study of seven LBVs (about 20% of the known Galactic population), some Wolf-Rayet stars, and massive binaries. We probe the environments that surround these massive stars with near-, mid-, and far-infrared images, investigating potential nebula/shells and the companion stars. Methods. To investigate large spatial scales, we used seeing-limited and near diffraction-limited adaptive optics images to obtain a differential diagnostic on the presence of circumstellar matter and to determine their extent. From those images, we also looked for the presence of binary companions on a wide orbit. Once a companion was detected, its gravitational binding to the central star was tested. Tests include the chance projection probability, the proper motion estimates with multi-epoch observations, flux ratio, and star separations. Results. We find that two out of seven of LBVs may have a wide orbit companion. Most of the LBVs display a large circumstellar envelope or several shells. In particular, HD168625, known for its rings, possesses several shells with possibly a large cold shell at the edge of which the rings are formed. For the first time, we have directly imaged the companion of LBV stars.

This work addresses a procedure to estimate fundamental stellar parameters such as T eff , logg, [Fe/H], and v sin i using a dimensionality reduction technique called Principal Component Analysis (PCA), applied to a large database of synthetic spectra. This technique shows promising results for inverting stellar parameters of observed targets from Gaia ESO Survey.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

This work describes a procedure to derive several fundamental parameters such as the effective temperature, surface gravity, equatorial rotational velocity and microturbulent velocity. In this work, we have written a numerical procedure in Python which finds the best fit between a grid of synthetic spectra and the observed spectra by minimizing a standard chi-square. LTE model atmospheres were calculated using the ATLAS9 code and were used as inputs to the spectrum synthesis code SYNSPEC48 in order to compute a large grid of synthetic Balmer line profiles. This new procedure has been applied to a large number of new observations (GIRAFFE spectra) of B and A stars members of the young open cluster NGC3293. These observations are part of the GAIA ESO Survey. Takeda's procedure was also used to derive rotational velocities and microturbulent velocities. The results have been compared to previous determinations by other authors and are found to agree with them. As a first result, we concluded that using this procedure, an accuracy of +-200 K could be achieved in effective temperature and +-0.2 dex in surface gravities.

Chemically inhomogeneous populations are observed in most globular clusters, but not in open clusters. Cluster mass seems to play a key role in the existence of multiple populations. Studying the chemical homogeneity of the most massive open clusters is necessary to better understand the mechanism of their formation and determine the mass limit under which clusters cannot host multiple populations. Here we studied NGC6705, that is a young and massive open cluster located towards the inner region of the Milky Way. This cluster is located inside the solar circle. This makes it an important tracer of the inner disk abundance gradient. This study makes use of BVI and ri photometry and comparisons with theoretical isochrones to derive the age of NGC6705. We study the density profile of the cluster and the mass function to infer the cluster mass. Based on abundances of the chemical elements distributed in the first internal data release of the Gaia-ESO Survey, we study elemental ratios and the chemical homogeneity of the red clump stars. Radial velocities enable us to study the rotation and internal kinematics of the cluster. The estimated ages range from 250 to 316Myr, depending on the adopted stellar model. Luminosity profiles and mass functions show strong signs of mass segregation. We derive the mass of the cluster from its luminosity function and from the kinematics, finding values between 3700 M$_{\odot}$ and 11 000 M$_{\odot}$. After selecting the cluster members from their radial velocities, we obtain a metallicity of [Fe/H]=0.10$\pm$0.06 based on 21 candidate members. Moreover, NGC6705 shows no sign of the typical correlations or anti-correlations between Al, Mg, Si, and Na, that are expected in multiple populations. This is consistent with our cluster mass estimate, which is lower than the required mass limit proposed in literature to develop multiple populations.

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

The derivation of radial velocities from large numbers of spectra that typically result from survey work, requires automation. However, except for the classical cases of slowly rotating late-type spectra, existing methods of measuring Doppler shifts require fine-tuning to avoid a loss of accuracy due to the idiosyncrasies of individual spectra. The radial velocity spectrometer (RVS) on the Gaia mission, which will start operating very soon, prompted a new attempt at creating a measurement pipeline to handle a wide variety of spectral types. The present paper describes the theoretical background on which this software is based. However, apart from the assumption that only synthetic templates are used, we do not rely on any of the characteristics of this instrument, so our results should be relevant for most telescope-detector combinations. We propose an approach based on the simultaneous use of several alternative measurement methods, each having its own merits and drawbacks, and conveying the spectral information in a different way, leading to different values for the measurement. A comparison or a combination of the various results either leads to a "best estimate" or indicates to the user that the observed spectrum is problematic and should be analysed manually. We selected three methods and analysed the relationships and differences between them from a unified point of view; with each method an appropriate estimator for the individual random error is chosen. We also develop a procedure for tackling the problem of template mismatch in a systematic way. Furthermore, we propose several tests for studying and comparing the performance of the various methods as a function of the atmospheric parameters of the observed objects. Finally, we describe a procedure for obtaining a knowledge-based combination of the various Doppler-shift measurements.

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 10^5 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium resolution (R~20,000) GIRAFFE spectrograph and the high resolution (R~47,000) UVES spectrograph. In this paper, we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, performs automatically sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is sigma~0.4 km s^-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s^-1) and to the radial velocities of the standard stars (~0.5 kms^-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the setups and instruments used for the survey will be established.

A programme to observe all A dwarfs in open clusters brighter than V=6.5 mag of various ages and in the field was initiated several years ago. In this work we present the current status of microturbulent velocity for A and F dwarfs. We have performed high resolution high signal-to-noise spectroscopy of stars well distributed in mass along the Main Sequence. Microturbulent velocities are derived iteratively by fitting grids of synthetic spectra calculated in LTE to observed spectra of 61 A field stars, 55 A and 58 F in open clusters (Pleiades, Coma Berenices, Hyades and the Ursa Major moving group). We compared our results to recent works and found a good agreement with their analytical formulation for the standard microturbulence. Our results show a broad maximum for microturbulent velocities in the range A5V to about A9V and a decrease (to $\sim 1$ km/s) for cooler and hotter stars as indicated in Smalley (2004).We also present a comparison to preliminary science results of Lobel et al. (2013) for the Gaia-ESO Public Spectroscopic Survey.

The colliding winds in a massive binary system generate synchrotron emission due to a fraction of electrons that have been accelerated to relativistic speeds around the shocks in the colliding-wind region. We studied the radio light curve of 9 Sgr = HD 164794, a massive O-type binary with a 9.1-yr period. We investigated whether the radio emission varies consistently with orbital phase and we determined some parameters of the colliding-wind region. We reduced a large set of archive data from the Very Large Array (VLA) to determine the radio light curve of 9 Sgr at 2, 3.6, 6 and 20 cm. We also constructed a simple model that solves the radiative transfer in the colliding-wind region and both stellar winds. The 2-cm radio flux shows clear phase-locked variability with the orbit. The behaviour at other wavelengths is less clear, mainly due to a lack of observations centred on 9 Sgr around periastron passage. The high fluxes and nearly flat spectral shape of the radio emission show that synchrotron radiation dominates the radio light curve at all orbital phases. The model provides a good fit to the 2-cm observations, allowing us to estimate that the brightness temperature of the synchrotron radiation emitted in the colliding-wind region at 2 cm is at least 4 x 10^8 K. The simple model used here already allows us to derive important information about the colliding-wind region. We propose that 9 Sgr is a good candidate for more detailed modelling, as the colliding-wind region remains adiabatic during the whole orbit thus simplifying the hydrodynamics.

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.

Context: Nonthermal radio emission in massive stars is expected to arise in wind-wind collisions occurring inside a binary system. One such case, the O-type star Cyg OB2 #9, was proven to be a binary only four years ago, but the orbital parameters remained uncertain. The periastron passage of 2011 was the first one to be observable under good conditions since the discovery of binarity. Aims: In this context, we have organized a large monitoring campaign to refine the orbital solution and to study the wind-wind collision. Methods: This paper presents the analysis of optical spectroscopic data, as well as of a dedicated X-ray monitoring performed with Swift and XMM. Results: In light of our refined orbital solution, Cyg OB2 #9 appears as a massive O+O binary with a long period and high eccentricity; its components (O5-5.5I for the primary and O3-4III for the secondary) have similar masses and similar luminosities. The new data also provide the first evidence that a wind-wind collision is present in the system. In the optical domain, the broad Ha line varies, displaying enhanced absorption and emission components at periastron. X-ray observations yield the unambiguous signature of an adiabatic collision because, as the stars approach periastron, the X-ray luminosity closely follows the 1/D variation expected in that case. The X-ray spectrum appears, however, slightly softer at periastron, which is probably related to winds colliding at slightly lower speeds at that time. Conclusions: It is the first time that such a variation has been detected in O+O systems, and the first case where the wind-wind collision is found to remain adiabatic even at periastron passage.

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.

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.

Six O-type stars were observed continuously by the CoRoT satellite during a 34.3-day run. The unprecedented quality of the data allows us to detect even low-amplitude stellar pulsations in some of these stars (HD 46202 and the binaries HD 46149 and Plaskett's star). These cover both opacity-driven modes and solar-like stochastic oscillations, both of importance to the asteroseismological modelling of O stars. Additional effects can be seen in the CoRoT light curves, such as binarity and rotational modulation. Some of the hottest O-type stars (HD 46223, HD 46150 and HD 46966) are dominated by the presence of red-noise: we speculate that this is related to a sub-surface convection zone.

The detection of pulsational frequencies in stellar photometry is required as input for asteroseismological modelling. The second short run (SRa02) of the CoRoT mission has provided photometric data of unprecedented quality and time-coverage for a number of O-type stars. We analyse the CoRoT data corresponding to three hot O-type stars, describing the properties of their light curves and we search for pulsational frequencies, which we then compare to theoretical model predictions. We determine the amplitude spectrum of the data, using the Lomb-Scargle and a multifrequency HMM-like technique. Frequencies are extracted by prewhitening, and their significance is evaluated under the assumption that the light curve is dominated by red noise. We search for harmonics, linear combinations and regular spacings among these frequencies. We use simulations with the same time sampling as the data as a powerful tool to judge the significance of our results. From the theoretical point of view, we use the MAD non-adiabatic pulsation code to determine the expected frequencies of excited modes. A substantial number of frequencies is listed, but none can be convincingly identified as being connected to pulsations. The amplitude spectrum is dominated by red noise. Theoretical modelling shows that all three O-type stars can have excited modes but the relation between the theoretical frequencies and the observed spectrum is not obvious. The dominant red noise component in the hot O-type stars studied here clearly points to a different origin than the pulsations seen in cooler O stars. The physical cause of this red noise is unclear, but we speculate on the possibility of sub-surface convection, granulation, or stellar wind inhomogeneities being responsible.

Detectable radio emission occurs during almost all phases of massive star evolution. I will concentrate on the thermal and non-thermal continuum emission from early-type stars. The thermal radio emission is due to free-free interactions in the ionized stellar wind material. Early ideas that this would lead to an easy and straightforward way of measuring the mass-loss rates were thwarted by the presence of clumping in the stellar wind. Multi-wavelength observations provide important constraints on this clumping, but do not allow its full determination. Non-thermal radio emission is associated with binarity. This conclusion was already known for some time for Wolf-Rayet stars and in recent years it has become clear that it is also true for O-type stars. In a massive-star binary, the two stellar winds collide and around the shocks a fraction of the electrons are accelerated to relativistic speeds. Spiralling in the magnetic field these electrons emit synchrotron radiation, which we detect as non-thermal radio emission. The many parameters that influence the resulting non-thermal radio fluxes make the modelling of these systems particularly challenging, but their study will provide interesting new insight into massive stars.

We develop 3-D models of the structured winds of massive hot stars with the Wind3D radiative transfer (RT) code. We investigate the physical properties of large-scale structures observed in the wind of the B-type supergiant HD 64760 with detailed line profile fits to Discrete Absorption Components (DACs) and rotational modulations observed with IUE in Si IV {\lambda}1395. We develop parameterized input models Wind3D with large-scale equatorial wind density- and velocity-structures, or so-called `Co-rotating Interaction Regions' (CIRs) and `Rotational Modulation Regions' (RMRs). The parameterized models offer important advantages for high-performance RT calculations over ab-initio hydrodynamic input models. The acceleration of the input model calculations permits us to simulate and investigate a wide variety of physical conditions in the extended winds of massive hot stars. The new modeling method is very flexible for constraining the dynamic and geometric wind properties of RMRs in HD 64760. We compute that the modulations are produced by a regular pattern of radial density enhancements that protrude almost linearly into the equatorial wind. We find that the modulations are caused by narrow `spoke-like' wind regions. We present a hydrodynamic model showing that the linearly shaped radial wind pattern can be caused by mechanical wave action at the base of the stellar wind from the blue supergiant.

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.

The SRa02 of the CoRoT space mission for Asteroseismology was partly devoted to stars belonging to the Mon OB2 association. An intense monitoring was performed on Plaskett's Star (HD47129) and the unprecedented quality of the light curve allows us to shed new light on this very massive, non-eclipsing binary system. We particularly aimed at detecting periodic variability which might be associated with pulsations or interactions between both components. We also searched for variations related to the orbital cycle which could help to constrain the inclination and the morphology of the binary system. A Fourier-based prewhitening and a multiperiodic fitting procedure were applied to analyse the time series and extract the frequencies of variations. We describe the noise properties to tentatively define an appropriate significance criterion, to only point out the peaks at a certain significance level. We also detect the variations related to the orbital motion and study them by using the NIGHTFALL program. The periodogram exhibits a majority of peaks at low frequencies. Among these peaks, we highlight a list of about 43 values, including notably two different sets of harmonic frequencies whose fundamental peaks are located at about 0.07 and 0.82d-1. The former represents the orbital frequency of the binary system whilst the latter could probably be associated with non-radial pulsations. The study of the 0.07d-1 variations reveals the presence of a hot spot most probably situated on the primary star and facing the secondary. The investigation of this unique dataset constitutes a further step in the understanding of Plaskett's Star. These results provide a first basis for future seismic modelling. The existence of a hot region between both components renders the determination of the inclination ambiguous.

We present multi-instruments and multi-wavelengths observations of the famous LBV star Pistol Star. These observations are part of a larger program about early O stars at different metallicities. The Pistol star has been claimed as the most massive star known, with 250 solar masses. We present the preliminary results based on X-Shooter spectra, as well as the observations performed with the VLTI-AMBER and the VLT-NACO adaptive optics. The X-shooter spectrograph allows to obtain simultaneously a spectrum from the UV to the K-band with a resolving power of $\sim$15000. The preliminary results obtained indicate that Pistol Star has similar properties of Eta Car, including shells of matter, but also the binarity. Other objects of the program, here briefly presented, were selected for their particular nature: early O stars with mass discrepancies between stellar evolution models and observations, discrepancies with the wind momentum luminosity relation.

First results of near-IR adaptive optics (AO)-assisted imaging, interferometry, and spectroscopy of this Luminous Blue Variable (LBV) are presented. They suggest that the Pistol Star is at least double. If the association is physical, it would reinforce questions concerning the importance of multiplicity for the formation and evolution of extremely massive stars.

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.

In colliding-wind binaries, shocks accelerate a fraction of the electrons up to relativistic speeds. These electrons then emit synchrotron radiation at radio wavelengths. Whether or not we detect this radiation depends on the size of the free-free absorption region in the stellar winds of both components. One expects long-period binaries to be detectable, but not the short-period ones. It was therefore surprising to find that Cyg OB2 No. 8A (P = 21.9 d) does show variability locked with orbital phase. To investigate this, we developed a model for the relativistic electron generation (including cooling and advection) and the radiative transfer of the synchrotron emission through the stellar wind. Using this model, we show that the synchrotron emitting region in Cyg OB2 No. 8A does extend far enough beyond the free-free absorption region to generate orbit-locked variability in the radio flux. This model can also be applied to other non-thermal emitters and will prove useful in interpreting observations from future surveys, such as COBRaS - the Cyg OB2 Radio Survey.

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.

Almost all early-type stars show Discrete Absorption Components (DACs) in their ultraviolet spectral lines. These can be attributed to Co-rotating Interaction Regions (CIRs): large-scale spiral-shaped structures that sweep through the stellar wind. We used the Zeus hydrodynamical code to model the CIRs. In the model, the CIRs are caused by "spots" on the stellar surface. Through the radiative acceleration these spots create fast streams in the stellar wind material. Where the fast and slow streams collide, a CIR is formed. By varying the parameters of the spots, we quantitatively fit the observed DACs in HD 64760. An important result from our work is that the spots do not rotate with the same velocity as the stellar surface. The fact that the cause of the CIRs is not fixed on the surface eliminates many potential explanations. The only remaining explanation is that the CIRs are due to the interference pattern of a number of non-radial pulsations.

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

This document is divided in two parts. The first part deals with the radial velocities (RV) distributions for B-type stars and nebulosities observed with the VLT-GIRAFFE in the Large and Small Magellanic Clouds towards the open clusters NGC2004 and NGC330. Thanks to the resolution of GIRAFFE spectra, we found that the RV distribution for the nebulosities in the LMC is bi-modal. This bi-modality can be interpreted, in term of dynamics, by the expansion of the LMC4 superbubble. The second part deals with the GAIA space mission and the determination of the radial velocities by using Radial Velocity Spectrometer (RVS) spectra. The methods to determine the radial velocities are presented as well as preliminary results on simulated RVS spectra.

The star Cyg OB2 No. 9 is a well-known non-thermal radio emitter. Recent theoretical work suggests that all such O-stars should be in a binary or a multiple system. However, there is no spectroscopic evidence of a binary component. Re-analysis of radio observations from the VLA of this system over 25 years has revealed that the non-thermal emission varies with a period of 2.35+-0.02 yr. This is interpreted as a strong suggestion of a binary system, with the non-thermal emission arising in a wind-collision region. We derived some preliminary orbital parameters for this putative binary and revised the mass-loss rate of the primary star downward from previous estimates.

We model the detailed time-evolution of Discrete Absorption Components (DACs) observed in P Cygni profiles of the Si IV lam1400 resonance doublet lines of the fast-rotating supergiant HD 64760 (B0.5 Ib). We adopt the common assumption that the DACs are caused by Co-rotating Interaction Regions (CIRs) in the stellar wind. We perform 3D radiative transfer calculations with hydrodynamic models of the stellar wind that incorporate these large-scale density- and velocity-structures. We develop the 3D transfer code Wind3D to investigate the physical properties of CIRs with detailed fits to the DAC shape and morphology. The CIRs are caused by irregularities on the stellar surface that change the radiative force in the stellar wind. In our hydrodynamic model we approximate these irregularities by circular symmetric spots on the stellar surface. We use the Zeus3D code to model the stellar wind and the CIRs, limited to the equatorial plane. We constrain the properties of large-scale wind structures with detailed fits to DACs observed in HD 64760. A model with two spots of unequal brightness and size on opposite sides of the equator, with opening angles of 20 +/- 5 degr and 30 +/- 5 degr diameter, and that are 20 +/- 5 % and 8 +/- 5 % brighter than the stellar surface, respectively, provides the best fit to the observed DACs. The recurrence time of the DACs compared to the estimated rotational period corresponds to spot velocities that are 5 times slower than the rotational velocity. The mass-loss rate of the structured wind model for HD 64760 does not exceed the rate of the spherically symmetric smooth wind model by more than 1 %. The fact that DACs are observed in a large number of hot stars constrains the clumping that can be present in their winds, as substantial amounts of clumping would tend to destroy the CIRs.

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 hydrodynamical models for Corotating Interaction Regions, which were used by Lobel (2007) to model the Discrete Absorption Components in HD 64760. We also discuss our failure to model the rotational modulations seen in the same star.

We discuss the development of the new radiative transfer code Wind3D. It solves the non-LTE radiative transport problem in moving stellar atmosphere models in three geometric dimensions. The code accepts arbitrary 3D velocity fields in Cartesian geometry without assumptions of axial symmetry. Wind3D is currently implemented as a fully parallelized (exact) accelerated lambda iteration scheme with a two level atom formulation. The numerical transfer scheme is efficient and very accurate to trace small variations of local velocity gradients on line opacity in strongly scattering dominated extended stellar winds. We investigate the detailed formation of P Cygni line profiles observed in ultraviolet spectra of massive stars. We compute the detailed shape of these resonance lines to model local enhancements of line opacity that can for instance be caused by clumping in supersonically expanding winds. Wind3D will be applied to hydrodynamic models to investigate physical properties of discrete absorption line components.

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.

We present a model for the non-thermal radio emission from presumably single O stars, in terms of synchrotron emission from relativistic electrons accelerated in wind-embedded shocks. These shocks are associated with an unstable, chaotic wind. The main improvement with respect to earlier models is the inclusion of the radial dependence of the shock velocity jump and compression ratio, based on 1D hydrodynamical simulations. The decrease of the velocity jump and the compression ratio as a function of radius produces a rapidly decreasing synchrotron emissivity. This effectively prohibits the models from reproducing the spectral shape of the observed non-thermal radio emission. We investigate a number of "escape routes" by which the hydrodynamical predictions might be reconciled with the radio observations. Although these escape routes reproduce the observed spectral shape, none of these escape routes are physically plausible. In particular, re-acceleration by feeding an electron distribution through a number of shocks, is in contradiction with current hydrodynamical simulations. These hydrodynamical simulations have their limitations, most notably the use of 1D. At present, it is not feasible to perform 2D simulations of the wind out to the distances required for synchrotron-emission models. Based on the current hydrodynamic models, we suspect that the observed non-thermal radio emission from O stars cannot be explained by wind-embedded shocks associated with the instability of the line-driving mechanism. The most likely alternative mechanism is synchrotron emission from colliding winds. That would imply that all O stars with non-thermal radio emission should be members of binary or multiple systems.

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 a model for the non-thermal radio emission from bright O stars, in terms of synchrotron emission from wind-embedded shocks. The model is an extension of an earlier one, with an improved treatment of the cooling of relativistic electrons. This improvement limits the synchrotron-emitting volume to a series of fairly narrow layers behind the shocks. We show that the width of these layers increases with increasing wavelength, which has important consequences for the shape of the spectrum. We also show that the strongest shocks produce the bulk of the emission, so that the emergent radio flux can be adequately described as coming from a small number of shocks, or even from a single shock. A single shock model is completely determined by four parameters: the position of the shock, the compression ratio and velocity jump of the shock, and the surface magnetic field. Applying a single shock model to the O5 If star Cyg OB2 No. 9 allows a good determination of the compression ratio and shock position and, to a lesser extent, the magnetic field and velocity jump. Our main conclusion is that strong shocks need to survive out to distances of a few hundred stellar radii. Even with multiple shocks, the shocks needed to explain the observed emission are stronger than predictions from time-dependent hydrodynamical simulations.

The solar photospheric carbon abundance has been determined from [C I], C I, CH vibration-rotation, CH A-X electronic and C2 Swan electronic lines by means of a time-dependent, 3D, hydrodynamical model of the solar atmosphere. Departures from LTE have been considered for the C I lines. These turned out to be of increasing importance for stronger lines and are crucial to remove a trend in LTE abundances with the strengths of the lines. Very gratifying agreement is found among all the atomic and molecular abundance diagnostics in spite of their widely different line formation sensitivities. The mean of the solar carbon abundance based on the four primary abundance indicators ([C I], C I, CH vibration-rotation, C_2 Swan) is log C = 8.39 +/- 0.05, including our best estimate of possible systematic errors. Consistent results also come from the CH electronic lines, which we have relegated to a supporting role due to their sensitivity to the line broadening. The new 3D based solar C abundance is significantly lower than previously estimated in studies using 1D model atmospheres.

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 present a theoretical model for the non-thermal radio emission from single hot stars, in terms of synchrotron radiation from electrons accelerated in wind-embedded shocks. The model is described by five independent parameters each with a straightforward physical interpretation. Applying the model to a high-quality observation of Cyg OB2 No. 9 (O5 If), we obtain meaningful constraints on most parameters. The most important result is that the outer boundary of the synchrotron emission region must lie between 500 and 2200 stellar radii. This means that shocks must persist up to that distance. We also find that relatively weak shocks (with a compression ratio < 3) are needed to produce the observed radio spectrum. These results are compatible with current hydrodynamical predictions. Most of our models also show a relativistic electron fraction that increases outwards. This points to an increasing efficiency of the acceleration mechanism, perhaps due to multiple acceleration, or an increase in the strength of the shocks. Implications of our results for non-thermal X-ray emission are discussed.

We present radio and submillimetre observations of the O4I(n)f star zeta Pup, and discuss structure in the outer region of its wind (~ 10-100 R*). The properties of bremsstrahlung, the dominant emission process at these wavelengths, make it sensitive to structure and allow us to study how the amount of structure changes in the wind by comparing the fluxes at different wavelengths. To look for variability, we acquired 3.6 and 6 cm observations with the Australia Telescope Compact Array (ATCA). We supplemented these with archive observations from the NRAO Very Large Array (VLA). We did not find variability at more than the +- 20 % level. The long integration time does allow an accurate determination of the fluxes at 3.6 and 6 cm. Converting these fluxes into a mass loss rate, we find Mdot = 3.5 x 10^{-6} Msun/yr. This value confirms the significant discrepancy with the mass loss rate derived from the Halpha profile, making zeta Pup an exception to the usually good agreement between Halpha and radio mass loss rates. We also observed zeta Pup at 850 mum with the James Clerk Maxwell Telescope (JCMT) and at 20 cm with the VLA. A smooth wind model shows that the millimetre fluxes are too high compared to the radio fluxes. While recombination of helium in the outer wind cannot be discounted as an explanation, the wealth of evidence for structure strongly suggests this as the explanation for the discrepancy. Model calculations show that the structure needs to be present in the inner ~ 70 R* of the wind, but that it decays significantly, or maybe even disappears, beyond that radius.