In the case of non-equilibrium structures, 100 non-equilibrium structures were created for each equilibrium structure. It was accomplished by displacing the equilibrium structures along linear combinations of normal mode coordinates computed at the DFTB3+MBD level within the harmonic approximation. It was ensured that the generated structures follow a Boltzmann distribution with a specific average energy difference from the corresponding equilibrium structure.