Spatial entanglement entropy in the ground state of the Lieb-Liniger model

We consider the entanglement between two spatial subregions in the Lieb-Liniger model of bosons in one spatial dimension interacting via a contact interaction. Using ground-state path integral quantum Monte Carlo we numerically compute the Renyi entropy of the reduced density matrix of the subsystem as a measure of entanglement. Our numerical algorithm is based on a replica method previously introduced by the authors, which we extend to efficiently study the entanglement of spatial subsystems of itinerant bosons. We confirm a logarithmic scaling of the Renyi entropy with subsystem size that is expected from conformal field theory, and compute the nonuniversal subleading constant for interaction strengths ranging over two orders of magnitude. In the strongly interacting limit, we find agreement with the known free fermion result.