To obtain representative structures, mixed morphologies have been simulated with ab-initio MD based on Density Functional Tight Binding theory using linear scaling self-consistent field calculations within the CP2K code. Initial configurations have been prepared using Packmol, targeting experimental values for densities and mole percentages. This structure is first equilibrated at 700 K in N pT (with velocity rescaling thermostat 41 at atmospheric pressure 42 ) for 7 ps (time step 1 fs), then annealed to 300 K within 10 ps. A final N pT equilibration followed for 5 ps.