We explore the luminosity L of magnetized white dwarfs and its effect on the mass-radius relation. We self-consistently obtain the interface between the electron degenerate gas dominated inner core and the outer ideal gas surface layer or envelope by incorporating both the components of gas throughout the model white dwarf. This is obtained by solving the set of magnetostatic equilibrium, photon diffusion and mass conservation equations in the Newtonian framework, for different sets of luminosity and magnetic field. We appropriately use magnetic opacity, instead of Kramer's opacity, wherever required. We show that the Chandrasekhar-limit is retained, even at high luminosity upto about 10^{-2} solar luminosity but without magnetic field, if the temperature is set constant inside the interface. However there is an increased mass for large-radius white dwarfs, an effect of photon diffusion. Nevertheless, in the presence of strong magnetic fields, with central strength of about 10^{14} G, super-Chandrasekhar white dwarfs, with masses of about 1.9 solar mass, are obtained even when the temperature inside the interface is kept constant. Most interestingly, small-radius magnetic white dwarfs remain super-Chandrasekhar even if their luminosity decreases to as low as about 10^{-20} solar luminosity. However, their large-radius counterparts in the same mass-radius relation merge with Chandrasekhar's result at low L. Hence, we argue for the possibility of highly magnetized, low luminous super-Chandrasekhar mass white dwarfs which, owing to their faintness, can be practically hidden.

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion.

We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we use an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds.

Asteroseismological analysis of NY Vir suggests that at least the outer 55 per cent of the star (in radius) rotates as a solid body and is tidally synchronized to the orbit. Detailed calculation of tidal dissipation rates in NY Vir fails to account for this synchronization. Recent observations of He core burning stars suggest that the extent of the convective core may be substantially larger than that predicted with theoretical models. We conduct a parametric investigation of sdB models generated with the Cambridge STARS code to artificially extend the radial extent of the convective core. These models with extended cores still fail to account for the synchronization. Tidal synchronization may be achievable with a non-MLT treatment of convection.

Convection in the cores of massive stars becomes anisotropic when they rotate. This anisotropy leads to a misalignment of the thermal gradient and the thermal flux, which in turn results in baroclinicity and circulation currents in the upper radiative zone. We show that this induces a much stronger meridional flow in the radiative zone than previously thought. This drives significantly enhanced mixing, though this mixing does not necessarily reach the surface. The extra mixing takes on a similar form to convective overshooting, and is relatively insensitive to the rotation rate above a threshold, and may help explain the large overshoot distances inferred from observations. This has significant consequences for the evolution of these stars by enhancing core-envelope mixing.

In a series of recent papers, it has been proposed that high field magnetic white dwarfs are the result of close binary interaction and merging. Population synthesis calculations have shown that the origin of isolated highly magnetic white dwarfs is consistent with the stellar merging hypothesis. In this picture, the observed fields are caused by an alpha-Omega dynamo driven by differential rotation. The strongest fields arise when the differential rotation equals the critical break-up velocity and result from the merging of two stars (one of which has a degenerate core) during common envelope evolution or from the merging of two white dwarfs. We now synthesise a population of binary systems to investigate the hypothesis that the magnetic fields in the magnetic cataclysmic variables also originate during stellar interaction in the common envelope phase. Those systems that emerge from common envelope more tightly bound form the cataclysmic variables with the strongest magnetic fields. We vary the common envelope efficiency parameter and compare the results of our population syntheses with observations of magnetic cataclysmic variables. We find that common envelope interaction can explain the observed characteristics of these magnetic systems if the envelope ejection efficiency is low.

The origin of magnetic fields in isolated and binary white dwarfs has been investigated in a series of recent papers. One proposal is that magnetic fields are generated through an alpha-omega dynamo during common envelope evolution. Here we present population synthesis calculations showing that this hypothesis is supported by observations of magnetic binaries.

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic and buoyancy effects exactly within the formalism of linear growth theories with nonlinear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic and shear instabilities as well as the magnetorotational instability (MRI). It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits.

We argue that an increased temperature in star-forming clouds alters the stellar initial mass function to be more bottom-light than in the Milky Way. At redshifts $z \gtrsim 6$, heating from the cosmic microwave background radiation produces this effect in all galaxies, and it is also present at lower redshifts in galaxies with very high star formation rates (SFRs). A failure to account for it means that at present, photometric template fitting likely overestimates stellar masses and star formation rates for the highest-redshift and highest-SFR galaxies. In addition this may resolve several outstanding problems in the chemical evolution of galactic halos.

We have simulated a common envelope interaction of a 0.88-M$_{\odot}$, 90-R$_{\odot}$, red giant branch star and a 0.6-M$_{\odot}$, compact companion with the smoothed particle hydrodynamics code, Phantom, from the beginning of the Roche lobe overflow phase to the beginning of the self-regulated inspiral, using three different resolutions. The duration of the Roche lobe overflow phase is resolution dependent and would lengthen with increased resolution beyond the $\sim$20 years observed, while the inspiral phase and the post-common envelope separation are largely independent of resolution. Mass transfer rates through the Lagrangian points drive the orbital evolution during the Roche lobe overflow phase, as predicted analytically. The absolute mass transfer rate is resolution dependent, but always within an order of magnitude of the analytical value. Similarly, the gravitational drag in the simulations is close to the analytical approximation. This gives us confidence that simulations approximate reality. The $L_2$ and $L_3$ outflow observed during Roche lobe overflow remains bound, forming a circumbinary disk that is largely disrupted by the common envelope ejection. However, a longer phase of Roche lobe overflow and weaker common envelope ejection typical of a more stable binary may result in a surviving circumbinary disk. Finally, we examine the density distribution resulting from the interaction for simulations that include or omit the phase of Roche lobe overflow. We conclude that the degree of stability of the Roche lobe phase may modulate the shape of the subsequent planetary nebula, explaining the wide range of post-common envelope planetary nebula shapes observed.

We evolve stellar models to study the common envelope (CE) interaction of an early asymptotic giant branch star of initial mass $5\,\rm M_{\odot}$ with a companion star of mass ranging from $0.1$ to $2\,\rm M_{\odot}$. We model the CE as a fast stripping phase in which the primary experiences rapid mass loss and loses about 80 per cent of its mass. The post-CE remnant is then allowed to thermally readjust during a Roche-lobe overflow (RLOF) phase and the final binary system and its orbital period are investigated. We find that the post-CE RLOF phase is long enough to allow nuclear burning to proceed in the helium shell. By the end of this phase, the donor is stripped of both its hydrogen and helium and ends up as carbon-oxygen white dwarf of mass about $0.8\,\rm M_{\odot}$. We study the sensitivity of our results to initial conditions of different companion masses and orbital separations at which the stripping phase begins. We find that the companion mass affects the final binary separation and that helium-shell burning causes the star to refill its Roche lobe leading to post-CE RLOF. Our results show that double mass transfer in such a binary interaction is able to strip the helium and hydrogen layers from the donor star without the need for any special conditions or fine tuning of the binary parameters.

Over half of all observed hot subdwarf B (sdB) stars are found in binaries, and over half of these are found in close configurations with orbital periods of 10$ \,\rm{d}$ or less. In order to estimate the companion masses in these predominantly single-lined systems, tidal locking has frequently been assumed for sdB binaries with periods less than half a day. Observed non-synchronicity of a number of close sdB binaries challenges that assumption and hence provides an ideal testbed for tidal theory. We solve the second-order differential equations for detailed 1D stellar models of sdB stars to obtain the tidal dissipation strength and hence to estimate the tidal synchronization time-scale owing to Zahn's dynamical tide. The results indicate synchronization time-scales longer than the sdB lifetime in all observed cases. Further, we examine the roles of convective overshooting and convective dissipation in the core of sdB stars and find no theoretical framework in which tidally-induced synchronization should occur.

Stars with unusual elemental abundances offer clues about rare astrophysical events or nucleosynthetic pathways. Stars with significantly depleted magnesium and enhanced potassium ([Mg/Fe] < -0.5; [K/Fe] > 1) have to date only been found in the massive globular cluster NGC 2419 and, to a lesser extent, NGC 2808. The origin of this abundance signature remains unknown, as does the reason for its apparent exclusivity to these two globular clusters. Here we present 112 field stars, identified from 454,180 LAMOST giants, that show significantly enhanced [K/Fe] and possibly depleted [Mg/Fe] abundance ratios. Our sample spans a wide range of metallicities (-1.5 < [Fe/H] < 0.3), yet none show abundance ratios of [K/Fe] or [Mg/Fe] that are as extreme as those observed in NGC 2419. If confirmed, the identified sample of stars represents evidence that the nucleosynthetic process producing the anomalous abundances ratios of [K/Fe] and [Mg/Fe] probably occurs at a wide range of metallicities. This would suggest that pollution scenarios that are limited to early epochs (such as Population III supernovae) are an unlikely explanation, although they cannot be ruled out entirely. This sample is expected to help guide modelling attempts to explain the origin of the Mg-K abundance signature.

A dynamo mechanism driven by differential rotation when stars merge has been proposed to explain the presence of strong fields in certain classes of magnetic stars. In the case of the high field magnetic white dwarfs (HFMWDs), the site of the differential rotation has been variously thought to be the common envelope, the hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently accreted by a degenerate core. We have shown previously that the observed incidence of magnetism and the mass distribution in HFMWDs are consistent with the hypothesis that they are the result of merging binaries during common envelope evolution. Here we calculate the magnetic field strengths generated by common envelope interactions for synthetic populations using a simple prescription for the generation of fields and find that the observed magnetic field distribution is also consistent with the stellar merging hypothesis. We use the Kolmogorov-Smirnov test to study the correlation between the calculated and the observed field strengths and find that it is consistent for low envelope ejection efficiency. We also suggest that field generation by the plunging of a giant gaseous planet on to a white dwarf may explain why magnetism among cool white dwarfs (including DZ white dwarfs) is higher than among hot white dwarfs. In this picture a super Jupiter residing in the outer regions of the planetary system of the white dwarf is perturbed into a highly eccentric orbit by a close stellar encounter and is later accreted by the white dwarf.

Rotation and magnetism are increasingly recognized as important phenomena in stellar evolution. Surface magnetic fields from a few to $20{,}000\,$G have been observed and models have suggested that magnetohydrodynamic transport of angular momentum and chemical composition could explain the peculiar composition of some stars. Stellar remnants such as white dwarfs have been observed with fields from a few to more than $10^{9}\,$G. We investigate the origin of and the evolution, on thermal and nuclear rather than dynamical time-scales, of an averaged large-scale magnetic field throughout a star's life and its coupling to stellar rotation. Large-scale magnetic fields sustained until late stages of stellar evolution with conservation of magnetic flux could explain the very high fields observed in white dwarfs. We include these effects in the Cambridge stellar evolution code using three time-dependant advection-diffusion equations coupled to the structural and composition equations of stars to model the evolution of angular momentum and the two components of the magnetic field. We present the evolution in various cases for a $3\rm\,M_{\odot}$ star from the beginning to the late stages of its life. Our particular model assumes that turbulent motions, including convection, favour small-scale field at the expense of large-scale field. As a result the large-scale field concentrates in radiative zones of the star and so is exchanged between the core and the envelope of the star as it evolves. The field is sustained until the end of the asymptotic giant branch, when it concentrates in the degenerate core.

