We present a model for the non-thermal emission from a colliding-wind binary. Relativistic protons and electrons are assumed to be accelerated through diffusive shock acceleration (DSA) at the global shocks bounding the wind-wind collision region. The non-linear effects of the back-reaction due to the cosmic ray pressure on the particle acceleration process and the cooling of the non-thermal particles as they flow downstream from the shocks are included. We explore how the non-thermal particle distribution and the keV-GeV emission changes with the stellar separation and the viewing angle of the system, and with the momentum ratio of the winds. We confirm earlier findings that DSA is very efficient when magnetic field amplification is not included, leading to significantly modified shocks. We also find that the non-thermal flux scales with the binary separation in a complicated way and that the anisotropic inverse Compton emission shows only a moderate variation with viewing angle due to the spatial extent of the wind-wind collision.

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.

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

Field's linear analysis of thermal instability is repeated using methods related to Whitham's theory of wave hierarchies, which brings out the physically relevant parameters in a much clearer way than in the original analysis. It is also used for the stability of non-equilibrium states and we show that for gas cooling behind a shock, the usual analysis is only quantitatively valid for shocks that are just able to trigger a transition to the cold phase. A magnetic field can readily be included and we show that this does not change the stability criteria. By considering steady shock solutions, we show that almost all plausible initial conditions lead to a magnetically dominated state on the unstable part of the equilibrium curve. These results are used to analyse numerical calculations of perturbed steady shock solutions and of shocks interacting with a warm cloud.

Studies of filamentary structures that are prevalent throughout the interstellar medium are of great significance to a number of astrophysical fields. Here, we present 3D hydrodynamic simulations of shock-filament interactions where the equation of state has been softened to become almost isothermal. We investigate the effect of such an isothermal regime on the interaction (where both the shock and filament are isothermal), and we examine how the nature of the interaction changes when the orientation of the filament, the shock Mach number, and the filament density contrast are varied. We find that only sideways-oriented filaments with a density contrast of $10^2$ form a three-rolled structure, dissimilar to the results of a previous study. Moreover, the angle of orientation of the filament plays a large role in the evolution of the filament morphology: the greater the angle of orientation, the longer and less turbulent the wake. Turbulent stripping of filament material leading to fragmentation of the core occurs in most filaments; however, filaments orientated at an angle of $85^{\circ}$ to the shock front do not fragment and are longer-lived. In addition, values of the drag time are influenced by the filament length, with longer filaments being accelerated faster than shorter ones. Furthermore, filaments in an isothermal regime exhibit faster acceleration than those struck by an adiabatic shock. Finally, we find that the drag and mixing times of the filament increase as the angle of orientation of the filament is increased.

We investigate the effect of mass-loading from embedded clouds on the evolution of supernova remnants and on the energy and momentum that they inject into an inhomogeneous interstellar medium. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit. The destruction of embedded clouds adds mass into the remnant, increasing its density and pressure, and decreasing its temperature. The remnant cools more quickly, is less able to do PdV work on the swept-up gas, and ultimately attains a lower final momentum (by up to a factor of two or more). We thus find that the injection of momentum is more sensitive to an inhomogeneous environment than previous work has suggested, and we provide fits to our results for the situation where the cloud mass is not limited. The behaviour of the remnant is more complex in situations where the cloud mass is finite and locally runs out. In the case of multiple supernovae in a clustered environment, later supernova explosions may encounter higher densities than previous explosions due to the prior liberation of mass from engulfed clouds. If the cloud mass is finite, later explosions may be able to create a sustained hot phase when earlier explosions have not been able to.

Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.

This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.1cm^{-3} and evolved with self-gravity included. No magnetic field was imposed. Resimulations at a higher resolution of a region extracted from this simulation (down to 0.039pc), show that the thermal instability forms sheets, then filaments and finally clumps. The width of the filaments increases over time, in one particular case from 0.26 to 0.56pc. Thereafter clumps with sizes of around 5pc grow at the intersections of filaments. 21 distinct clumps, with properties similar to those observed in molecular clouds, are found by using the FellWalker algorithm to find minima in the gravitational potential. Not all of these are gravitationally bound, but the convergent nature of the flow and increasing central density suggest they are likely to form stars. Further simulation of the most massive clump shows the gravitational collapse to a density >10^6 cm^{-3}. These results provide realistic initial conditions that can be used to study feedback in individual clumps, interacting clumps and the entire molecular cloud complex.

The similarities, or otherwise, of a shock or wind interacting with a cloud of density contras t$ \chi = 10$ were explored in a previous paper. Here, we investigate such interactions with clouds of higher density contrast. We compare the adiabatic hydrodynamic interaction of a Mach 10 shock with a spherical cloud of $\chi = 10^{3}$ with that of a cloud embedded in a wind with identical parameters to the post-shock flow. We find that initially there are only minor morphological differences between the shock-cloud and wind-cloud interactions, compared to when $\chi = 10$. However, once the transmitted shock exits the cloud, the development of a turbulent wake and fragmentation of the cloud differs between the two simulations. On increasing the wind Mach number, we note the development of a thin, smooth tail of cloud material, which is then disrupted by the fragmentation of the cloud core and subsequent `mass-loading' of the flow. We find that the normalized cloud mixing time ($t_{mix}$) is shorter at higher $\chi$. However, a strong Mach number dependence on tmix and the normalized cloud drag time, $t'_{drag}$, is not observed. Mach-number-dependent values of $t_{mix}$ and $t'_{drag}$ from comparable shock-cloud interactions converge towards the Mach-number-independent time-scales of the wind-cloud simulations. We find that high $\chi$ clouds can be accelerated up to 80-90 per cent of the wind velocity and travel large distances before being significantly mixed. However, complete mixing is not achieved in our simulations and at late times the flow remains perturbed.

We investigate the structure and X-ray emission from the colliding stellar winds in massive star binaries. We find that the opening angle of the contact discontinuity (CD) is overestimated by several formulae in the literature at very small values of the wind momentum ratio. We find also that the shocks in the primary (dominant) and secondary winds flare by approx 20 degrees compared to the CD, and that the entire secondary wind is shocked when the wind momentum ratio < 0.02. Analytical expressions for the opening angles of the shocks, and the fraction of each wind that is shocked, are provided. We find that the X-ray luminosity scales with the wind momentum ratio, and that the spectrum softens slightly as the wind momentum ratio decreases.

The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD46223, from NGC 2244 occurred 1.73 (+0.34,-0.25)Myr (1$\sigma$ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD46150, the age is set for the famous ionised region at more than ten times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40M$_{\odot}$ star in a thin sheet-like molecular cloud. We form the 135,000M$_{\odot}$ cloud from thermally-unstable diffuse interstellar medium under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5Myr, including cavity size, stellar age, magnetic field and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback.

We report on the first detection of a global change in the X-ray emitting properties of a wind-wind collision, thanks to XMM-Newton observations of the massive SMC system HD5980. While its lightcurve had remained unchanged between 2000 and 2005, the X-ray flux has now increased by a factor of ~2.5, and slightly hardened. The new observations also extend the observational coverage over the entire orbit, pinpointing the lightcurve shape. It has not varied much despite the large overall brightening, and a tight correlation of fluxes with orbital separation is found, without any hysteresis effect. Moreover, the absence of eclipses and of absorption effects related to orientation suggests a large size for the X-ray emitting region. Simple analytical models of the wind-wind collision, considering the varying wind properties of the eruptive component in HD5980, are able to reproduce the recent hardening and the flux-separation relationship, at least qualitatively, but they predict a hardening at apastron and little change in mean flux, contrary to observations. The brightness change could then possibly be related to a recently theorized phenomenon linked to the varying strength of thin-shell instabilities in shocked wind regions.

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 interaction of a hot, high-velocity wind with a cold, dense molecular cloud has often been assumed to resemble the evolution of a cloud embedded in a post-shock flow. However, no direct comparative study of these two processes currently exists in the literature. We present 2D adiabatic hydrodynamical simulations of the interaction of a Mach 10 shock with a cloud of density contrast $\chi = 10$ and compare our results with those of a commensurate wind-cloud simulation. We then investigate the effect of varying the wind velocity, effectively altering the wind Mach number $M_{wind}$, on the cloud's evolution. We find that there are significant differences between the two processes: 1) the transmitted shock is much flatter in the shock-cloud interaction; 2) a low-pressure region in the wind-cloud case deflects the flow around the edge of the cloud in a different manner to the shock-cloud case; 3) there is far more axial compression of the cloud in the case of the shock. As $M_{wind}$ increases, the normalised rate of mixing is reduced. Clouds in winds with higher $M_{wind}$ also do not experience a transmitted shock through the cloud's rear and are more compressed axially. In contrast with shock-cloud simulations, the cloud mixing time normalised by the cloud-crushing time-scale $t_{cc}$ increases for increasing $M_{wind}$ until it plateaus (at $t_{mix} \simeq 25 \, t_{cc}$) at high $M_{wind}$, thus demonstrating the expected Mach scaling. In addition, clouds in high Mach number winds are able to survive for long durations and are capable of being moved considerable distances.

We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with randomly distributed high density condensations. The stellar winds are introduced using appropriate non-rotating Geneva stellar evolution models. In the 15 Msun star case, the wind has very little effect, spreading around a few neighbouring clumps before becoming overwhelmed by the cloud collapse. In contrast, in the 40 Msun, 60 Msun and 120 Msun star cases, the more powerful stellar winds create large cavities and carve channels through the cloud, breaking out into the surrounding tenuous medium during the wind phase and considerably altering the cloud structure. After 4.97 Myrs, 3.97 Myrs and 3.01 Myrs respectively, the massive stars explode as supernovae (SNe). The wind-sculpted surroundings considerably affect the evolution of these SN events as they both escape the cloud along wind-carved channels and sweep up remaining clumps of cloud/wind material. The `cloud' as a coherent structure does not survive the SN from any of these stars, but only in the 120 Msun case is the cold molecular material completely destabilised and returned to the unstable thermal phase. In the 40 Msun and 60 Msun cases, coherent clumps of cold material are ejected from the cloud by the SN, potentially capable of further star formation.

