OH maser mapping of the evolved star HD 179821: evidence for interacting outflows

The evolved star HD 179821 continues to be the subject of much debate as to whether it is a nearby     post-asymptotic giant branch (post-AGB) star or a distant     high initial mass     post-red supergiant. We have mapped the OH maser emission around HD 179821 in the 1612- and 1667-MHz lines with the MERLIN interferometer array at a resolution of 0.4 arcsec and 0.35 km s⁻¹. The OH emission lies in a thick shell with inner and outer radii of 1.3 and         and expansion velocity of 30 km s⁻¹. Although we find some evidence for acceleration and for deviations from spherical symmetry, the bulk of the maser emission is consistent with a constant-velocity spherical shell. The extent of the shell agrees with H₂O and OH dissociation models and supports a distance estimate of 6 kpc. However, the shell is incomplete and appears to have been disrupted by more recent collimated outflow activity within the last 1500 yr. We suggest that this activity is also responsible for the active envelope chemistry (in particular the presence of HCO⁺) and for the apparent offset of the star from the centre of the shell. The luminous yellow hypergiant star IRC +10420 also shows signs of recent outflows, and HD 179821 may be at a similar, perhaps slightly earlier, phase of evolution. We suggest that the SiO thermal emission arises from the same detached envelope as the OH maser emission as in IRC +10420. If so then this would strengthen the connection between these two stars and probably rule out a post-AGB status for HD 179821.     

A detailed abundance analysis of the hot post-AGB star ZNG-1 in M10

We present a model-atmosphere analysis for the bright ( V ∼13) star ZNG-1, in the globular cluster M10. From high-resolution ( R ∼40 000) optical spectra we confirm ZNG-1 to be a post-asymptotic giant branch (post-AGB) star. The derived atmospheric parameters are T _eff=26 500±1000 K and log  g =3.6±0.2 dex . A differential abundance analysis reveals a chemical composition typical of hot post-AGB objects, with ZNG-1 being generally metal poor, although helium is approximately solar. The most interesting feature is the large carbon underabundance of more than 1.3 dex. This carbon deficiency, along with an observed nitrogen enhancement relative to other elements, may suggest that ZNG-1 evolved off the AGB before the third dredge-up occurred. Also, iron depletions observed in other similar stars suggest that gas–dust fractionation in the AGB progenitor could be responsible for the observed composition of these objects. However, we need not invoke either scenario since the chemical composition of ZNG-1 is in good agreement with abundances found for a Population II star of the same metallicity.     

Protodiscs around hot magnetic rotator stars

We develop equations and obtain solutions for the structure and evolution of a protodisc region that is initially formed with no radial motion and super-Keplerian rotation speed when wind material from a hot rotating star is channelled towards its equatorial plane by a dipole-type magnetic field. Its temperature is around 10⁷ K because of shock heating and the inflow of wind material causes its equatorial density to increase with time. The centrifugal force and thermal pressure increase relative to the magnetic force and material escapes at its outer edge. The protodisc region of a uniformly rotating star has almost uniform rotation and will shrink radially unless some instability intervenes. In a star with angular velocity increasing along its surface towards the equator, the angular velocity of the protodisc region decreases radially outwards and magnetorotational instability (MRI) can occur within a few hours or days. Viscosity resulting from MRI will readjust the angular velocity distribution of the protodisc material and may assist in the formation of a quasi-steady disc. Thus, the centrifugal breakout found in numerical simulations for uniformly rotating stars does not imply that quasi-steady discs with slow outflow cannot form around magnetic rotator stars with solar-type differential rotation.     

Modelling the spectra of colliding winds in the Wolf–Rayet WC7+O binaries WR 42 and WR 79

We have obtained complete phase coverage of the WC7+O binaries WR 42 = HD 97152 and WR 79 = HD 152270 with high signal-to-noise ratio (S/N), moderate-resolution spectra. Remarkable orbital phase-locked profile variations of the C  iii λ 5696 line are observed and interpreted as arising from colliding wind effects. Within this scenario, we have modelled the spectra using a purely geometrical model that assumes a cone-shaped wind–wind interaction region which partially wraps around the O star. Such modelling holds the exciting promise of revealing a number of interesting parameters for WR+O binaries, such as the orbital inclination, the streaming velocity of material in the interaction region and the ratio of wind momentum flux. Knowledge of these parameters in turn leads to the possibility of a better understanding of WR star masses, mass-loss rates and wind region characteristics.     

The interaction of planetary nebulae and their asymptotic giant branch progenitors with the interstellar medium

