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.

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.

The asteroseismic and planetary studies, like all research related to stars, need precise and accurate stellar atmospheric parameters as input. We aim at deriving the effective temperature (Teff), the surface gravity (log g), the metallicity ([Fe/H]), the projected rotational velocity (v sin i) and the MK type for 169 F, G, K, and M-type Kepler targets which were observed spectroscopically from the ground with five different instruments. We use two different spectroscopic methods to analyse 189 high-resolution, high-signal-to-noise spectra acquired for the 169 stars. For 67 stars, the spectroscopic atmospheric parameters are derived for the first time. KIC 9693187 and 11179629 are discovered to be double-lined spectroscopic binary systems. The results obtained for those stars for which independent determinations of the atmospheric parameters are available in the literature are used for a comparative analysis. As a result, we show that for solar-type stars the accuracy of present determinations of atmospheric parameters is +/- 150 K in Teff, +/- 0.15 dex in [Fe/H], and +/-? 0.3 dex in log g. Finally, we confirm that the curve-of-growth analysis and the method of spectral synthesis yield systematically different atmospheric parameters when they are applied to stars hotter than 6,000 K.

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.

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.

Many early-type stars are in binary systems. A number of them shows radio emissivity with periodic variability. This variability is associated with non-thermal synchrotron radiation emitted by relativistic electrons. The strong shocks necessary to accelerate the electrons up to high energies are produced by the collision of the radiatively-driven stellar winds. A study of the non-thermal emission is necessary in order to investigate O-star colliding wind binaries. Here preliminary results of our modeling of the colliding winds in Cyg OB2 No.9 are presented.

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.

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.

The main goal of our present work is to provide, for the first time, a simple computational tool that can be used to compute the brightness, the spectral index, the polarization, the time variability and the spectrum of the non-thermal light (both synchrotron and inverse Compton, IC) associated with the plasma dynamics resulting from given relativistic magnetohydrodynamics (RMHD) simulations. The proposed method is quite general, and can be applied to any scheme for RMHD and to all non-thermal emitting sources, e.g. pulsar wind nebulae (PWNe), and in particular to the Crab Nebula (CN) as in the present proceeding. Here only the linear optical and X-ray polarization that characterizes the PWNe synchrotron emission is analyzed in order to infer information on the inner bulk flow structure, to provide a direct investigation of the magnetic field configuration, in particular the presence and the strength of a poloidal component, and to understand the origin of some emitting features, such as the knot, whose origins are still uncertain. The inverse Compton radiation is examined to disentangle the different contributions to radiation from the magnetic field and the particle energy distribution function, and to search for a possible hadronic component in the emitting PWN, and thus for the presence of ions in the wind. If hadronic radiation was found in a PWN, young supernova remnants would provide a natural birth-place of the cosmic-rays (CRs) up to the so-called knee in the CR spectrum.

In this paper we complete the set of diagnostic tools for synchrotron emitting sources presented by Del Zanna et al. (Astron. Astrophys. 453, 621, 2006) with the computation of inverse Compton radiation from the same relativistic particles. Moreover we investigate, for the first time, the gamma-ray emission properties of Pulsar Wind Nebulae in the light of the axisymmetric jet-torus scenario. The method consists in evolving the relativistic MHD equations and the maximum energy of the emitting particles. The particle energy distribution function is split in two components: the radio one connected to a relic population born at the outburst of the supernova and the other associated to the wind population continuously accelerated at the termination shock and emitting up to the gamma-ray band. We consider the general Klein-Nishina cross section and three different photon targets: the nebular synchrotron photons, far-infrared thermal ones and the cosmic microwave background. The overall synchrotron spectrum is fitted assuming an excess of injected particles and a steeper power law with respect to previous models. The TeV emission has the correct shape but is in excess of the data. This is due to the nebular magnetic field structure as obtained by the simulations. The jet-torus morphology is visible in high-resolution gamma-ray synthetic maps too. We present a preliminary exploration of time variability in the X and gamma-ray bands.

We present a complete set of diagnostic tools aimed at reproducing synthetic non-thermal (synchrotron and/or Inverse Compton, IC) emissivity, integrated flux energy, polarization and spectral index simulated maps in comparison to observations. The time dependent relativistic magnetohydrodynamic (RMHD) equations are solved with a shock capturing code together with the evolution of the maximum particles energy. Applications to Pulsar Wind Nebulae (PWNe) are shown.

A complete set of diagnostic tools aimed at producing synthetic synchrotron emissivity, polarization, and spectral index maps from relativistic MHD simulations is presented. As a first application we consider here the case of the emission from Pulsar Wind Nebulae (PWNe). The proposed method is based on the addition, on top of the basic set of MHD equations, of an extra equation describing the evolution of the maximum energy of the emitting particles. This equation takes into account adiabatic and synchrotron losses along streamlines for the distribution of emitting particles and its formulation is such that it is easily implemented in any numerical scheme for relativistic MHD. Application to the axisymmetric simulations of PWNe, analogous to those described by Del Zanna et al. (2004, A&A, 421, 1063), allows direct comparison between the numerical results and observations of the inner structure of the Crab Nebula, and similar objects, in the optical and X-ray bands. We are able to match most of the observed features typical of PWNe, like the equatorial torus and the polar jets, with velocities in the correct range, as well as finer emission details, like arcs, rings and the bright knot, that turn out to arise mainly from Doppler boosting effects. Spectral properties appear to be well reproduced too: detailed spectral index maps are produced for the first time and show softening towards the PWN outer borders, whereas spectral breaks appear in integrated spectra. The emission details are found to strongly depend on both the average wind magnetization (here approximately 2%), and on the magnetic field shape.

We present here the first effort to compute synthetic synchrotron polarization maps of Pulsar Wind Nebulae (PWNe). Our goal is to highlight how polarization can be used as an additional diagnostic tool for the flow structure in the inner regions of these nebulae. Recent numerical simulations suggest the presence of flow velocities ~0.5 c in the surroundings of the termination shock, where most of the high energy emission comes from. We construct polarization maps taking into account relativistic effects like Doppler boosting and position angle swing. The effect of different bulk velocities is clarified with the help of a toy-model consisting of a uniformly emitting torus. We also present a map based on recent numerical simulations of the entire nebula and compare it with presently available data. The comparison with upcoming high resolution observations could provide new insight into the inner structure of the nebula and put constraints on the geometrical properties of the magnetic field.