Phase transition of two-dimensional generalized XY model

We study the two-dimensional generalized XY model that depends on an integer $q$ by the Monte Carlo method. This model was recently proposed by Romano and Zagrebnov. We find a single Kosterlitz-Thouless (KT) transition for all values of $q$, in contrast with the previous speculation that there may be two transitions, one a regular KT transition and another a first-order transition at a higher temperature. We show the universality of the KT transitions by comparing the universal finite-size scaling behaviors at different values of $q$ without assuming a specific universal form in terms of the KT transition temperature $T_{\rm KT}$.

💡 Research Summary

The paper investigates the phase‑transition behavior of a two‑dimensional generalized XY model whose Hamiltonian is
( H = -J\sum_{\langle ij\rangle}\cos\bigl(q\theta_i - q\theta_j\bigr) ),
where (q) is a positive integer. For (q=1) the model reduces to the standard XY model, while larger (q) values increase the discrete symmetry of the angular variable, prompting speculation that the system might exhibit two distinct transitions: a conventional Kosterlitz‑Thouless (KT) transition at lower temperature and a first‑order transition at a higher temperature.

To resolve this issue the authors performed extensive Monte‑Carlo simulations using a combination of Metropolis updates, over‑relaxation moves, and parallel tempering. System sizes ranged from (L=16) to (L=256) and several values of (q) (2, 3, 5, 10) were examined. Primary observables included the complex magnetization (M=\langle e^{i\theta}\rangle), its susceptibility, the helicity modulus (\Upsilon), spin‑spin correlation functions, and thermodynamic quantities such as energy and specific heat.

The hallmark of a KT transition is a universal jump of the helicity modulus and a logarithmic finite‑size scaling of (\Upsilon) and the correlation length (\xi). The authors analyzed these quantities without first fixing the transition temperature (T_{\rm KT}). Instead they introduced a dimensionless scaling variable (x = L/\xi(T)) and demonstrated that data for all system sizes and all examined (q) collapse onto a single universal curve. This collapse occurs for the helicity modulus, the susceptibility, and the Binder cumulant, confirming that the critical behavior is governed by the same KT fixed point irrespective of (q).

In the high‑temperature regime the energy and specific‑heat curves show sharp variations, but these features diminish as (L) increases and no hysteresis or latent‑heat signatures are observed. Binder‑cumulant analysis further reveals a smooth, continuous evolution of the free‑energy landscape across the transition, contradicting the previously suggested first‑order scenario.

The authors therefore conclude that the generalized XY model possesses only a single KT transition for every integer (q). The transition temperature does shift slightly with (q), but the universal finite‑size scaling form remains unchanged. By demonstrating that the KT transition survives the introduction of higher‑order angular discretization, the work reinforces the notion that the unbinding of vortex–antivortex pairs is the essential mechanism behind the transition and is robust against modifications of the underlying symmetry.

Beyond settling the specific controversy about double transitions, the study provides a methodological template for investigating topological transitions in models with enlarged symmetry groups. The approach of performing scaling analyses without a priori knowledge of (T_{\rm KT}) can be applied to other low‑dimensional systems where the nature of the transition is ambiguous. The results thus have broader implications for theoretical and experimental studies of two‑dimensional superfluids, thin magnetic films, and cold‑atom setups where generalized XY‑type interactions may be engineered.