Isolated versus common envelope dynamos in planetary nebula progenitors

The origin, evolution and role of magnetic fields in the production and shaping of proto-planetary nebulae (PPNe) and planetary nebulae (PNe) are a subject of active research. Most PNe and PPNe are axisymmetric with many exhibiting highly collimated outflows; however, it is important to understand whether such structures can be generated by isolated stars or require the presence of a binary companion. Towards this end, we study a dynamical, large-scale α−Ω interface dynamo operating in a 3.0 M_⊙ Asymptotic Giant Branch (AGB) star in both an isolated setting and a setting in which a low-mass companion is embedded inside the envelope. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For the isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. Furthermore, we investigate the energy requirements that convection must satisfy to accomplish this by analogy to the Sun. For the common envelope case, a robust dynamo results, unbinding the envelope under a range of conditions. Two qualitatively different types of explosion may arise: (i) magnetically induced, possibly resulting in collimated bipolar outflows and (ii) thermally induced from turbulent dissipation, possibly resulting in quasi-spherical outflows. A range of models is presented for a variety of companion masses.     

Extrasolar planets and the rotation and axisymmetric mass-loss of evolved stars

I examine the implications of the recently found extrasolar planets on the planet-induced axisymmetric mass-loss model for the formation of elliptical planetary nebulae (PNe). This model attributes the low departure from spherical mass-loss of upper asymptotic giant branch (AGB) stars to envelope rotation which results from deposition of orbital angular momentum of the planets. Since about half of all PNe are elliptical, i.e., have low equatorial to polar density contrast, it was predicted that about 50 per cent of all Sun-like stars have Jupiter-like planets around them, i.e., a mass about equal to that of Jupiter, M _J, or more massive. In the light of the new findings that only 5 per cent of Sun-like stars have such planets, and a newly proposed mechanism for axisymmetric mass-loss, the cool magnetic spots model, I revise this prediction. I predict that indeed ∼50 per cent of PN progenitors do have close planets around them, but the planets can have much lower masses, as low as ∼0.01 M _J, in order to spin-up the envelopes of AGB stars efficiently. To support this claim, I follow the angular momentum evolution of single stars with main-sequence mass in the range of 1.3–2.4 M_⊙ , as they evolve to the post-AGB phase. I find that single stars rotate much too slowly to possess any significant non-spherical mass-loss as they reach the upper AGB. It seems, therefore, that planets, in some cases even Earth-like planets, are sufficient to spin-up the envelope of these AGB stars for them to form elliptical PNe. The prediction that on average several such planets orbit each star, as in the Solar system, still holds.     

Turbulent dynamo in asymptotic giant branch stars

Using recent results on the operation of turbulent dynamos, we show that a turbulent dynamo may amplify a large-scale magnetic field in the envelopes of asymptotic giant branch (AGB) stars. We propose that a slow rotation of the AGB envelope can fix the symmetry axis, leading to the formation of an axisymmetric magnetic field structure. Unlike solar-type αω dynamos, the rotation has only a small role in amplifying the toroidal component of the magnetic field; instead of an αω dynamo we propose an α ² ω . The magnetic field may reach a value of     , where B _e is the equipartition (between the turbulent and magnetic energy densities) magnetic field. The large-scale magnetic field is strong enough for the formation of magnetic cool spots on the AGB stellar surface. The spots may regulate dust formation, and hence the mass-loss rate, leading to axisymmetric mass loss and the formation of elliptical planetary nebulae (PNe). Despite its role in forming cool spots, the large-scale magnetic field is too weak to play a dynamic role and directly influence the wind from the AGB star, as required by some models. We discuss other possible problems in models where the magnetic field plays a dynamic role in shaping the AGB winds, and argue that they cannot explain the formation of non-spherical PNe.     

An explosive end to intermediate-mass zero-metallicity stars and early Universe nucleosynthesis

We use the Cambridge stellar evolution code stars to model the evolution of 5 and  7 M_⊙  zero-metallicity stars. With enhanced resolution at the hydrogen- and helium-burning shell in the asymptotic giant branch (AGB) phases, we are able to model the entire thermally pulsing AGB (TP-AGB) phase. The helium luminosities of the thermal pulses are significantly lower than in higher metallicity stars so there is no third dredge-up. The envelope is enriched in nitrogen by hot-bottom burning of carbon that was previously mixed in during second dredge-up. There is no s -process enrichment owing to the lack of third dredge-up. The thermal pulses grow weaker as the core mass increases and they eventually cease. From then on the star enters a quiescent burning phase which lasts until carbon ignites at the centre of the star when the CO core mass is  1.36 M_⊙  . With such a high degeneracy and a core mass so close to the Chandrasekhar mass, we expect these stars to explode as type 1.5 supernovae, very similar to type Ia supernovae but inside a hydrogen-rich envelope.     

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.     

