Abstract
We develop a random sampling method for calculating the time evolution of the R\'{e}nyi 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\'{e}nyi 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_s× N_s matrix corresponding to N_s sites at half filling, the sampling cost becomes O(2^α N_s) with a constant α≪ 1, in contrast to the conventional algorithm with the O(2^N_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.