i.e. the convolution of g t (t) and e -αλ 2 mn t . g v (r, z) is the product of the optical absorption coefficient and the optical fluence, as defined in Eq. S2. By utilizing Eq. 6, we can generate spatial temperature distributions within the absorber and matrix at various time points (Figure ). In frequency domain T is the product of the Fourier transforms of g t (t) and e -αλ 2 mn t , allowing an easy way to study location dependent amplitudes of temperature changes (see Figure )