An extremely slow nova?

The binary companion to the peculiar F supergiant HD 172481 is shown to be a Mira variable with a pulsation period of 312 d. Its characteristics are within the normal range found for solitary Miras of that period, although its pulsation amplitude and mass-loss rate ̇ ∼3×10⁻⁶ M_⊙ yr⁻¹ are higher than average. Reasons are given for suspecting that the F supergiant, which has L ∼10⁴ L_⊙, is a white dwarf burning hydrogen accreted from its companion.     

Sakurai's object: the ionized nebula at radio wavelengths

Sakurai's object (V4334 Sgr) is a planetary nebula nucleus which is undergoing its final helium shell flash. This is the first of these rare and important events to be observable with non-optical instruments. We report the first radio detection, using a short (2-h) observation with the Very Large Array (VLA) at 4.86 GHz. The radio emission structure is coincident with the 34-arcsec diameter planetary nebula seen in optical emission lines. We find a statistical distance ∼ 3.8 ± 0.6 kpc, with a range of 1.9 <  D  < 5.3 kpc, depending on the planetary nebula (PN) mass. While we have no direct evidence for a new (post-flash) stellar wind, we estimate an upper limit to the mass-loss rate due to any such wind of 1.7 × 10⁻⁷ M⊙ yr⁻¹. The number of emitting knots in the radio-visible nebula indicates an electron density of ∼ 2 × 10⁸ m⁻³ in those knots, and a total emitting ionized mass of ∼ 0.15 M⊙, at an assumed distance of 3.8 kpc. The radio flux density indicates an Hβ flux of ∼ 6 × 10⁻¹⁶ W m⁻², suggesting an extinction E ( B  −  V ) ∼ 1.15, comparable with reddening estimates in the direction of V4334 Sgr.     

Cygnus X-3 with ISO: investigating the wind

We have observed the energetic binary Cygnus X-3 in both quiescent and flaring states between 4 and 16 μm using the ISO satellite. We find that the quiescent source shows the thermal free–free spectrum typical of a hot, fast stellar wind, such as from a massive helium star. The quiescent mass-loss rate arising from a spherically symmetric, non-accelerating wind is found to be in the range (0.4–2.9)×10⁻⁴ M_⊙ yr⁻¹, consistent with other infrared and radio observations, but considerably larger than the 10⁻⁵ M_⊙ yr⁻¹ deduced from both the orbital change and the X-ray column density. There is rapid, large-amplitude flaring at 4.5 and 11.5 μm at the same time as enhanced radio and X-ray activity, with the infrared spectrum apparently becoming flatter in the flaring state. We believe that non-thermal processes are operating, perhaps along with enhanced thermal emission.     

Near-infrared polarimetry and modelling of the dusty young planetary nebula IRAS 19306+1407

We present near-infrared polarimetric images of the dusty circumstellar envelope (CSE) of IRAS 19306+1407, acquired at the United Kingdom Infrared Telescope (UKIRT) using the UKIRT 1–5 μm Imager Spectrometer (UIST) in conjunction with the half-waveplate module IRPOL2. We present additional 450- and 850-μm photometry data obtained with the Submillimetre Common-User Bolometer Array (SCUBA) at the James Clerk Maxwell Telescope (JCMT), as well as archived Hubble Space Telescope ( HST ) F606W - and F814W -filter images. The CSE structure in polarized flux at J and K bands shows an elongation north of north-east and south of south-west with two bright scattering shoulders north-west and south-east. These features are not perpendicular to each other and could signify a recent 'twist' in the outflow axis. We model the CSE using an axisymmetric light scattering ( als ) code to investigate the polarization produced by the CSE, and an axisymmetric radiation transport ( dart ) code to fit the spectral energy distribution. A good fit was achieved with the als and dart models using silicate grains, 0.1–0.4 μm with a power-law size distribution of a ^−3.5, and an axisymmetric shell geometry with an equator-to-pole ratio of 7:1. The spectral type of the central star is determined to be B1 i supporting previous suggestions that the object is an early planetary nebula. We have constrained the CSE and interstellar extinction as 2.0 and 4.2 mag, respectively, and have estimated a distance of 2.7 kpc. At this distance, the stellar luminosity is ∼4500 L_⊙ and the mass of the CSE is ∼0.2 M_⊙. We also determine that the mass loss lasted for ∼5300 yr with a mass-loss rate of ∼3.4 × 10⁻⁵ M_⊙ yr⁻¹.     

Type IIP supernova SN 2004et: a multiwavelength study in X-ray, optical and radio

