Anomalous finite-size scaling in higher-order processes with absorbing states

We study standard and higher-order birth-death processes on fully connected networks, within the perspective of large-deviation theory (also referred to as Wentzel-Kramers-Brillouin (WKB) method in some contexts). We obtain a general expression for the leading and next-to-leading terms of the stationary probability distribution of the fraction of "active" sites, as a function of parameters and network size $N$. We reproduce several results from the literature and, in particular, we derive all the moments of the stationary distribution for the $q$-susceptible-infected-susceptible ($q-SIS$) model, i.e., a high-order epidemic model requiring of $q$ active ("infected") sites to activate an additional one. We uncover a very rich scenario for the fluctuations of the fraction of active sites, with non-trivial finite-size-scaling properties. In particular, we show that the variance-to-mean ratio diverges at criticality for $[1 \leq q\leq 3]$, with a maximal variability at $q=2$, confirming that complex-contagion processes can exhibit peculiar scaling features including wild variability and that the leading-order in a large-deviation approach does not suffice to describe them: next-to-leading terms are essential to capture the intrinsic singularity at the origin of systems with absorbing states.