Multi-discontinuity algorithm for world-line Monte Carlo simulations

We introduce a novel multi-discontinuity algorithm for efficient global update of world-line configurations in Monte Carlo simulations of interacting quantum systems. This new algorithm is a generalization of the two-discontinuity algorithms introduced in Refs. [N. Prokof'ev, B. Svistunov, and I. Tupitsyn, Phys. Lett. A {\bf 238}, 253 (1998)] and [O. Sylju{\aa}sen and A. Sandvik, Phys. Rev. E {\bf 66}, 046701 (2002)] . This generalization is particularly effective for studying Bose-Einstein condensates (BEC) of composite particles. In particular, we demonstrate the utility of the generalized algorithm by simulating a Hamiltonian for an S=1 anti-ferromagnet with strong uniaxial single-ion anisotropy. The multi-discontinuity algorithm not only solves the freezing problem that arises in this limit, but also allows for efficiently computing the off-diagonal correlator that characterizes a BEC of composite particles.