Bosonic reaction-diffusion processes on scale-free networks
  Reaction-diffusion processes can be adopted to model a large number of dynamics on complex networks, such as transport processes or epidemic outbreaks. In most cases, however, they have been studied from a fermionic perspective, in which each vertex can be occupied by at most one particle. While still useful, this approach suffers from some drawbacks, the most important probably being the difficulty to implement reactions involving more than two particles simultaneously. Here we introduce a general framework for the study of bosonic reaction-diffusion processes on complex networks, in which there is no restriction on the number of interacting particles that a vertex can host. We describe these processes theoretically by means of continuous time heterogeneous mean-field theory and divide them into two main classes: steady state and monotonously decaying processes. We analyze specific examples of both behaviors within the class of one-species process, comparing the results (whenever possible) with the corresponding fermionic counterparts. We find that the time evolution and critical properties of the particle density are independent of the fermionic or bosonic nature of the process, while differences exist in the functional form of the density of occupied vertices in a given degree class k. We implement a continuous time Monte Carlo algorithm, well suited for general bosonic simulations, which allow us to confirm the analytical predictions formulated within mean-field theory. Our results, both at the theoretical and numerical level, can be easily generalized to tackle more complex, multi-species, reaction-diffusion processes, and open a promising path for a general study and classification of this kind of dynamical systems on complex networks. 
Andrea Baronchelli, Michele Catanzaro, Romualdo Pastor-Satorras