High resolution radio observations of the colliding-wind binary WR140

Milli-arcsecond resolution Very Long Baseline Array (VLBA) observations of the archetype WR+O star colliding-wind binary (CWB) system WR140 are presented for 23 epochs between orbital phases 0.74 and 0.97. At 8.4 GHz, the emission in the wind-collision region (WCR) is clearly resolved as a bow-shaped arc that rotates as the orbit progresses. We interpret this rotation as due to the O star moving from SE to approximately E of the WR star, which leads to solutions for the orbit inclination of 122+/-5 deg, the longitude of the ascending node of 353+/-3 deg, and an orbit semi-major axis of 9.0+/-0.5 mas. The distance to WR140 is determined to be 1.85+/-0.16 kpc, which requires the O star to be a supergiant. The inclination implies the mass of the WR and O star to be 20+/-4 and 54+/-10 solar masses respectively. We determine a wind-momentum ratio of 0.22, with an expected half-opening angle for the WCR of 63 deg, consistent with 65+/-10 deg derived from the VLBA observations. Total flux measurements from Very Large Array (VLA) observations show the radio emission from WR140 is very closely the same from one orbit to the next, pointing strongly toward emission, absorption and cooling mechanism(s) that are controlled largely by the orbital motion. The synchrotron spectra evolve dramatically through the orbital phases observed, exhibiting both optically thin and optically thick emission. We discuss a number of absorption and cooling mechanisms that may determine the evolution of the synchrotron spectrum with orbital phase.Comment: Accepted by ApJ, to appear in v623, April 20, 2005. 14 pages, 13 figs, requires emulateapj.cls. A version with full resolution figs can be obtained from http://www.drao.nrc.ca/~smd/preprint/wr140_data.pd

ISO spectroscopy of circumstellar dust in the Herbig Ae systems AB Aur and HD 163296

Using both the Short- and Long-wavelength Spectrometers on board the Infrared Space Observatory (ISO), we have obtained infrared spectra of the Herbig Ae systems AB Aur and HD 163296. In addition, we obtained ground-based N band images of HD 163296. Our results can be summarized as follows: (1) The main dust components in AB Aur are amorphous silicates, iron oxide and PAHs; (2) The circumstellar dust in HD 163296 consists of amorphous silicates, iron oxide, water ice and a small fraction of crystalline silicates; (3) The infrared fluxes of HD 163296 are dominated by solid state features; (4) The colour temperature of the underlying continuum is much cooler in HD 163296 than in AB Aur, pointing to the existence of a population of very large (mm sized) dust grains in HD 163296; (5) The composition and degree of crystallization of circumstellar dust are poorly correlated with the age of the central star. The processes of crystallization and grain growth are also not necessarily coupled. This means that either the evolution of circumstellar dust in protoplanetary disks happens very rapidly (within a few Myr), or that this evolution is governed by factors other than stellar mass and age.Comment: 6 pages, 2 figures, accepted for publication in Astronomy & Astrophysic

Discovery of extended radio emission in the young cluster Wd1

We present 10 micron, ISO-SWS and Australia Telescope Compact Array observations of the region in the cluster Wd1 in Ara centred on the B[e] star Ara C. An ISO-SWS spectrum reveals emission from highly ionised species in the vicinity of the star, suggesting a secondary source of excitation in the region. We find strong radio emission at both 3.5cm and 6.3cm, with a total spatial extent of over 20 arcsec. The emission is found to be concentrated in two discrete structures, separated by 14''. The westerly source is resolved, with a spectral index indicative of thermal emission. The easterly source is clearly extended and nonthermal (synchrotron) in nature. Positionally, the B[e] star is found to coincide with the more compact radio source, while the southerly lobe of the extended source is coincident with Ara A, an M2 I star. Observation of the region at 10micron reveals strong emission with an almost identical spatial distribution to the radio emission. Ara C is found to have an extreme radio luminosity in comparison to prior radio observations of hot stars such as O and B supergiants and Wolf-Rayet stars, given the estimated distance to the cluster. An origin in a detatched shell of material around the central star is therefore suggested; however given the spatial extent of the emission, such a shell must be relatively young (10^3 yrs). The extended non thermal emission associated with the M star Ara A is unexpected; to the best of our knowledge this is a unique phenomenon. SAX (2-10keV) observations show no evidence of X-ray emission, which might be expected if a compact companion were present.Comment: 5 pages including encapsulated figures, figure 3 separate. Accepted for MNRAS pink page

