Simulation trajectory files were analyzed using the ummap package, which was modified for this work to enable determination of branching in micelle structures following the method presented in Conchuir et al. (see below). The trajectory data comprise simulation snapshots taken every 2000 DPD time steps (i. e. 20 DPD time units ≈ 1 ns). The N agg (number average) values are extracted as previously described in Anderson et al. after partitioning the surfactants into sub-micellar and micellar aggregates. N agg was plotted as a function of time to understand self-assembly behavior. Shape metrics determined to understand the average morphology of micelles were calculated as reported in the original ummap paper using cut off values ϵ = 0.333, ϵ rod = 0.5, ϵ disc = 1.0, N rod = N disc = 100. These values were determined by trial and error to ensure that observed structures (i. e. from visualizing simulation outputs) corresponded to the shapes described through automatic identification using the ummap package. In this work we combine spherical, oblate and prolate spheroids resulting from the shape analysis into one 'spheroid' class as the distinction between these structures does not not provide detail relating to viscosity (all structures in the spheroid class would result in low viscosity). Micelles with gyration tensors that reported as rod or ellipsoid by ummap are classed as 'rod-like'. This was decided upon after visualising the resulting ellipsoids resulting from the gyration tensor which were found in reality to be exclusively 'bent' rods. All other classifications remain the same as the original citation except lamellar which refers to a structure with two principle extents (defined using the gyration tensor) longer than the bounding box size. In this work, these structures are generally branched micelles (determined via visual inspection) and we refer to these as 2D extended. Those with three extents larger than the box dimensions are referred to as 3D extended. It should be noted that for aggregates that self-interact across the boundaries only the unreplicated structure (suitably unwrapped) is used for shape analysis (and aggregation size) which limits their total extents. Hence care must be taken in interpreting the 1D to 3D extended structures where replicas can be staggered and stacked within the box. For these a clearer interpretation is given using the below 'skeleton' method which also indicates connections between the replicated structures.