We present three-dimensional simulations of a protoplanetary disk subject to the effect of a nearby (0.3pc distant) supernova, using a time-dependent flow from a one dimensional numerical model of the supernova remnant (SNR), in addition to constant peak ram pressure simulations. Simulations are performed for a variety of disk masses and inclination angles. We find disk mass-loss rates that are typically 1e-7 to 1e-6 Msol/yr (but peak near 1e-5 Msol/yr during the "instantaneous" stripping phase) and are sustained for around 200 yr. Inclination angle has little effect on the mass loss unless the disk is close to edge-on. Inclined disks also strip asymmetrically with the trailing edge ablating more easily. Since the interaction lasts less than one outer rotation period, there is not enough time for the disk to restore its symmetry, leaving the disk asymmetrical after the flow has passed. Of the low-mass disks considered, only the edge-on disk is able to survive interaction with the SNR (with 50% of its initial mass remaining). At the end of the simulations, disks that survive contain fractional masses of SN material up to 5e-6. This is too low to explain the abundance of short-lived radionuclides in the early solar system, but a larger disk and the inclusion of radiative cooling might allow the disk to capture a higher fraction of SN material.

We have modelled the evolution of cometary HII regions produced by zero-age main-sequence stars of O and B spectral types, which are driving strong winds and are born off-centre from spherically symmetric cores with power-law ($\alpha = 2$) density slopes. A model parameter grid was produced that spans stellar mass, age and core density. Exploring this parameter space we investigated limb-brightening, a feature commonly seen in cometary HII regions. We found that stars with mass $M_\star \geq 12\, \mathrm{M}_\odot$ produce this feature. Our models have a cavity bounded by a contact discontinuity separating hot shocked wind and ionised ambient gas that is similar in size to the surrounding HII region. Due to early pressure confinement we did not see shocks outside of the contact discontinuity for stars with $M_\star \leq 40\, \mathrm{M}_\odot$, but the cavities were found to continue to grow. The cavity size in each model plateaus as the HII region stagnates. The spectral energy distributions of our models are similar to those from identical stars evolving in uniform density fields. The turn-over frequency is slightly lower in our power-law models due to a higher proportion of low density gas covered by the HII regions.

We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morphology that in projection appears filamentary. The stellar winds are introduced using appropriate Geneva stellar evolution models. In the 15 solar mass star case, the wind forms a narrow bipolar cavity with minimal effect on the parent cloud. In the 40 solar mass star case, the more powerful stellar wind creates a large cylindrical cavity through the centre of the cloud. After 12.5 Myrs and 4.97 Myrs respectively, the massive stars explode as supernovae (SNe). In the 15 solar mass star case, the SN material and energy is primarily deposited into the molecular cloud surroundings over ~10^5 years before the SN remnant escapes the cloud. In the 40 solar mass star case, a significant fraction of the SN material and energy rapidly escapes the molecular cloud along the wind cavity in a few tens of kiloyears. Both SN events compress the molecular cloud material around them to higher densities (so may trigger further star formation), and strengthen the magnetic field, typically by factors of 2-3 but up to a factor of 10. Our simulations are relevant to observations of bubbles in flattened ring-like molecular clouds and bipolar HII regions.

We present 3D magnetohydrodynamic numerical simulations of the adiabatic interaction of a shock with a dense, filamentary cloud. We investigate the effects of various filament lengths and orientations on the interaction using different orientations of the magnetic field, and vary the Mach number of the shock, the density contrast of the filament, and the plasma beta, in order to determine their effect on the evolution and lifetime of the filament. We find that in a parallel magnetic field filaments have longer lifetimes if they are orientated more 'broadside' to the shock front, and that an increase in the density contrast hastens the destruction of the cloud, in terms of the modified cloud-crushing time-scale, tcs. The combination of a mild shock and a perpendicular or oblique field provides the best condition for extending the life of the filament, with some filaments able to survive almost indefinitely since they are cocooned by the magnetic field. A high value for the density contrast does not initiate large turbulent instabilities in either the perpendicular or oblique field cases but rather draws the filament out into long tendrils which may eventually fragment. In addition, flux ropes are only formed in parallel magnetic fields. The length of the filament is, however, not as important for the evolution and destruction of a filament.

We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure.

We present 3D hydrodynamic adiabatic simulations of a shock interacting with a dense, elongated cloud. We compare how the nature of the interaction changes with the filament's length and its orientation to the shock, and with the shock Mach number and the density contrast of the filament. We then examine the differences with respect to 3D spherical-cloud calculations. We find significant differences in the morphology of the interaction when M=10 and chi=100: in many cases 3 parallel rolls are formed, and spread further apart with time, and periodic vortex shedding can occur off the ends of oblique filaments. Sideways-on filaments are accelerated more quickly, and initially lose mass more quickly than spherical clouds due to their greater surface area to volume ratio. However, at late stages they lose mass more slowly, due to the reduced relative speed between the filament and the postshock flow. The acceleration and mixing timescales can vary by a factor of 2 as the filament orientation changes. Oblique filaments can achieve transverse velocities up to 10% of the shock speed. Some aspects of our simulations are compared against experimental and numerical work on rigid cylinders.

The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.

We present 3D hydrodynamic simulations of the adiabatic interaction of a shock with a dense, spherical cloud. We compare how the nature of the interaction changes with the Mach number of the shock, $M$, and the density contrast of the cloud, $\chi$. We examine the differences with 2D axisymmetric calculations, perform detailed resolution tests, and compare "inviscid" results to those obtained with the inclusion of a $k$-$\epsilon$ subgrid turbulence model. Resolutions of 32-64 cells per cloud radius are the minimum necessary to capture the dominant dynamical processes in 3D simulations, while the 3D inviscid and $k$-$\epsilon$ simulations typically show very good agreement. Clouds accelerate and mix up to 5 times faster when they are poorly resolved. The interaction proceeds very similarly in 2D and 3D - although non-azimuthal modes lead to different behaviour, there is very little effect on key global quantities such as the lifetime of the cloud and its acceleration. In particular, we do not find significant differences in the hollowing or "voiding" of the cloud between 2D and 3D simulations with $M=10$ and $\chi=10$, which contradicts previous work in the literature.

Suzaku observations of the Wolf-Rayet binary WR 140 (WC7pd+O5.5fc) were made at four different times around periastron passage in 2009 January. The spectra changed in shape and flux with the phase. As periastron approached, the column density of the low-energy absorption increased, which indicates that the emission from the wind-wind collision plasma was absorbed by the dense W-R wind. The spectra can be mostly fitted with two different components: a warm component with kT=0.3--0.6 keV and a dominant hot component with kT~3 keV. The emission measure of the dominant, hot component is not inversely proportional to the distance between the two stars. This can be explained by the O star wind colliding before it has reached its terminal velocity, leading to a reduction in its wind momentum flux. At phases closer to periastron, we discovered a cool plasma component in a recombining phase, which is less absorbed. This component may be a relic of the wind-wind collision plasma, which was cooled down by radiation, and may represent a transitional stage in dust formation.

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

