Lastly, we note that the FF parameters for the CGMD model presented in this work are not fully validated to experimental viscosity curves (for the slurry), or porosity and density measurements (for the dry electrode). Instead, here we aim at qualitatively studying the effect of Si volume expansion on a given graphite/Si microstructure, for which obtaining highly experimentally validated microstructures is not of paramount importance. Nevertheless, using the above presented method to generate microstructures has various advantages when compared to cheaper stochastic microstructure generation. For instance, the model has the potential to be improved by fitting FF parameters to reproduce experimental physical properties. Additionally, DEM methods to simulate electrode calendering could be easily included in the workflow . Summing up, Figure shows the overall workflow used to generate the graphite/Si composite microstructures. The structures were equilibrated using a Nose-Hoover barostat and thermostat at 300 K and 1 atm. The simulations were performed using LAMMPS molecular dynamics software [57] using periodic boundary conditions in all directions. The FF parameters used for this work are included in Table . All simulations were performed in the TGCC Joliot-Curie (CEA) supercomputer, in 64 core AMD Rome (Epyc) processors. More details on the generation of CGMD electrode structures can be found in Ref .