A Very Large Array 3.6cm continuum survey of Galactic Wolf-Rayet stars

We report the results of a survey of radio continuum emission of Galactic Wolf-Rayet stars north of declination -46 degrees. The observations were obtained at 8.46 GHz (3.6cm) using the Very Large Array (VLA), with an angular resolution of about 6 x 9 arcsec and typical rms noise of 0.04 mJy/beam. Our survey of 34 WR stars resulted in 15 definite and 5 probable detections, 13 of these for the first time at radio wavelengths. All detections are unresolved. Time variations in flux are confirmed in the cases of WR98a, WR104, WR105 and WR125. WR79a and WR89 are also variable in flux and we suspect they are also non-thermal emitters. Thus, of our sample 20-30% of the detected stars are non-thermal emiters. Average mass loss rates determinations obtained excluding definite and suspected non-thermal cases give similar values for WN (all subtypes) and WC5-7 stars, while a lower value was obtained for WC8-9 stars. Uncertainties in stellar distances largely contribute to the observed scatter in mass loss rates. Upper limits to the mass loss rates were obtained in cases of undetected sources or for sources which probably show additional non-thermal emission.Comment: 21 pages, Latex, 21 postscript figures, to be published in The Astronomical Journal, May 200

Non-thermal high-energy emission from colliding winds of massive stars

Colliding winds of massive star binary systems are considered as potential sites of non-thermal high-energy photon production. This is motivated merely by the detection of synchrotron radio emission from the expected colliding wind location. Here we investigate the properties of high-energy photon production in colliding winds of long-period WR+OB-systems. We found that in the dominating leptonic radiation process anisotropy and Klein-Nishina effects may yield spectral and variability signatures in the gamma-ray domain at or above the sensitivity of current or upcoming gamma-ray telescopes. Analytical formulae for the steady-state particle spectra are derived assuming diffusive particle acceleration out of a pool of thermal wind particles, and taking into account adiabatic and all relevant radiative losses. For the first time we include their advection/convection in the wind collision zone, and distinguish two regions within this extended region: the acceleration region where spatial diffusion is superior to convective/advective motion, and the convection region defined by the convection time shorter than the diffusion time scale. The calculation of the Inverse Compton radiation uses the full Klein-Nishina cross section, and takes into account the anisotropic nature of the scattering process. This leads to orbital flux variations by up to several orders of magnitude which may, however, be blurred by the geometry of the system. The calculations are applied to the typical WR+OB-systems WR 140 and WR 147 to yield predictions of their expected spectral and temporal characteristica and to evaluate chances to detect high-energy emission with the current and upcoming gamma-ray experiments. (abridged)Comment: 67 pages, 24 figures, submitted to Ap

Radio emission models of Colliding-Wind Binary Systems

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.Comment: 18 pages, 18 figures; Accepted by Astronomy and Astrophysics, July 8, 200

Computational core and fixed-point organisation in Boolean networks

In this paper, we analyse large random Boolean networks in terms of a constraint satisfaction problem. We first develop an algorithmic scheme which allows to prune simple logical cascades and under-determined variables, returning thereby the computational core of the network. Second we apply the cavity method to analyse number and organisation of fixed points. We find in particular a phase transition between an easy and a complex regulatory phase, the latter one being characterised by the existence of an exponential number of macroscopically separated fixed-point clusters. The different techniques developed are reinterpreted as algorithms for the analysis of single Boolean networks, and they are applied to analysis and in silico experiments on the gene-regulatory networks of baker's yeast (saccaromices cerevisiae) and the segment-polarity genes of the fruit-fly drosophila melanogaster.Comment: 29 pages, 18 figures, version accepted for publication in JSTA

Implications of Cognitive Load for Hypothesis Generation and Probability Judgment

We tested the predictions of HyGene (Thomas et al., 2008) that both divided attention at encoding and judgment should affect the degree to which participants’ probability judgments violate the principle of additivity. In two experiments, we showed that divided attention during judgment leads to an increase in subadditivity, suggesting that the comparison process for probability judgments is capacity limited. Contrary to the predictions of HyGene, a third experiment revealed that divided attention during encoding leads to an increase in later probability judgment made under full attention. The effect of divided attention during encoding on judgment was completely mediated by the number of hypotheses participants generated, indicating that limitations in both encoding and recall can cascade into biases in judgments