The Suzaku X-ray observatory monitored the supermassive binary system Eta Carinae 10 times during the whole 5.5 year orbital cycle between 2005-2011. This series of observations presents the first long-term monitoring of this enigmatic system in the extremely hard X-ray band between 15-40 keV. During most of the orbit, the 15-25 keV emission varied similarly to the 2-10 keV emission, indicating an origin in the hard energy tail of the kT ~4 keV wind-wind collision (WWC) plasma. However, the 15-25 keV emission declined only by a factor of 3 around periastron when the 2-10 keV emission dropped by two orders of magnitude due probably to an eclipse of the WWC plasma. The observed minimum in the 15-25 keV emission occurred after the 2-10 keV flux had already recovered by a factor of ~3. This may mean that the WWC activity was strong, but hidden behind the thick primary stellar wind during the eclipse. The 25-40 keV flux was rather constant through the orbital cycle, at the level measured with INTEGRAL in 2004. This result may suggest a connection of this flux component to the gamma-ray source detected in this field. The Helium-like Fe Kalpha line complex at ~6.7 keV became strongly distorted toward periastron as seen in the previous cycle. The 5-9 keV spectra can be reproduced well with a two-component spectral model, which includes plasma in collision equilibrium (CE) and a plasma in non-equilibrium ionization (NEI) with tau ~1e11 cm-3 s-1. The NEI plasma increases in importance toward periastron.

We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $\beta$ in cloud material can vary rapidly as clouds collide with one another, but low values of $\beta$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds.

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 X-ray emission from a simulated massive stellar cluster is investigated. The emission is calculated from a 3D hydrodynamical model which incorporates the mechanical feedback from the stellar winds of 3 O-stars embedded in a giant molecular cloud clump containing 3240 M$_{\odot}$ of molecular material within a 4 pc radius. A simple prescription for the evolution of the stars is used, with the first supernova explosion at t=4.4 Myrs. We find that the presence of the GMC clump causes short-lived attenuation effects on the X-ray emission of the cluster. However, once most of the material has been ablated away by the winds the remaining dense clumps do not have a noticable effect on the attenuation compared with the assumed interstellar medium column. We determine the evolution of the cluster X-ray luminosity, L$_X$, and spectra, and generate synthetic images. The intrinsic X-ray luminosity drops from nearly 10$^{34}$ ergs s$^{-1}$ while the winds are `bottled up', to a near constant value of 1.7$\times 10^{32}\rm ergs s^{-1}$ between t=1-4 Myrs. L$_X$ reduces slightly during each star's red supergiant stage due to the depressurization of the hot gas. However, L$_X$ increases to $\approx 10^{34}\rm\,ergs s^{-1}$ during each star's Wolf-Rayet stage. The X-ray luminosity is enhanced by 2-3 orders of magnitude to $\sim 10^{37}\rm ergs s^{-1}$ for at least 4600 yrs after each supernova, at which time the blast wave leaves the grid and the X-ray luminosity drops. The X-ray luminosity of our simulation is generally considerably fainter than predicted from spherically-symmetric bubble models, due to the leakage of hot gas material through gaps in the outer shell. This process reduces the pressure within our simulation and thus the X-ray emission. However, the X-ray luminosities and temperatures which we obtain are comparable to similarly powerful massive young clusters.

X-ray emission from the supermassive binary system Eta Carinae declines sharply around periastron. This X-ray minimum has two distinct phases - the lowest flux phase in the first ~3 weeks and a brighter phase thereafter. In 2009, the Chandra X-ray Observatory monitored the first phase five times and found the lowest observed flux at ~1.9e-12 ergs cm-2 s-1 (3-8 keV). The spectral shape changed such that the hard band above ~4 keV dropped quickly at the beginning and the soft band flux gradually decreased to its lowest observed value in ~2 weeks. The hard band spectrum had begun to recover by that time. This spectral variation suggests that the shocked gas producing the hottest X-ray gas near the apex of the wind-wind collision (WWC) is blocked behind the dense inner wind of the primary star, which later occults slightly cooler gas downstream. Shocked gas previously produced by the system at earlier orbital phases is suggested to produce the faint residual X-ray emission seen when the emission near the apex is completely blocked by the primary wind. The brighter phase is probably caused by the re-appearance of the WWC plasma, whose emissivity significantly declined during the occultation. We interpret this to mean that the X-ray minimum is produced by a hybrid mechanism of an occultation and a decline in emissivity of the WWC shock.

We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

Blowouts can occur when a dense shell confining hot, high pressure, gas ruptures. The venting gas inflates a blister on the surface of the shell. Here we examine the growth of such blisters on the surfaces of wind-blown-bubbles (WBBs) and supernova remnants (SNRs) due to shell rupture caused by the Vishniac instability. On WBBs the maximum relative size of the blister (R_bstall/R) is found to grow linearly with time, but in many cases the blister radius will not exceed 20 per cent of the bubble radius. Thus blowouts initiated by the Vishniac instability are unlikely to have a major effect on the global dynamics and properties of the bubble. The relative size of blisters on SNRs is even smaller than on WBBs, with blisters only growing to a radius comparable to the thickness of the cold shell of SNRs. The small size of the SNR blowouts is, however, in good agreement with observations of blisters in the Vela SNR. The difference in relative size between WBB and SNR blisters is due to the much higher speed at which gas vents out of WBBs, which translates into a greater energy flux through a rupture of a given size from interior gas of a given pressure. Larger blisters are possible if shell ruptures are bigger than expected. We expect the observed velocity structure of SNR shells to be affected by the presence of blisters until the shell is no longer susceptible to ruptures, since the initial expansion of blisters is faster than the ongoing expansion of the shell.

We use 3D hydrodynamical models to investigate the effects of massive star feedback from winds and supernovae on inhomogeneous molecular material left over from the formation of a massive stellar cluster. We simulate the interaction of the mechanical energy input from a cluster with 3 O-stars into a giant molecular cloud (GMC) clump containing 3240 solar masses of molecular material within a 4 pc radius. The cluster wind blows out of the molecular clump along low-density channels, into which denser clump material is entrained. We find that the densest molecular regions are surprisingly resistant to ablation by the cluster wind, in part due to shielding by other dense regions closer to the cluster. Nonetheless, molecular material is gradually removed by the cluster wind during which mass-loading factors in excess of several 100 are obtained. Because the clump is very porous, 60-75 per cent of the injected wind energy escapes the simulation domain, with the difference being radiated. After 4.4 Myr, the massive stars in our simulation begin to explode as supernovae. The highly structured environment into which the SN energy is released allows even weaker coupling to the remaining dense material and practically all of the SN energy reaches the wider environment. The molecular material is almost completely dispersed and destroyed after 6 Myr. The escape fraction of ionizing radiation is estimated to be about 50 per cent during the first 4 Myr of the cluster's life. A similar model with a larger and more massive GMC clump reveals the same general picture, though more time is needed for it to be destroyed.

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.

We simulate the effects of massive star feedback, via winds and SNe, on inhomogeneous molecular material left over from the formation of a massive stellar cluster. We use 3D hydrodynamic models with a temperature dependent average particle mass to model the separate molecular, atomic, and ionized phases. We find that the winds blow out of the molecular clump along low-density channels, and gradually ablate denser material into these. However, the dense molecular gas is surprisingly long-lived and is not immediately affected by the first star in the cluster exploding.

We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution.

In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent.

X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (O9.7 I+b2 v, PA = 21 d) and B (O8 III+o9 v, PB = 6 d) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three temperature thermal plasma model, characterised by cool, moderate, and hot plasma components at kT ~ 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of ~ 0.2 x 10^22 cm-2. Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of ~ 7 x 10^-13 erg s-1 cm-2, do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer the model and observations are in far better agreement, which lends support to the previous suggestion of mass-transfer in the O8 III + o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.

The Chandra Carina Complex Project (CCCP) has shown that the Carina Nebula displays bright, spatially-complex soft diffuse X-ray emission. Here we `sum up' the CCCP diffuse emission work by comparing the global morphology and spectrum of Carina's diffuse X-ray emission to other famous sites of massive star formation with pronounced diffuse X-ray emission: M17, NGC 3576, NGC 3603, and 30 Doradus. All spectral models require at least two diffuse thermal plasma components to achieve adequate spectral fits, a softer component with kT = 0.2--0.6 keV and a harder component with kT = 0.5--0.9 keV. In several cases these hot plasmas appear to be in a state of non-equilibrium ionization that may indicate recent and current strong shocks. A cavity north of the embedded giant HII region NGC 3576 is the only region studied here that exhibits hard diffuse X-ray emission; this emission appears to be nonthermal and is likely due to a recent cavity supernova, as evidenced by a previously-known pulsar and newly-discovered pulsar wind nebula also seen in this cavity. All of these targets exhibit X-ray emission lines that are not well-modeled by variable-abundance thermal plasmas and that might be attributed to charge exchange at the shock between the hot, tenuous, X-ray-emitting plasma and cold, dense molecular material; this is likely evidence for dust destruction at the many hot/cold interfaces that characterize massive star-forming regions.

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

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

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

I discuss some of the important aspects of the phenomena of colliding winds in massive binary star systems with a particular focus on WR140.

We present the results from the spectroscopic follow-up of WR140 (WC7 + O4-5) during its last periastron passage in January 2009. This object is known as the archetype of colliding wind binaries and has a relatively large period (~ 8 years) and eccentricity (~ 0.89). We provide updated values for the orbital parameters, new estimates for the WR and O star masses and new constraints on the mass-loss rates.

