The results presented in the previous section show behaviour consistent with the formation of extended, multi-chain aggregates in poor solvents, as observed experimentally. Although the increase in R g is relatively limited due to system-size and time-scale limitations, the steady growth of the RMS R g with aggregate size suggests that large aggregates are feasible, with values already reaching multiple times that of a single chain. Previous Monte Carlo simulations of a generic bead-spring model of a semiflexible polymer indicated that the thermodynamically favoured aggregate in a poor solvent is one in which all monomers between chains overlap to give a fully stacked bundle of chains. However, the formation of fully overlapping aggregates cannot explain the large values of R g observed experimentally or in our simulations in a poor solvent. Instead, to explain the observed forma-