The depletion of carbon by extra mixing in metal-poor giants

There is an apparent dichotomy between the metal-poor  ([Fe/H]≤−2)  yet carbon-normal giants and their carbon-rich counterparts. The former undergo significant depletion of carbon on the red giant branch after they have undergone first dredge-up, whereas the latter do not appear to experience significant depletion. We investigate this in the context that the extra mixing occurs via the thermohaline instability that arises due to the burning of  ³He  . We present the evolution of [C/Fe], [N/Fe] and  ¹²C/¹³C  for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a  1.5 M_⊙  AGB companion, one having received  0.01 M_⊙  of material and the other having received  0.1 M_⊙  . We find the behaviour of the carbon-normal metal-poor stars is well reproduced by this mechanism. In addition, our models also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich stars. This extra-mixing mechanism is able to reproduce the observed properties of both carbon-normal and carbon-rich stars.     

A new population of planetary nebulae discovered in the Large Magellanic Cloud – I. Preliminary sample

We report our initial discovery of 73 new planetary nebulae (PNe) in the Large Magellanic Cloud (LMC) following confirmatory 2dF spectroscopy on the Anglo-Australian Telescope. Preliminary candidate sources come from a 10 per cent sub-area of our new deep, high-resolution Hα map of the central 25 deg² of the LMC obtained with the UK Schmidt Telescope. The depth of the high-resolution map was extended to   R _equiv∼ 22  for  Hα (4.5 × 10⁻¹⁷ erg cm⁻² s⁻¹Å⁻¹)  by a process of multi-exposure median co-addition of a dozen 2-h Hα exposures. The resulting map is at least 1-mag deeper than the best wide-field narrow-band LMC images currently available. This depth, combined with our selection technique, has also led to the discovery of extended asymptotic giant branch (AGB) haloes around many new and previously known LMC PNe for the first time. Once complete, our new survey is expected to triple the LMC PN population and have significant implications for the LMC PN luminosity function, kinematics, abundance gradients, chemical evolution and, via study of the AGB haloes, the initial to final mass relation for low- to intermediate-mass stars.     

Magnetic field, dust and axisymmetrical mass loss on the asymptotic giant branch

I propose a mechanism for axisymmetrical mass loss on the asymptotic giant branch (AGB) that may account for the axially symmetric structure of elliptical planetary nebulae. The proposed model operates for slowly rotating AGB stars, having angular velocities in the range of 10⁻⁴ω _Kep  ω  10⁻² ω_Kep, where ω_Kep is the equatorial Keplerian angular velocity. Such angular velocities could be gained from a planet companion of mass  0.1  M _Jupiter, which deposits its orbital angular momentum to the envelope at late stages, or even from single stars that are fast rotators on the main sequence. The model assumes that dynamo magnetic activity results in the formation of cool spots, above which dust forms much more easily. The enhanced magnetic activity towards the equator results in a higher dust formation rate there, and hence higher mass-loss rate. As the star ascends the AGB, both the mass-loss rate and magnetic activity increase rapidly, and hence the mass loss becomes more asymmetrical, with higher mass-loss rate closer to the equatorial plane.     

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.     

Cooling curves and initial models for low-mass white dwarfs (<0.25 M⊙) with helium cores

We present a detailed calculation of the evolution of low-mass (<0.25 M_⊙) helium white dwarfs. These white dwarfs (the optical companions to binary millisecond pulsars) are formed via long-term, low-mass binary evolution. After detachment from the Roche lobe, the hot helium cores have a rather thick hydrogen layer with mass between 0.01 and 0.06 M_⊙. As a result of mixing between the core and outer envelope, the surface hydrogen content ( X _surf) is 0.5–0.35 , depending on the initial value of the heavy element Z and the initial secondary mass. We found that the majority of our computed models experience one or two hydrogen shell flashes. We found that the mass of the helium dwarf in which the hydrogen shell flash occurs depends on the chemical composition. The minimum helium white dwarf mass in which a hydrogen flash takes place is 0.213 M_⊙ ( Z =0.003), 0.198 M_⊙ ( Z =0.01), 0.192 M_⊙ ( Z =0.02) or 0.183 M_⊙ ( Z =0.03). The duration of the flashes (independent of chemical composition) is between a few ×10⁶ and a few ×10⁷ yr. In several flashes the white dwarf radius will increase so much that it forces the model to fill its Roche lobe again. Our calculations show that the cooling history of the helium white dwarf depends dramatically on the thickness of the hydrogen layer. We show that the transition from a cooling white dwarf with a temporarily stable hydrogen-burning shell to a cooling white dwarf in which almost all residual hydrogen is lost in a few thermal flashes (via Roche lobe overflow) occurs between 0.183 and 0.213 M_⊙ (depending on the heavy element value).     

