Tracking the variation of entanglement R\'enyi negativity: an efficient quantum Monte Carlo method

Although the entanglement entropy probing novel phases and phase transitions numerically via quantum Monte Carlo (QMC) has achieved huge success in pure ground states of quantum many-body systems, numerical explorations on mixed states remain limited, despite the fact that most real-world systems are non-isolated. Meanwhile, entanglement negativity, as a rarely computable entanglement monotone for mixed states, is significant in characterizing mixed-state entanglement, such as in systems with two disconnected regions, dissipation or at finite temperature. However, efficient numerical approaches are scarce to calculate this quantity in large-scale and high-dimensional systems, especially when we need to access how it varies with certain parameters to study critical behaviors. Within the reweight-annealing frame, we present an accessible and efficient QMC algorithm, which is able to achieve the values as well as tracking the variation of the R\'enyi version of entanglement negativity on some specified parameter path. Our algorithm makes it feasible to directly study the role that entanglement plays at the critical point and in different phases for mixed states in high dimensions numerically. In addition, this method is accessible and easy to parallelize on computers. Through this method, different intrinsic mechanisms in quantum and thermal criticalities with the same universal class have been revealed clearly through the numerical calculations on R\'enyi negativity.
Comment: 10 pages, 4 figures