Diffusion of elements in the atmosphere and envelope of a star can drastically alter its surface composition, leading to extreme chemical peculiarities. We consider the case of hot subdwarfs, where surface helium abundances range from practically zero to almost 100 percent. Since hot subdwarfs can form via a number of different evolution channels, a key question concerns how the formation mechanism is connected to the present surface chemistry. A sequence of extreme horizontal branch star models was generated by producing post-common envelope stars from red giants. Evolution was computed with MESA from envelope ejection up to core-helium ignition. Surface abundances were calculated at the zero-age horizontal branch for models with and without diffusion. A number of simulations also included radiative levitation. The goal was to study surface chemistry during evolution from cool giant to hot subdwarf and determine when the characteristic subdwarf surface is established. Only stars leaving the giant branch close to core-helium ignition become hydrogen-rich subdwarfs at the zero-age horizontal branch. Diffusion, including radiative levitation, depletes the initial surface helium in all cases. All subdwarf models rapidly become more depleted than observations allow. Surface abundances of other elements follow observed trends in general, but not in detail. Additional physics is required.

Wickramasinghe et al. (2014) and Briggs et al. (2015) have proposed that the strong magnetic fields observed in some single white dwarfs (MWDs) are formed by a dynamo driven by differential rotation when two stars, the more massive one with a degenerate core, merge during common envelope (CE) evolution (Ferrario et al., 2015b). We synthesize a population of binaries to investigate if fields in the magnetic cataclysmic variables (MCVs) may also originate during stellar interaction in the CE phase.

The most frequently seen binary companions to millisecond pulsars (MSPs) are helium white dwarfs (He WDs). The standard rejuvenation mechanism, in which a low- to intermediate-mass companion to a neutron star fills its Roche lobe between central hydrogen exhaustion and core helium ignition, is the most plausible formation mechanism. We have investigated whether the observed population can realistically be formed via this mechanism. We used the Cambridge STARS code to make models of Case B RLOF with Reimers' mass loss from the donor. We find that the range of initial orbital periods required to produce the currently observed range of orbital periods of MSPs is extremely narrow. To reduce this fine tuning, we introduce a companion enhanced wind (CEW) that strips the donor of its envelope more quickly so that systems can detach at shorter periods. Our models indicate that the fine tuning can be significantly reduced if a CEW is active. Because significant mass is lost owing to a CEW we expect some binary pulsars to accrete less than the 0.1Msun needed to spin them up to millisecond periods. This can account for mildly recycled pulsars present along the entire Mc-Porb relation. Systems with Pspin > 30ms are consistent with this but too few of these mildly recycled pulsars have yet been observed to make a significant comparison.

The combination of asteroseismologically-measured masses with abundances from detailed analyses of stellar atmospheres challenges our fundamental knowledge of stars and our ability to model them. Ancient red-giant stars in the Galactic thick disc are proving to be most troublesome in this regard. They are older than 5 Gyr, a lifetime corresponding to an initial stellar mass of about $1.2{\mathrm{M}_{\odot}}$. So why do the masses of a sizeable fraction of thick-disc stars exceed $1.3{\mathrm{M}_{\odot}}$, with some as massive as $2.3{\mathrm{M}_{\odot}}$ ? We answer this question by considering duplicity in the thick-disc stellar population using a binary population-nucleosynthesis model. We examine how mass transfer and merging affect the stellar mass distribution and surface abundances of carbon and nitrogen. We show that a few per cent of thick-disc stars can interact in binary star systems and become more massive than $1.3{\mathrm{M}_{\odot}}$. Of these stars, most are single because they are merged binaries. Some stars more massive than $1.3{\mathrm{M}_{\odot}}$ form in binaries by wind mass transfer. We compare our results to a sample of the APOKASC data set and find reasonable agreement except in the number of these thick-disc stars more massive than $1.3{\mathrm{M}_{\odot}}$. This problem is resolved by the use of a logarithmically-flat orbital-period distribution and a large binary fraction.

We estimate the extent of overshooting inwards from the bottom of the intershell convective zone in thermal pulses in (S)AGB stars. We find that the buoyancy is so strong that any overshooting should be negligible. The temperature inversion at the bottom of the convective zone adds to the stability of the region. Any mixing that occurs in this region is highly unlikely to be due to convective overshooting, and so must be due to another process.

We study the interaction between stellar irradiation and tidal heating in gaseous planets with short orbital periods. The intentionally simplified atmospheric model we employ makes the problem analytically tractable and permits the derivation of useful scaling relations. We show that many tidal models provide thermal feedback, producing interior radiative zones and leading to enhanced g-mode dissipation with a wide spectrum of resonances. These resonances are dynamically tuned by the thermal feedback, and so represent a novel form of thermomechanical feedback, coupling vibrational modes to the very slow thermal evolution of the planet. We then show that stellar irradiation allows the heat produced by these modes to be trapped at depth with high efficiency, leading to entropy increase in the central convective region, as well as expansion of the planet's radius sufficient to match observed swelling. We find that thermally driven winds play an essential role in this process by making the thermal structure of the atmosphere spherically symmetric within a few scale heights of the photosphere. We characterise the relationship between the swelling factor, the orbital period and the host star and determine the timescale for swelling. We show that these g-modes suffice to produce bloating on the order of the radius of the planet over $\mathrm{Gyr}$ timescales when combined with significant insolation and we provide analytic relations for the relative magnitudes of tidal heating and insolation.

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

The candidate Thorne-\.Zytkow object (T\.ZO), HV2112, is becoming a well-studied if enigmatic object. A key point of its candidacy as a T\.ZO is whether or not it resides in the Small Magellanic Cloud (SMC). HV2112 has detections in a series of photometric catalogues which have resulted in contradictory estimates of its proper motion and, therefore, its membership within the SMC. This letter seeks to resolve the issue of the SMC membership of HV2112 through a reanalysis of extant photometric data. We also demonstrate the difficulties and downfalls inherent in considering a range of catalogue proper motions. We conclude that the proper motion, and associated ancillary radial velocity, positional and photometric properties, are fully consistent with HV2112 being within the SMC and thus it remains a candidate T\.ZO.

We present the Kepler photometric light-variation analysis of the late-type double-lined binary system V568 Lyr that is in the field of the high metallicity old open cluster NGC 6791. The radial velocity and the high-quality short-cadence light curve of the system are analysed simultaneously. The masses, radii and luminosities of the component stars are $M_1 = 1.0886\pm0.0031\, M{\odot}$, $M_2 = 0.8292 \pm 0.0026\, M{\odot}$, $R_1 = 1.4203\pm 0.0058\, R{\odot}$, $R_2 = 0.7997 \pm 0.0015\, R{\odot}$, $L_1 = 1.85\pm 0.15\, L{\odot}$, $L_2 = 0.292 \pm 0.018\, L{\odot}$ and their separation is $a = 31.060 \pm 0.002\, R{\odot}$. The distance to NGC 6791 is determined to be $4.260\pm 0.290\,$kpc by analysis of this binary system. We fit the components of this well-detached binary system with evolution models made with the Cambridge STARS and TWIN codes to test low-mass binary star evolution. We find a good fit with a metallicity of $Z = 0.04$ and an age of $7.704\,$Gyr. The standard tidal dissipation, included in TWIN is insufficient to arrive at the observed circular orbit unless it formed rather circular to begin with.

We consider the blue loops in the Hertzsprung-Russell diagram that occur when intermediate-mass stars begin core helium burning. It has long been known that the excess of helium above the burning shell, the result of the contraction of the convective core during core hydrogen burning, has the effect of making such stars redder and larger than they would be otherwise. The outward motion of the burning shell in mass removes this excess and triggers the loop. Hitherto nobody has attempted to demonstrate why the excess helium has this effect. We consider the effect of the local opacity, which is reduced by excess helium, the shell fuel supply, which is also reduced, and the local mean molecular weight, which is increased. We demonstrate that the mean molecular weight is the decisive reddening factor. The opacity has a much smaller effect and a reduced fuel supply actually favours blueward motion.

A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently incorporated into a degenerate core. We synthesize a population of binary systems to investigate the stellar merging hypothesis for observed single HFMWDs. Our calculations provide mass distribution and the fractions of white dwarfs that merge during a common envelope phase or as double degenerate systems in a post common envelope phase. We vary the common envelope efficiency parameter alpha and compare with observations. We find that this hypothesis can explain both the observed incidence of magnetism and the mass distribution of HFMWDs for a wide range of alpha. In this model, the majority of the HFMWDs are of the Carbon Oxygen type and merge within a common envelope. Less than about a quarter of a per cent of HFMWDs originate from double degenerate stars that merge after common envelope evolution and these populate the high-mass tail of the HFMWD mass distribution.

We examine the growing class of binary millisecond pulsars known as redbacks. In these systems the pulsar's companion has a mass between 0.1 and about 0.5 solar masses in an orbital period of less than 1.5 days. All show extended radio eclipses associated with circumbinary material. They do not lie on the period-companion mass relation expected from the canonical intermediate-mass X-ray binary evolution in which the companion filled its Roche lobe as a red giant and has now lost its envelope and cooled as a white dwarf. The redbacks lie closer to, but usually at higher period than, the period-companion mass relation followed by cataclysmic variables and low-mass X-ray binaries. In order to turn on as a pulsar mass accretion on to a neutron star must be sufficiently weak, considerably weaker than expected in systems with low-mass main-sequence companions driven together by magnetic braking or gravitational radiation. If a neutron star is formed by accretion induced collapse of a white dwarf as it approaches the Chandrasekhar limit some baryonic mass is abruptly lost to its binding energy so that its effective gravitational mass falls. We propose that redbacks form when accretion induced collapse of a white dwarf takes place during cataclysmic variable binary evolution because the loss of gravitational mass makes the orbit expand suddenly so that the companion no longer fills its Roche lobe. Once activated, the pulsar can ablate its companion and so further expand the orbit and also account for the extended eclipses in the radio emission of the pulsar that are characteristic of these systems. The whole period-companion mass space occupied by the redbacks can be populated in this way.

Many classes of objects and events are thought to form in binary star systems after a phase in which a core and companion spiral to smaller separation inside a common envelope (CE).Such a phase can end with the merging of the two stars or with the ejection of the envelope to leave a surviving binary system.The outcome is usually predicted by calculating the separation to which the stars must spiral to eject the envelope, assuming that the ratio of the core--envelope binding energy to the change in orbital energy is equal to a constant efficiency factor $\alpha$. If either object would overfill its Roche lobe at this end-of-CE separation, then the stars are assumed to merge. It is unclear what critical radius should be compared to the end-of-CE Roche lobe for stars which have developed cores before the start of a CE phase. After improving the core radius formulae in the widely used BSE rapid evolution code, we compare the properties of populations in which the critical radius is chosen to be the pre-CE core radius or the post-CE stripped remnant radius. Our improvements to the core radius formulae and the uncertainty in the critical radius significantly affect the rates of merging in CE phases of most types. We find the types of systems for which these changes are most important.

