Study of the fully frustrated clock model using the Wang–Landau algorithm

Antiferromagnetic triangular Blume-Capel model with hard-core exclusions

Using Monte Carlo simulation, we analyze phase transitions of two antiferromagnetic (AFM) triangular Blume-Capel (BC) models with AFM interactions between third-nearest neighbors. One model has hard-core exclusions between the nearest-neighbor (1NN) particles (3NN1 model) and the other has them between the nearest-neighbor and next-nearest-neighbor particles (3NN12 model). Finite-size scaling analysis reveals that in these models, the transition from the paramagnetic to long-range order (LRO) AFM phase is either of the first order or goes through an intermediate phase which might be attributed to the Berezinskii-Kosterlitz-Thouless (BKT) type. The properties of the low-temperature phase transition to the AFM phase of the 1NN, 3NN1, and 3NN12 models are found to be very similar for almost all values of a normalized single-ion anisotropy parameter, 0 < δ < 1.5. Higher temperature behavior of the 3NN12 and 3NN1 models is rather different from that of the 1NN model. Three phase transitions are observed for the 3NN12 model: from the paramagnetic phase to the phase with domains of the LRO AFM phase at T(c), from this structure to the diluted frustrated BKT-type phase at T(2), and from the frustrated phase to the AFM LRO phase at T(1). For the 3NN12 model, T(c) > T(2) > T(1) at 0 < δ < 1.15 (range I), T(c) ≈ T(2) > T(1) at 1.15 < δ < 1.3 (range II), and T(c) = T(2) = T(1) at 1.3 < δ < 1.5 (range III). For the 3NN1 model, T(c) ≈ T(2) > T(1) at 0 < δ < 1.2 (range II) and T(c) = T(2) = T(1) at 1.2 < δ < 1.5 (range III). There is only one first-order phase transition in range III. The transition at T(c) is of the first order in range II and either of a weak first order or a second order in range I.

Nonequilibrium critical relaxation of the order parameter and energy in the two-dimensional ferromagnetic Potts model

The static and dynamic critical properties of the ferromagnetic q -state Potts models on a square lattice with q=2 and 3 are numerically studied via the nonequilibrium relaxation method. The relaxation behavior of both the order parameter and energy as well as that of the second moments are investigated, from which static and dynamic critical exponents can be obtained. We find that the static exponents thus obtained from the relaxation of the order parameter and energy together with the second moments of the order parameter exhibit a close agreement with the exact exponents, especially for the case of the q=2 (Ising) model, when care is taken in the choice of the initial states for the relaxation of the second moments. As for the case of q=3 , the estimates for the static exponents become less accurate, but still exhibit reasonable agreement with the exactly known static exponents. The dynamic critical exponent for the q=2 (Ising) model is estimated from the relaxation of the second moments of the order parameter with mixed initial conditions to give z(q=2) approximately 2.1668(19) .

Fast flat-histogram method for generalized spin models

We present a Monte Carlo method that efficiently computes the density of states for spin models having any number of interaction per spin. By combining a random walk in the energy space with collective updates controlled by the microcanonical temperature, our method yields dynamic exponents close to their ideal random-walk values, reduced equilibrium times, and very low statistical error in the density of states. The method can host any density of states estimation scheme, including the Wang-Landau algorithm and the transition matrix method. Our approach proves remarkably powerful in the numerical study of models governed by long-range interactions, where it is shown to reduce the algorithm complexity to that of a short-range model with the same number of spins. We apply the method to the -state Potts chains with power-law decaying interactions in their first-order regime; we find that conventional local-update algorithms are outperformed already for sizes above a few hundred spins. By considering chains containing up to spins, which we simulated in fairly reasonable time, we obtain estimates of transition temperatures correct to five-figure accuracy. Finally, we propose several efficient schemes aimed at estimating the microcanonical temperature.