The Binary Population and Spectral Synthesis (BPASS) suite of binary stellar evolution models and synthetic stellar populations provides a framework for the physically motivated analysis of both the integrated light from distant stellar populations and the detailed properties of those nearby. We present a new version 2.1 data release of these models, detailing the methodology by which BPASS incorporates binary mass transfer and its effect on stellar evolution pathways, as well as the construction of simple stellar populations. We demonstrate key tests of the latest BPASS model suite demonstrating its ability to reproduce the colours and derived properties of resolved stellar populations, including well- constrained eclipsing binaries. We consider observational constraints on the ratio of massive star types and the distribution of stellar remnant masses. We describe the identification of supernova progenitors in our models, and demonstrate a good agreement to the properties of observed progenitors. We also test our models against photometric and spectroscopic observations of unresolved stellar populations, both in the local and distant Universe, finding that binary models provide a self-consistent explanation for observed galaxy properties across a broad redshift range. Finally, we carefully describe the limitations of our models, and areas where we expect to see significant improvement in future versions.

Wolf-Rayet (WR) stars are massive stars that have lost most or all of their hydrogen via powerful stellar winds. Recent observations have indicated that hydrogen-free WR stars have cooler temperatures than those predicted by current evolutionary models. To investigate how varying mass-loss rate affects WR evolution, we have created a grid of pure helium star models. We compare our results with Galactic and LMC WR observations and show that the temperature ranges of observed WR stars can be reproduced by varying the mass-loss rate, which effectively determines the size of the helium envelope around the core. We also find that WN and WO stars arise from more massive stars, whereas WC stars come from lower masses. This contradicts the standard Conti scenario by which WN and WC stars evolve in an age sequence. We also predict the magnitudes of our models at core-collapse and compare with observations of nearby progenitors of Type Ib/c supernovae. We confirm the findings of previous studies that suggest WR stars are the progenitors of core-collapse supernovae: the progenitors would remain unobserved except in the cases where the progenitor is a low-mass helium giant, as is the case of iPTF13bvn.