Metal-rich carbon stars in the Sagittarius dwarf spheroidal galaxy

We present spectroscopic observations from the Spitzer Space Telescope of six carbon-rich asymptotic giant branch (AGB) stars in the Sagittarius dwarf spheroidal galaxy (Sgr dSph) and two foreground Galactic carbon stars. The band strengths of the observed C₂H₂ and SiC features are very similar to those observed in Galactic AGB stars. The metallicities are estimated from an empirical relation between the acetylene optical depth and the strength of the SiC feature. The metallicities are higher than those of the Large Magellanic Cloud, and close to Galactic values. While the high metallicity could imply an age of around 1 Gyr, for the dusty AGB stars, the pulsation periods suggest ages in excess of 2 or 3 Gyr. We fit the spectra of the observed stars using the dusty radiative transfer model and determine their dust mass-loss rates to be in the range  1.0–3.3 × 10⁻⁸ M_⊙ yr⁻¹  . The two Galactic foreground carbon-rich AGB stars are located at the far side of the solar circle, beyond the Galactic Centre. One of these two stars shows the strongest SiC feature in our present Local Group sample.     

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.     

Red giant stellar collisions in the Galactic Centre

We show that collisions with stellar-mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al. We show that the missing objects are low-mass giants and asymptotic giant branch stars in the range  1–3 M_⊙  . Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of smoothed particle hydrodynamic calculations, restricted three-body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main-sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the  10.5 < K < 12  band. However, depleting the brighter  ( K < 10.5)  objects out to the required radius would require a large population of black hole impactors which would in turn deplete the  10.5 < K < 12  giants in a region much larger than is observed. We conclude that collisions with stellar-mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar-mass black holes in the Galactic Centre.     

Full evolution of low-mass white dwarfs with helium and oxygen cores

We study the full evolution of low-mass white dwarfs with helium and oxygen cores. We revisit the age dichotomy observed in many white dwarf companions to millisecond pulsar on the basis of white dwarf configurations derived from binary evolution computations. We evolve 11 dwarf sequences for helium cores with final masses of 0.1604, 0.1869, 0.2026, 0.2495, 0.3056, 0.3333, 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . In addition, we compute the evolution of five sequences for oxygen cores with final masses of 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . A metallicity of   Z = 0.02  is assumed. Gravitational settling, chemical and thermal diffusion are accounted for during the white dwarf regime. Our study reinforces the result that diffusion processes are a key ingredient in explaining the observed age and envelope dichotomy in low-mass helium-core white dwarfs, a conclusion we arrived at earlier on the basis of a simplified treatment for the binary evolution of progenitor stars. We determine the mass threshold where the age dichotomy occurs. For the oxygen white dwarf sequences, we report the occurrence of diffusion-induced, hydrogen-shell flashes, which, as in the case of their helium counterparts, strongly influence the late stages of white dwarf cooling. Finally, we present our results as a set of white dwarf mass–radius relations for helium and oxygen cores.     

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.     

Merged binary white dwarf evolution: rapidly accreting carbon–oxygen white dwarfs and the progeny of extreme helium stars

We have examined the evolution of merged low-mass double white dwarfs that become luminous helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a carbon–oxygen (CO) white dwarf. After a certain mass is accumulated, a helium shell flash occurs, the radius and luminosity increase and the star becomes a yellow giant. Mass accretion is stopped artificially when the total mass reaches a pre-determined value. When the mass above the helium-burning shell becomes small enough, the star evolves blueward almost horizontally in the Hertzsprung–Russell diagram. The theoretical models for the merger of a 0.6-M_⊙ CO white dwarf with a 0.3-M_⊙ He white dwarf agree very well with the observed locations of extreme helium stars in the  log  T _eff–log  g   diagram, with their observed rates of blueward evolution, and with luminosities and masses obtained from their pulsations. Together with predicted merger rates for  CO+He  white dwarf pairs, the evolutionary time-scales are roughly consistent with the observed numbers of extreme helium stars. Predicted surface carbon and oxygen abundances can be consistent with the observed values if carbon and oxygen produced in the helium shell during a previous asymptotic giant branch phase are assumed to exist in the helium zone of the initial CO white dwarfs. These results establish the  CO+He  white dwarf merger as the best, if not only, viable model for the creation of extreme helium stars and, by association, the majority of R Coronae Borealis stars.     

