Helicity modulus in the bilayer XY model by worm algorithm

Abstract: The behaviour of the helicity modulus has been frequently employed to investigate the onset of the topological order characterizing the low-temperature phase of the two-dimensional XY-model. We here present how the analysis based on the use of this key quantity can be applied to the study of the properties of coupled layers. To this aim, we first discuss how to extend the popular worm algorithm to a layered sample, and in particular to the evaluation of the longitudinal helicity, that we introduce taking care of the fact that the virtual twist representing the elastic deformation one applies to properly define the helicity modulus can act on a single layer or on all of them. We then apply the method to investigate the bilayer XY-model, showing how the helicity modulus can be used to determine the phase diagram of the model as a function of temperature and inter-layer coupling strength