The two Berezinskii–Kosterlitz–Thouless phase transitions of the two-dimensional 5-state clock model are studied on infinite strips using the DMRG algorithm. Because of the open boundary conditions, the helicity modulus ϒ 2 is computed by imposing twisted magnetic fields at the two boundaries. Its scaling behavior is in good agreement with the existence of essential singularities with σ = 1/2 at the two transitions. The predicted universal values of ϒ 2 are shown to be reached in the thermodynamic limit. The fourth-order helicity modulus is observed to display a dip at the high-temperature BKT transition, like the XY model, and shown to take a new universal value at the low-temperature BKT transition. Finally, the scaling behavior of magnetization at the low-temperature transition is compatible with η = 1/4.
The two Berezinskii–Kosterlitz–Thouless phase transitions of the two-dimensional 5-state clock model are studied on infinite strips using the DMRG algorithm. Because of the open boundary conditions, the helicity modulus ϒ 2 is computed by imposing twisted magnetic fields at the two boundaries. Its scaling behavior is in good agreement with the existence of essential singularities with σ = 1/2 at the two transitions. The predicted universal values of ϒ 2 are shown to be reached in the thermodynamic limit. The fourth-order helicity modulus is observed to display a dip at the high-temperature BKT transition, like the XY model, and shown to take a new universal value at the low-temperature BKT transition. Finally, the scaling behavior of magnetization at the low-temperature transition is compatible with η = 1/4.
Abstract In a two-dimensional (2D) spin system, the XY model, characterized by planar rotational symmetry, exhibits a unique phenomenon known as the Berezinskii–Kosterlitz–Thouless (BKT) transition. In contrast, the clock model, which introduces discrete rotational symmetry, exhibits the BKT transition at two different temperatures due to this discreteness. In this study, we numerically investigate the BKT transition for XY and six-state clock models over various 2D lattices. We employ two primary methods: the Monte Carlo method, which analyzes the size dependence of the ratio of the correlation functions for two different distances, and a machine-learning approach to classify the different phases—namely, the low-temperature ordered phase, the intermediate BKT phase, and the high-temperature disordered phase. We identify the BKT transition temperatures for the XY and six-state models on honeycomb, kagome, and diced lattices. Combined with the previously calculated data for the triangular lattice, we then compare these values with the second-order phase transition temperatures of the 2D Ising model, for which exact solutions are known. Our results indicate that the ratio of the BKT transition temperatures for each lattice relative to the Ising model transition temperatures are close, although the values are not universal.
Abstract In a two-dimensional (2D) spin system, the XY model, characterized by planar rotational symmetry, exhibits a unique phenomenon known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. In contrast, the clock model, which introduces discrete rotational symmetry, exhibits the BKT transition at two different temperatures due to this discreteness. In this study, we numerically investigate the BKT transition for XY and six-state clock models over various two-dimensional lattices. We employ two primary methods: the Monte Carlo method, which analyzes the size dependence of the ratio of the correlation functions for two different distances, and a machine-learning approach to classify the different phases --- namely, the low-temperature ordered phase, the intermediate BKT phase, and the high-temperature disordered phase. We identify the BKT transition temperatures for the XY and six-state models on honeycomb, kagome, and diced lattices. Combined with the previously calculated data for the triangular lattice, we then compare these values with the second-order phase transition temperatures of the 2D Ising model, for which exact solutions are known. Our results indicate that the ratio of the BKT transition temperatures for each lattice relative to the Ising model transition temperatures are close, although the values are not universal.
Abstract In this paper, we thoroughly examined the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in the two-dimensional XY model on the honeycomb lattice. To address its thermodynamical behavior, we combined standard numerical Monte Carlo simulations with the simulated annealing (SA) protocol and entropic simulations based on the Wang-Landau (WL) algorithm. The transition temperature was determined using the second (ϒ) and fourth-order (ϒ4) helicity modulus as the order parameter. Our best finite-size scaling analysis suggests T BKT = 0.575(8) from SA and T BKT = 0.576(3) from WL. These values deviate significantly from the expected theoretical value of 1 / 2 . We believe that this discrepancy suggests that the theoretical assumptions regarding the analytical calculation may need to be revisited. Additionally, we calculated the vortex density and the formation energy of the vortex-antivortex pairs, where the obtained vortex formation energy is 2μ = 5.80(12). Upon comparison with the square lattice, our results support the notion of instability of the honeycomb lattice to support the spin long-range order and give additional backing to the critical behavior we found.
The q-state clock model, sometimes called the discrete XY model, is known to show a second-order (symmetry-breaking) phase transition in two dimensions for q≤4 (q=2 corresponds to the Ising model). On the other hand, the q→∞ limit of the model corresponds to the XY model, which shows the infinite order (non-symmetry-breaking) Berezinskii-Kosterlitz-Thouless (BKT) phase transition in two dimensions. Interestingly, the two-dimensional clock model with q≥5 is predicted to show three different phases and two associated phase transitions. There are varying opinions about the actual characters of phases and the associated transitions. In this work, we develop the basic and higher-order mean-field (MF) theories to study the q-state clock model systematically. Our MF calculations reaffirm that, for large q, there are three phases: (broken) Z_{q} symmetric ferromagnetic phase at the low temperature, emergent U(1) symmetric BKT phase at the intermediate temperature, and paramagnetic (disordered) phase at the high temperature. The phase transition at the higher temperature is found to be of the BKT type, and the other transition at the lower temperature is argued to be a large-order spontaneous symmetry-breaking type (the largeness of transition order yields the possibility of having some of the numerical characteristics of a BKT transition). The higher-order MF theory developed here better characterizes phases by estimating the spin-spin correlation between two neighbors.
The 2D XY model exhibits an essential phase transition, which was predicted long ago—by Berezinskii, Kosterlitz and Thouless (BKT)—to be driven by the (un)binding of vortex–anti-vortex pairs. This transition has been confirmed for the standard lattice action, and for actions with distinct couplings, in agreement with universality. Here we study a highly unconventional formulation of this model, which belongs to the class of topological lattice actions: it does not have any couplings at all, but just a constraint for the relative angles between nearest neighbour spins. By means of dynamical boundary conditions we measure the helicity modulus ϒ, which shows that this formulation performs a BKT phase transition as well. Its finite size effects are amazingly mild, in contrast to other lattice actions. This provides one of the most precise numerical confirmations ever of a BKT transition in this model. On the other hand, up to the lattice sizes that we explored, there are deviations from the spin wave approximation, for instance for the Binder cumulant U 4 and for the leading finite size correction to ϒ. Finally we observe that the (un)binding mechanism follows the usual pattern, although free vortices do not require any energy in this formulation. Due to this observation, one should reconsider an aspect of the established picture, which estimates the critical temperature based on this energy requirement.
Abstract The Berezinskii–Kosterlitz–Thouless (BKT) transition of the two-dimensional XY model on the honeycomb lattice is investigated using both the techniques of Neural network (NN) and Monte Carlo simulations. It is demonstrated in the literature that, with certain plausible assumptions, the associated critical temperature TBKT,H is found to be 1/2 exactly. Surprisingly, the value of TBKT,H obtained from our NN calculations is 0.572(3), which deviates significantly from 1/2. In addition, based on the helicity modulus, the TBKT,H determined is 0.576(4), agreeing well with that resulting from the NN estimation. It will be interesting to carry out a more detailed analytic calculation to obtain a theoretical value consistent with the numerical result reached here.