The determination of progenitor masses for planetary nebulae

We have evaluated the likely progenitor masses M _PG of nebulae having elliptical, circular and bipolar morphologies, using observed ratios between the populations of these sources, and deduced central star mass functions. We find that most bipolar nebulae (BPNe) are likely to arise from progenitors having mass M _PG>2.3 M_⊙ and spectral types earlier than A3.2, whilst circular sources are associated with progenitors of mass 1.0 M_⊙< M _PG<1.2 M_⊙ and spectral range G1.9–F7.8 . Elliptical sources arise from intermediate-mass progenitors. The procedures employed to determine these values are relatively insensitive to uncertainties in scaleheights and population ratios, and completely insensitive to uncertainties in the distance scale. They are, however, dependent upon the precise forms adopted for the initial–final and central star mass functions, and we discuss the sensitivity of M _PG to uncertainties in these functions.     

The case and fate of HD 75767 – neutron star or supernova?

We report the discovery of the nearby  ( d = 24 pc)  HD 75767 as an eight billion year old quadruple system consisting of a distant M dwarf pair, HD 75767 C–D, in orbit around the known short-period   P = 10.25 d  single-lined binary HD 75767 A–B, the primary of which is a solar-like G star. On the reasonable assumption of synchronous orbital rotation as well as rotational and orbital coplanarity for the inner pair, we get   M _B= 0.96 M_⊙  for the unseen HD 75767 B, that is, the case of a massive white dwarf. Upon future evolution, mass transfer towards HD 75767 B will render the   M _A= 0.96 M_⊙  G-type primary, now a turnoff star, to become a helium white dwarf of   M _A∼ 0.33 M_⊙  . Depending on the mass accretion rate, accretion efficiency and composition of the massive white dwarf, this in turn may result in a collapse of HD 75767 B with the formation of a millisecond pulsar, i.e. the creation of a low-mass binary pulsar (LMBP), or, instead, a Type Ia supernova explosion and the complete disruption of HD 75767 B. Irrespective of which scenario applies, we point to the importance of the distant M dwarfs as the likely agents for the formation of the inner, short-period HD 75767 A–B pair, and hence a path that particularly avoids preceding phases of common envelope evolution.     

Continuum emission around AGB stars at 1.2 mm

It is generally acknowledged that the mass-loss of asymptotic giant branch (AGB) stars undergoes variations on different time-scales. We constructed models for the dust envelopes for a sample of AGB stars to assess whether mass-loss variations influence the spectral energy distribution. To constrain the variability, extra observations at millimetre wavelengths (1.2 mm) were acquired. From the analysis of the dust models, two indications for the presence of mass-loss variations can be found, being (1) a dust temperature at the inner boundary of the dust envelope that is far below the dust condensation temperature and (2) an altered density distribution with respect to  ρ( r ) ∝ r ⁻²  resulting from a constant mass-loss rate. For five out of the 18 studied sources a two-component model of the envelope is required, consisting of an inner region with a constant mass-loss rate and an outer region with a less steep density distribution. For one source an outer region with a steeper density distribution was found. Moreover, in a search for time variability in our data set at 1.2 mm, we found that WX Psc shows a large relative time variation of 34 per cent which might partially be caused by variable molecular line emission.     

The origin of the lead-rich stars in the Galactic halo: investigation of model parameters for the s-process

Several stars at the low-metallicity extreme of the Galactic halo show large spreads of lead and associated 'heavy' s-process elements ([Pb/hs]). Theoretically, an s-process pattern should be obtained from an AGB star with a fixed metallicity and initial mass. For the third dredge-up and the s-process model, several important properties depend primarily on the core mass of AGB stars. Zijlstra reported that the initial-to-final mass relation steepens at low metallicity, due to low mass-loss efficiency. This might affect the model parameters of the AGB stars, e.g. the overlap factor and the neutron irradiation time, in particular at low metallicity. The calculated results do indeed show that the overlap factor and the neutron irradiation time are significantly small at low metallicities, especially for  3.0 M_⊙ AGB  stars. The scatter of [Pb/hs] found in low metallicities can therefore be explained naturally when varying the initial mass of the low-mass AGB stars.