The very bright red star HV2112 in the Small Magellanic Cloud could be a massive Thorne-Zytkow Object, a supergiant-like star with a degenerate neutron core. With its luminosity of over $10^5\,\rm L_\odot$, it could also be a super asymptotic giant branch star, a star with an oxygen/neon core supported by electron degeneracy and undergoing thermal pulses with third dredge up. Both TZOs and SAGB stars are expected to be rare. Abundances of heavy elements in HV2112's atmosphere, as observed to date, do not allow us to distinguish between the two possibilities based on the latest models. Molybdenum and rubidium can be enhanced by both the irp-process in a TZO or by the s-process in SAGB stars. Lithium can be generated by hot bottom burning at the base of the convective envelope in either. HV2112's enhanced calcium could thus be the key determinant. A SAGB star is not able to synthesise its own calcium but it may be possible to produce this in the final stages of the process that forms a TZO, when the degenerate electron core of a giant star is tidally disrupted by a neutron star. Hence our calculations indicate that HV2112 is most likely a genuine TZO.

We provide a set of analytic fits to the radii of pre-mainsequence stars in the mass range 0.1 < M/Msun < 8.0. We incorporate the formulae in N-body cluster models for evolution from the beginning of pre-main sequence. In models with 1,000 stars and high initial cluster densities, pre-mainsequence evolution causes roughly twice the number of collisions between stars than in similar models with evolution begun only from the zero-age main sequence. The collisions are often all part of a runaway sequence that creates one relatively massive star.

We examine the growing data set of binary millisecond pulsars that are thought to have a helium white dwarf companion. These systems are believed to form when a low- to intermediate-mass companion to a neutron star fills its Roche lobe between central hydrogen exhaustion and core helium ignition. We confirm that our own stellar models reproduce a well-defined period-companion mass relation irrespective of the details of the mass transfer process. With magnetic braking this relation extends to periods of less than 1d for a 1Msun giant donor. With this and the measured binary mass functions we calculate the orbital inclination of each system for a given pulsar mass. We expect these inclinations to be randomly oriented in space. If the masses of the pulsars were typically 1.35Msun then there would appear to be a distinct dearth of high-inclination systems. However if the pulsar masses are more typically 1.55 to 1.65Msun then the distribution of inclinations is indeed indistinguishable from random. If it were as much as 1.75Msun then there would appear to be an excess of high-inclination systems. Thus with the available data we can argue that the neutron star masses in binary millisecond pulsars recycled by mass transfer from a red giant typically lie around 1.6Msun and that there is no preferred inclination at which these systems are observed. Hence there is reason to believe that pulsar beams are either sufficiently broad or show no preferred direction relative to the pulsar's spin axis which is aligned with the binary orbit. This is contrary to some previous claims, based on a subset of the data available today, that there might be a tendency for the pulsar beams to be perpendicular to their spin.

It is likely that at least some planetary nebulae are composed of matter which was ejected from a binary star system during common-envelope (CE) evolution. For these planetary nebulae the ionizing component is the hot and luminous remnant of a giant which had its envelope ejected by a companion in the process of spiralling-in to its current short-period orbit. A large fraction of CE phases which end with ejection of the envelope are thought to be initiated by low-mass red giants, giants with inert, degenerate helium cores. We discuss the possible end-of-CE structures of such stars and their subsequent evolution to investigate for which structures planetary nebulae are formed. We assume that a planetary nebula forms if the remnant reaches an effective temperature greater than 30 kK within 10^4 yr of ejecting its envelope. We assume that the composition profile is unchanged during the CE phase so that possible remnant structures are parametrized by the end-of-CE core mass, envelope mass and entropy profile. We find that planetary nebulae are expected in post-CE systems with core masses greater than about 0.3 solar masses if remnants end the CE phase in thermal equilibrium. We show that whether the remnant undergoes a pre-white dwarf plateau phase depends on the prescribed end-of-CE envelope mass. Thus, observing a young post-CE system would constrain the end-of CE envelope mass and post-CE evolution.

Observations of magnetic A, B and O stars show that the poloidal magnetic flux per unit mass has an upper bound of 10^-6.5 G cm^2/g. A similar upper bound is found for magnetic white dwarfs even though the highest magnetic field strengths at their surfaces are much larger. For magnetic A and B stars there also appears to be a well defined lower bound below which the incidence of magnetism declines rapidly. According to recent hypotheses, both groups of stars may result from merging stars and owe their strong magnetism to fields generated by a dynamo mechanism as they merge. We postulate a simple dynamo that generates magnetic field from differential rotation. The growth of magnetic fields is limited by the requirement that the poloidal field stabilizes the toroidal and vice versa. While magnetic torques dissipate the differential rotation, toroidal field is generated from poloidal by an Omega dynamo. We further suppose that mechanisms that lead to the decay of toroidal field lead to the generation of poloidal. Both poloidal and toroidal fields reach a stable configuration which is independent of the size of small initial seed fields but proportional to the initial differential rotation. We pose the hypothesis that strongly magnetic stars form from the merging of two stellar objects. The highest fields are generated when the merge introduces differential rotation that amounts to critical break up velocity within the condensed object. Calibration of a simplistic dynamo model with the observed maximum flux per unit mass for main-sequence stars and white dwarfs indicates that about 1.5x10^-4 of the decaying toroidal flux must appear as poloidal. The highest fields in single white dwarfs are generated when two degenerate cores merge inside a common envelope or when two white dwarfs merge by gravitational-radiation angular momentum loss.

Our improving understanding of the life cycle of planetary systems prompts investigations of the role of the Galactic environment before, during and after Asymptotic Giant Branch (AGB) stellar evolution. Here, we investigate the interplay between stellar mass loss, Galactic tidal perturbations, and stellar flybys for evolving stars which host one planet, smaller body or stellar binary companion and reside in the Milky Way's bulge or disc. We find that the potential evolutionary pathways from a main sequence (MS) to a white dwarf (WD) planetary system are a strong function of Galactocentric distance only with respect to the prevalence of stellar flybys. Planetary ejection and collision with the parent star should be more common towards the bulge. At a given location anywhere in the Galaxy, if the mass loss is adiabatic, then the secondary is likely to avoid close flybys during AGB evolution, and cannot eventually escape the resulting WD because of Galactic tides alone. Partly because AGB mass loss will shrink a planetary system's Hill ellipsoid axes by about 20 to 40 per cent, Oort clouds orbiting WDs are likely to be more depleted and dynamically excited than on the MS.

We propose a new method to infer the star formation histories of resolved stellar populations. With photometry one may plot observed stars on a colour-magnitude diagram (CMD) and then compare with synthetic CMDs representing different star formation histories. This has been accomplished hitherto by parametrising the model star formation history as a histogram, usually with the bin widths set by fixed increases in the logarithm of time. A best fit is then found with maximum likelihood methods and we consider the different means by which a likelihood can be calculated. We then apply Bayesian methods by parametrising the star formation history as an unknown number of Gaussian bursts with unknown parameters. This parametrisation automatically provides a smooth function of time. A Reversal Jump Markov Chain Monte Carlo method is then used to find both the most appropriate number of Gaussians, thus avoiding avoid overfitting, and the posterior probability distribution of the star formation rate. We apply our method to artificial populations and to observed data. We discuss the other advantages of the method: direct comparison of different parametrisations and the ability to calculate the probability that a given star is from a given Gaussian. This allows the investigation of possible sub-populations.

The birth and death of planets may be affected by mass outflows from their parent stars during the T-Tauri or post-main-sequence phases of stellar evolution. These outflows are often modelled to be isotropic, but this assumption is not realistic for fast rotators, bipolar jets and supernovae. Here we derive the general equations of motion for the time evolution of a single planet, brown dwarf, comet or asteroid perturbed by anisotropic mass loss in terms of a complete set of planetary orbital elements, the ejecta velocity, and the parent star's co-latitude and longitude. We restrict our application of these equations to 1) rapidly rotating giant stars, and 2) arbitrarily-directed jet outflows. We conclude that the isotropic mass-loss assumption can safely be used to model planetary motion during giant branch phases of stellar evolution within distances of hundreds of au. In fact, latitudinal mass loss variations anisotropically affect planetary motion only if the mass loss is asymmetric about the stellar equator. Also, we demonstrate how constant-velocity, asymmetric bipolar outflows in young systems incite orbital inclination changes. Consequently, this phenomenon readily tilts exoplanetary orbits external to a nascent disc on the order of degrees.

The helioseismic observations of the internal rotation profile of the Sun raise questions about the two-dimensional (2D) nature of the transport of angular momentum in stars. Here we derive a convective prescription for axisymmetric (2D) stellar evolution models. We describe the small scale motions by a spectrum of unstable linear modes in a Boussinesq fluid. Our saturation prescription makes use of the angular dependence of the linear dispersion relation to estimate the anisotropy of convective velocities. We are then able to provide closed form expressions for the thermal and angular momentum fluxes with only one free parameter, the mixing length. We illustrate our prescription for slow rotation, to first order in the rotation rate. In this limit, the thermodynamical variables are spherically symetric, while the angular momentum depends both on radius and latitude. We obtain a closed set of equations for stellar evolution, with a self-consistent description for the transport of angular momentum in convective regions. We derive the linear coefficients which link the angular momentum flux to the rotation rate ($\Lambda$- effect) and its gradient ($\alpha$-effect). We compare our results to former relevant numerical work.

Models of rotationally-driven dynamos in stellar radiative zones have suggested that magnetohydrodynamic transport of angular momentum and chemical composition can dominate over the otherwise purely hydrodynamic processes. A proper consideration of the interaction between rotation and magnetic fields is therefore essential. Previous studies have focused on a magnetic model where the magnetic field strength is derived as a function of the stellar structure and angular momentum distribution. We have adapted our one-dimensional stellar rotation code, RoSE, to model the poloidal and toroidal magnetic field strengths with a pair of time-dependent advection-diffusion equations coupled to the equations for the evolution of the angular momentum distribution and stellar structure. This produces a much more complete, though still reasonably simple, model for the magnetic field evolution. Our model reproduces well observed surface nitrogen enrichment of massive stars in the Large Magellanic Cloud. In particular it reproduces a population of slowly-rotating nitrogen-enriched stars that cannot be explained by rotational mixing alone alongside the traditional rotationlly-enriched stars. The model further predicts a strong mass-dependency for the dynamo-driven field. Above a threshold mass, the strength of the magnetic dynamo decreases abruptly and so we predict that more massive stars are much less likely to support a dynamo-driven field than less massive stars.

Rotation has a number of important effects on the evolution of stars. Apart from structural changes because of the centrifugal force, turbulent mixing and meridional circulation caused by rotation can dramatically affect a star's chemical evolution. This leads to changes in the surface temperature and luminosity as well as modifying its lifetime. Observationally rotation decreases the surface gravity, causes enhanced mass loss and leads to surface abundance anomalies of various chemical isotopes. The replication of these physical effects with simple stellar evolution models is very difficult and has resulted in the use of numerous different formulations to describe the physics. Using stellar evolution calculations based on several physical models we discuss the features of the resulting simulated stellar populations which can help to distinguish between the models.

The existence of bright quasars at high redshifts implies that supermassive black holes were able to form in the early Universe. Though a number of mechanisms to achieve this have been proposed, none yet stands out. A recent suggestion is the formation of quasi-stars, initially stellar-mass black holes accreting from hydrostatic giant-like envelopes of gas, formed from the monolithic collapse of pre-galactic gas clouds. In this work, we modify the Cambridge STARS stellar evolution package to construct detailed models of the evolution of these objects. We find that, in all of our models, the black hole inside the envelope is able to reach slightly more than one-tenth of the total mass of the system before hydrostatic equilibrium breaks down. This breakdown occurs after a few million years of evolution. We show that the mechanism which causes the hydrostatic evolution to end is present in polytropic models. We also show that the solutions are highly sensitive to the size of the inner boundary radius and that no physical solutions exist if the inner boundary is chosen to be less than about 0.3 of the Bondi radius.

