performed through the open source molecular dynamics software LAMMPS took ~6 hours and ~6.5 days by using one and two node(s) (128 GB of RAM) composed of 2 processors (Intel® Xeon® CPU E5-2680 v4 @ 2.40GHz, 14 cores) on the MatriCs platform (Université de Picardie-Jules Verne), respectively. The initial structure (prior the slurry equilibration) was generated to match the experimental composition and AM particles size distribution by locating randomly the AM and CBD particles in a simulation box big enough (200 × 200 × 500 μm, x × y × z) to avoid significant particles' overlap. All the details about the random structure generation can be found in the Supporting Information of Ref. The slurry simulation was performed in a NPT environment, at 298 K and 1 atm, while the drying simulation was performed in a NVT environment at 353 K, while the pressure was considered as described in the Computational section. The boundary conditions of the slurry were considered as periodic in all the directions (x, y and z), while for the z direction two repulsive planes have been added at the bottom and top of the slurry structures to constrain the z direction, for which not periodic boundary conditions were applied. The slurry and drying simulations used a timestep of 1 and 0.1 ns, while the number of timestep was 13×10 6 and 100×10 , respectively.