Interaction with the interstellar medium (ISM) cannot be ignored in understanding planetary nebula (PN) evolution and shaping. In an effort to understand the range of shapes observed in the outer envelopes of PNe, we have run a comprehensive set of three-dimensional hydrodynamic simulations, from the beginning of the asymptotic giant branch (AGB) superwind phase until the end of the post-AGB/PN phase. A 'triple-wind' model is used, including a slow AGB wind, fast post-AGB wind and third wind reflecting the linear movement through the ISM. A wide range of stellar velocities, mass-loss rates and ISM densities have been considered.
We find that ISM interaction strongly affects outer PN structures, with the dominant shaping occurring during the AGB phase. The simulations predict four stages of PN–ISM interaction whereby (i) the PN is initially unaffected, (ii) then limb-brightened in the direction of motion, (iii) then distorted with the star moving away from the geometric centre, and (iv) finally so distorted that the object is no longer recognizable as a PN and may not be classed as such. Parsec-size shells around PNe are predicted to be common. The structure and brightness of ancient PNe are largely determined by the ISM interaction, caused by rebrightening during the second stage; this effect may address the current discrepancies in Galactic PN abundance. The majority of PNe will have tail structures. Evidence for strong interaction is found for all known PNe in globular clusters.     

Dust formation and inhomogeneous mass-loss from asymptotic giant branch stars

We examine the flow from asymptotic giant branch (AGB) stars when along a small solid angle the optical depth resulting from dust is very large. We consider two types of flows. In the first, small cool spots are formed on the surface of slowly rotating AGB stars. Large quantities of dust are expected to be formed above the surface of these cool spots. We propose that if the dust formation occurs during the last AGB phase when the mass-loss rate is high, the dust shields the region above it from the stellar radiation. This leads to both further dust formation in the shaded region and, owing to lower temperature and pressure, the convergence of the stream toward the shaded region, and the formation of a flow having a higher density than its surroundings. This density contrast can be as high as ∼4. A concentration of magnetic cool spots toward the equator will lead to a density contrast of up to a few between the equatorial and polar directions. This process can explain the positive correlation between high mass-loss rate and a larger departure from sphericity in progenitors of elliptical planetary nebulae. In the second type of flow, the high density in the equatorial plane is formed by a binary interaction, where the secondary star is close to, but outside the AGB envelope. The shielding of the radiation by dust results in a very slow and dense flow in the equatorial plane. We suggest this flow as an alternative explanation for the equatorial dense matter found at several hundred astronomical units around several post-AGB binary systems.     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

The enigma of the oldest 'nova': the central star and nebula of CK Vul

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova or a very late thermal pulse object. Following its eruption in ad 1670, the star remained visible for 2 yr. A 15-arcsec nebula was discovered in the 1980s, but the star itself has not been detected since the eruption. We here present radio images which reveal a 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Hα images show a bipolar nebula with a longest extension of 70 arcsec, with the previously known compact nebula at its waist. The emission-line ratios show that the gas is shock-ionized, at velocities  >100 km s⁻¹  . Dust emission yields an envelope mass of  ∼5 × 10⁻² M_⊙  . Echelle spectra indicate outflow velocities up to 360 km s⁻¹. From a comparison of images obtained in 1991 and 2004 we find evidence for expansion of the nebula, consistent with an origin in the 1670 explosion; the measured expansion is centred on the radio source. No optical or infrared counterpart is found at the position of the radio source. The radio emission is interpreted as thermal free–free emission from gas with   T _e∼ 10⁴ K  . The radio source may be due to a remnant circumbinary disc, similar to those seen in some binary post-AGB stars. We discuss possible classifications of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova eruption on a cool, low-mass white dwarf or a thermal pulse induced by accretion from a circumbinary disc.     

Magnetized relativistic jets and long-duration GRBs from magnetar spin-down during core-collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond,   B ∼ 10¹⁵ G  , magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers  (∼10⁵¹–10⁵² erg s⁻¹)  , the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers  (∼10⁵⁰ erg s⁻¹)  , the wind is subfast but still super-Alfvénic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≲30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated  (∼5–10°)  relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.     

X-ray emission by a shocked fast wind from the central stars of planetary nebulae

We calculate the X-ray emission from the shocked fast wind blown by the central stars of planetary nebulae (PNe) and compare with observations. Using spherically symmetric self-similar solutions, we calculate the flow structure and X-ray temperature for a fast wind slamming into a previously ejected slow wind. We find that the observed X-ray emission of six PNe can be accounted for by shocked wind segments that were expelled during the early-PN phase, if the fast wind speed is moderate,   v ₂∼ 400–600 km s⁻¹  , and the mass-loss rate is a few times  10⁻⁷ M_⊙ yr⁻¹  . We find, as proposed previously, that the morphology of the X-ray emission is in the form of a narrow ring inner to the optical bright part of the nebula. The bipolar X-ray morphology of several observed PNe, which indicates an important role of jets, rather than a spherical fast wind, cannot be explained by the flow studied here.     

