Spectroscopy of the γ-ray burst GRB 021004: a structured jet ploughing through a massive stellar wind

We present spectra of the afterglow of the γ-ray burst GRB 021004 taken with the ISIS spectrograph on the William Herschel Telescope (WHT) and with the Focal Reducer/Low Dispersion Spectrograph 1 (FORS1) on the Very Large Telescope (VLT) at three epochs spanning 0.49–6.62 d after the burst. We observe strong absorption probably coming from the host galaxy, alongside absorption in H  i , Si  iv and C  iv with blueshifts of up to 2900 km s⁻¹ from the explosion centre, which we assume originates close to the progenitor. We find no significant variability of these spectral features. We investigate the origin of the outflowing material and evaluate various possible progenitor models. The most plausible explanation is that these result in the fossil stellar wind of a highly evolved Wolf–Rayet (WR) star. However, ionization from the burst itself prevents the existence of H  i , Si  iv and C  iv close to the afterglow surface where the fast stellar wind should dominate, and large amounts of blueshifted hydrogen are not expected in a WR star wind. We propose that the WR star wind is enriched by a hydrogen-rich companion, and that the GRB has a structured jet geometry in which the γ-rays emerge in a small opening angle within the wider opening angle of the cone of the afterglow. This scenario is able to explain both the spectral-line features and the irregular light curve of this afterglow.     

Shock breakout in Type Ibc supernovae and application to GRB 060218/SN 2006aj

Recently, a soft blackbody component was observed in the early X-ray afterglow of GRB 060218, which was interpreted as shock breakout from the thick wind of the progenitor Wolf–Rayet (WR) star of the underlying Type Ic supernova 2006aj. In this paper, we present a simple model for computing the characteristic quantities (including energy, temperature and time duration) for the transient event from the shock breakout in Type Ibc supernovae produced by the core-collapse of WR stars surrounded by dense winds. In contrast to the case of a star without a strong wind, the shock breakout occurs in the wind region rather than inside the star, caused by the large optical depth in the wind. We find that, for the case of a WR star with a dense wind, the total energy of the radiation generated by the supernova shock breakout is larger than that in the case of the same star without a wind by a factor of >10. The temperature can be either hotter or colder, depending on the wind parameters. The time duration is larger caused by the increase in the effective radius of the star due to the presence of a thick wind. Then, we apply the model to GRB 060218/SN 2006aj. We show that, to explain both the temperature and the total energy of the blackbody component observed in GRB 060218 by the shock breakout, the progenitor WR star has to have an unrealistically large core radius (the radius at optical depth of 20), larger than 100 R_⊙. In spite of this disappointing result, our model is expected to have important applications to the observations on Type Ibc supernovae in which the detection of shock breakout will provide important clues to the progenitors of Type Ibc supernovae.     

Exospheric models for the X-ray emission from single Wolf–Rayet stars

We review existing ROSAT detections of single Galactic Wolf–Rayet (WR) stars and develop wind models to interpret the X-ray emission. The ROSAT data, consisting of bandpass detections from the ROSAT All-Sky Survey (RASS) and some pointed observations, exhibit no correlations of the WR X-ray luminosity ( L _X) with any star or wind parameters of interest (e.g. bolometric luminosity, mass-loss rate or wind kinetic energy), although the dispersion in the measurements is quite large. The lack of correlation between X-ray luminosity and wind parameters among the WR stars is unlike that of their progenitors, the O stars, which show trends with such parameters. In this paper we seek to (i) test by how much the X-ray properties of the WR stars differ from the O stars and (ii) place limits on the temperature T _X and filling factor f _X of the X-ray-emitting gas in the WR winds. Adopting empirically derived relationships for T _X and f _X from O-star winds, the predicted X-ray emission from WR stars is much smaller than observed with ROSAT . Abandoning the T _X relation from O stars, we maximize the cooling from a single-temperature hot gas to derive lower limits for the filling factors in WR winds. Although these filling factors are consistently found to be an order of magnitude greater than those for O stars, we find that the data are consistent (albeit the data are noisy) with a trend of in WR stars, as is also the case for O stars.     

The Galactic hybrid Wolf–Rayet WN 7o/CE + O7V((f)) binary system WR 145

A spectroscopic study of the binary Wolf–Rayet (WR)+O system WR 145 is performed, in order to determine the radial velocity orbits of the individual stars, the angle of orbital inclination and the stellar masses. The emission and absorption components are separated from the original spectra, allowing us to confirm the spectral classification WN 7o/CE of the hybrid WR component and to derive a spectral classification O7V((f)) for the O star. A study of the wind-collision properties is performed. Fitting the radial velocity and full width at half-maximum of the excess emission with Lührs' model results in an inclination angle of   i = 63°  , leading to estimates of the stellar masses:   M _WR= 18 M_⊙  and   M _O= 31 M_⊙  . Both of these masses are compatible with those of other stars of similar types.     

