The ZIF-8 simulation showed imidazolate linker decomposition as well as limited H 2 and N 2 formation at 1750 K on this timescale (Supplementary Figures ). There was significant structural change as many zinc imidazolate bonds were broken at 1750 K, forming a trigonal planar zinc environment. This is reflected in the average zinc coordination number dropping to three (Figure ). With quenching, this coordination number recovers to almost four, forming a zeolitic-like structure, which is distinct from the initial crystalline MOF structure (Supplementary Figure ). The stability of ZIF-8 linkers relative to CALF-20 linkers decomposing more readily is reflected in experimental results which show that CALF-20 and ZIF-8 decompose at approximately 650 K and 870 K respectively. Initially this quench was conducted at 2000 K for CALF-20 and ZIF-8, however we found instabilities that led to close contact between atoms, where atoms drifted to within less than 1 Å separation, causing unstable trajectories. While diverse environments were seen over the 15 ps simulation, it appears that the 1 ns simulation created new environments which the MLP was unable to treat due to the limited timescale of the DFT training data. However, the foundational model was stable for 2000 K simulations while the finetuned MACE-MP-0 model was stable up to 1000 K for CALF-20. The stabilities we observed are reported in Supplementary Table . This suggests that the regions of the PES explored in large spatial and temporal scale simulations are poorly parametrized for the CALF-20 and ZIF-8 specific trained models. This is despite the "from scratch" and finetuned models having significantly lower error than MACE-MP-0 at 2000 K. Additionally, the transferability of these models was investigated by validating them against 1000 frames randomly sampled from a 15 ps trajectory of MOF-5 and MOF-10 at 300 K. These validations show that the foundational model most accurately predicts MOF-5 and MOF-10 followed by the "from scratch" model and then the finetuned (Supplementary Table ). This reflects the order of high temperature stability for each model. Overall, this suggests the regions of the PES absent in the MOF-temperature datasets are poorly captured in the "from scratch" model and finetuned models compared to the foundational model. This is likely due to the structural diversity in the foundational training datasets creating a well parametrized PES, leading to higher model stability and greater transferability.