A new mechanism for the formation of cometary tails behind dense clouds or globules is discussed. Numerical hydrodynamical models show that when a dense shell of swept-up matter overruns a cloud, material in the shell is focussed behind the cloud to form a tail. This mode of tail formation is completely distinct from other methods, which involve either the removal of material from the cloud, or shadowing from a strong, nearby source of ionization. This mechanism is relevant to the cometary tails seen in planetary nebulae and to the interaction of superbubble shells with dense clouds.

Three dimensional (3D) adaptive-mesh refinement (AMR) hydrodynamical simulations of the wind-wind collision between the enigmatic super-massive star \etacar and its mysterious companion star are presented which include radiative driving of the stellar winds, gravity, optically-thin radiative cooling, and orbital motion. Simulations with static stars with a periastron passage separation reveal that the preshock companion star's wind speed is sufficiently reduced that radiative cooling in the postshock gas becomes important, permitting the runaway growth of non-linear thin shell (NTSI) instabilities which massively distort the WCR. However, large-scale simulations which include the orbital motion of the stars, show that orbital motion reduces the impact of radiative inhibition, and thus increases the acquired preshock velocities. As such, the postshock gas temperature and cooling time see a commensurate increase, and sufficient gas pressure is preserved to stabilize the WCR against catastrophic instability growth. We then compute synthetic X-ray spectra and lightcurves and find that, compared to previous models, the X-ray spectra agree much better with {\it XMM-Newton} observations just prior to periastron. The narrow width of the 2009 X-ray minimum can also be reproduced. However, the models fail to reproduce the extended X-ray mimimum from previous cycles. We conclude that the key to explaining the extended X-ray minimum is the rate of cooling of the companion star's postshock wind. If cooling is rapid then powerful NTSIs will heavily disrupt the WCR. Radiative inhibition of the companion star's preshock wind, albeit with a stronger radiation-wind coupling than explored in this work, could be an effective trigger.

The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.

We present preliminary results of a 3-month campaign carried out in the framework of the Mons project, where time-resolved Halpha observations are used to study the wind and circumstellar properties of a number of OB stars.

Hydrodynamical models of colliding hypersonic flows are presented which explore the dependence of the resulting dynamics and the characteristics of the derived X-ray emission on numerical conduction and viscosity. For the purpose of our investigation we present models of colliding flow with plane-parallel and cylindrical divergence. Numerical conduction causes erroneous heating of gas across the contact discontinuity which has implications for the rate at which the gas cools. We find that the dynamics of the shocked gas and the resulting X-ray emission are strongly dependent on the contrast in the density and temperature either side of the contact discontinuity, these effects being strongest where the postshock gas of one flow behaves quasi-adiabatically while the postshock gas of the other flow is strongly radiative. Introducing additional numerical viscosity into the simulations has the effect of damping the growth of instabilities, which in some cases act to increase the volume of shocked gas and can re-heat gas via sub-shocks as it flows downstream. The resulting reduction in the surface area between adjacent flows, and therefore of the amount of numerical conduction, leads to a commensurate reduction in spurious X-ray emission, though the dynamics of the collision are compromised. The simulation resolution also affects the degree of numerical conduction. A finer resolution better resolves the interfaces of high density and temperature contrast and although numerical conduction still exists the volume of affected gas is considerably reduced. However, since it is not always practical to increase the resolution, it is imperative that the degree of numerical conduction is understood so that inaccurate interpretations can be avoided. This work has implications for the dynamics and emission from astrophysical phenomena which involve high Mach number shocks.

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3.

The X-ray emission from the wind-wind collision in short-period massive O+O-star binaries is investigated. The emission is calculated from three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, orbital motion of the stars, and radiative cooling of the shocked plasma. Changes in the amount of stellar occultation and circumstellar attenuation introduce phase-dependent X-ray variability in systems with circular orbits, while strong variations in the intrinsic emission also occur in systems with eccentric orbits. The X-ray emission in eccentric systems can display strong hysteresis, with the emission softer after periastron than at corresponding orbital phases prior to periastron, reflecting the physical state of the shocked plasma at these times. Furthermore, the rise of the luminosity to maximum does not necessarily follow a 1/D law. Our models further demonstrate that the effective circumstellar column can be highly energy dependent. We simulate Chandra and Suzaku observations, and fit these using standard XSPEC models. We find that the recovered temperatures from two or three-temperature mekal fits are comparable to those from fits to the emission from real systems with similar stellar and orbital parameters/nature. We also find that when the global abundance is thawed in the spectral fits, sub-solar values are exclusively returned, despite the calculations using solar values as input (abridged).

In this work the thermal emission over cm to sub-mm wavelengths from the winds in short-period O+O-star binaries is investigated (potential non-thermal emission is presently ignored). The calculations are based on three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, orbital motion of the stars, and radiative cooling of the shocked plasma. The thermal emission arises from the stellar winds and the region where they collide. We investigate the flux and spectrum from a variety of models as a function of orbital phase and orientation of the observer, and compare to the single star case. The emission from the wind-wind collision region (WCR) is strongly dependent on its density and temperature, being optically thick in radiative systems, and optically thin in adiabatic systems. The flux from systems where the WCR is highly radiative, as investigated for the first time in this work, can be over an order of magnitude greater than the combined flux from identically typed single stars. This excess occurs over a broad range of wavelengths from cm to sub-mm. In contrast, when the WCR is largely adiabatic, a significant excess in the thermal flux occurs only below 100 GHz. Eccentric systems may show order of magnitude or greater flux variability, especially if the plasma in the WCR transitions from an adiabatic to a radiative state and vice-versa - in such cases the flux can display significant hysteresis with stellar separation. We further demonstrate that clumping can affect the variability of radio lightcurves. We investigate the spectral index of the emission, and often find indices steeper than +0.6. Our predictions are of interest to future observations with the next generation of radio and sub-mm telescopes (abridged).

In the young stellar object (YSO) phase of their lives, massive stars drive bi-polar molecular outflows. These outflows produce beautiful, often hourglass shaped, cavities. The central star possesses a powerful stellar wind (v ~ 2000 km s^-1), and possibly a dense equatorial disk wind (v ~ 400 km s^-1), which collide with the inner surface of the bi-polar cavity and produces hot (T ~ 10^5 - 10^8 K) shocked plasma. A reverse shock is formed at the point where the ram pressure between the preshock flow balances the thermal pressure of the postshock flow and provides a site for the acceleration of non-thermal particles to relativistic energies. Hydrodynamical models of the wind interaction, coupled with calculations of the non-thermal energy spectrum, are used to explore the observable synchrotron and gamma-ray emission from these objects.

2D axis-symmetric hydrodynamical simulations are presented which explore the interaction of stellar and disk winds with surrounding infalling cloud material. The star, and its accompanying disk, blow winds inside a cavity cleared out by an earlier jet. The collision of the winds with their surroundings generates shock heated plasma which reaches temperatures up to ~10^8 K. Attenuated X-ray spectra are calculated from solving the equation of radiative transfer along lines-of-sight. This process is repeated at various epochs throughout the simulations to examine the evolution of the intrinsic and attenuated flux. We find that the dynamic nature of the wind-cavity interaction fuels intrinsic variability in the observed emission on timescales of several hundred years. This is principally due to variations in the position of the reverse shock which is influenced by changes in the shape of the cavity wall. The collision of the winds with the cavity wall can cause clumps of cloud material to be stripped away. Mixing of these clumps into the winds mass-loads the flow and enhances the X-ray emission measure. The position and shape of the reverse shock plays a key role in determining the strength and hardness of the X-ray emission. In some models the reverse shock is oblique to much of the stellar and disk outflows, whereas in others it is closely normal over a wide range of polar angles. For reasonable stellar and disk wind parameters the integrated count rate and spatial extent of the intensity peak for X-ray emission agree with \textit{Chandra} observations of the deeply embedded MYSOs S106 IRS4, Mon R2 IRS3 A, and AFGL 2591.(abridged)

The powerful wind-wind collision in massive star binaries creates a region of high temperature plasma and accelerates particles to relativistic energies. I briefly summarize the hydrodynamics of the wind-wind interaction and the observational evidence, including recent $\gamma$-ray detections, of non-thermal emission from such systems. I then discuss existing models of the non-thermal emission and their application to date, before concluding with some future prospects.

The dynamics of the wind-wind collision in massive stellar binaries is investigated using three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, the orbital motion of the stars, and radiative cooling of the shocked plasma. In this first paper we restrict our study to main-sequence O+O binaries. The nature of the wind-wind collision region is highly dependent on the degree of cooling of the shocked plasma, and the ratio of the flow timescale of the shocked plasma to the orbital timescale. The pre-shock wind speeds are lower in close systems as the winds collide prior to their acceleration to terminal speeds. Radiative inhibition may also reduce the pre-shock wind speeds. Together, these effects can lead to rapid cooling of the post-shock gas. Radiative inhibition is less important in wider systems, where the winds are accelerated to higher speeds before they collide, and the resulting collision region can be largely adiabatic. In systems with eccentric orbits, cold gas formed during periastron passage can persist even at apastron, before being ablated and mixed into its surroundings and/or accelerated out of the system.

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

