We also used the cluster models of the zeolite frameworks (See SI Section 2) for computing the free energy for the proton exchange reactions. These calculations were performed using a static approach, in which the cluster models for the reactant and the transition state, were first optimized. We performed the calculations at the level of PBE with Grimme's-D3 dispersion correction 73 and 6-311+G(d,p) basis set using the Gaussian software. The initial configuration of the transition state structure was taken from the corresponding TASS trajectory. The transition state optimization was then performed to locate the stationary point with the imaginary frequency for the required normal mode. Then, intrinsic reaction coordinate was computed to confirm the transition state and to determine the reactant and the product states. Free energies of the reactant and the transition states were calculated by adding entropic contributions at 300 K computed using the model partition functions: