Large deviations of density fluctuations in the boundary-driven quantum symmetric simple inclusion process

We consider the boundary driven Quantum Symmetric Simple Inclusion Process (QSSIP) which describes a one-dimensional system of bosonic particles with stochastic nearest-neighbor hopping, modeled as a Brownian motion, with gain/loss processes at the endpoints of the chain driving the system out-of-equilibrium. Although the averaged QSSIP dynamics differs from that of the Quantum Symmetric Simple Exclusion Process (QSSEP)-the analogous system where bosons are replaced by fermions-we show that, paradoxically, the dynamics of their matrices of two-point functions, along with all their fluctuations, coincide. In contrary, the underlying classical models differs significantly, as the bosonic statistics allow the inclusion of multiple particles at the same site, in contrast to (symmetric) simple exclusion processes (SSEP). We provide an exact derivation of the large deviation function of density fluctuations in QSSIP and, as a consequence, in the classical inclusion process (SSIP) by exploiting its quantum formulation. Remarkably, our study highlights that, both in QSSEP and QSSIP, fluctuations of the local densities are typically classical, i.e., the cumulant generating functions of the local densities are asymptotically self-averaging and converge toward those of the classical SSEP and SSIP, realization-wise. This provides a test of the conjectured almost sure classical behavior of transport fluctuations, at leading order in the system size, in noisy diffusive quantum many-body systems.