We present X-ray, broad-band optical and low-frequency radio observations of the bright type IIP supernova SN 2004et. The Chandra X-ray Observatory observed the supernova at three epochs, and the optical coverage spans a period of ∼470 d since explosion. The X-ray emission softens with time, and we characterize the X-ray luminosity evolution as   L _X∝ t ^−0.4  . We use the observed X-ray luminosity to estimate a mass-loss rate for the progenitor star of  ∼2 × 10⁻⁶ M_⊙ yr⁻¹  . The optical light curve shows a pronounced plateau lasting for about 110 d. Temporal evolution of photospheric radius and colour temperature during the plateau phase is determined by making blackbody fits. We estimate the ejected mass of ⁵⁶Ni to be  0.06 ± 0.03 M_⊙  . Using the expressions of Litvinova & Nadëzhin we estimate an explosion energy of  (0.98 ± 0.25) × 10⁵¹ erg  . We also present a single epoch radio observation of SN 2004et. We compare this with the predictions of the model proposed by Chevalier, Fransson & Nymark. These multiwavelength studies suggest a main-sequence progenitor mass of  ∼20 M_⊙  for SN 2004et.     

The expansion of the giant filamentary shell, DEM 171, in the Magellanic Bridge

The giant filamentary shell, DEM 171, is found to be expanding at approximately 37 km s⁻¹. A supernova and stellar wind origin are both explored as possible causes for the expanding shell. A stellar wind origin would imply a mass-loss rate of the order of 10⁻⁵ M_⊙ yr⁻¹, indicating that it could be caused by a Wolf–Rayet star. A number of blue stars are found to lie within the shell and one is identified as a Wolf–Rayet candidate.     

Supernova 1994aj: a probe for pre-supernova evolution and mass-loss from the progenitor

Extensive photometric and spectroscopic observations of SN 1994aj until 540 d after maximum light have been obtained. The photometry around maximum suggests that the SN belongs to the Type II Linear class, with a peak absolute magnitude of M _V∼−17.8 (assuming H ₀=75 km s⁻¹ Mpc⁻¹). The spectra of SN 1994aj were unusual, with the presence of a narrow line with a P Cygni profile on top of the broad Balmer line emission. This narrow feature is attributed to the presence of a dense superwind surrounding the SN. At 100–120 d after maximum light the SN ejecta start to interact with this circumstellar material. The SN luminosity decline rates slowed down [γ _R =0.46 mag (100 d)⁻¹], becoming less steep than the average late luminosity decline of normal SN II [∼1 mag (100 d)⁻¹]. This dense ( ˙M / u _W∼10¹⁵ g cm⁻¹) wind was confined to a short distance from the progenitor ( R _out=∼5×10¹⁶ cm), and results from a very strong mass-loss episode ( ˙M =10⁻³ M⊙ yr⁻¹), which terminated shortly before explosion (∼5–10 yr).     

Peering through the stellar wind of IGR J19140+0951 with simultaneous INTEGRAL/RXTE observations

We have used the RXTE and INTEGRAL satellites simultaneously to observe the high-mass X-ray binary (HMXB) IGR J19140+0951. The spectra obtained in the 3–80 keV range have allowed us to perform a precise spectral analysis of the system along its binary orbit. The spectral evolution confirms the supergiant nature of the companion star and the neutron star nature of the compact object. Using a simple stellar wind model to describe the evolution of the photoelectric absorption, we were able to restrict the orbital inclination angle in the range 38°–75°. This analysis leads to a wind mass-loss rate from the companion star of  ∼5 × 10⁻⁸ M_⊙ yr⁻¹  , consistent with an OB I spectral type. We have detected a soft excess in at least four observations, for the first time for this source. Such soft excesses have been reported in several HMXBs in the past. We discuss the possible origin of this excess, and suggest, based on its spectral properties and occurrences around the superior conjunction, that it may be explained as the reprocessing of the X-ray emission originating from the neutron star by the surrounding ionized gas.     

Theoretical X-ray line profiles from colliding wind binaries

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 _∞)  and definite blueshifts 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 _∞)  , 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⁻⁶ M_⊙ 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 linewidths and shifts is reduced in a low-inclination binary. The extreme case is a binary with   i = 0°  , for which we would expect no line profile variation.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

Magnetically collimated jets with high mass flux

We present HST /WFPC2 observations of UGC 4483, an irregular galaxy in the M81/NGC 2403 complex. Stellar photometry was carried out with HSTphot, and is complete to V ≃26.0 and I ≃24.7. We measure the red giant branch tip at I =23.56±0.10, and calculate a distance modulus of μ ₀=27.53±0.12 (corresponding to a distance of 3.2±0.2 Mpc), placing UGC 4483 within the NGC 2403 subgroup. We were able to measure properties of a previously known young star cluster in UGC 4483, finding integrated magnitudes of V =18.66±0.21 and I =18.54±0.10 for the stellar contribution (integrated light minus H α and [O  iii ] contribution), corresponding to an age of ∼10–15 Myr and an initial mass of ∼10⁴ M_⊙. This is consistent with the properties of the cluster's brightest stars, which were resolved in the data for the first time. Finally, a numerical analysis of the galaxy's stellar content yields a roughly constant star formation rate of 1.3×10⁻³ M_⊙ yr⁻¹ and mean metallicity of [Fe/H]=−1.3 dex from 15 Gyr ago to the present.     

Silicate and hydrocarbon emission from Galactic M supergiants

