To understand the time dependent temperature distribution within the ZSM-5 pellets during and following the T-jump laser pulses, the spatial temperature profile and its evolution over time was simulated in three dimensions. Assuming radial symmetry about the excitation laser pulse centre, the model was based on the 2D axisymmetric heat equation in cylindrical coordinates solved numerically using finite element analysis in Mathematica Version 13.2. The model was solved with logarithmic time steps, to allow high time resolution immediately after the laser pulse, while still following the cooling profile out as far as 1 second. Further details are provided in the ESI, Section 8.