The circumstellar environment of Wolf–Rayet stars and gamma-ray burst afterglows

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long gamma-ray bursts (GRBs). We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind/stalled-wind interface, rebrightening occurs and a bump is seen in the afterglow light curve. However, our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final wind density,   A _*  , from our stellar models are of the same order (≲1) as some of the values inferred from observed afterglow light curves. We do not reproduce the lowest   A _*  values below 0.5 inferred from afterglow observations. For these cases, we suggest that the progenitors could have been a WO-type Wolf–Rayet (WR) star or a very low-metallicity star. Finally, we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two WR stellar types, WC and WO, in the afterglow light curve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.     

Stellar wind accretion in GX 301−2: evidence for a high-density stream

The X-ray binary system GX 301−2 consists of a neutron star in an eccentric orbit accreting from the massive early-type star Wray 977. It has previously been shown that the X-ray orbital light curve is consistent with the existence of a gas stream flowing out from Wray 977 in addition to its strong stellar wind. Here, X-ray monitoring observations by the Rossi X-ray Timing Explorer ( RXTE )/All-Sky Monitor and pointed observations by the RXTE /Proportional Counter Array over the past decade are analysed. We analyse both the flux and column density dependence on orbital phase. The wind and stream dynamics are calculated for various system inclinations, companion rotation rates and wind velocities, as well as parametrized by the stream width and density. These calculations are used as inputs to determine both the expected accretion luminosity and the column density along the line-of-sight to the neutron star. The model luminosity and column density are compared to observed flux and column density versus orbital phase, to constrain the properties of the stellar wind and the gas stream. We find that the change between bright and medium intensity levels is primarily due to decreased mass loss in the stellar wind, but the change between medium and dim intensity levels is primarily due to decreased stream density. The mass-loss rate in the stream exceeds that in the stellar wind by a factor of ∼2.5. The quality of the model fits is better for lower inclinations, favouring a higher mass for Wray 977 in its allowed range of  40–60 M_⊙  .     

Violent stellar merger model for transient events

We derive the constraints on the mass ratio for a binary system to merge in a violent process. We find that the secondary-to-primary stellar mass ratio should be  0.003 ≲ ( M ₂/ M ₁) ≲ 0.15  . A more massive secondary star will keep the primary stellar envelope in synchronized rotation with the orbital motion until merger occurs. This implies a very small relative velocity between the secondary star and the primary stellar envelope at the moment of merger, and therefore very weak shock waves, and low-flash luminosity. A too low-mass secondary will release small amount of energy, and will expel small amount of mass, which is unable to form an inflated envelope. It can, however, produce a quite luminous but short flash when colliding with a low-mass main-sequence star.
Violent and luminous mergers, which we term mergebursts, can be observed as V838 Monocerotis-type events, where a star undergoes a fast brightening lasting days to months, with a peak luminosity of up to  ∼10⁶ L_⊙  followed by a slow decline at very low effective temperatures.     

X-ray emission from planetary nebulae calculated by 1D spherical numerical simulations

We calculate the X-ray emission from both constant and time-evolving shocked fast winds blown by the central stars of planetary nebulae (PNe) and compare our calculations with observations. Using spherically symmetric numerical simulations with radiative cooling, we calculate the flow structure and the X-ray temperature and luminosity of the hot bubble formed by the shocked fast wind. We find that a constant fast wind gives results that are very close to those obtained from the self-similar solution. We show that in order for a fast shocked wind to explain the observed X-ray properties of PNe, rapid evolution of the wind is essential. More specifically, the mass-loss rate of the fast wind should be high early on when the speed is  ∼300–700 km s⁻¹  , and then it needs to drop drastically by the time the PN age reaches ∼1000 yr. This implies that the central star has a very short pre-PN (post-asymptotic giant branch) phase.     

UBV photometry, UV spectroscopy and radio observations of the peculiar binary V Sagittae