The X-ray emission from the super-massive star Eta Carinae is simulated using a three dimensional model of the wind-wind collision. In the model the intrinsic X-ray emission is spatially extended and energy dependent. Absorption due to the unshocked stellar winds and the cooled postshock material from the primary LBV star is calculated as the intrinsic emission is ray-traced along multiple sightlines through the 3D spiral structure of the circumstellar environment. The observable emission is then compared to available X-ray data, including the lightcurve observed by the Rossi X-ray Timing Explorer (RXTE) and spectra observed by XMM-Newton. The orientation and eccentricity of the orbit are explored, as are the wind parameters of the stars and the nature and physics of their close approach. Our modelling supports a viewing angle with an inclination of ~ 42 degrees, consistent with the polar axis of the Homunculus nebula (Smith 2006), and the projection of the observer's line-of-sight onto the orbital plane has an angle of ~ 0 - 30 degrees in the prograde direction on the apastron side of the semi-major axis. However, there are significant discrepancies between the observed and model lightcurves and spectra through the X-ray minimum. In particular, the hard flux in our synthetic spectra is an order of magnitude greater than observed. Further calculations reveal that radiative inhibition significantly reduces the preshock velocity of the companion wind. As a consequence the hard X-ray emission is quenched, but it is unclear whether the long duration of the minimum is due solely to this mechanism alone. Models incorporating a collapse/disruption of the WCR and/or reduced preshock companion wind velocities bring the predicted emission and the observations into much better agreement (abridged).

The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous single-cloud studies. In reality, the flow behind the shock is also likely to be turbulent, with non-uniform density, pressure and velocity structure created as the shock sweeps over inhomogenities upstream of the cloud. To address these twin issues we use a sub-grid compressible k-epsilon turbulence model to estimate the properties of the turbulence generated in shock-cloud interactions and the resulting increase in the transport coefficients that the turbulence brings. A detailed comparison with the output from an inviscid hydrodynamical code puts these new results into context. We find that cloud destruction in inviscid and k-epsilon models occurs at roughly the same speed when the post-shock flow is smooth and when the density contrast between the cloud and inter-cloud medium is less than 100. However, there are increasing and significant differences as this contrast increases. Clouds subjected to strong ``buffeting'' by a highly turbulent post-shock environment are destroyed significantly quicker. Additional calculations with an inviscid code where the post-shock flow is given random, grid-scale, motions confirms the more rapid destruction of the cloud. Our results clearly show that turbulence plays an important role in shock-cloud interactions, and that environmental turbulence adds a new dimension to the parameter space which has hitherto been studied (abridged).

We show that the mechanical energy injection rate generated as the intra-cluster medium (ICM) flows around cold clouds may be sufficient to power the optical and near infra-red emission of nebulae observed in the central regions of a sample of seven galaxy clusters. The energy injection rate is extremely sensitive to the velocity difference between the ICM and cold clouds, which may help to explain why optical and infra-red luminosity is often larger than expected in systems containing AGNs. We also find that mass recycling is likely to be important for the dynamics of the ICM. This effect will be strongest in the central regions of clusters where there is more than enough cold gas for its evaporation to contribute significantly to the density of the hot phase.

We present a 3D dynamical model of the orbital induced curvature of the wind-wind collision region in binary star systems. Momentum balance equations are used to determine the position and shape of the contact discontinuity between the stars, while further downstream the gas is assumed to behave ballistically. An archimedean spiral structure is formed by the motion of the stars, with clear resemblance to high resolution images of the so-called ``pinwheel nebulae''. A key advantage of this approach over grid or smoothed particle hydrodynamic models is its significantly reduced computational cost, while it also allows the study of the structure obtained in an eccentric orbit. The model is relevant to symbiotic systems and Gamma-ray binaries, as well as systems with O-type and Wolf-Rayet stars. As an example application, we simulate the X-ray emission from hypothetical O+O and WR+O star binaries, and describe a method of ray tracing through the 3D spiral structure to account for absorption by the circumstellar material in the system. Such calculations may be easily adapted to study observations at wavelengths ranging from the radio to Gamma-ray.

We report on variations in important X-ray emission lines in a series of Chandra grating spectra of the supermassive colliding wind binary star Eta Carinae, including key phases around the X-ray minimum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of Eta Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission lines provide the only direct measure of the flow dynamics of the companion's wind along the wind-wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with symmetric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant velocity shifts and correlated increases in line widths through the observations. The R = forbidden-to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit in all observations, suggesting that the line-forming region is >1.6 stellar radii from the companion star. We show that simple geometrical models cannot simultaneously fit both the observed centroid variations and changes in line width as a function of phase. We show that the observed profiles can be fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision boundary. We use this analysis to help constrain the line formation region as a function of orbital phase, and the orbital geometry.

Recent data support the idea that the filaments observed in H_alpha emission near the centres of some galaxy clusters were shaped by bulk flows within their intracluster media. We present numerical simulations of evaporated clump material interacting with impinging winds to investigate this possibility. In each simulation, a clump falls due to gravity while the drag of a wind retards the fall of evaporated material leading to elongation of the tail. However, we find that long filaments can only form if the outflowing wind velocity is sufficiently large, 10^8 cm s^-1. Otherwise, the tail material sinks almost as quickly as the cloud. For reasonable values of parameters, the morphological structure of a tail is qualitatively similar to those observed in clusters. Under certain conditions, the kinematics of the tail resemble those reported in Hatch et al.(2006). A comparison of the observations with the numerical results indicates that the filaments are likely to be a few tens of Myrs old. We also present arguments which suggest that the momentum transfer, from an outflowing wind, in the formation of these filaments is probably significant. As a result, tail formation could play a role in dissipating some of the energy injected by a central AGN close to the cluster centre where it is needed most. The trapping of energy by the cold gas may provide an additional feedback mechanism that helps to regulate the heating of the central regions of galaxy clusters and couple the AGN to the ICM.

We present the results of optical wavelength observations of the unusual SMC eclipsing binary system HD 5980 obtained in 1999 and 2004--2005. Radial velocity curves for the erupting LBV/WR object (star A) and its close WR-like companion (star B) are obtained by deblending the variable emission-line profiles of N IV and N V lines under the simplistic assumption that these lines originate primarily in the winds of star A and star B. The derived masses M_A=58--79 Mo and M_B=51--67 Mo, are more consistent with the stars' location near the top of the HRD than previous estimates. The presence of a wind-wind interaction region is inferred from the orbital phase-dependent behavior of He I P Cygni absorption components. The emission-line intensities continued with the declining trend previously seen in UV spectra. The behavior of the photospheric absorption lines is consistent with the results of Schweickhardt (2002) who concludes that the third object in the combined spectrum, star C, is also a binary system with P(starC)~96.5 days, e=0.83. The data used in this paper will be made publicly available for further analysis.

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.

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

Massive WR+O star systems produce high-temperature, shock-heated plasma where the wind of the WR star and that of its binary companion collide - the wind-collision region (WCR). The WCR is a source of thermal (e.g. hard X-rays) and non-thermal (e.g. synchrotron) emission, the latter arising from electrons and ions accelerated to relativistic energies. These colliding wind binaries provide an excellent laboratory for the study of particle acceleration at higher mass, photon and magnetic energy densities than exist in SNRs. Recent models of the non-thermal emission from WR 140 have provided insight into this process.

Clumping in hot star winds can significantly affect estimates of mass-loss rates, the inferred evolution of the star and the environmental impact of the wind. A hydrodynamical simulation of a colliding winds binary (CWB) with clumpy winds reveals that the clumps are rapidly destroyed after passing through the confining shocks of the wind-wind collision region (WCR) for reasonable parameters of the clumps if the flow in the WCR is adiabatic. Despite large density and temperature fluctuations in the post-shock gas, the overall effect of the interaction is to smooth the existing structure in the winds. Averaged over the entire interaction region, the resulting X-ray emission is very similar to that from the collision of smooth winds. The insensitivity of the X-ray emission to clumping suggests it is an excellent diagnostic of the stellar mass-loss rates in wide CWBs, and may prove to be a useful addition to existing techniques for deriving mass-loss rates, many of which are extremely sensitive to clumping. Clumpy winds also have implications for a variety of phenomena at the WCR: particle acceleration may occur throughout the WCR due to supersonic MHD turbulence, re-acceleration at multiple shocks, and re-connection; a statistical description of the properties of the WCR may be required for studies of non-equilibrium ionization and the rate of electron heating; and the physical mixing of the two winds will be enhanced, as seems necessary to trigger dust formation.

