All training was originally performed with the main (stable) development branch of NeqiIP, which does not include the computation of the stress, and the potentials resulting from this branch were unable to describe the volume variation of the zeolites, and could not be used in constant-pressure MD simulations. To overcome this shortcoming, the ML potential obtained from the main branch was updated by performing a single step of a training adopting the stress-branch version of NequIP, which encoded the evaluation of the virial on the final form of the deployed potential. To check that this update did not mess up the potential, the errors on the training set were computed with both the versions, which resulted in negligible differences on all the metrics considered (below the third decimal point).