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.

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.

We investigate the kinematic evolution of planetesimals in self-gravitating discs, combining Smoothed Particle Hydrodynamical (SPH) simulations of the disc gas with a gravitationally coupled population of test particle planetesimals. We find that at radii of 10s of au (which is where we expect planetesimals to be possibly formed in such discs) the planetesimals' eccentricities are rapidly pumped to values $>$ 0.1 within the timescales for which the disc is in the self-gravitating regime. The high resulting velocity dispersion and the lack of planetesimal concentration in the spiral arms means that the collision timescale is very long and that the effect of those collisions that do occur is destructive rather than leading to further planetesimal growth. We also use the SPH simulations to calibrate Monte Carlo dynamical experiments: these can be used to evolve the system over long timescales and can be compared with analytical solutions of the diffusion equation in particle angular momentum space. We find that if planetesimals are only formed in a belt at large radius then there is significant scattering of objects to small radii; nevertheless the majority of planetesimals remain at large radii. If planetesimals indeed form at early evolutionary stages, when the disc is strongly self-gravitating, then the results of this study constrain their spatial and kinematic distribution at the end of the self-gravitating phase.

We investigate the red supergiant problem: the apparent dearth of Type IIP supernova progenitors with masses between 16 and 30 M_sun. Although red supergiants with masses in this range have been observed, none have been identified as progenitors in pre-explosion images. We show that by failing to take into account the additional extinction resulting from the dust produced in the red supergiant winds, we risk underestimating the luminosity of the most massive red supergiants at the end of their lives. We estimate the initial masses of all Type IIP progenitors for which observations exist and analyse the resulting population. We find that the most likely maximum mass for a Type IIP progenitor is 21^{+2}_{-1} M_sun. This is in closer agreement with the limit predicted from single star evolution models.