In order to generate a diverse and physically representative set of zeolite configurations, we performed ab initio molecular dynamics (AIMD) simulations, coupling Born-Oppenheimer dynamics with evaluation of the energy and atomic forces at the DFT level. AIMD simulations were run on the 4 strained structures of each of the 17 zeolite frameworks considered. QUICKSTEP/CP2K with the same computational setting adopted in the relaxation of the distorted lattices was employed.