We now return to the bi-exponential cooling of the zeolite samples that was observed experimentally. The simulated sample cooling was also found to be bi-exponential regardless of the value of thermal conductivity used, with only the magnitude of the time constants varying with κ. Fit parameters are reported in the ESI, Section 9. The τ1:τ2 ratio and A1 / A2 for a given start temperature were the same, whether averaged over the FWHM diameter of the probe beam, or at different single-positions in the sample. The simulations therefore indicate that the bi-exponential cooling is not a consequence of averaging the temperature over the probed region of the sample. Although both cooling processes are slowed by the presence of water clusters, τ2 is more strongly affected, leading to the τ2/τ1 ratio being proportional to the level of water clusters at low temperatures. The neglect of the properties of water from the simulation explains the inability to match the simulated cooling to experimental results at temperatures ≤ 125 °C.