where β m , η n and λ mn are eigenvalues obtained according to boundary conditions (see Supplementary section S1). While G contains an infinite sum of modes, the decay time of each mode is 1 /αλ 2 mn , which decreases as integer values m and n increase. Hence, the sum can be truncated once sufficiently short time scales have been reached.