Scaling behavior in the number theoretic division model of self-organized criticality

We revisit the number theoretic division model of self-organized criticality [Phys. Rev. Lett. 101, 158702 (2008)]. The model consists of a pool of $M-1$ ordered integers $\{2,3,\cdots ,M\}$, and the aim is to dynamically form a primitive set of integers, where no number can be divided or divisible by others. Using intensive simulation studies and finite-size scaling method, we find the primitive set size fluctuations in the division model to show power spectral density of the form $1\mathrm{/}{f}^{\alpha }$ in the frequency regime $1\mathrm{/}M\ll f\ll 1\mathrm{/}2$ with $\alpha \approx 2$ (different from $\alpha \approx 1.80(1)$ as reported previously) along with an additional scaling in terms of the system size $\sim M^b$. We also show similar power spectra properties for a class of random walks with a power-law distributed jump size (Lévy flights).