A fast solver for systems of reaction-diffusion equations

In this paper we present a fast algorithm for the numerical solution of systems of reaction-diffusion equations, $\partial_t u + a \cdot \nabla u = \Delta u + F (x, t, u)$, $x \in \Omega \subset \mathbf{R}^3$, $t > 0$. Here, $u$ is a vector-valued function, $u \equiv u(x, t) \in \mathbf{R}^m$, $m$ is large, and the corresponding system of ODEs, $\partial_t u = F(x, t, u)$, is stiff. Typical examples arise in air pollution studies, where $a$ is the given wind field and the nonlinear function $F$ models the atmospheric chemistry.Comment: 8 pages, 3 figures, to appear in Proc. 13th Domain Decomposition Conference, Lyon, October 200

The radial velocity curve of HD153919 (4U1700-37) revisited

We have re-analysed all available high-resolution ultraviolet IUE spectra of the high-mass X-ray binary HD153919/4U1700-37. The radial velocity semi-amplitude of 20.6 +/- 1.0 km/s and orbital eccentricity of 0.22 +/- 0.04 agree very well with the values obtained earlier from optical spectra. They disagree with earlier conclusions for the same data reduced by Heap & Corcoran (1992) and by Stickland & Lloyd (1993).Comment: 6 pages, latex, figure included, Astronomy & Astrophysics, in pres

The Primordial Binary Population in OB Associations

For understanding the process of star formation it is essential to know how many stars are formed as singles or in multiple systems, as a function of environment and binary parameters. This requires a characterization of the primordial binary population, which we define as the population of binaries that is present just after star formation has ceased, but before dynamical and stellar evolution have significantly altered its characteristics. In this article we present the first results of our adaptive optics survey of 200 (mainly) A-type stars in the nearby OB association Sco OB2. We report the discovery of 47 new candidate companions of Sco OB2 members. The next step will be to combine these observations with detailed simulations of young star clusters, in order to find the primordial binary population.Comment: 2 pages, 1 figure, poster paper to appear in proceedings of IAU Coll. 191 "The environments and evolution of binary and multiple stars

The properties of single WO stars

The enigmatic oxygen sequence Wolf-Rayet (WO) stars represent a very late stage in massive star evolution, although their exact nature is still under debate. The spectra of most of the WO stars have never been analysed through detailed modelling with a non-local thermodynamic equilibrium expanding atmosphere code. Here we present preliminary results of the first homogeneous analysis of the (apparently) single WOs.Comment: 2 pages, 1 figure, Proceedings IAU Symposium No. 307, 2014, 'New windows on massive stars: asteroseismology, interferometry, and spectropolarimetry

The origin of the runaway high-mass X-ray binary HD153919/4U1700-37

Based on its Hipparcos proper motion, we propose that the high-mass X-ray binary HD153919/4U1700-37 originates in the OB association Sco OB1. At a distance of 1.9 kpc the space velocity of 4U1700-37 with respect to Sco OB1 is 75 km/s. This runaway velocity indicates that the progenitor of the compact X-ray source lost about 7 Msun during the (assumed symmetric) supernova explosion. The system's kinematical age is about 2 +/- 0.5 million years which marks the date of the supernova explosion forming the compact object. The present age of Sco OB1 is <8 Myr; its suggested core, NGC 6231, seems to be somewhat younger (~5 Myr). If HD153919/4U1700-37 was born as a member of Sco OB1, this implies that the initially most massive star in the system terminated its evolution within 30 Msun. With these parameters the evolution of the binary system can be constrained.Comment: 6 pages, latex, 3 embedded ps figures, to appear in A&

The N Enrichment and Supernova Ejection of the Runaway Microquasar LS 5039

We present an investigation of new optical and ultraviolet spectra of the mass donor star in the massive X-ray binary LS 5039. The optical band spectral line strengths indicate that the atmosphere is N-rich and C-poor, and we classify the stellar spectrum as type ON6.5 V((f)). The N-strong and C-weak pattern is also found in the stellar wind P Cygni lines of N V 1240 and C IV 1550. We suggest that the N-enrichment may result from internal mixing if the O-star was born as a rapid rotator, or the O-star may have accreted N-rich gas prior to a common-envelope interaction with the progenitor of the supernova. We re-evaluated the orbital elements to find an orbital period of P=4.4267 +/- 0.0010 d. We compared the spectral line profiles with new non-LTE, line-blanketed model spectra, from which we derive an effective temperature T_eff = 37.5 +/- 1.7 kK, gravity log g = 4.0 +/- 0.1, and projected rotational velocity V sin i = 140 +/- 8 km/s. We fit the UV, optical, and IR flux distribution using a model spectrum and extinction law with parameters E(B-V)= 1.28 +/- 0.02 and R= 3.18 +/- 0.07. We confirm the co-variability of the observed X-ray flux and stellar wind mass loss rate derived from the H-alpha profile, which supports the wind accretion scenario for the X-ray production in LS 5039. Wind accretion models indicate that the compact companion has a mass M_X/M_sun = 1.4 +/- 0.4, consistent with its identification as a neutron star. The observed eccentricity and runaway velocity of the binary can only be reconciled if the neutron star received a modest kick velocity due to a slight asymmetry in the supernova explosion (during which >5 solar masses was ejected).Comment: 38 pages, 9 figures; 2004, ApJ, 600, Jan. 10 issue, in press Discussion revised thanks to comments from P. Podsiadlowsk

The Rich Mid-Infrared Environments of Two Highly-Obscured X-ray Binaries: Spitzer Observations of IGR J16318-4848 and GX 301-2

We present the results of Spitzer mid-infrared spectroscopic observations of two highly-obscured massive X-ray binaries: IGR J16318-4848 and GX301-2. Our observations reveal for the first time the extremely rich mid-infrared environments of this type of source, including multiple continuum emission components (a hot component with T > 700 K and a warm component with T ~ 180 K) with apparent silicate absorption features, numerous HI recombination lines, many forbidden ionic lines of low ionization potentials, and pure rotational H2 lines. This indicates that both sources have hot and warm circumstellar dust, ionized stellar winds, extended low-density ionized regions, and photo-dissociated regions. It appears difficult to attribute the total optical extinction of both sources to the hot and warm dust components, which suggests that there could be an otherwise observable colder dust component responsible for the most of the optical extinction and silicate absorption features. The observed mid-infrared spectra are similar to those from Luminous Blue Variables, indicating that the highly-obscured massive X-ray binaries may represent a previously unknown evolutionary phase of X-ray binaries with early-type optical companions. Our results highlight the importance and utility of mid-infrared spectroscopy to investigate highly-obscured X-ray binaries.Comment: To appear in ApJ Letter