The kinetic curves' molar balances (Supplementary Information) were implemented in MATLAB R2022b. A non-linear least squares method was used to fit the effective optical path length (or photon flux for the batch system) in the balance. The differential equations were solved with a variable-step, variable-order (VSVO) solver for stiff differential equations.