We have performed high-speed UBV photometric observations on the peculiar binary V Sagittae. Using three new eclipse timings we update the orbital ephemeris and convert it to a dynamical time-scale (TDB). We also searched for quasi-periodic oscillations but did not detect them. Using the Wilson–Devinney algorithm we have modelled the light curve to find the stellar parameters of V Sge. We find that the system is a detached binary but that the primary star is very close to filling its Roche lobe, while the secondary star fills 90 per cent of its Roche lobe volume. We find temperatures of the primary and the secondary star to be T ₁=41 000 K and T ₂=22 000 K. We find i =72° and masses of 0.8 M_⊙ and 3.3 M_⊙ for the primary and secondary stars respectively. De-archived Hubble Space Telescope ( HST ) spectroscopy of V Sge shows evidence of mass loss via a wind or winds. In addition we report radio observations of V Sge during an optical high state at 2 cm, 3.6 cm and 6 cm wavelengths. The 3.6 cm emission is increased by a factor of more than six compared with an earlier detection in a previous optical high state.     

A first detailed study of the colliding wind WR+O binary WR 30a

We present a detailed, extensive investigation of the photometric and spectroscopic behaviour of WR 30a. This star is definitely a binary system with a period around 4.6 d. We propose the value         . The identification of the components as WO4+O5((f)) indicates a massive evolved binary system; the O5 component is a main-sequence or, more likely, a giant star. The radial velocities of the O star yield a circular orbit with an amplitude         and a mass function of 0.013     . The spectrum of WR 30a exhibits strong profile variations of the broad emission lines that are phase-locked with the orbital period. We report the detection of the orbital motion of the WO component with     , but this should be confirmed by further observations. If correct, it implies a mass ratio     . The star exhibits sinusoidal light variations of amplitude 0.024 mag peak-to-peak with the minimum of light occurring slightly after the conjunction with the O star in front. On the basis of the phase-locked profile variations of the C  iv λ 4658 blend in the spectrum of the WO, we conclude that a wind–wind collision phenomenon is present in the system. We discuss some possibilities for the geometry of the interaction region.     

ROSAT observations of V471 Tauri, showing that stellar activity is determined by rotation, not age

I present pointed ROSAT PSPC observations of the pre-cataclysmic binary V471 Tauri. The hard X-ray emission (>0.4 keV) is not eclipsed by the K star, demonstrating conclusively that this component cannot be emitted by the white dwarf. Instead I show that its spectrum and luminosity are consistent with coronal emission from the tidally spun-up K star. The star is more active than other K stars in the Hyades, but equally active as K stars in the Pleiades with the same rotation periods, demonstrating that rotation — and not age — is the key parameter in determining the level of stellar activity.   The soft X-ray emission (<0.4 keV) is emitted predominately by the white dwarf and is modulated on its spin period. I find that the pulse profile is stable on time-scales of hours and years, supporting the idea that it is caused by the opacity of accreted material. The profile itself shows that the magnetic field configuration of the white dwarf is dipolar and that the magnetic axis passes through the centre of the star.   There is an absorption feature in the light curve of the white dwarf, which occurs at a time when our line of sight passes within a stellar radius of the K star. The column density and duration of this feature imply a volume and mass for the absorber that are similar to those of coronal mass ejections of the Sun.   Finally I suggest that the spin–orbit beat period detected in the optical by Clemens et al. may be the result of the interaction of the K-star wind with the magnetic field of the white dwarf.     

On the first generation of stars

We argue that the first stars may have spanned the conventional mass range rather than be identified with the very massive objects  (∼100–10³ M_⊙)  favoured by numerical simulations. Specifically, we find that magnetic field generation processes acting in the first protostellar systems suffice to produce fields that exceed the threshold for magneto-rotational instability (MRI) to operate, and thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on protostellar mass scales below  ∼50 M_⊙  . Such fields allow primordial star formation to occur at essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion and feedback in much the same way as occurs in conventional molecular clouds.     

High-resolution spectroscopic observations of post-asymptotic giant branch candidates from the Edinburgh–Cape Survey

High spectral resolution  ( R ∼ 40 000)  and signal-to-noise optical spectra, obtained at the Very Large Telescope (VLT), are presented for three post–asymptotic giant branch (AGB) candidates selected from the Edinburgh–Cape (EC) Faint Blue Object Survey. The stellar atmospheric parameters and chemical compositions, derived using sophisticated non-local thermodynamic equilibrium calculations, reveal that EC 14102-1337 and EC 20068-7324 are both in an evolved post–horizontal branch (HB) evolutionary state. However, EC 11507-2253 is most likely a post-AGB star.     

The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities   Z = 10⁻⁴  and 10⁻⁸. We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of   Z = 10⁻⁵  for any mass of star and that stars above  2 M_⊙  do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen-burning region, occurs in the mass range  2 M_⊙  to around  4 M_⊙  . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of Carbon-Enhanced Metal-Poor (CEMP) stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge-up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.