Extrasolar planets and belts of debris orbiting post-main-sequence single stars may become unbound as the evolving star loses mass. In multiple star systems, the presence or co-evolution of the additional stars can significantly complicate the prospects for orbital excitation and escape. Here, we investigate the dynamical consequences of multi-phasic, nonlinear mass loss and establish a criterion for a system of any stellar multiplicity to retain a planet whose orbit surrounds all of the parent stars. For single stars which become white dwarfs, this criterion can be combined with the Chandrasekhar Limit to establish the maximum allowable mass loss rate for planet retention. We then apply the criterion to circumbinary planets in evolving binary systems over the entire stellar mass phase space. Through about 10^5 stellar evolutionary track realizations, we characterize planetary ejection prospects as a function of binary separation, stellar mass and metallicity. This investigation reveals that planets residing at just a few tens of AU from a central concentration of stars are susceptible to escape in a wide variety of multiple systems. Further, planets are significantly more susceptible to ejection from multiple star systems than from single star systems for a given system mass. For system masses greater than about 2 Solar masses, multiple star systems represent the greater source of free-floating planets.

The Sch\"onberg-Chandrasekhar (SC) limit is a well-established result in the understanding of stellar evolution. It provides an estimate of the point at which an evolved isothermal core embedded in an extended envelope begins to contract. We investigate contours of constant fractional mass in terms of homology invariant variables U and V and find that the SC limit exists because the isothermal core solution does not intersect all the contours for an envelope with polytropic index 3. We find that this analysis also applies to similar limits in the literature including the inner mass limit for polytropic models of quasi-stars. Consequently, any core solution that does not intersect all the fractional mass contours exhibits an associated limit and we identify several relevant cases where this is so. We show that a composite polytrope is at a fractional core mass limit when its core solution touches but does not cross the contour of the corresponding fractional core mass. We apply this test to realistic models of helium stars and find that stars typically expand when their cores are near a mass limit. Furthermore, it appears that stars that evolve into giants have always first exceeded an SC-like limit.

We discuss mass transfer in eccentric binaries containing a white dwarf and a neutron star (WD--NS binaries). We show that such binaries are produced from field binaries following a series of mass transfer episodes that allow the white dwarf to form before the neutron star. We predict the orbital properties of binaries similar to the observed WD--NS binary J1141+6545, and show that they will undergo episodic mass transfer from the white dwarf to the neutron star. Furthermore, we describe oil-on-water, a two-phase SPH formalism that we have developed in order to model mass transfer in such binaries.

The effects of rapid rotation on stellar evolution can be profound. We are now beginning to gather enough data to allow a realistic comparison between different physical models. Two key tests for any theory of stellar rotation are first whether it can match observations of the enrichment of nitrogen, and potentially other elements, in clusters containing rapid rotators and secondly whether it can reproduce the observed broadening of the main sequence in the Hertzsprung-Russel diagram. Models of stellar rotation have been steadily increasing in number and complexity over the past two decades but the lack of data makes it difficult to determine whether such additions actually give a closer reflection of reality. One of the most poorly explored features of stellar rotation models is the treatment of angular momentum transport within convective zones. If we treat the core as having uniform specific angular momentum the angular momentum distribution in the star, for a given surface rotation, is dramatically different from what it is when we assume the star rotates as a solid body. The uniform specific angular momentum also generates strong shears which can drive additional transport of chemical elements close to the boundary of a convection zone. A comparison of different models and their reproduction of observable properties with otherwise identical input physics is essential to properly distinguish between them. We compare detailed grids of stellar evolution tracks of intermediate and high-mass stars produced using several models for rotation generated with our new stellar rotation code.

Subdwarf B stars show chemical peculiarities that cannot be explained by diffusion theory alone. Both mass loss and turbulence have been invoked to slow down atomic diffusion in order to match observed abundances. The fact that some sdB stars show pulsations gives upper limits on the amount of mass loss and turbulent mixing allowed. Consequently, non-adiabatic asteroseismology has the potential to decide which process is responsible for the abundance anomalies. We compute for the first time seismic properties of sdB models with atomic diffusion included consistently during the stellar evolution. The diffusion equations with radiative forces are solved for H, He, C, N, O, Ne, Mg, Fe and Ni. We examine the effects of various mass-loss rates and mixed surface masses on the abundances and mode stability. It is shown that the mass-loss rates needed to simulate the observed He abundances (10^{-14}<=Mdot [Msun/yr]<=10^{-13}) are not consistent with observed pulsations. We find that for pulsations to be driven the rates should be Mdot<=10^{-15} Msun/yr. On the other hand, weak turbulent mixing of the outer 10^{-6} Msun can explain the He abundance anomalies while still allowing pulsations to be driven. The origin of the turbulence remains unknown but the presence of pulsations gives tight constraints on the underlying turbulence model.

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed edge on. Transitions between the two implies that the discs may be warped and precessing. Type II X-ray outbursts are thought to occur when the warped disc interacts with the fast stellar wind. We suggest that tides from a misaligned companion neutron star can cause the observed effects. We make numerical models of a Be star decretion disc in which the spin of the Be star is misaligned with the orbital axis of a neutron star companion. Tidal torques from the neutron star truncate the disc at a radius small enough that the neutron star orbit does not intersect the disc unless the eccentricity or misalignment is very large. A magnetic torque from the Be star that is largest at the equator, where the rotation is fastest, is approximated by an inner boundary condition. There are large oscillations in the mass and inclination of the disc as it moves towards a steady state. These large variations may explain the observed changes from Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find the tidal timescale on which the disc warps, precesses and reaches a steady state to be around a year up to a few hundred years. If present, the oscillations in mass and disc inclination occur on a fraction of this timescale depending on the orbital parameters of the binary. The timescales associated with the tidal torque for observed Be star binaries suggest that these effects are important in all but the longest period binaries.

In the symbiotic (or WD+RG) channel of the single degenerate scenario for type Ia supernovae (SNe Ia) the explosions occur a relatively long time after star formation. The birthrate from this channel would be too low to account for all observed SNe Ia were it not for some mechanism to enhance the rate of accretion on to the white dwarf. A tidally enhanced stellar wind, of the type which has been postulated to explain many phenomena related to giant star evolution in binary systems, can do this. Compared to mass-stripping this model extends the space of SNe Ia progenitors to longer orbital periods and hence increases the birthrate to about $0.0069\,{\rm yr^{-1}}$ for the symbiotic channel. Two symbiotic stars, T CrB and RS Oph, considered to be the most likely progenitors of SNe Ia through the symbiotic channel, are well inside the period--companion mass space predicted by our models.

Models of rapidly rotating massive stars at low metallicities show significantly different evolution and higher metal yields compared to non-rotating stars. We estimate the spin-down time-scale of rapid rotating non-convective stars supporting an alpha-Omega dynamo. The magnetic dynamo gives rise to mass loss in a magnetically controlled stellar wind and hence stellar spin down owing to loss of angular momentum. The dynamo is maintained by strong horizontal rotation-driven turbulence which dominates over the Parker instability. We calculate the spin-down time-scale and find that it could be relatively short, a small fraction of the main-sequence lifetime. The spin-down time-scale decreases dramatically for higher surface rotations suggesting that rapid rotators may only exhibit such high surface velocities for a short time, only a small fraction of their main-sequence lifetime.

When a core collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the SN, may reach a substantial space velocity. With binary population synthesis modelling at solar and one fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf-Rayet stars. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc SN or a long or short gamma-ray burst occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed WR runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic SN progenitors travel shorter distances than the progenitors of other SN types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP SN can fly farthest about 48 pc on average at solar metallicity. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway WR stars and that some 30 per cent of GRBs could occur a hundred parsecs or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star--black hole mergers that are roughly ten times more common than neutron star--neutron star mergers. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick.

The stellar origin of gamma-ray bursts can be explained by the rapid release of energy in a highly collimated, extremely relativistic jet. This in turn appears to require a rapidly spinning highly magnetised stellar core that collapses into a magnetic neutron star or a black hole within a relatively massive envelope. They appear to be associated with type Ib/c supernovae but, with a birthrate of around 10^{-6}-10^{-5} per year per galaxy, they are considerably rarer than such supernovae in general. To satisfy all these requirements we hypothesize a binary star model that ends with the merging of an oxygen neon white dwarf with the carbon-oxygen core of a naked helium star during a common envelope phase of evolution. The rapid spin and high magnetic field are natural consequences of such a merging. The evolution that leads to these progenitors is convoluted and so naturally occurs only very rarely. To test the hypothesis we evolve a population of progenitors and find that the rate is as required. At low metallicity we calculate that a similar fraction of stars evolve to this point and so would expect the gamma-ray burst rate to correlate with the star formation rate in any galaxy. This too is consistent with observations. These progenitors, being of intermediate mass, differ radically from the usually postulated high-mass stars. Thus we can reconcile observations that the bursts occur close to but not within massive star associations.

We consider the spin angular momentum evolution of the accreting components of Algol-type binary stars. In wider Algols the accretion is through a disc so that the accreted material can transfer enough angular momentum to the gainer that material at its equator should be spinning at break-up. We demonstrate that even a small amount of mass transfer, much less than required to produce today's mass ratios, transfers enough angular momentum to spin the gainer up to this critical rotation velocity. However the accretors in these systems have spins typically between 10 and $40\,$per cent of the critical rate. So some mechanism for angular momentum loss from the gainers is required. We consider generation of magnetic fields in the radiative atmospheres in a differentially rotating star and the possibility of angular momentum loss driven by strong stellar winds in the intermediate mass stars, such as the primaries of the Algols. Differential rotation, induced by the accretion itself, may produce such winds which carry away enough angular momentum to reduce their rotational velocities to the today's observed values. We apply this model to two systems with initial periods of 5\,d, one with initial masses 5 and $3\,\rm{M}_{\odot}$ and the other with 3.2 and $2\,\rm{M}_{\odot}$. Our calculations show that, if the mass outflow rate in the stellar wind is about $10\,$per cent of the accretion rate and the dipole magnetic field is stronger than about $1\,$kG, the spin rate of the gainer is reduced to below break-up velocity even in the fast phase of mass transfer. Larger mass loss is needed for smaller magnetic fields. The slow rotation of the gainers in the classical Algol systems is explained by a balance between the spin-up by mass accretion and spin-down by a stellar wind linked to a magnetic field.

GD 356 is unique among magnetic white dwarfs because it shows Zeeman-split Balmer lines in pure emission. The lines originate from a region of nearly uniform field strength (delta B/B is approximately 0.1) that covers 10 per cent of the stellar surface in which there is a temperature inversion. The energy source that heats the photosphere remains a mystery but it is likely to be associated with the presence of a companion. Based on current models we use archival Spitzer IRAC observations to place a new and stringent upper limit of 12 Jupiter masses for the mass of such a companion. In the light of this result and the recent discovery of a 115 min photometric period for GD 356, we exclude previous models that invoke accretion and revisit the unipolar inductor model that has been proposed for this system. In this model a highly conducting planet with a metallic core orbits the magnetic white dwarf and, as it cuts through field lines, a current is set flowing between the two bodies. This current dissipates in the photosphere of the white dwarf and causes a temperature inversion. Such a planet is unlikely to have survived the RGB/AGB phases of evolution so we argue that it may have formed from the circumstellar disc of a disrupted He or CO core during a rare merger of two white dwarfs. GD 356 would then be a white dwarf counterpart of the millisecond binary pulsar PSR 1257+12 which is known to host a planetary system.

