The characterization of the thermal behaviour of zeolites is made experimentally difficult by the combination of several factors. First, the analysis of thermal expansion requires the presence of single crystal data, which are difficult to obtain on fully evacuated/activated samples; in some cases in the literature, organic templates or solvent molecules are still present in the pores, which affects the thermal behaviour. Secondly, several zeolitic frameworks undergo phase transitions upon heating, which makes their thermal expansion harder to analyse. Therefore, the computational study of the thermal behaviour of zeolites is an area of active research. In a previous work in our group, the thermal expansion coefficient of the pure-silica zeolites were systematically computed through DFT in the quasi-harmonic approximation (QHA) , i.e., in a low-temperature regime, finding that all the investigated frameworks displayed NTE within the temperature range of 10-300 K. The exploitation of DFT and QHA to investigate the thermal properties of the whole zeolitic database is a very demanding task from a computational point of view, therefore an alternative cheaper and as accurate approach is desired. The MLIP was trained on trajectories at different volumes and degrees of strain, therefore the rotation of the RUMs should be captured by the model.