### 7.1  Regarding the Mathematical Definition of SD

To implement synaptic depression in our model, we introduced an adaptation variable $A$ to our model, governed by equations 3.2 to 3.4. These equations represent synaptic depression as a convolution with an alpha kernel with rate $α$ and timescale $τA$. We solved the convolution integral to express it as a set of differential equations (see equations 3.5 and 3.6). However, equation 3.6 contains an erroneous scaling of the input firing rate $r$ with the capitation time constant $τA$. The convolution integral in equation 3.2 can be expressed as a second-order differential equation:
$τAA¨=-2A˙-AτA+αr.$
(7.1)
The equation can be transformed into a set of two coupled first-order differential equations under a simple change of variables:
$τAA˙=B,$
(7.2)
$τAB˙=-2B-A+ατAr.$
(7.3)
Equation 7.3 differs from equation 3.6 merely by a constant scaling of $r$ with $τA$. Since we have never varied $τA$ in this study, this has no impact on our results. Still, the reported values of $α$ for which we find synchronized bursting hold only for synaptic depression given by equations 7.2 and 7.3.