The surface magnetic field strength of white dwarfs is observed to vary from very little to around 10^9 G. Here we examine the proposal that the strongest fields are generated by dynamo action during the common envelope phase of strongly interacting stars that leads to binary systems containing at least one white dwarf. The resulting magnetic field depends strongly on the electrical conductivity of the white dwarf, the lifetime of the convective envelope and the variability of the magnetic dynamo. We assess the various energy sources available and estimate necessary lifetimes of the common envelope. In the case of a dynamo that leads a randomly oriented magnetic field we find that the induced field is confined to a thin boundary layer at the surface of the white dwarf. This then decays away rapidly upon dispersal of the common envelope. The residual field is typically less than 10^-8 times the strength of the external field. Only in the case where there is some preferential direction to the dynamo-generated field can an induced field, that avoids rapid decay, be produced. We show that a surface field of magnitude a few per cent of the external field may be produced after a few Myr. In this case the residual field strength is roughly proportional to the lifetime of the dynamo activity.

The low-mass X-ray binary microquasar GRO J1655-40 is observed to have a misalignment between the jets and the binary orbital plane. Since the current black hole spin axis is likely to be parallel to the jets, this implies a misalignment between the spin axis of the black hole and the binary orbital plane. It is likely the black holes formed with an asymmetric supernova which caused the orbital axis to misalign with the spin of the stars. We ask whether the null hypothesis that the supernova explosion did not affect the spin axis of the black hole can be ruled out by what can be deduced about the properties of the explosion from the known system parameters. We find that this null hypothesis cannot be disproved but we find that the most likely requirements to form the system include a small natal black hole kick (of a few tens of km/s) and a relatively wide pre-supernova binary. In such cases the observed close binary system could have formed by tidal circularisation without a common envelope phase.

We study the formation of ultra-compact binaries (AM CVn stars and ultra-compact X-ray binaries) with emphasis on the surface chemical abundances of the donors in these systems. Hydrogen is not convincingly detected in the spectra of these systems. Three different proposed formation scenarios involve different donor stars, white dwarfs, helium stars or evolved main-sequence stars. Using detailed evolutionary calculations we show that the abundances of helium WD donors and evolved main-sequence stars are close to equilibrium CNO-processed material, and the detailed abundances correlate with the core temperature and thus mass of the MS progenitors. Evolved MS donors typically have traces of H left. For hybrid or CO white dwarf donors, the carbon and oxygen abundances depend on the temperature of the helium burning and thus on the helium core mass of the progenitors. For helium star donors in addition to their mass, the abundances depend strongly on the amount of helium burnt before mass transfer starts and can range from unprocessed and thus almost equal to CNO-processed matter, to strongly processed and thus C/O rich and N-deficient. We briefly discuss the relative frequency of these cases for helium star donors, based on population synthesis results. Finally we give diagnostics for applying our results to observed systems and find that the most important test is the N/C ratio, which can indicate the formation scenario as well as, in some cases, the mass of the progenitor of the donor. In addition, if observed, the N/O, O/He and O/C ratios can distinguish between helium star and WD donors. Applied to the known systems we find evidence for WD donors in the AM CVn systems GP Com, CE 315 and SDSS J0804+16 and evidence for hybrid WD or very evolved helium star donors in the UCXBs 4U 1626-67 and 4U 0614+09. [Abridged]

We model the overall shape of an accretion disc in a semi-detached binary system in which mass is transfered on to a spinning black hole the spin axis of which is misaligned with the orbital rotation axis. We assume the disc is in a steady state. Its outer regions are subject to differential precession caused by tidal torques of the companion star. These tend to align the outer parts of the disc with the orbital plane. Its inner regions are subject to differential precession caused by the Lense-Thirring effect. These tend to align the inner parts of the disc with the spin of the black hole. We give full numerical solutions for the shape of the disc for some particular disc parameters. We then show how an analytic approximation to these solutions can be obtained for the case when the disc surface density varies as a power law with radius. These analytic solutions for the shape of the disc are reasonably accurate even for large misalignments and can be simply applied for general disc parameters. They are particularly useful when the numerical solutions would be slow.

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are thought to be caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion, $\sigma_k \approx 15\rm km s^{-1}$, is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically $i \ge 25^\circ$. Alternatively a higher velocity kick distribution, $\sigma_k = 265 \rm km s^{-1}$, is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularisation of the binary orbits has taken place. We discuss briefly how this might be achieved.

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities Z = 10-4 and Z = 10-8 . 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-5 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 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.

To measure the onset of mass transfer in eccentric binaries we have developed a two-phase SPH technique. Mass transfer is important in the evolution of close binaries, and a key issue is to determine the separation at which mass transfer begins. The circular case is well understood and can be treated through the use of the Roche formalism. To treat the eccentric case we use a newly-developed two phase system. The body of the donor star is made up from high-mass "water" particles, whilst the atmosphere is modelled with low-mass "oil" particles. Both sets of particles take part fully in SPH interactions. To test the technique we model circular mass-transfer binaries containing a 0.6 Msun donor star and a 1 Msun white dwarf; such binaries are thought to form cataclysmic variable (CV) systems. We find that we can reproduce a reasonable CV mass-transfer rate, and that our extended atmosphere gives a separation that is too large by aproximately 16%, although its pressure scale height is considerably exaggerated. We use the technique to measure the semi-major axis required for the onset of mass transfer in binaries with a mass ratio of q=0.6 and a range of eccentricities. Comparing to the value obtained by considering the instantaneous Roche lobe at pericentre we find that the radius of the star required for mass transfer to begin decreases systematically with increasing eccentricity.

An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic stellar evolution code SSE. We find that, providing the physics package within STARS is set up correctly to match the parameters of the models used to construct SSE, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in SSE would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.

Angular momentum loss (AML) mechanisms and dynamical evolution owing to magnetic braking and gravitational radiation in relativistic binary stars (RBS) are studied with use of physical parameters collected from the literature. We have calculated and compared AML time scales for the RBS with non-degenerate components and double degenerate (DD) systems.

White dwarfs with surface magnetic fields in excess of $1 $MG are found as isolated single stars and relatively more often in magnetic cataclysmic variables. Some 1,253 white dwarfs with a detached low-mass main-sequence companion are identified in the Sloan Digital Sky Survey but none of these is observed to show evidence for Zeeman splitting of hydrogen lines associated with a magnetic field in excess of 1MG. If such high magnetic fields on white dwarfs result from the isolated evolution of a single star then there should be the same fraction of high field white dwarfs among this SDSS binary sample as among single stars. Thus we deduce that the origin of such high magnetic fields must be intimately tied to the formation of cataclysmic variables. CVs emerge from common envelope evolution as very close but detached binary stars that are then brought together by magnetic braking or gravitational radiation. We propose that the smaller the orbital separation at the end of the common envelope phase, the stronger the magnetic field. The magnetic cataclysmic variables originate from those common envelope systems that almost merge. We propose further that those common envelope systems that do merge are the progenitors of the single high field white dwarfs. Thus all highly magnetic white dwarfs, be they single stars or the components of MCVs, have a binary origin. This hypothesis also accounts for the relative dearth of single white dwarfs with fields of 10,000 - 1,000,000G. Such intermediate-field white dwarfs are found preferentially in cataclysmic variables. In addition the bias towards higher masses for highly magnetic white dwarfs is expected if a fraction of these form when two degenerate cores merge in a common envelope. Similar scenarios may account for very high field neutron stars.

The microquasar GRO J1655-40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady state disc warped by Lense-Thirring precession and find the timescale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in the latter case but we find that, in both cases, the lifetime of the mass transfer state is at most a few times the alignment timescale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter-aligned after its supernova kick but that it is most likely to be close to alignment rather than counteralignment now.

We investigate the evolution of accretion luminosity $L_{\rm acc}$ and stellar luminosity ${L_\ast}$ in pre-mainsequence stars. We make the assumption that when the star appears as a Class II object, the major phase of accretion is long past, and the accretion disc has entered its asymptotic phase. We use an approximate stellar evolution scheme for accreting pre-mainsequence stars based on Hartmann, Cassen & Kenyon, 1997. We show that the observed range of values $k = L_{\rm acc}/L_\ast$ between 0.01 and 1 can be reproduced if the values of the disc mass fraction $M_{\rm disc}/M_*$ at the start of the T Tauri phase lie in the range 0.01 -- 0.2, independent of stellar mass. We also show that the observed upper bound of $L_{\rm acc} \sim L_\ast$ is a generic feature of such disc accretion. We conclude that as long as the data uniformly fills the region between this upper bound and observational detection thresholds, then the degeneracies between age, mass and accretion history severely limit the use of this data for constraining possible scalings between disc properties and stellar mass.

