Search

The wind collision region (WCR) in a colliding wind binary (CWB) is a particularly violent place, as such, it is surprising that it is also a region where significant quantities of interstellar dust can form. In extreme cases, approximately 30% of the total mass loss rate of a system can be converted into dust. These regions are poorly understood, as observation and simulation of these systems are difficult. In our previous paper we simulated dust growth in CWB systems using an advected scalar model and found our model to be suitable for qualitative study. For this paper we simulated the periodic dust forming CWB (WCd) system WR140 with our dust model, to determine how dust growth changes over the systems periastron passage. We found that dust production increases significantly at periastron passage, which is consistent with IR emission of the surrounding dusty shell. We also find that the dust production rate of the system decreases rapidly as the stars recede from each other, though the rate of decrease is significantly lower than the rate of increase during periastron passage. This was found to be due to strong cooling and its associated thermal instabilities, resulting in cool, high-density pockets of gas in the WCR where dust forms. The WCR also shows a degree of hysteresis, resulting in a radiative post-shock flow even when the stars are separated enough for the region to behave adiabatically., Comment: 9 pages, 9 figures, 5 tables

Dust production is one of the more curious phenomena observed in massive binary systems with interacting winds. The high temperatures, UV photon flux and violent shocks should destroy any dust grains that condense. However, in some extreme cases dust production yields of approximately 30% of the total mass loss rate of the stellar winds have been observed. In order to better understand this phenomenon a parameter space exploration was performed using a series of numerical models of dust producing carbon phase Wolf-Rayet (WCd) systems. These models incorporated a passive scalar dust model simulating dust growth, destruction and radiative cooling. We find that reasonable dust yields were produced by these simulations. Significant changes in the dust yield were caused by changing the mass loss rates of the stars, with a greater mass loss rate contributing to increased dust yields. Similarly, a close orbit between the stars also resulted in higher dust yields. Finally, a high velocity wind shear, which induces Kelvin-Helmholtz (KH) instabilities and wind mixing, drastically increases the dust yields., Comment: 13 pages, 18 figures, 10 tables

[Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$\mu$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($\theta_{B}$). We derive the power spectrum of the intensity and $\theta_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $\theta_{B}$ varies coherently along the filaments. The observed power spectrum of $\theta_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $\theta_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $\theta_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars., Comment: Accepted for publication in Astronomy & Astrophysics

We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic flux is redistributed within the filament from the centre towards the envelope. The rate of the collapse is inversely proportional to the fractional ionisation and two gravitationally-driven ambipolar diffusion regimes for the collapse are observed as predicted in Mouschovias & Morton (1991). For high values of the ionisation coefficient, that is $X \geq 10^{-7}$, the gas is strongly coupled to the magnetic field and the Jeans length is larger than the ambipolar diffusion length scale. Then the collapse is governed by magnetically-regulated ambipolar diffusion. For $X \lesssim 10^{-8}$, the gas is weakly coupled to the magnetic field and the magnetic support is removed by gravitationally-dominated ambipolar diffusion. Here, neutrals and ions only collide sporadically, that is the ambipolar diffusion length scale is larger than the Jeans length. When decaying turbulence is included, additional support is provided to the filament. This slows down the collapse of the filament even in the absence of a magnetic field. When a magnetic field is present, the collapse rate increases by a ratio smaller than for the non-magnetic case. This is because of a speed-up of the ambipolar diffusion due to larger magnetic field gradients generated by the turbulence and because the ambipolar diffusion aids the dissipation of turbulence below the ambipolar diffusion length scale. The highest increase in the rate is observed for the lowest ionisation coefficient and the highest turbulent intensity., Comment: 10 pages, 11 figures, accepted for publication by A&A

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., Comment: 17 pages, 10 figures. Accepted to MNRAS, minor changes and additions to text, primarily in Section 6

We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced., Comment: 8 pages, 4 figures and 3 tables. Accepted by MNRAS

Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problems that have received attention in the last few years. The coverage is not comprehensive, and the authors have drawn on their own results more and touched more briefly on those of others than they would in a normal review. Theirs is a personal contribution to the issue honouring Leon, one of the truly great gentlemen, wits, and most insightful of astrophysicists., Comment: 13 pages, no figures, to be published in Astronomy and Geophysics

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., Comment: 5 pages, to appear in the Proceedings of Cosmic Rays & Their Interstellar Medium Environment CRISM-2011 (Montpellier, 26th June - 1st July 2011), to be published in an issue of MEMORIE della Societ\`a Astronomica Italiana, eds. A. Marcowith, A.Bykov, K.Ferriere, Y.Gallant, P.Hennebelle, T.Montmerle

To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the shocked cloud with a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 3D, we investigate the effect of different magnetic field orientations including parallel, perpendicular and oblique to the shock normal. We find that the angle between the shock normal and the magnetic field must be small to produce clouds with properties similar to observed molecular clouds., Comment: 10 pages, 12 figures, accepted for publication in MNRAS

We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one bright and compact arising from three condensations (N, E and S), and the other weak and extended along the filament. While the first component shows broad lines (linewidths of ~4-7 kms-1) in both SiO J=2-1 and SiO J=3-2, the second one is only detected in SiO J=2-1 and has narrow lines (~0.8 kms-1). The maps of CO and its isotopologues show that low-density filaments are intersecting the IRDC and appear to merge toward the densest portion of the cloud. This resembles the molecular structures predicted by flow-driven, shock-induced and magnetically-regulated cloud formation models. As in outflows associated with low-mass star formation, the excitation temperatures and fractional abundances of SiO toward N, E and S, increase with velocity from ~6 to 40 K, and from ~1E-10 to >1E-8 respectively, over a velocity range of ~7 kms-1. Since 8 micron sources, 24 micron sources and/or extended 4.5 micron emission are detected in N, E and S, broad SiO is likely produced in outflows associated with high-mass protostars. The excitation temperatures and fractional abundances of the narrow SiO lines, however, are very low (~9 K and ~1E-11, respectively), and consistent with the processing of interstellar grains by the passage of a shock with vs~12 kms-1. This emission could be generated i) by a large-scale shock, perhaps remnant of the IRDC formation process; ii) by decelerated or recently processed gas in large-scale outflows driven by 8 micron and 24 micron sources; or iii) by an undetected and widespread population of lower mass protostars. High-angular resolution observations are needed to disentangle between these three scenarios., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

Star formation occurs in dark molecular regions where the number density of hydrogen nuclei, nH, exceeds 10^4 cc and the fractional ionization is 10^-7 or less. Dust grains with sizes ranging up to tenths of microns and perhaps down to tens of nanometers contain just under one percent of the mass. Recombination on grains is important for the removal of gas phase ions, which are produced by cosmic rays penetrating the dark regions. Collisions of neutrals with charged grains contribute significantly to the coupling of the magnetic field to the neutral gas. Consequently, the dynamics of the grains must be included in the magnetohydrodynamic models of large scale collapse, the evolution of waves and the structures of shocks important in star formation., Comment: 17 pages, 3 figures, accepted for publication in the Journal of Plasma Physics special issue in honour of Prof. Padma K. Shukla

Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-type shocks interacting with density perturbations. We include a self-consistent calculation of the thermal and ionisation balances and a fluid treatment of grains. We identify features that develop when a multifluid shock encounters a density inhomogeneity to investigate whether any part of the precursor region ever behaves in a quasi-steady fashion. If it does the shock may be modelled approximately without solving the time-dependent hydromagnetic equations. Simulations were made for initially steady oblique C-type shocks encountering density inhomogeneities. For a semi-finite inhomogeneity with a density larger than the surrounding medium, a transmitted shock evolves from being J-type to a steady C-type shock on a timescale comparable to the ion-flow time through it. A sufficiently upstream part of the precursor of an evolving J-type shock is quasi-steady. The ion-flow timescale is also relevant for the evolution of a shock moving into a region of decreasing density. The models for shocks propagating into regions in which the density increases and then decreases to its initial value cannot be entirely described in terms of the results obtained for monotonically increasing and decreasing densities. For the latter model, the long-term evolution to a C-type shock cannot be approximated by quasi-steady models., Comment: 11 pages, 9 figures

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging and studied the effect of single-sized and multiple sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionisation structures including appropriate microphysics. At low densities the grains do not play any significant role in the shock dynamics. At high densities, the ionisation fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allows a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena., Comment: 9 pages, 4 figures, accepted for publication in MNRAS

In 1986 Alex Dalgarno published a paper entitled "Is Interstellar Chemistry Useful?" By the middle 1970s, and perhaps even earlier, Alex had hoped that astronomical molecules would prove to: possess significant diagnostic utility; control many of the environments in which they exist; stimulate a wide variety of physicists and chemists who are at least as fascinated by the mechanisms forming and removing the molecules as by astronomy. His own research efforts have contributed greatly to the realization of that hope. This paper contains a few examples of: how molecules are used to diagnose large-scale dynamics in astronomical sources including star forming regions and supernovae; the ways in which molecular processes control the evolution of astronomical objects such as dense cores destined to become stars and very evolved giant stars; theoretical and laboratory investigations that elucidate the processes producing and removing astronomical molecules and allow their detection., Comment: 15 pages, 1 figure, to be published in Proceedings of the Dalgarno Celebratory Symposium:Contributions to Atomic, Molecular, and Optical Physics, Astrophysics and Atmospheric Physics, Eds. James Babb, Kate Kirby and Hossein Sadeghpour

Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest., Comment: 4 pages, to be published in the proceedings of Fifth International Conference on Physics of Dusty Plasmas

We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar diffusion then plays an important role in the formation process. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D numerical simulations, we follow the evolution of the fast-mode wave. The simulations are done with a multifluid, adaptive mesh refinement MHD code. Initial perturbations with wavelengths that are 2 orders of magnitude larger than the dissipation length are strongly affected by the ion-neutral drift. Only in situations where there are large variations in the magnetic field corresponding to a highly turbulent gas can fast-mode waves generate dense cores. This means that, in most cores, no substructure can be produced. However, Core D of TMC-1 is an exception to this case. Due to its atypically high ionisation fraction, waves with wavelengths up to 3 orders of magnitude greater than the dissipation length can be present. Such waves are only weakly affected by ambipolar diffusion and can produce dense substructure without large wave-amplitudes. Our results also explain the observed transition from Alfvenic turbulent motion at large scales to subsonic motions at the level of dense cores., Comment: 7 pages; 7 figures; accepted by A&A

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

To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud, we follow the dynamical evolution of the cloud with a time-dependent axisymmetric magnetohydrodynamic code. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium interacts with the cloud. As a fast-mode shock propagates through the cloud, the gas behind it becomes thermally unstable. The $\beta$ value of the gas also becomes much smaller than the initial value of order unity. These conditions are ideal for magnetohydrodynamic waves to produce high-density clumps embedded in a rarefied warm medium. A slow-mode shock follows the fast-mode shock. Behind this shock a dense shell forms, which subsequently fragments. This is a primary region for the formation of massive stars. Our simulations show that only weak and moderate-strength shocks can form cold clouds which have properties typical of giant molecular clouds., Comment: 7 pages, 6 figures, accepted by Astronomy and Astrophysics

We have performed 2D hydromagnetic simulations with an adaptive mesh refinement code to examine the response of a pre-existing initially spherical dense core to a non-linear fast-mode wave. One key parameter is the ratio of the wavelength to the initial core radius. If that ratio is large and the wave amplitude is sufficient, significant compression of the core occurs, as envisaged by Myers & Lazarian (1998) in their "turbulent cooling flow" picture. For smaller values of that ratio, an initial value of the ratio of the thermal pressure to magnetic pressure of 0.2, and sufficiently large wave amplitude, the scattering induces the production of dense substructure in the core. This substructure may be related to that detected in the dense core associated with the cyanopolyyne peak in TMC-1. Our simulations also show that short-wavelength waves, contrary to large-wavelength waves, do not confine dense cores., Comment: 7 pages, 6 figures, accepted by MNRAS

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

Using two dimensional simulations, we study the formation of structures with a high-density contrast by magnetohydrodynamic waves in regions in which the ratio of thermal to magnetic pressure is small. The initial state is a uniform background perturbed by fast-mode wave. Our most significant result is that dense structures persist for far longer in a two-dimensional simulation than in the one-dimensional case. Once formed, these structures persist as long as the fast-mode amplitude remains high., Comment: 6 pages, 7 figures, accepted by MNRAS

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

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

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

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., Comment: 17 pages, 11 figures (10 postscript + 1 gif)

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