Following our discovery of unidentified infrared (UIR) band emission in a number of M supergiants in h and χ Per, we have obtained 10-μm spectra of a sample of 60 galactic M supergiants. Only three new sources, V1749 Cyg, UW Aql and IRC+40 427, appear to show the UIR bands; the others show the expected silicate emission or a featureless continuum. The occurrence of UIR-band emission in M supergiants is therefore much higher in the h and χ Per cluster than in the Galaxy as a whole. Possible explanations for the origin and distribution of UIR bands in oxygen-rich supergiants are discussed. We use our spectra to derive mass-loss rates ranging from 10⁻⁸ to 10⁻⁴ M_⊙ yr⁻¹ for the new sample, based on the power emitted in the silicate feature. The relationship between mass-loss rate and luminosity for M supergiants is discussed, and correlations are explored between their mid-infrared emission properties.     

A binary merger model for the formation of the Supernova 1987A triple-ring nebula

We examine a binary merger model for the formation of the mysterious triple-ring nebula surrounding Supernova 1987A, which still has not been convincingly explained in detailed hydrodynamical calculations. During the merger of 15 and  5 M_⊙  binary systems, mass is ejected primarily at mid-latitudes for a sufficiently evolved primary, as demonstrated by Morris & Podsiadlowski. This material is swept up by the fast wind of the central star during its post-merger blue supergiant phase, leading to a density contrast of ∼150 in the outer rings at the time of the supernova. The equatorial ring probably formed later when the star contracted to become a blue supergiant. The asymmetry between the northern and southern outer rings can be explained by a 10 per cent asymmetry during the merger, perhaps due to a pulsational instability in the common envelope.
We present a parameter study from which we determine a mass-loss rate in the blue supergiant wind in the range  1.5–3 × 10⁻⁷ M_⊙ yr⁻¹  in agreement with previous estimates. The morphology of the best model is consistent with the well-known Hubble Space Telescope image at better than 5 per cent and is also in broad agreement with light-echo observations. The circumstellar environment on larger scales (up to 3 pc) is also investigated. We conclude with a brief discussion of the bipolar nebulae surrounding the Galactic stars, Sheridan 25, HD 168625 and η Carinae.     

Spitzer observations of acetylene bands in carbon-rich asymptotic giant branch stars in the Large Magellanic Cloud

We investigate the molecular bands in carbon-rich asymptotic giant branch (AGB) stars in the Large Magellanic Cloud (LMC), using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope ( SST ) over the 5–38 μm range. All 26 low-resolution spectra show acetylene (C₂H₂) bands at 7 and 14 μm. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 μm C₂H₂  band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 μm C₂H₂  band by circumstellar dust emission. This 14-μm band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C₂H₂ column density, derived from the 13.7 μm band, shows a gas mass-loss rate in the range 3 × 10⁻⁶ to 5 × 10⁻⁵ M_⊙ yr⁻¹. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 μm C₂H₂  band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C₂H₂  abundance among LMC and Galactic stars. This reflects the effect of the third dredge-up bringing self-produced carbon to the surface, leading to high carbon-to-oxygen ratio at low metallicity.     

On the late spectral types of cataclysmic variable secondaries

We investigate why the spectral type of most cataclysmic variable (CV) secondaries is significantly later than that of a zero-age main-sequence (ZAMS) star with the same mean density. Using improved stellar input physics, tested against observations of low-mass stars at the bottom of the main sequence, we calculate the secular evolution of CVs with low-mass donors. We consider sequences with different mass transfer rates and with a different degree of nuclear evolution of the donor prior to mass transfer.
Systems near the upper edge of the gap ( P ∼3–6 h) can be reproduced by models with a wide range of mass transfer rates from 1.5×10⁻⁹ M_⊙ yr⁻¹ to 10⁻⁸ M_⊙ yr⁻¹. Evolutionary sequences with a small transfer rate and donors that are substantially evolved off the ZAMS (central hydrogen content 0.05–0.5) reproduce CVs with late spectral types above P ≳6 h. Systems with the most discrepant (late) spectral type should have the smallest donor mass at any given P .
Consistency with the period gap suggests that the mass transfer rate increases with decreasing donor mass for evolved sequences above the period gap. In this case, a single-parameter family of sequences with varying X _c and increasing mass transfer rate reproduces the full range of observed spectral types. This would imply that CVs with such evolved secondaries dominate the CV population.     

The local star formation rate and radio luminosity density

We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×10²² W Hz⁻¹  and density   φ ∗=(4.8±1.1)×10⁻⁴ Mpc⁻³.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×10¹⁹ W Hz⁻¹ Mpc⁻³  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ _1.4 GHz=(2.10±0.45±0.04)×10⁻² M_⊙ yr⁻¹ Mpc⁻³  for a Salpeter initial mass function from  0.1–125 M_⊙  and Hubble constant of 50 km s⁻¹ Mpc⁻¹. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

X-ray analysis of the cluster 3C 129

The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×10⁴⁴ erg s⁻¹ in 0.2–2.4 keV and 2.7×10⁴⁴ erg s⁻¹ in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M_⊙ yr⁻¹. The central gas density is 6×10⁻³ cm⁻³, and the ICM mass is 3.6×10¹³ M_⊙. The total cluster mass is ∼5×10¹⁴ M_⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.