We report the results of an X-ray observing campaign on the massive, evolved star Eta Carinae, concentrating on the 2003 X-ray minimum as seen by the XMM-Newton observatory. These are the first spatially-resolved X-ray monitoring observations of the stellar X-ray spectrum during the minimum. The hard X-ray emission, believed to be associated with the collision of Eta Carinae's wind with the wind from a massive companion star, varied strongly in flux on timescales of days, but not significantly on timescales of hours. The lowest X-ray flux in the 2-10 keV band seen by XMM-Newton was only 0.7% of the maximum seen by RXTE just before the X-ray minimum. The slope of the X-ray continuum above 5 keV did not vary in any observation, which suggests that the electron temperature of the hottest plasma associated with the stellar source did not vary significantly at any phase. Through the minimum, the absorption to the stellar source increased by a factor of 5-10 to NH ~3-4E23 cm-2. The thermal Fe XXV emission line showed significant excesses on both the red and blue sides of the line outside the minimum and exhibited an extreme red excess during the minimum. The Fe fluorescence line at 6.4 keV increased in equivalent width from 100 eV outside the minimum to 200 eV during the minimum. From these observed features, we discuss two possible causes of the X-ray minimum; the eclipse of the X-ray plasma and an intrinsic fading of the X-ray emissivity. The drop in the colliding wind X-ray emission also revealed the presence of an additional X-ray component which exhibited no variation on timescales of weeks to years. This component may be produced by the collision of high speed outflows at v \~1000-2000 km s-1 from Eta Carinae with ambient gas within a few thousand AU from the star.

High-resolution radio observations have revealed that non-thermal radio emission in WR stars arises where the stellar wind of the WR star collides with that of a binary companion. These colliding-wind binary (CWB) systems offer an important laboratory for investigating the underlying physics of particle acceleration. Hydrodynamic models of the binary stellar winds and the wind-collision region (WCR) that account for the evolution of the electron energy spectrum, largely due to inverse Compton cooling, are now available. Radiometry and imaging obtained with the VLA, MERLIN, EVN and VLBA provide essential constraints to these models. Models of the radio emission from WR146 and WR147 are shown, though these very wide systems do not have defined orbits and hence lack a number of important model parameters. Multi-epoch VLBI imaging of the archetype WR+O star binary WR140 through a part of its 7.9-year orbit has been used to define the orbit inclination, distance and the luminosity of the companion star to enable the best constraints for any radio emitting CWB system. Models of the spatial distribution of relativistic electrons and ions, and the magnetic energy density are used to model the radio emission, and also to predict the high energy emission at X-ray and gamma-ray energies. It is clear that high-energy facilities e.g. GLAST and VERITAS, will be important for constraining particle acceleration parameters such as the spectral index of the energy spectrum and the acceleration efficiency of both ions and electrons, and in turn, identify unique models for the radio spectra. This will be especially important in future attempts to model the spectra of WR140 throughout its complete orbit. A WCR origin for the synchrotron emission in O-stars, the progenitors of WR stars, is illustrated by observations of Cyg OB2 No. 9.

We use hydrodynamical models of the wind-collision region (WCR) in the archetype colliding-wind system WR140 to determine the spatial and spectral distribution of the radio, X-ray and gamma-ray emission from shock accelerated electrons. Our calculations are for orbital phase 0.837 when the observed radio emission is close to maximum. Using the observed thermal X-ray emission together with the radio emission to constrain the mass-loss rates, we find that the O-star mass-loss rate is consistent with recent reductions for O4-5 supergiants. We demonstrate that radio VLBI observations of the WCR fail to constrain the opening angle. The observed low frequency turnover at ~3 GHz in the radio emission is due to free-free absorption, since models based on the Razin effect have an unacceptably large fraction of energy in non-thermal particles. The index of the non-thermal electron energy distribution is flatter than the canonical value for diffusive shock acceleration, namely p<2. Several mechanisms are discussed that could lead to such an index. Tighter constraints on p and the nature of the shocks in WR140 will be obtained from future observations at MeV and GeV energies, for which we generally predict lower fluxes than previous work. Since the high stellar photon fluxes prevent the acceleration of electrons beyond gamma > 1e5-1e6, TeV emission from CWB systems will provide unambiguous evidence of pion-decay emission from accelerated ions. We finish by commenting on the emission and physics of the multiple wind collisions in dense stellar clusters, paying particular attention to the Galactic Centre (abridged).

A key process within astronomy is the exchange of mass, momentum, and energy between diffuse plasmas in many types of astronomical sources (including planetary nebulae, wind-blown bubbles, supernova remnants, starburst superwinds, and the intracluster medium) and dense, embedded clouds or clumps. This transfer affects the large scale flows of the diffuse plasmas as well as the evolution of the clumps. I review our current understanding of mass-injection processes, and examine intermediate-scale structure and the global effect of mass-loading on a flow. I then discuss mass-loading in a variety of diffuse sources.

We have examined a stream-source model for the production of the cometary tails observed in the Helix Nebula NGC 7293 in which a transonic or moderately supersonic stream of ionized gas overruns a source of ionized gas. Hydrodynamic calculations reveal velocity structures which are in good agreement with the observational data on tail velocities and are consistent with observations of the nebular structure. The results also are indicative of a stellar atmosphere origin for the cometary globules. Tail remnants persist for timescales long enough for their identification with faint striations visible in the nebula gas to be plausible.

WR+O star binary systems exhibit synchrotron emission arising from relativistic electrons accelerated where the wind of the WR star and that of its massive binary companion collide - the wind-collision region (WCR). These ``colliding-wind'' binaries (CWB), provide an excellent laboratory for the study of particle acceleration, with the same physical processes as observed in SNRs, but at much higher mass, photon and magnetic energy densities. WR140 is the best studied CWB, and high resolution radio observations permit a determination of several system parameters, particularly orbit inclination and distance, that are essential constraints for newly developed models of CWBs. We show a model fit to the radio data at orbital phase 0.9, and show how these models may be used to predict the high energy emission from WR140.

We have developed radiative transfer models of the radio emission from colliding-wind binaries (CWB) based on a hydrodynamical treatment of the wind-collision region (WCR). The archetype of CWB systems is the 7.9-yr period binary WR140, which exhibits dramatic variations at radio wavelengths. High-resolution radio observations of WR140 permit a determination of several system parameters, particularly orbit inclination and distance, that are essential for any models of this system. A model fit to data at orbital phase 0.9 is shown, and some short comings of our model described.

Radio emission models of colliding wind binaries (CWBs) have been discussed by Dougherty et al. (2003). We extend these models by considering the temporal and spatial evolution of the energy distribution of relativistic electrons as they advect downstream from their shock acceleration site. The energy spectrum evolves significantly due to the strength of inverse-Compton (IC) cooling in these systems, and a full numerical evaluation of the synchrotron emission and absorption coefficients is made. We have demonstrated that the geometry of the WCR and the streamlines of the flow within it lead to a spatially dependent break frequency in the synchrotron emission. We therefore do not observe a single, sharp break in the synchrotron spectrum integrated over the WCR, but rather a steepening of the synchrotron spectrum towards higher frequencies. We also observe that emission from the wind-collision region (WCR) may appear brightest near the shocks, since the impact of IC cooling on the non-thermal electron distribution is greatest near the contact discontinuity (CD). We show that the thermal emission from a CWB can mimic a thermal plus non-thermal composite spectrum if the thermal emission from the WCR becomes comparable to that from the unshocked winds, and obtain a number of insights relevant to models of closer systems such as WR140. Finally, we apply our new models to the very wide system WR147. The acceleration of non-thermal electrons appears to be very efficient in our models of WR147, and we suggest that the shock structure may be modified by feedback from the accelerated particles. (Abridged).

WR140 is the archetype long-period colliding wind binary (CWB) system, and is well known for dramatic variations in its synchrotron emission during its 7.9-yr, highly eccentric orbit. This emission is thought to arise from relativistic electrons accelerated at the global shocks bounding the wind-collision region (WCR). The presence of non-thermal electrons and ions should also give rise to X-ray and gamma-ray emission from several separate mechanisms, including inverse-Compton cooling, relativistic bremsstrahlung, and pion decay. We describe new calculations of this emission and make some preliminary predictions for the new generation of gamma-ray observatories. We determine that WR140 will likely require several Megaseconds of observation before detection with INTEGRAL, but should be a reasonably strong source for GLAST.

WR146 is a WC6+O8 colliding-wind binary (CWB) system with thermal emission from the stellar winds of the two stars, and bright non-thermal emission from the wind-collision region (WCR) where the winds collide. We present high resolution radio observations from 1.4 to 43 GHz that give one of the best quality radio spectra of any CWB to date. Observations at 22 GHz now span 8 years, and reveal the proper motion of the system, allowing comparison of multi-epoch data. VLBI observations show the location of the WCR relative to the stellar components, from which the wind momentum ratio can be shown to be 0.06+/-0.15. The radio spectrum and the spatial distribution of emission are modelled, and we determine the contribution of both stellar winds and the WCR to the observed emission. We show that our current models fail to account for the high frequency spectrum of WR146, and also produce too much emission far from the stagnation point of the wind collision.

