$XY$ model with higher-order exchange

An $XY$ model, generalized by inclusion of up to an infinite number of higher-order pairwise interactions with an exponentially decreasing strength, is studied by spin-wave theory and Monte Carlo simulations. At low temperatures the model displays a quasi-long-range order phase characterized by an algebraically decaying correlation function with the exponent $\eta = T/[2 \pi J(p,\alpha)]$, nonlinearly dependent on the parameters $p$ and $\alpha$ that control the number of the higher-order terms and and the decay rate of their intensity, respectively. At higher temperatures the system shows a crossover from the continuous Berezinskii-Kosterlitz-Thouless to the first-order transition for the parameter values corresponding to a highly nonlinear shape of the potential well. The role of of topological excitations (vortices) in changing the nature of the transition is discussed.