Many functional properties of particle system rely on collective behaviour and the type of superstructures formed when thousands of particles come together. Self-assembly, agglomeration and aggregation depend sensitively on the size and shape of particles present, as well as the size distribution and the mixtures of shapes within a given sample, which makes simulation of these superstructures and their properties challenging. Here we present a new, flexible, software package for the simulations of ordered and disorder aggregates of faceted polyhedral particle from the nanoscale to the micronscale, which is capable of including size distributions and mixtures of multiple particle shapes defined by the User, subject to additional User-defined interactions. Following relaxation using molecular dynamics a number of characterisation tools are provided, including interfacial probabilities and distribution functions. The software is applicable to a range of problems from nanoparticle assembly to additive manufacturing.Aggregates of nanoparticles can provide a range of unique properties that are different from isolated, individual entities. Depending on the particle size(s) and shape(s), the overall porosity, strength and complexity of the aggregates the reactivity, transport and efficiency of entire samples can be controlled [1][2][3][4][5][6][7]. Examples include the tribological properties of powders [8], the toxicological properties of colloids [9], the mechanical properties of composites [10], the collective electromagnetic properties exhibited in surface plasmon resonances and nanoporous materials formed via aggregation or self-assembly [11].Characterising the superstructure of aggregates is much more difficult than ordered self-assemblies of nanoparticles, due to the variations in particles sizes, shapes, interaction points, types of interactions and the inherent randomness. Probing with experimental methods can be challenging due to interference from the probe and the collection of either averaged results that represent entire samples, or local results (limited by the probe size) that may not be representative at all. Both local interactions and long-ranged statistics can be obtained from computational methods, but these have their challenges too. Simulations without periodic boundary conditions (PBCs) lack the ability to control the number density, while simulations with PBCs impose an unnatural translational symmetry and necessitates larger systems or repetition to develop the right statistics. The use of classical atomistic simulations methods limits the number of nanoparticles that can be practically handled, while the use of electronic structure methods limits the size and number of nanoparticles that can be handled, leading to an artificial bias toward minority features such as edges and corners. The aggregation of most realistic systems are driven by formation of either facet-facet interfaces (depending on the size and type of the facets) and the formation of voids. This cannot be achieved if we are limite...