We use the Cambridge stellar evolution code STARS to model the evolution of 5-7 solar mass zero-metallicity stars. With enhanced resolution at the hydrogen and helium burning shell in the AGB phases, we are able to model the entire thermally pulsing asymptotic giant branch (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 solar mass. 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.

We compare our latest single and binary stellar model results from the Cambridge STARS code to several sets of observations. We examine four stellar population ratios, the number of blue to red supergiants, the number of Wolf-Rayet stars to O supergiants, the number of red supergiants to Wolf-Rayet stars and the relative number of Wolf-Rayet subtypes, WC to WN stars. These four ratios provide a quantitative measure of nuclear burning lifetimes and the importance of mass loss during various stages of the stars' lifetimes. In addition we compare our models to the relative rate of type Ib/c to type II supernovae to measure the amount of mass lost over the entire lives of all stars. We find reasonable agreement between the observationally inferred values and our predicted values by mixing single and binary star populations. However there is evidence that extra mass loss is required to improve the agreement further, to reduce the number of red supergiants and increase the number of Wolf-Rayet stars.

We consider the shape of an accretion disc whose outer regions are misaligned with the spin axis of a central black hole and calculate the steady state form of the warped disc in the case where the viscosity and surface densities are power laws in the distance from the central black hole. We discuss the shape of the resulting disc in both the frame of the black hole and that of the outer disc. We note that some parts of the disc and also any companion star maybe shadowed from the central regions by the warp. We compute the torque on the black hole caused by the Lense-Thirring precession and hence compute the alignment and precession timescales. We generalise the case with viscosity and hence surface density independent of radius to more realistic density distributions for which the surface density is a decreasing function of radius. We find that the alignment timescale does not change greatly but the precession timescale is more sensitive. We also determine the effect on this timescale if we truncate the disc. For a given truncation radius, the the timescales are less affected for more sharply falling density distributions.

We use the Cambridge stellar evolution code STARS to model the evolution and nucleosynthesis of zero-metallicity intermediate-mass stars. We investigate the effect of duplicity on the nucleosynthesis output of these systems and the potential abundances of the secondaries. The surfaces of zero-metallicity stars are enriched in CNO elements after second dredge up. During binary interaction, such as Roche lobe overflow or wind accretion, metals can be released from these stars and the secondaries enriched in CNO isotopes. We investigate the formation of the two most metal poor stars known, HE 0107-5240 and HE 1327-2326. The observed carbon and nitrogen abundances of HE 0107-5240 can be reproduced by accretion of material from the companion-enhanced wind of a seven solar star after second dredge-up, though oxygen and sodium are underproduced. We speculate that HE 1327-2326, which is richer in nitrogen and strontium, may similarly be formed by wind accretion in a later AGB phase after third dredge-up.

We introduce a set of stellar models for massive stars whose evolution has been affected by mass transfer in a binary system, at a range of metallicities. As noted by other authors, the effect of such mass transfer is frequently more than just rejuvenation. We find that, whilst stars with convective cores which have accreted only H-rich matter rejuvenate as expected, those stars which have accreted He-rich matter (for example at the end stages of conservative mass transfer) evolve in a way that is qualitatively similar to rejuvenated stars of much higher metallicity. Thus the effects of non-conservative evolution depend strongly on whether He-rich matter is amongst the portion accreted or ejected. This may lead to a significant divergence in binary evolution paths with only a small difference in initial assumptions. We compare our models to observed systems and find approximate formulae for the effect of mass accretion on the effective age and metallicity of the resulting star.

Based on the work of OPAL (Iglesias & Rogers 1996) and Alexander & Ferguson (1994), we construct a series of opacity tables for various metallicities Z=0, 0.00001, 0.00003, 0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08 and 0.1. These tables can be easily used in Eggleton's stellar evolution code in place of the old tables without changing the code. The OPAL tables are used for ${\log}_{10}(T/{\rm K})> 3.95$ and Alexander's for ${\log}_{10}(T/{\rm K})< 3.95$. At ${\rm log}_{10}(T/{\rm K})=3.95$, the two groups' data fit well for all hydrogen mass fractions. Conductive opacities are included by reciprocal addition according to the formulae of Yakovlev & Urpin (1980). For different metallicities $1M_\odot$ models are constructed with the new opacity tables. Comparison of 1 and $5M\odot$ models constructed with the older OPAL tables (Iglesias & Rogers 1992) shows that the new opacities have most effect in the late stages of evolution, the extension of the blue loop during helium burning for intermediate-mass and massive stars.

Doubly-degenerate binary systems consisting of two white dwarfs both composed of carbon and oxygen and close enough that mass is transferred from the less massive to the more massive are possible progenitors of type Ia supernovae. If the mass transfer rate is slow enough that the accreting white dwarf can reach 1.38 solar masses then it can ignite carbon degenerately at its centre. This can lead to a thermonuclear runaway and thence a supernova explosion. However if the accretion rate is too high the outer layers of the white dwarf heat up too much and carbon ignites there non-degenerately. A series of mild carbon flashes can then propagate inwards and convert the carbon to neon relatively gently. There is no thermonuclear runaway and no supernova. We have determined the critical rate for accretion and find it to be generally about two-fifths of the Eddington rate. In the case of actual bianry stars the mass transfer rate falls off as mass is transferred so we find that this rate can be very much exceeded at the onset of mass transfer. We determine combinations of white dwarf masses that can lead to central degenerate carbon ignition.

Cataclysmic variables undergo periodic nova explosions during which a finite mass of material is expelled on a short timescale. The system widens and, as a result, the mass-transfer rate drops. This state of hibernation may account for the variety of cataclysmic variable types observed in systems of similar mass and period. In the light of recent changes to the theory of nova ignition and magnetic braking we investigate whether hibernation remains a viable mechanism for creating cataclysmic variable diversity. We model the ratio of time spent as dwarf novae (DNe) to nova-like systems (NLs). Above a critical mass-transfer rate the system is NL and below it a DN. The dominant loss of angular momentum is by magnetic braking but the rate is uncertain. It is also uncertain what fraction of the mass accreted is expelled during the novae. We compare the models of the ratios against the period of the system for different magnetic braking rates and different ejected masses with the ratio of the number of observed NLs to DNe. We deduce that a rate of angular momentum loss a factor of ten smaller than that traditionally assumed is necessary if hibernation is to account for the observed ratios.

Thanks to a stellar evolution code able to compute through the C-flash we link the binary population synthesis of single degenerate progenitors of type Ia supernovae (SNe Ia) to their physical condition at the time of ignition. We show that there is a large range of possible ignition densities and we detail how their probability distribution depends on the accretion properties. The low density peak of this distribution qualitatively reminds of the clustering of the luminosities of Branch-normal SNe Ia. We tighten the possible range of initial physical conditions for explosion models: they form a one-parameter family, independent of the metallicity. We discuss how these results may be modified if we were to relax our hypothesis of a permanent Hachisu wind or if we were to include electron captures.

The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and binary evolution. Our preferred model starts with 12000 single stars and 12000 binaries placed in a Galactic tidal field at 8.0 kpc from the Galactic Centre. Our choices for the initial conditions and for the primordial binary population are explained in detail. At 4 Gyr, the age of M67, the total mass has reduced by 90% as a result of mass loss and stellar escapes. The mass and half-mass radius of luminous stars in the cluster are a good match to observations although the model is more centrally concentrated than observations indicate. The stellar mass and luminosity functions are significantly flattened by preferential escape of low-mass stars. We find that M67 is dynamically old enough that information about the initial mass function is lost, both from the current luminosity function and from the current mass fraction in white dwarfs. The model contains 20 blue stragglers at 4 Gyr which is slightly less than the 28 observed in M67. Nine are in binaries. The blue stragglers were formed by a variety of means and we find formation paths for the whole variety observed in M67. Both the primordial binary population and the dynamical cluster environment play an essential role in shaping the population. A substantial population of short-period primordial binaries (with periods less than a few days) is needed to explain the observed number of blue stragglers in M67.

In the light of recent recalculations of the \el{19}{F}(alpha,p)\el{22}{Ne} reaction rate we present results of the expected yield of \el{19}{F} from Wolf-Rayet (WR) stars. We have computed models using the upper and lower limits for the rate in addition to the recommended rate and hence we constrain the uncertainty in the yield with respect to this reaction. We find a yield of 3.1x10^{-4}M_sol of \el{19}{F} with our recommended rate and a difference of a factor of two between the yields computed with the upper and lower limits. In comparison with previous work we find a difference in the yield of approximately a factor of 4, connected to a different choice of mass loss. Model uncertainties must be carefully evaluated in order to obtain a reliable estimate of the yield of fluorine from WR stars together with its uncertainties.

Evidence of high-velocity features such as those seen in the near-maximum spectra of some Type Ia Supernovae (eg SN 2000cx) has been searched for in the available SNIa spectra observed earlier than one week before B maximum. Recent observational efforts have doubled the number of SNeIa with very early spectra. Remarkably, all SNeIa with early data (7 in our RTN sample and 10 from other programmes) show signs of such features, to a greater or lesser degree, in CaII IR, and some also in SiII 6255A line. High-velocity features may be interpreted as abundance or density enhancements. Abundance enhancements would imply an outer region dominated by Si and Ca. Density enhancements may result from the sweeping up of circumstellar material by the highest velocity SN ejecta. In this scenario, the high incidence of HVFs suggests that a thick disc and/or a high-density companion wind surrounds the exploding white dwarf, as may be the case in Single Degenerate systems. Large-scale angular fluctuations in the radial density and abundance distribution may also be responsible: this could originate in the explosion, and would suggest a deflagration as the more likely explosion mechanism. CSM-interaction and surface fluctuations may coexist, possibly leaving different signatures on the spectrum. In some SNe the HVFs are narrowly confined in velocity, suggesting the ejection of blobs of burned material.

One possible fate of an accreting white dwarf is explosion in a type Ia supernova. However, the route to the thermonuclear runaway has always been uncertain owing to the lack of a convective model consistent with the Urca process. We derive a formalism for convective motions involving two radial flows. This formalism provides a framework for convective models that guarantees self-consistency for chemistry and energy budget, allows time-dependence and describes the interaction of convective motions with the global contraction or expansion of the star. In the one-stream limit, we reproduce several already existing convective models and allow them to treat chemistry. We also suggest as a model easy to implement in a stellar evolution code. We apply this formalism to convective Urca cores in Chandrasekhar mass white dwarfs. We stress that in degenerate matter, nuclear reactions that change the number of electrons strongly influence the convective velocities. We point out the sensitivity of the energy budget on the mixing. We illustrate our model by computing {\it stationary} convective cores with Urca nuclei. We show that even a very small mass fraction of Urca nuclei ($10^{-8}$) strongly influences the convective velocities. Finally, we present preliminary computations of the late evolution of a close to Chandrasekhar mass C+O white dwarf including the convective Urca process.

We derive a new formalism for convective motions involving two radial flows. This formalism provides a framework for convective models that guarantees consistency for the chemistry and the energy budget in the flows, allows time-dependence and accounts for the interaction of the convective motions with the global contraction or expansion of the star. In the one-stream limit the formalism reproduces several existing convective models and allows them to treat the chemistry in the flows. We suggest a version of the formalism that can be implemented easily in a stellar evolution code. We then apply the formalism to convective Urca cores in Chandrasekhar mass white dwarfs and compare it to previous studies. We demonstrate that, in degenerate matter, nuclear reactions that change the number of electrons strongly influence the convective velocities and we show that the net energy budget is sensitive to the mixing. We illustrate our model by computing stationary convective cores with Urca nuclei. Even a very small mass fraction of Urca nuclei (as little as $10^{-8}$) strongly influences the convective velocities. We conclude that the proper modelling of the Urca process is essential for determining the ignition conditions for the thermonuclear runaway in Chandrasekhar-mass white dwarfs.

We determine the possible masses and radii of the progenitors of white dwarfs in binaries from fits to detailed stellar evolution models and use these to reconstruct the mass-transfer phase in which the white dwarf was formed. We confirm the earlier finding that in the first phase of mass transfer in the binary evolution leading to a close pair of white dwarfs, the standard common-envelope formalism equating the energy balance in the system, does not work. An algorithm equating the angular momentum balance can explain the observations. This conclusion is now based on ten observed systems rather than three. With the latter algorithm the separation does not change much for approximately equal mass binaries. Assuming constant efficiency in the standard \alpha-formalism and a constant value of \gamma, we investigate the effect of both methods on the change in separation in general and conclude that when there is observational evidence for strong shrinkage of the orbit, the \gamma-algorithm also leads to this. We then extend our analysis to all close binaries with at least one white dwarf component and reconstruct the mass transfer phases that lead to these binaries. We find that all observations can be explained with a single value of \gamma, making the \gamma-algorithm a useful tool to predict the outcome of common-envelope evolution. We discuss the consequences of our findings for different binary populations in the Galaxy, including massive binaries, for which the reconstruction method cannot be used (abbriged).

A long standing problem with asymptotic giant branch (AGB) star models has been their inability to produce the low-luminosity carbon stars in the Large and Small Magellanic Clouds. Dredge-up must begin earlier and extend deeper. We find this for the first time in our models of LMC metallicity. Such features are not found in our models of SMC metallicity. The fully implicit and simultaneous stellar evolution code STARS has been used to calculate the evolution of AGB stars with metallicities of Z=0.008 and Z=0.004, corresponding to the observed metallicities of the Large and Small Magellanic Clouds, respecitively. Third dredge-up occurs in stars of 1Msol and above and carbon stars were found for models between 1Msol and 3Msol. We use the detailed models as input physics for a population synthesis code and generate carbon star luminosity functions. We now find that we are able to reproduce the carbon star luminosity function of the LMC without any manipulation of our models. The SMC carbon star luminosity function still cannot be produced from our detailed models unless the minimum core mass for third dredge-up is reduced by 0.06Msol.

We discuss the mass ranges over which we find AGB and super-AGB stars. The most massive super-AGB stars are candidate progenitors for type II core-collapse SNe. We discuss the two supernovae, SN1980K and SN2003gd that provide some restrictions on the upper mass limit of super-AGB stars.

We present details of our investigation of the progenitors to core-collapse supernovae. We discuss observations and the theory of the lowest-mass stars to explode as supernovae.

The post-Asymptotic Giant Branch (AGB) star HR4049 is in an eccentric binary system with a relatively short period probably surrounded by a dusty circumbinary disk. Extremely anomalous oxygen isotopic ratios, O16/O17 ~ O16/O18 ~ 7, have been measured from CO_2 molecules likely residing in the disk. Such a composition cannot be explained in the framework of AGB and post-AGB evolution while it can be qualitatively associated with the nucleosynthesis occurring during nova outbursts. We discuss nova models, the presence of a white dwarf companion to HR4049 and possible scenarios for the dynamical evolution of this binary system. Circumbinary disks in which mixing occurs between red-giant and nova material may also be invoked as the site of formation of some rare types of meteoritic presolar grains.

ABRIDGED: To compare the chemistries of stars in the Milky Way dSph galaxies with stars in the Galaxy, we have compiled a large sample of Galactic stellar abundances from the literature. As found in previous studies, the [alpha/Fe] ratios of most stars in the dSph galaxies are generally lower than similar metallicity Galactic stars in this extended sample. Our kinematically selected stars confirm that this is true for all major stellar components of the Galaxy, including the halo, thin disk, and thick disk. There is marginal overlap in the low [alpha/Fe] ratios between dSph stars and Galactic halo stars on extreme retrograde orbits, but this is not supported by other element ratios, like [Ba/Y]. Thus, the full chemical signature of most of the dSph stars are distinct from the main components of the Galaxy. This result rules out continuous merging of low mass galaxies similar to these dSph satellites during the formation of the Galaxy, which is discuss further. The new environments offered by the dSph galaxies also allow us to examine fundamental assumptions related to nucleosynthesis. The metal-poor stars ([Fe/H] < -1.8) in the dSph galaxies have [alpha/Fe] ratio that are consistent with predictions from the alpha-process (alpha-rich freeze out). This may also explain the very low Y abundances, and could suggest a separate r-process site for this light (first peak r-process) element. In stars with higher metallicities ([Fe/H] > -1.8), contributions from the s-process are expected; [Ba/Y] is still much higher in the dSph stars than similar metallicity Galactic stars though, consistent with s-process contributions from only the low metallicity AGB stars in dSph galaxies. Finally, the Na-Ni trend in Galactic halo stars is confirmed, but discuss this in terms of nucleosynthesis, and not the accretion of dSphs.

We present maps of the nature of single star progenitors of supernovae and their remnants in mass and metallicity space. We find our results are similar to others but we have gone further in varying the amount of mixing and using various mass-loss schemes to see how the maps change. We find that extra-mixing, in the form of convective overshooting, moves boundaries such as the minimum mass for a supernova or WR star to lower masses. We also find that the pre-WR mass-loss determines the shape of our maps. We find that different mass-loss rates lead to quite different results. We find that the rise in luminosity at 2nd dredge-up places quite tight constraints on the masses of some progenitors and in particular the progenitor of supernova 2003gd.

We present results of the evolution of AGB stars of 3M_sol and 5M_sol with solar metallicity calculated with the Eggleton stellar evolution code (STARS), which has a fully implicit and simultaneous method for solving for the stellar structure, convective mixing and nuclear burning. We introduce the concept of a viscous mesh in order to improve the numerical stability of the calculations. For the 5M_sol star, we evolve through 25 thermal pulses and their associated third dredge-up events. We obtain a maximum helium luminosity of 1.7x10^9 L_sol and significantly deep dredge-up after the second pulse. Strong hot-bottom burning is observed after the 5th pulse. The 3M_sol model is evolved through 20 thermal pulse events and we find third dredge-up after the 7th pulse. During the 14th pulse sufficient carbon has been brought to the surface to produce a carbon star. We find that dredge-up and the transformation into a carbon star occur at significantly smaller core masses (0.584M_sol and 0.608M_sol, respectively) than in previous calculations for 3M_sol.

We present the discovery of a red supergiant star that exploded as supernova 2003gd in the nearby spiral galaxy M74. The Hubble Space Telescope (HST) and the Gemini Telescope imaged this galaxy 6 to 9 months before the supernova explosion and subsequent HST images confirm the positional coincidence of the supernova with a single,resolved star which is an 8 +4/-2 solar mass red supergiant. This confirms both stellar evolution models and supernova theories which predict that type II-Plateau supernovae have cool red supergiants as their immediate progenitor stars.

We study the evolution of stars that may be the progenitors of long-soft gamma-ray burst (GRBs) -- rotating naked helium stars presumed to have lost their envelopes to winds or companions. Our aim is to investigate the formation and development of single and binary systems and from this population evaluate the rates of interesting individual species. Using a rapid binary evolution algorithm, that enables us to model the most complex binary systems and to explore the effect of metallicity on GRB production, we draw the following conclusions. First we find that, if we include an approximate treatment of angular momentum transport by mass loss, the resulting spin rates for single stars become too low to form a centrifugally supported disc that can drive a GRB engine. Second massive stars in binaries result in enough angular momentum -- due to spin-orbit tidal interactions -- to form a centrifugally supported disc and are thus capable of supplying a sufficient number of progenitors. This holds true even if only a small fraction of bursts are visible to a given observer and the GRB rate is several hundred times larger than the observed rate. Third low-metallicity stars aid the formation of a rapidly rotating, massive helium cores at collapse and so their evolution is likely to be affected by the local properties of the ISM. This effect could increase the GRB formation rate by a factor of 5-7 at Z=Z_solar/200. Finally we quantify the effects of mass loss, common-envelope evolution and black-hole formation and show that more stringent constraints to many of these evolution parameters are needed in order to draw quantitative conclusions from population synthesis work.

We have included opacity tables in our stellar evolution code that enable us to accurately model the structure of stars composed of mixtures with carbon and oxygen independently enhanced relative to solar. We present tests to demonstrate the effects of the new tables. Two of these are practical examples, the effect on the evolution of a thermally pulsing asymptotic giant branch star and a Wolf-Rayet Star. The changes are small but perceptible.

We use N-body simulations of star clusters to investigate the possible dynamical origins of the observed spread in core radius among intermediate-age and old star clusters in the Large Magellanic Cloud (LMC). Two effects are considered, a time-varying external tidal field and variations in primordial hard binary fraction. Simulations of clusters orbiting a point-mass galaxy show similar core radius evolution for clusters on both circular and elliptical orbits and we therefore conclude that the tidal field of the LMC has not yet significantly influenced the evolution of the intermediate-age clusters. The presence of large numbers of hard primordial binaries in a cluster leads to core radius expansion; however, the magnitude of the effect is insufficient to explain the observations. Further, the range of binary fractions required to produce significant core radius growth is inconsistent with the observational evidence that all the LMC clusters have similar stellar luminosity functions.

The supernova SN2001du in the galaxy NGC1365 (19+/-2Mpc), is a core-collapse event of type II-P. Images of this galaxy, have been taken with HST approximately 6.6 years before discovery and include the supernova position on the WFPC2 field of view. We have observed the supernova with the WFPC2 to allow accurate differential astrometry of SN2001du on the pre-explosion frames. There is a marginal detection (3-sigma) of a source close to the supernova position on the prediscovery V-band frame, however it is not precisely coincident and we do not believe it to be a robust detection of a point source. We conclude that there is no stellar progenitor at the supernova position and derive sensitivity limits of the prediscovery images which provide an upper mass limit for the progenitor star. We estimate that the progenitor had a mass of less than 15M_sol. We revisit two other nearby SNe II-P which have high quality pre-explosion images, and refine the upper mass limits for the progenitor stars. Finally we compile all the direct information available for the progenitors of eight nearby core-collapse supernovae and compare their mass estimates. These are compared with the latest stellar evolutionary models of pre-supernova evolution which have attempted to relate metallicity and mass to the supernovae type. Reasonable agreement is found for the lower mass events (generally the II-P), but some discrepancies appear at higher masses. (abridged).

The progenitors of Type Ia supernovae (SNe Ia) have not been identified. Though they are no longer fashionable we investigate the consequences if a significant number of SNe Ia were edge-lit detonations (ELDs) of carbon/oxygen white dwarfs that have accreted a critical mass of helium. Our best understanding of the Phillips relation between light curve speed and peak luminosity assigns both these phenomena to the amount of $^{56}$Ni produced. In ELDs there are two sites of $^{56}$Ni synthesis. If the peak luminosity is determined primarily by the C/O ratio in the core it is primarily a function of its progenitor's initial mass. If the light curve decay speed is determined by the total mass of iron group elements ejected this is a function of the total mass of the ELD at the time of explosion. In general these two masses are correlated and an empirical relation between peak luminosity and light curve shape can be expected. However when we perform population synthesis for progenitors of different metallicities we find a systematic shift in this relation that would make distant SNe Ia fainter than those nearby. The abundances of alpha-rich isotopes, such as $^{44}$Ca, in the solar system indicate that only about 40 per cent of SNe Ia are edge-lit so any systematic effect that could be present would be correspondingly diluted. If we examine only the small subset of ELDs that accrete from a naked helium star, rather than a He white dwarf, the systematic effect disappears.

We have developed a numerical code with which we study the effects of 2D perturbations on stellar structure. We present new numerical and analytical results on the heating of a main-sequence star in a binary system by its companion.

We present a rapid binary evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but because orbits generally circularize before Roche-lobe overflow the outcome of the interactions of systems with the same semi-latus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries, however, the initial separations should be distributed according to the observed distribution of semi-latera recta rather than periods or semi-major axes.

We present high quality ground-based VRI images of the site of the Type II-P SN1999em (in NGC1637) taken before explosion, which were extracted from the CFHT archive. We determine a precise position of the SN on these images to an accuracy of 0.17''. The host galaxy is close enough (7.5 +/- 0.5 Mpc) that the bright supergiants are resolved as individual objects, however we show that there is no detection of an object at the SN position before explosion that could be interpreted as the progenitor star. By determining the sensitivity limits of the VRI data, we derive bolometric luminosity limits for the progenitor. Comparing these to standard stellar evolutionary tracks which trace evolution up to the point of core carbon ignition, we initially derive an upper mass limit of approximately 12M_sol. However we present evolutionary calculations that follow 7-12M_sol stars throughout their C-burning lifetime and show that we can restrict the mass of the progenitor even further. Our calculations indicate that progenitors initially of 8-10M_sol, undergoing expected mass loss, can also be excluded because a second dredge up sends them to somewhat higher luminosities than a star of initially 12M_sol. These results limit the progenitor's initial main-sequence mass to a very narrow range of 12 +/- 1 M_sol. We discuss the similarities between the Type II-P SNe 1999em and 1999gi and their progenitor mass limits, and suggest that SN Type II-P originate only in intermediate mass stars of 8-12M_sol, which are in the red supergiant region and that higher mass stars produce the other Type II sub-types. (Abridged).

Recent observations suggest that long-duration gamma-ray bursts (GRBs) and their afterglows are produced by highly relativistic jets emitted in core-collapse explosions. The pre-explosive ambient medium provides a natural test for the most likely progenitors of GRBs. Those stars that shed their envelopes most readily have short jet crossing times and are more likely to produce a GRB. We construct a simple computational scheme to explore the expected contribution of the presupernova ejecta of single Wolf-Rayet (WR) stars to the circumstellar environment. Using detailed stellar tracks for the evolution of massive stars, we discuss the effects that the initial main sequence mass, metallicity, rotation and membership in a binary system have on the ambient medium. We extend the theory of GRB afterglows in winds to consider the effect of the relativistic fireball propagating through the WR ejecta. We make specific predictions for the interaction of the relativistic blast wave with the density bumps that arise when the progenitor star rapidly loses a large fraction of its initial mass or when the ejected wind interacts with the external medium and decelerates. A re-brightening of the afterglow with a spectrum redder than the typical synchrotron spectrum (as seen in GRB 970508, GRB 980326 and GRB 000911) is predicted. We also calculate the luminosity of the reflected echo that arises when circumstellar material Compton-scatters the prompt radiation and examine the spectral signatures expected from the interaction of the GRB afterglow with the ejected medium.

Masses and progenitor evolutionary states of TypeII supernovae remain almost unconstrained by direct observations. Only one robust observation of a progenitor (SN1987A) and one plausible observation (SN1993J) are available. Neither matched theoretical predictions and in this Letter we report limits on a third progenitor (SN1999gi). The Hubble Space Telescope has imaged the site of the TypeII-P supernova SN1999gi with the WFPC2 in two filters (F606W and F300W) prior to explosion. The distance to the host galaxy (NGC3184) of 7.9Mpc means that the most luminous, massive stars are resolved as single objects in the archive images. The supernova occurred in a resolved, young OB association 2.3kpc from the centre of NGC3184 with an association age of about 4Myrs. Follow-up images of SN1999gi with WFPC2 taken 14 months after discovery determine the precise position of the SN on the pre-explosion frames. An upper limit of the absolute magnitude of the progenitor is estimated (M_v >= -5.1). By comparison with stellar evolutionary tracks this can be interpreted as a stellar mass, and we determine an upper mass limit of 9(+3/-2)M_solar. We discuss the possibility of determining the masses or mass limits for numerous nearby core-collapse supernovae using the HST archive enhanced by our current SNAP programme.

Optically thin cooling gas at most temperatures above 30K will make condensations by pressure pushing material into cool dense regions. This works without gravity. Cooling condensations will flatten and become planar/similarity solutions. Most star formation may start from cooling condensations - with gravity only important in the later stages. The idea that some of the dark matter could be pristine white dwarfs that condensed slowly on to planetary sized seeds without firing nuclear reactions is found lacking. However, recent observations indicate fifty times more halo white dwarfs than have been previously acknowledged; enough to make the halo fraction observed as MACHOS. A cosmological census shows that only 1% of the mass of the Universe is of known constitution.

We present a state-of-the-art N-body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N-body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N-body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N = 10000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.

Using Eggleton's stellar evolution code, we carry out 150 runs of Pop I binary evolution calculations, with the initial primary mass between 1 and 8 solar masses the initial mass ratio between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle, or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass of the primary, the initial mass ratio, and the radius of the primary at the onset of RLOF. We also find that a carbon-oxygen white dwarf with mass as low as 0.33 solar masses may be formed if the primary's initial mass is around 2.5 solar masses.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

We present analytic formulae that approximate the evolution of stars for a wide range of mass and metallicity. Stellar luminosity, radius and core mass are given as a function of age, M and Z, for all phases from the zero-age main-sequence up to, and including, the remnant stages. For the most part we find continuous formulae accurate to within 5% of detailed models. These formulae are useful for purposes such as population synthesis that require very rapid but accurate evaluation of stellar properties, and in particular for use in combination with N-body codes. We describe a mass loss prescription that can be used with these formulae and investigate the resulting stellar remnant distribution.

We investigate the gravonuclear instabilities reported by Bono et al. (1997a,b) during the onset of helium-shell burning at the end of the horizontal-branch (HB) phase. These instabilities are characterized by relaxation oscillations within the helium shell which lead to loops in the evolutionary tracks. We find the occurrence of these instabilities depends critically on how the breathing pulses are suppressed near the end of the HB phase. If they are suppressed by omitting the gravitational energy term in the stellar structure equations, then the helium profile within the core at the end of the HB phase will contain a broad region of varying helium abundance. The helium-burning shell which forms in this region is too thick to be unstable, and gravonuclear instabilities do not occur. If, on the other hand, the breathing pulses are suppressed by prohibiting any increase in the central helium abundance, then the final helium profile can exhibit a large discontinuity at the edge of the helium-exhausted core. The helium shell which forms just exterior to this discontinuity is then much thinner and can be thermally unstable. Even in this case, however, the gravonuclear instabilities disappear as soon as the nuclear burning broadens the helium shell into its characteristic S-shape. We conclude that the gravonuclear instabilities found by Bono et al. are a consequence of the ad hoc procedure used to suppress the breathing pulses.

The position of pre-main-sequence or protostars in the Hertzsprung--Russell diagram is often used to determine their mass and age by comparison with pre-main-sequence evolution tracks. On the assumption that the stellar models are accurate, we demonstrate that, if the metallicity is known, the mass obtained is a good estimate. However, the age determination can be very misleading because it is significantly (generally different by a factor of two to five) dependent on the accretion rate and, for ages less than about one million years, the initial state of the star. We present a number of accreting protostellar tracks that can be used to determine age if the initial conditions can be determined and the underlying accretion rate has been constant in the past. Because of the balance established between the Kelvin-Helmholtz, contraction timescale and the accretion timescale a pre-main-sequence star remembers its accretion history. Knowledge of the current accretion rate, together with an H--R-diagram position gives information about the rate of accretion in the past but does not necessarily improve any age estimate. We do not claim that ages obtained by comparison with these particular accreting tracks are likely to be any more reliable than those from comparisons with non-accreting tracks. Instead we stress the unreliability of any such comparisons and use the disparities between various tracks to estimate the likely errors in age and mass estimates. We also show how a set of coeval accreting objects do not appear coeval when compared with non-accreting tracks. Instead accreting pre-main-sequence stars of around a solar mass are likely to appear older than those of either smaller or larger mass.

We present theoretical models for the evolution of T Tauri stars surrounded by circumstellar discs. The models include the effects of pre-main-sequence stellar and time dependent disc evolution, and incorporate the effects of stellar magnetic fields acting on the inner disc. For single stars, consistency with observations in Taurus-Auriga demands that disc dispersal occurs rapidly, on much less than the viscous timescale of the disc, at roughly the epoch when heating by stellar radiation first dominates over internal viscous dissipation. Applying the models to close binaries, we find that because the initial conditions for discs in binaries are uncertain, studies of extreme mass ratio systems are required to provide a stringent test of theoretical disc evolution models. We also note that no correlation of the infra-red colours of T Tauri stars with their rotation rate is observed, in apparent contradiction to the predictions of simple magnetospheric accretion models.

We show how the second sequence seen lying above the main sequence in cluster colour magnitude diagrams results from binaries with a large range of mass ratios and not just from those with equal masses. We conclude that the presence of a densely populated second sequence, with only sparse filling in between it and the single star main sequence, does not necessarily imply that binary mass ratios are close to unity.

We present a method of determining lower limits on the masses of pre-main-sequence (PMS) stars and so constraining the PMS evolutionary tracks. This method uses the red-shifted absorption feature observed in some emission-line profiles of T Tauri stars indicative of infall. The maximum velocity of the accreting material measures the potential energy at the stellar surface, which, combined with an observational determination of the stellar radius, yields the stellar mass. This estimate is a lower limit owing to uncertainties in the geometry and projection effects. Using available data, we show that the computed lower limits can be larger than the masses derived from PMS evolutionary tracks for M < 0.5 M_\odot. Our analysis also supports the notion that accretion streams do not impact near the stellar poles but probably hit the stellar surface at moderate latitudes.

This paper describes the first optical spectroscopic survey of class I sources (also known as embedded sources and protostars) in the Taurus-Auriga dark cloud. We detect 10 of the 24 known class I sources in the cloud at 5500-9000 A. All detected class I sources have strong H-alpha emission; most also have strong [O I] and [S II] emission. These data - together with high quality optical spectra of T Tauri stars in the Taurus-Auriga cloud - demonstrate that forbidden emission lines are stronger and more common in class I sources than in T Tauri stars. Our results also provide a clear discriminant in the frequency of forbidden line emission between weak-emission and classical T Tauri stars. In addition to strong emission lines, three class I sources have prominent TiO absorption bands. The M-type central stars of these sources mingle with optically visible T Tauri stars in the HR diagram and lie somewhat below both the birthline for spherical accretion and the deuterium burning sequence for disc accretion.

Dynamical models of star clusters are maturing in the sense that effects other than simple point particle dynamics are taken into account. We summarize the relevance of and prospects for this new generation of N-body models.

We investigate the dependence of theoretically generated mass - (absolute magnitude) relations on stellar models. Using up to date physics we compute models in the mass range 0.1 < m < 1M_sun. We compare the solar-metallicity models with our older models, with recent models computed by others, and also with an empirical mass - (absolute magnitude) relation that best fits the observed data. At a given mass below 0.6M_sun the effective temperatures differ substantially from model to model. However taken individually each set of models is in good agreement with observations in the mass - luminosity plane. A minimum in the derivative dm/dM_V at M_V = 11.5, which is due to H_2 formation and establishment of a fully convective stellar interior, is present in all photometric bands, for all models. This minimum leads to a maximum in the stellar luminosity function for Galactic disk stars at M_V = 11.5, M_bol = 9.8. Stellar models should locate this maximum in the stellar luminosity function at the same magnitude as observations. Models which incorporate the most realistic theoretical atmospheres and the most recent equation of state and opacities can satisfy this constraint. These models are also in best agreement with the most recent luminosity - (effective temperature) and mass-luminosity data. Each set of our models of a given metallicity (with 0.2 > [Fe/H] > -2.3) shows a maximum in -dm/dM_bol, which moves to brighter bolometric magnitudes with decreasing metallicity. The change in location of the maximum, as a function of [Fe/H], follows the location of structure in luminosity functions for stellar populations with different metal abundances. This structure seen in all observed stellar populations can be accounted for by the mass--luminosity relation.

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiter's mass around a solar-like star could have dynamically survived in an orbit with a period as short as $\sim10\,$hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun.

We present a simple and efficient, yet reasonably accurate, equation of state, which at the moderately low temperatures and high densities found in the interiors of stars less massive than the Sun is substantially more accurate than its predecessor by Eggleton, Faulkner & Flannery. Along with the most recently available values in tabular form of opacities, neutrino loss rates, and nuclear reaction rates for a selection of the most important reactions, this provides a convenient package of input physics for stellar modelling. We briefly discuss a few results obtained with the updated stellar evolution code.