Metastable phases of SiO 2 , also known as siliceous zeolites, are typically formed by crystallisation in the presence of templating agents. This produces nanoporous structures that have applications in adsorption, separation and gas storage. More than 200 unique zeolite topologies have been observed from simple tetrahedral SiO 4 building blocks connected by their edges, though much less have been obtained as pure silica compounds. This structural diversity of zeolitic frameworks induces a wide variety of their properties, including geometric characteristics (e.g. channel dimensionality, connectivity and pore space) which are directly relevant for adsorption-related applications, but also in physical properties such as mechanical behaviour. Indeed, among the large array of topologies observed in known zeolites, there are unique structures that display mechanical properties unusual in the inorganic world, such as large-scale elastic anisotropy. This also includes a structure that displays complete auxeticity, a rare property demonstrated by few materials, and even fewer crystalline structures. Moreover, it has been proposed that mechanical stability, reflected in the elastic properties of the frameworks, plays a key role in the selection of feasible zeolite structures, i.e. experimentally synthesised structures are among the most mechanically stable potential frameworks. Commonly, the computational determination of mechanical properties of crystalline materials in the elastic regime involves the optimisation and energy determination of many deformation modes, requiring greater computational cost compared to routine optimisations. This can be performed using first-principles methods for the determination of energies at the quantum chemical level, and efficient computing techniques have allowed for the calculation of extensive databases of material properties. Recently de Jong et. al. reported the elastic properties for 1181 inorganic solids (most of them nonporous structures), and this database is being continually updated. However, even with current advanced computing resources, high-throughput methods can only be applied to a small fraction of realised crystal structures as even for simple quaternary inorganic compounds there is estimated 10 10 combinations. As an alternative, computationally cheap classical interatomic pair potentials can also be employed 10 however, the resulting accuracy is expected to be significantly less.