As a characteristic property of all quantum systems, entanglement participates in many important quantum phenomena. In this proceeding, we employ it in the study of quantum field theories at finite density. We incorporate evaluations of entanglement entropy using the replica trick into MC simulations of $O(N)$ models at finite density with the worm algorithm and present some initial results for the nonlinear $O(4)$ model in 3 dimensions.