Non-emission-line young stars of intermediate mass

We present optical spectra of four intermediate-mass candidate young stellar objects that have often been classified as Herbig Ae/Be stars. Typical Herbig Ae/Be emission features are not present in the spectra of these stars. Three of them, HD 36917, HD 36982 and HD 37062, are members of the young Orion nebula cluster (ONC). This association constrains their ages to be ≲1 Myr. The lack of appreciable near-infrared excess in them suggests the absence of hot dust close to the central star. However, they do possess significant amounts of cold and extended dust as revealed by the large excess emission observed at far-infrared wavelengths. The fractional infrared luminosities  ( L _ir/ L _★)  and the dust masses computed from IRAS fluxes are systematically lower than those found for Herbig Ae/Be stars but higher than those for Vega-like stars. These stars may thus represent the youngest examples of the Vega phenomenon known so far. In contrast, the other star in our sample, HD 58647, is more likely to be a classical Be star, as is evident from the low   L _ir/ L _★  , the scarcity of circumstellar dust, the low polarization, the presence of H α emission and near-infrared excess, and the far-infrared spectral energy distribution consistent with free–free emission similar to other well-known classical Be stars.     

Single and binary evolution of Population III stars and their supernova explosions

We present stellar evolution calculations for Population III stars for both single- and binary-star evolutions. Our models include 10- and  16.5-M_⊙  single stars and a  10-M_⊙  model star that undergoes an episode of accretion resulting in a final mass of  16.1 M_⊙  . For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late-stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single-star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near-solar heavy metal abundance stars and not including both the evolution and the supernova explosions in one work.     

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.     

Optical spectroscopy of the quiescent counterpart to EXO 0748−676

We present phase resolved optical spectroscopy and X-ray timing of the neutron star X-ray binary EXO 0748−676 after the source returned to quiescence in the autumn of 2008. The X-ray light curve displays eclipses consistent in orbital period, orbital phase and duration with the predictions and measurements before the return to quiescence. Hα and He  i emission lines are present in the optical spectra and show the signature of the orbit of the binary companion, placing a lower limit on the radial velocity semi-amplitude of   K ₂ > 405 km s⁻¹  . Both the flux in the continuum and the emission lines show orbital modulations, indicating that we observe the hemisphere of the binary companion that is being irradiated by the neutron star. Effects due to this irradiation preclude a direct measurement of the radial velocity semi-amplitude of the binary companion; in fact, no stellar absorption lines are seen in the spectrum. Nevertheless, our observations place a stringent lower limit on the neutron star mass of   M ₁ > 1.27 M_⊙  . For the canonical neutron star mass of   M ₁= 1.4 M_⊙  , the mass ratio is constrained to  0.075 < q < 0.105  .     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Bulk viscosity of mixed nucleon–hyperon–quark matter in neutron stars

We calculate the coefficient of bulk viscosity by considering the non-leptonic weak interactions in the cores of hybrid stars with both hyperons and quarks. We first determine the dependence of the production rate of neutrons on the reaction rate of quarks in the non-leptonic processes, that is,  Γ _n = K _s Γ _s +Γ_Λ+ 2Γ_Σ−  . The conversion rate,   K _s   , in our scenario is a complicated function of baryon number density. We also consider the medium effect of quark matter on bulk viscosity. Using these results, we estimate the limiting rotation of the hybrid stars, which may suppress the r-mode instability more effectively. Hybrid stars should be the candidates for the extremely rapid rotators.     

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.     

Methanol masers and their environment at high resolution

We have used the Australia Telescope Compact Array (ATCA) to make high-resolution images of the 6.7-GHz 5₁ → 6₀A⁺ maser transition of methanol towards 33 sources in the Galactic plane. Including the results from 12 methanol sources in the literature, we find that 17 out of 45 sources have curved or linear morphology. Most of the 17 have a velocity gradient along the line, which is consistent with masers lying in an edge-on circumstellar disc surrounding a massive star. We also made simultaneous continuum observations of the sources at 8.6 GHz, in order to image any associated H  ii region. 25 of the sources are associated with an ultracompact H  ii region, with a detection limit of ∼0.5 mJy beam⁻¹. We argue that the methanol sources without an associated H  ii region represent less massive embedded stars, not an earlier stage in the lifetime of the star, as previously suggested.     

