We adopt the neural network (NN) flow method to study the Berezinskii-Kosterlitz-Thouless (BKT) phase transitions of the two-dimensional q-state clock model with q≥4. The NN flow consists of a sequence of the same units that proceed with the flow. This unit is a variational autoencoder trained by the data of Monte Carlo configurations in unsupervised learning. To gauge the difference among the ensembles of Monte Carlo configurations at different temperatures and the uniqueness of the ensemble of NN-flow states, we adopt the Jensen-Shannon divergence (JSD) as the information-distance measure "thermometer." This JSD thermometer compares the probability distribution functions of the mean spin value of two ensembles of states. Our results show that the NN flow will flow an arbitrary spin state to some state in a fixed-point ensemble of states. The corresponding JSD of the fixed-point ensemble takes a unique profile with peculiar features, which can help to identify the critical temperatures of BKT phase transitions of the underlying Monte Carlo configurations.
We adopt the neural network (NN) flow method to study the Berezinskii-Kosterlitz-Thouless (BKT) phase transitions of the two-dimensional q-state clock model with q≥4. The NN flow consists of a sequence of the same units that proceed with the flow. This unit is a variational autoencoder trained by the data of Monte Carlo configurations in unsupervised learning. To gauge the difference among the ensembles of Monte Carlo configurations at different temperatures and the uniqueness of the ensemble of NN-flow states, we adopt the Jensen-Shannon divergence (JSD) as the information-distance measure "thermometer." This JSD thermometer compares the probability distribution functions of the mean spin value of two ensembles of states. Our results show that the NN flow will flow an arbitrary spin state to some state in a fixed-point ensemble of states. The corresponding JSD of the fixed-point ensemble takes a unique profile with peculiar features, which can help to identify the critical temperatures of BKT phase transitions of the underlying Monte Carlo configurations.
We demonstrate that a machine learning technique with a simple feedforward neural network can sensitively detect two successive phase transitions associated with the Berezinskii-Kosterlitz-Thouless (BKT) phase in q-state clock models simultaneously by analyzing the weight matrix components connecting the hidden and output layers. We find that the method requires only a data set of the raw spatial spin configurations for the learning procedure. This data set is generated by Monte Carlo thermalizations at selected temperatures. Neither prior knowledge of, for example, the transition temperatures, number of phases, and order parameters nor processed data sets of, for example, the vortex configurations, histograms of spin orientations, and correlation functions produced from the original spin-configuration data are needed, in contrast with most of previously proposed machine learning methods based on supervised learning. Our neural network evaluates the transition temperatures as T2/J=0.921 and T1/J=0.410 for the paramagnetic-to-BKT transition and BKT-to-ferromagnetic transition in the eight-state clock model on a square lattice. Both critical temperatures agree well with those evaluated in the previous numerical studies.
We propose two machine-learning methods based on neural networks, which we respectively call the phase-classification method and the temperature-identification method, for detecting different types of phase transitions in the XXZ models without prior knowledge of their critical temperatures. The XXZ models have exchange couplings which are anisotropic in the spin space where the strength is represented by a parameter $\mathrm{\ensuremath{\Delta}}(>0)$. The models exhibit second-order phase transition when $\mathrm{\ensuremath{\Delta}}>1$, but the Berezinskii-Kosterlitz-Thouless (BKT) phase transition when $\mathrm{\ensuremath{\Delta}}<1$. In the phase-classification method, the neural network is trained using spin or vortex configurations of well-known classical spin models other than the XXZ models, e.g., the Ising models and the XY models, to classify those of the XXZ models to corresponding phases. We demonstrate that the trained neural network successfully detects the phase transitions for both $\mathrm{\ensuremath{\Delta}}>1$ and $\mathrm{\ensuremath{\Delta}}<1$, and the evaluated critical temperatures coincide well with those evaluated by conventional numerical calculations. In the temperature-identification method, on the other hand, the neural network is trained so as to identify temperatures at which the input spin or vortex configurations are generated by the Monte Carlo thermalization. The critical temperatures are evaluated by analyzing the optimized weight matrix, which coincides with the result of numerical calculation for the second-order phase transition in the Ising-like XXZ model with $\mathrm{\ensuremath{\Delta}}=1.05$ but cannot be determined uniquely for the BKT transition in the XY-like XXZ model with $\mathrm{\ensuremath{\Delta}}=0.95$.
Berezinskii--Kosterlitz--Thouless transition of the two-dimensional XY model on the honeycomb lattice
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 21 exactly. Surprisingly, the value of TBKT,H obtained from our NN calculations is 0.560(9) which deviates significantly from 21. In addition, based on the helicity modulus, the TBKT,H determined is 0.571(8) agreeing well with that resulting from the NN estimation. The outcomes presented in this study indicate that a detailed analytic calculation is desirable to solve the found discrepancy.
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.
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.