Because of its key role in feedback in star formation and galaxy formation, we examine the nature of the interaction of a flow with discrete sources of mass injection. We show the results of two-dimensional numerical simulations in which we explore a range of configurations for the mass sources and study the effects of their proximity on the downstream flow. The mass sources act effectively as a single source of mass injection if they are so close together that the ratio of their combined mass injection rate is comparable to or exceeds the mass flux of the incident flow into the volume that they occupy. The simulations are relevant to many diffuse sources, such as planetary nebulae and starburst superwinds, in which a global flow interacts with material evaporating or being ablated from the surface of globules of cool, dense gas.

We investigate the Mach number dependence of the radiative overstability of shocks. We find that the stability of radiative shocks increases with decreasing Mach number, with the result that M=2 shocks require cooling exponents < -1.2 to be overstable. The frequency of oscillation of the fundamental mode has a strong Mach number dependence. We find that feedback between the cooling region and the cold dense layer (CDL) further downstream is a function of Mach number, with stronger feedback and oscillation of the boundary between the CDL and the cooling region occuring at lower Mach numbers. This feedback can be quantified in terms of the reflection coefficient of sound waves. An interesting finding is that the stability properties of low Mach number shocks can be dramatically altered if the shocked gas is able to cool to temperatures less than the pre-shock value, and it is probable that low Mach number astrophysical shocks will be overstable in a variety of situations. The results of this work are relevant to astrophysical shocks with low Mach numbers, such as supernova remnants (SNRs) immersed in a hot interstellar medium (e.g., within a starburst region), and shocks in molecular clouds, where time-dependent chemistry can lead to overstability. (Abridged)

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 new analysis of an archived Chandra HETGS X-ray spectrum of the WR+O colliding wind binary Gamma Velorum. The spectrum is dominated by emission lines from astrophysically abundant elements: Ne, Mg, Si, S and Fe. From a combination of broad-band spectral analysis and an analysis of line flux ratios we infer a wide range of temperatures in the X-ray emitting plasma (~4-40 MK). As in the previously published analysis, we find the X-ray emission lines are essentially unshifted, with a mean FWHM of 1240 +/- 30 km/s. Calculations of line profiles based on hydrodynamical simulations of the wind-wind collision predict lines that are blueshifted by a few hundred km/s. The lack of any observed shift in the lines may be evidence of a large shock-cone opening half-angle (> 85 degrees), and we suggest this may be evidence of sudden radiative braking. From the R and G ratios measured from He-like forbidden-intercombination-resonance triplets we find evidence that the Mg XI emission originates from hotter gas closer to the O star than the Si XIII emission, which suggests that non-equilibrium ionization may be present.

The collision of the hypersonic winds in early-type binaries produces shock heated gas, which radiates thermal X-ray emission, and relativistic electrons, which emit nonthermal radio emission. We review our current understanding of the emission in these spectral regions and discuss models which have been developed for the interpretation of this emission. Physical processes which affect the resulting emission are reviewed and ideas for the future are noted.

X-ray line profiles represent a new way of studying the winds of massive stars. In particular, they enable us to probe in detail the wind-wind collision in colliding wind binaries, providing new insights into the structure and dynamics of the X-ray-emitting regions. We present the key results of new analyses of high-resolution Chandra X-ray spectra of two important colliding wind systems, Gamma Velorum and WR140. The lines of Gamma Vel are essentially unshifted from their rest wavelengths, which we suggest is evidence of a wide shock opening angle, indicative of sudden radiative braking. The widths of the lines of WR140 are correlated with ionization potential, implying non-equilibrium ionization. The implications of these results for the radio emission from these systems are discussed, as are some of the future directions for X-ray line profile modelling of colliding wind binaries.

Massive colliding wind binary stars serve as laboratories for the study of strong-shock physics. In these systems fundamental flow parameters (velocities, densities, directions) are often well known, and photon fields important for Compton scattering and ionization are usually defined or constrained as well. In this paper we review observed X-ray and radio characteristics of some important colliding wind systems. We discuss how X-ray and radio observations of these systems can be used to test models of the collision of massive flows in the presence of intense UV radiation and magnetic fields.

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 report the first detection of X-ray emission associated with the Homunculus Nebula which surrounds the supermassive star Eta Carinae. The emission is characterized by a temperature in excess of 100 MK, and is consistent with scattering of the time-delayed X-ray flux associated with the star. The nebular emission is bright in the northwestern lobe and near the central regions of the Homunculus, and fainter in the southeastern lobe. We also report the detection of an unusually broad Fe K fluorescent line, which may indicate fluorescent scattering off the wind of a companion star or some other high velocity outflow. The X-ray Homunculus is the nearest member of the small class of Galactic X-ray reflection nebulae, and the only one in which both the emitting and reflecting sources are distinguishable.

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

We present a linear stability analysis of a flow undergoing conductively-driven mass-loading from embedded clouds. We find that mass-loading damps isobaric and isentropic perturbations, and in this regard is similar to the effect of thermal conduction, but is much more pronounced where many embedded clumps exist. The stabilizing influence of mass-loading is wavelength independent against isobaric (condensing) perturbations, but wavelength dependent against isentropic (wave-like) perturbations. We derive equations for the degree of mass-loading needed to stabilize such perturbations. We have also made 1D numerical simulations of a mass-loaded radiative shock and demonstrated the damping of the overstability when mass-loading is rapid enough.

One aspect of supernova remnant evolution that is relatively unstudied is the influence of an AGN environment. A high density ambient medium and a nearby powerful continuum source will assist the cooling of shocked ejecta and swept-up gas. Motion of the surrounding medium relative to the remnant will also affect the remnant morphology. In an extension to previous work we have performed 2D hydrodynamical calculations of SNR evolution in an AGN environment, and have determined the evolutionary behaviour of cold gas in the remnant. The cold gas will contribute to the observed broad line emission in AGNs, and we present preliminary theoretical line profiles from our calculations. A more detailed comparison with observations will be performed in future work. The SNR-AGN interaction may be useful as a diagnostic of the AGN wind.

We present calculations of the spatial and spectral distribution of the radio emission from a wide WR+OB colliding-wind binary system based on high-resolution hydrodynamical simulations and solutions to the radiative transfer equation. We account for both thermal and synchrotron radio emission, free-free absorption in both the unshocked stellar wind envelopes and the shocked gas, synchrotron self-absorption, and the Razin effect. The applicability of these calculations to modelling radio images and spectra of colliding-wind systems is demonstrated with models of the radio emission from the wide WR+OB binary WR147. Its synchrotron spectrum follows a power-law between 5 and 15 GHz but turns down to below this at lower and higher frequencies. We find that while free-free opacity from the circum-binary stellar winds can potentially account for the low-frequency turnover, models that also include a combination of synchrotron self-absorption and Razin effect are favoured. We argue that the high-frequency turn down is a consequence of inverse-Compton cooling. We present our resulting spectra and intensity distributions, along with simulated MERLIN observations of these intensity distributions. From these we argue that the inclination of the WR147 system to the plane of the sky is low. We summarise by considering extensions of the current model that are important for models of the emission from closer colliding wind binaries, in particular the dramatically varying radio emission of WR140.

We present theoretical X-ray line profiles from a range of model colliding wind systems. In particular, we investigate the effects of varying the stellar mass-loss rates, the wind speeds, and the viewing orientation. We find that a wide range of theoretical line profile shapes is possible, varying with orbital inclination and phase. At or near conjunction, the lines have approximately Gaussian profiles, with small widths (HWHM ~ 0.1 v_\infty) and definite blue- or redshifts (depending on whether the star with the weaker wind is in front or behind). When the system is viewed at quadrature, the lines are generally much broader (HWHM ~ v_\infty), flat-topped and unshifted. Local absorption can have a major effect on the observed profiles - in systems with mass-loss rates of a few times 10^{-6} Msol/yr the lower energy lines (E <~ 1 kev) are particularly affected. This generally results in blueward-skewed profiles, especially when the system is viewed through the dense wind of the primary. The orbital variation of the line widths and shifts is reduced in a low inclination binary. The extreme case is a binary with i = 0 degrees, for which we would expect no line profile variation.

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.

We investigate which shocked wind is responsible for the majority of the X-ray emission in colliding wind binaries, an issue where there is some confusion in the literature, and which we show is more complicated than has been assumed. We find that where both winds rapidly cool (typically close binaries), the ratio of the wind speeds is often more important than the momentum ratio, because it controls the energy flux ratio, and the faster wind is generally the dominant emitter. When both winds are largely adiabatic (typically long-period binaries), the slower and denser wind will cool faster and the stronger wind generally dominates the X-ray luminosity.