Photometric modelling of starspots – I. A Barnes–Evans-like surface brightness–colour relation using (Ic−K)

In the first part of this work, the empirical correlation of stellar surface brightness F _V with ( I _c− K ) broad-band colour is investigated by using a sample of stars cooler than the Sun. A bilinear correlation is found to represent well the brightness of G, K and M giant stars. The change in slope occurs at ( I _c− K )∼2.1 or at about the transition from K to M spectral types. The same relationship is also investigated for dwarf stars and found to be distinctly different from that of the giants. The dwarf star correlation differs by an average of −0.4 in ( I _c− K ) or by a maximum in F _V of ∼−0.1, positioning it below that of the giants, with both trends tending towards convergence for the hotter stars in our sample. The flux distribution derived from the F _V −( I _c− K ) relationship for the giant stars, together with that derived from an F _V −( V − K ) relationship and the blackbody flux distribution, is then utilized to compute synthetic light V and colour ( V − R )_c, ( V − I )_c and ( V − K ) curves of cool spotted stars. We investigate the effects on the amplitudes of the curves by using these F _V –colour relations and by assuming the effective gravity of the spots to be lower than the gravity of the unspotted photosphere. We find that the amplitudes produced by using the F _V −( I _c− K ) relationship are larger than those produced by the other two brightness correlations, meaning smaller and/or warmer spots.     

On the magnetic structure and wind parameter profiles of Alfvén wave driven winds in late-type supergiant stars

Cool stars at giant and supergiant evolutionary phases present low-velocity and high-density winds, responsible for the observed high mass-loss rates. Although presenting high luminosities, radiation pressure on dust particles is not sufficient to explain the wind acceleration process. Among the possible solutions to this still unsolved problem, Alfvén waves are, probably, the most interesting for their high efficiency in transfering energy and momentum to the wind. Typically, models of Alfvén wave driven winds result in high-velocity winds if they are not highly damped. In this work, we determine self-consistently the magnetic field geometry and solve the momentum, energy and mass conservation equations, to demonstrate that even a low-damped Alfvén wave flux is able to reproduce the low-velocity wind. We show that the magnetic flux tubes expand with a super-radial factor of S > 30 near the stellar surface, larger than that used in previous semi-empirical models. The rapid expansion results in a strong spatial dilution of the wave flux. We obtained the wind parameter profiles for a typical supergiant star of  16 M_⊙  . The wind is accelerated in a narrow region, coincident with the region of high divergence of the magnetic field lines, up to 100 km s⁻¹. For the temperature, we obtained a slight decrease near the surface for low-damped waves, because the wave heating mechanism is less effective than the radiative losses. The peak temperature occurs at   r ≃ 1.5  r ₀  reaching 6000 K. Propagating outwards, the wind cools down mainly due to adiabatic expansion.     

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.     

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.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

An asteroseismic signature of helium ionization

We investigate the influence of the ionization of helium on the low-degree acoustic oscillation frequencies in model solar-type stars. The signature in the oscillation frequencies characterizing the ionization-induced depression of the first adiabatic exponent γ is a superposition of two decaying periodic functions of frequency ν, with 'frequencies' that are approximately twice the acoustic depths of the centres of the He  i and He  ii ionization regions. That variation is probably best exhibited in the second frequency difference  Δ₂ν _{n ,  l} ≡ν _{n −1,  l} − 2ν _{n ,  l} +ν _{n +1,  l}   . We show how an analytic approximation to the variation of γ leads to a simple representation of this oscillatory contribution to Δ₂ν which can be used to characterize the γ variation, our intention being to use it as a seismic diagnostic of the helium abundance of the star. We emphasize that the objective is to characterize γ, not merely to find a formula for Δ₂ν that reproduces the data.     

On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

Doppler images and chromospheric variability of TWA 17

We present Doppler imaging and a Balmer line analysis of the weak-line T Tauri star TWA 17. Spectra were taken in 2006 with the University College London Echelle Spectrograph on the Anglo-Australian Telescope. Using least-squares deconvolution to improve the effective signal-to-noise ratio, we produced a Doppler map of the surface spot distribution. This shows similar features to maps of other rapidly rotating T Tauri stars, i.e. a polar spot with more spots extending out of it down to the equator.
In addition to the photospheric variability, the chromospheric variability was studied using the Balmer emission. The mean Hα profile has a narrow component consistent with rotational broadening and a broad component extending out to 220 km s⁻¹. The variability in Hα suggests that the chromosphere has at least one slingshot prominence  3 R _*  above the surface.