We investigate the partition function zeros of the two-dimensional p-state clock model in the complex temperature plane by using the Wang-Landau method. For p=5, 6, 8, and 10, we propose a modified energy representation to enumerate exact irregular energy levels for the density of states without any binning artifacts. Comparing the leading zeros between different p's, we provide strong evidence that the upper transition at p=6 is indeed of the Berezinskii-Kosterlitz-Thouless (BKT) type in contrast to the claim of the previous Fisher zero study [Phys. Rev. E 80, 042103 (2009)10.1103/PhysRevE.80.042103]. We find that the leading zeros of p=6 at the upper transition collapse onto the zero trajectories of the larger p's including the XY limit while the finite-size behavior of p=5 differs from the converged behavior of p≥6 within the system sizes examined. In addition, we argue that the nondivergent specific heat in the BKT transition is responsible for the small partition function magnitude that decreases exponentially with increasing system size near the leading zero, fundamentally limiting access to large systems in search for zeros with an estimator under finite statistical fluctuations.
We investigate the partition function zeros of the two-dimensional p-state clock model in the complex temperature plane by using the Wang-Landau method. For p=5, 6, 8, and 10, we propose a modified energy representation to enumerate exact irregular energy levels for the density of states without any binning artifacts. Comparing the leading zeros between different p's, we provide strong evidence that the upper transition at p=6 is indeed of the Berezinskii-Kosterlitz-Thouless (BKT) type in contrast to the claim of the previous Fisher zero study [Phys. Rev. E 80, 042103 (2009)10.1103/PhysRevE.80.042103]. We find that the leading zeros of p=6 at the upper transition collapse onto the zero trajectories of the larger p's including the XY limit while the finite-size behavior of p=5 differs from the converged behavior of p≥6 within the system sizes examined. In addition, we argue that the nondivergent specific heat in the BKT transition is responsible for the small partition function magnitude that decreases exponentially with increasing system size near the leading zero, fundamentally limiting access to large systems in search for zeros with an estimator under finite statistical fluctuations.
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 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.
We investigate the six-state clock model with nearest-neighbor interactions on the square lattice. We obtain the density of states of the finite systems up to L=28 using the Wang-Landau sampling. With the density of states and the Fisher zero approach, we successfully find two different critical temperatures 0.632(2) and 0.997(2) for the clock model. It seems that this study supports the recent results by [Lapilli Phys. Rev. Lett. 96, 140603 (2006)] that the transitions are not of Kosterlitz and Thouless type.
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.
In this paper, we revisit the q-states clock model for small systems. We present results for the thermodynamics of the q-states clock model from q=2 to q=20 for small square lattices L×L, with L ranging from L=3 to L=64 with free-boundary conditions. Energy, specific heat, entropy and magnetization are measured. We found that the Berezinskii-Kosterlitz-Thouless (BKT)-like transition appears for q>5 regardless of lattice size, while the transition at q=5 is lost for L<10; for q≤4 the BKT transition is never present. We report the phase diagram in terms of q showing the transition from the ferromagnetic (FM) to the paramagnetic (PM) phases at a critical temperature T1 for small systems which turns into a transition from the FM to the BKT phase for larger systems, while a second phase transition between the BKT and the PM phases occurs at T2. We also show that the magnetic phases are well characterized by the two dimensional (2D) distribution of the magnetization values. We make use of this opportunity to do an information theory analysis of the time series obtained from the Monte Carlo simulations. In particular, we calculate the phenomenological mutability and diversity functions. Diversity characterizes the phase transitions, but the phases are less detectable as q increases. Free boundary conditions are used to better mimic the reality of small systems (far from any thermodynamic limit). The role of size is discussed.
The two-dimensional q-state clock model for $q\ensuremath{\ge}5$ undergoes two Berezinskii-Kosterlitz-Thouless (BKT) phase transitions as temperature decreases. Here we report an extensive worm-type simulation of the square-lattice clock model for q=5--9 in a pair of flow representations, from high- and low-temperature expansions, respectively. By finite-size scaling analysis of susceptibilitylike quantities, we determine the critical points with a precision improving over the existing results. Due to the dual flow representations, each point in the critical region is observed to simultaneously exhibit a pair of anomalous dimensions, which are ${\ensuremath{\eta}}_{1}=1/4$ and ${\ensuremath{\eta}}_{2}=4/{q}^{2}$ at the two BKT transitions. Further, the approximate self-dual points ${\ensuremath{\beta}}_{\mathrm{sd}}(L)$, defined by the stringent condition that the susceptibilitylike quantities in both flow representations are identical, are found to be nearly independent of system size L and behave as ${\ensuremath{\beta}}_{\mathrm{sd}}\ensuremath{\simeq}q/2\ensuremath{\pi}$ asymptotically at the large-q limit. The exponent $\ensuremath{\eta}$ at ${\ensuremath{\beta}}_{\mathrm{sd}}$ is consistent with 1/q within statistical error as long as $q\ensuremath{\ge}5$. Based on this, we further conjecture that $\ensuremath{\eta}({\ensuremath{\beta}}_{\mathrm{sd}})=1/q$ holds exactly and is universal for systems in the q-state clock universality class. Our work provides a vivid demonstration of rich phenomena associated with the duality and self-duality of the clock model in two dimensions.