Abstract
Entanglement entropy (EE) plays a central role in the intersection of quantum information science and condensed matter physics. However, scanning the EE for two-dimensional and higher-dimensional systems still remains challenging. To address this challenge, we propose a quantum Monte Carlo scheme capable of extracting large-scale data of Rényi EE with high precision and low technical barrier. Its advantages lie in the following aspects: a single simulation can obtain the continuous variation curve of EE with respect to parameters, greatly reducing the computational cost; the algorithm implementation is simplified, and there is no need to alter the spacetime manifold during the simulation, making the code easily extendable to various many-body models. Additionally, we introduce a formula to calculate the derivative of EE without resorting to numerical differentiation from dense EE data. We then demonstrate the feasibility of using EE and its derivative to find phase transition points, critical exponents, and various phases.