Entanglement entropy dynamics of non-Gaussian states in free boson systems: Random sampling approach

We develop a random sampling method for calculating the time evolution of the R\'enyi entanglement entropy after a quantum quench from an insulating state in free boson systems. Because of the non-Gaussian nature of the initial state, calculating the R\'enyi entanglement entropy calls for the exponential cost of computing a matrix permanent. We numerically demonstrate that a simple random sampling method reduces the computational cost of a permanent; for an ${N}_{\mathrm{s}}\ifmmode\times\else\texttimes\fi{}{N}_{\mathrm{s}}$ matrix corresponding to ${N}_{\mathrm{s}}$ sites at half filling, the sampling cost becomes $O({2}^{\ensuremath{\alpha}{N}_{\mathrm{s}}})$ with a constant $\ensuremath{\alpha}\ensuremath{\ll}1$, in contrast to the conventional algorithm with the $O({2}^{{N}_{\mathrm{s}}})$ number of summations requiring the exponential time cost. Although the computational cost is still exponential, this improvement allows us to obtain the entanglement entropy dynamics in free boson systems for more than 100 sites. We present several examples of the entanglement entropy dynamics in low-dimensional free boson systems.