We analyze new high-resolution Chandra X-ray images of the Wolf-Rayet binary system WR147. This system contains a WN8 star with an early-type companion located 0.6'' to its north, and is the only known early-type binary with a separation on the sky large enough for the wind-wind collision between the stars to currently be resolved at X-ray energies. The 5 ksec Chandra HRC-I image provides the first direct evidence for spatially extended X-ray emission in an early-type binary system. The X-ray emission peaks close to the position of the radio bow shock and north of the WN8 star. A deeper X-ray image is needed to accurately determine the degree of spatial extension, to exactly align the X-ray and optical/radio frames, and to determine whether part of the detected X-ray emission arises in the individual stellar winds. Simulated X-ray images of the wind-wind collision have a FWHM consistent with the data, and maximum likelihood fits suggest that a deeper observation may also constrain the inclination and wind momentum ratio of this system. However, as the WR wind dominates the colliding wind X-ray emission it appears unlikely that the mass-loss rate and the terminal velocity of the companion wind can be separately determined from X-ray observations. We also note an inconsistency between numerical and analytical estimates of the X-ray luminosity ratio of the stronger and weaker wind components, and conclude that the analytical results are in error.

A newly developed 3-D Monte Carlo model is used, in conjunction with a multi-line non-LTE radiative transfer model, to determine the mass-loss rate of the Wolf-Rayet (W-R) star in the massive binary \object{V444 Cyg} (WN5+O6). This independent estimate of mass-loss rate is attained by fitting the observed \HeI (5876) \AA and \HeII (5412) \AA line profiles, and the continuum light curves of three Stokes parameters ((I, Q, U)) in the (V) band simultaneously. The high accuracy of our determination arises from the use of many observational constraints, and the sensitivity of the continuum polarization to the mass-loss rate. Our best fit model suggests that the mass-loss rate of the system is (\dot{M}_{\WR}=0.6(\pm 0.2) \times 10^{-5} M_{\sun} \mathrm{yr}^{-1} ), and is independent of the assumed distance to \object{V444 Cyg}. The fits did not allow a unique value for the radius of the W-R star to be derived. The range of the volume filling factor for the W-R star atmosphere is estimated to be in the range of 0.050 (for $R_{\WR}=5.0 R_{\sun}$) to 0.075 (for $R_{\WR}=2.5 R_{\sun}$). We also found that the blue-side of \HeI (5876 ) \AA and \HeII (5412) \AA lines at phase 0.8 is relatively unaffected by the emission from the wind-wind interaction zone and the absorption by the O-star atmosphere; hence, the profiles at this phase are suitable for spectral line fittings using a spherical radiative transfer model.

We present X-ray spectral fits to a recently obtained Chandra grating spectrum of Eta Carinae, one of the most massive and powerful stars in the Galaxy and which is strongly suspected to be a colliding wind binary system. Hydrodynamic models of colliding winds are used to generate synthetic X-ray spectra for a range of mass-loss rates and wind velocities. They are then fitted against newly acquired Chandra grating data. We find that due to the low velocity of the primary wind (~500 km/s), most of the observed X-ray emission appears to arise from the shocked wind of the companion star. We use the duration of the lightcurve minimum to fix the wind momentum ratio at 0.2. We are then able to obtain a good fit to the data by varying the mass-loss rate of the companion and the terminal velocity of its wind. We find that Mdot ~ 1e-5 Msol/yr and v ~ 3000 km/s. With observationally determined values of ~500-700 km/s for the velocity of the primary wind, our fit implies a primary mass-loss rate of Mdot ~ 2.5e-4 Msol/yr. This value is smaller than commonly inferred, although we note that a lower mass-loss rate can reduce some of the problems noted by Hillier et al. (2001) when a value as high as 1e-3 Msol/yr is used. The wind parameters of the companion are indicative of a massive star which may or may not be evolved. The line strengths appear to show slightly sub-solar abundances, although this needs further confirmation. Based on the over-estimation of the X-ray line strengths in our model, and re-interpretation of the HST/FOS results, it appears that the homunculus nebula was produced by the primary star.

IRS 13E is an infrared, mm and X-ray source in the Galactic Centre. We present the first Chandra X-ray spectrum for IRS 13E and show that it is consistent with a luminous and highly absorbed X-ray binary system. Since the X-ray luminosity is too large for a solitary star, our interpretation is that of an early-type long-period binary with strong colliding winds emission. This naturally explains the observed X-ray spectrum and count rate as well as its lack of significant short term variability. Due to the short lifetime of any nebula 0.2 pc from the putative central super-massive black-hole, we argue that the primary of IRS 13E has exited the LBV phase in the last few thousand years.

We present similarity solutions for adiabatic bubbles that are blown by winds having time independent mechanical luminosities and that are each mass-loaded by the hydrodynamic ablation of distributed clumps. The mass loading is `switched-on' at a specified radius (with free-expansion of the wind interior to this point) and injects mass at a rate per unit volume proportional to M^delta r^lambda where delta = 4/3 (1) if the flow is subsonic (supersonic) with respect to the clumps. In the limit of negligible mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. The presence of mass loading heats the flow, which leads to a reduction in the Mach number of the supersonic freely-expanding flow, and weaker jump conditions across the inner shock. In solutions with large mass loading, it is possible for the wind to connect directly to the contact discontinuity without first passing through an inner shock. For a solution that gives the mass of swept-up ambient gas to be less than the sum of the masses of the wind and ablated material, lambda < -2. Maximum possible values for the ratio of ablated mass to wind mass occur when mass loading starts very close to the bubble center and when the flow is supersonic with respect to the clumps over the entire bubble radius. The maximum temperature in the bubble often occurs near the onset of mass loading, and in some cases can be many times greater than the post-inner-shock temperature. Our solutions are relevant to eg stellar wind-blown bubbles, galactic winds, etc. This work complements Pittard et al (2001) where it was assumed that clumps were evaporated through conductive energy transport.

We show that a cooled region of shocked supernova ejecta forms in a type II supernova-QSO wind interaction, and has a density, an ionization parameter, and a column density compatible with those inferred for the high ionization component of the broad emission line regions in QSOs. The calculations are based on the assumption that the ejecta flow is described initially by a similarity solution investigated by Chevalier (1982) and Nadyozhin (1985) and is spherically symmetric. Heating and cooling appropriate for gas irradiated by a nearby powerful continuum source is included in our model, together with reasonable assumptions for the properties of the QSO wind. The model results are also in agreement with observational correlations and imply reasonable supernova rates.

We present similarity solutions for mass-loaded adiabatic bubbles that are blown by winds having time dependent mechanical luminosities. We consider mass-loading through conductive evaporation of clumps. In the limit of little mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. We find that the Mach number in the shocked mass-loaded wind shows a radial dependence that varies qualitatively from solution to solution. In some cases it is everywhere less than unity in the frame of the clumps being evaporated, while in others it is everywhere greater than unity. In some solutions the mass-loaded shocked wind undergoes one or two sonic transitions in the clump frame. Maximum possible values of the ratio of evaporated mass to stellar wind mass are found as a consequence of the evaporation rates dependence on temperature and the lowering of the temperature by mass-loading. Mass-loading tends to reduce the emissivity in the interior of the bubble relative to its limb, whilst simultaneously increasing the central temperature relative to the limb temperature.

We present several models for IRS 13E, an infrared, mm and X-ray source in the Galactic Center. Our favored interpretation is that of an early-type binary with strong colliding winds emission. This naturally explains the observed X-ray count rate and the strong infrared emission lines, and has a distinct advantage over competing hypotheses based upon a single star or black-hole system.

We analyse two ASCA observations of the highly eccentric O9III + B1III binary Iota Orionis obtained at periastron and apastron. Based on the assumption of a strong colliding winds shock between the stellar components, we expected to see significant variation in the X-ray emission between these phases. The observations proved otherwise: the X-ray luminosities and spectral distributions were remarkably similar. The only noteworthy feature was the hint of a proximity effect during periastron passage, supported also in the optical. We discuss the accuracy of our results, and also analyse archival ROSAT observations. We investigate why we do not see a clear colliding winds signature. A simple model shows that the wind attenuation to the expected position of the shock apex is negligible throughout the orbit, which poses the puzzling question of why the expected 1/D variation (ie. a factor of 7.5) in the intrinsic luminosity is not seen in the data. Two scenarios are proposed: either the colliding winds emission is unexpectedly weak such that intrinsic shocks in the winds dominate the emission, or, alternatively, that the emission observed is colliding winds emission but in a more complex form than we would naively expect. Complex hydrodynamical models are then analyzed. Despite strongly phase-variable emission from the models, both were consistent with the observations. We find that if the mass-loss rates of the stars are low then intrinsic wind shocks could dominate the emission. However, when we assume higher mass-loss rates of the stars, we find that the observed emission could also be consistent with a purely colliding winds origin. To distinguish between the different models X-ray observations with improved phase coverage will be necessary.