We study the phase diagram of the two-dimensional fully frustrated XY model (FFXY) and of two related models, a lattice discretization of the Landau-Ginzburg-Wilson Hamiltonian for the critical modes of the FFXY model, and a coupled Ising-XY model. We present Monte Carlo simulations on square lattices L × L , L ≲ 10 3 . We show that the low-temperature phase of these models is controlled by the same line of Gaussian fixed points as in the standard XY model. We find that, if a model undergoes a unique transition by varying temperature, then the transition is of first order. In the opposite case we observe two very close transitions: a transition associated with the spin degrees of freedom and, as temperature increases, a transition where chiral modes become critical. If they are continuous, they belong to the Kosterlitz-Thouless and to the Ising universality class, respectively. Ising and Kosterlitz-Thouless behavior is observed only after a preasymptotic regime, which is universal to some extent. In the chiral case, the approach is nonmonotonic for most observables, and there is a wide region in which finite-size scaling is controlled by an effective exponent ν eff ≈ 0.8. This explains the result ν ≈ 0.8 of many previous studies using smaller lattices.
The fully frustrated XY model with Villain interaction on a square lattice is studied by means of Monte Carlo simulations. On the basis of the universal jump condition it is argued that there are two distinct transitions in the model, corresponding to the loss of XY order and
Z2 order, respectively. The Kosterlitz-Thouless (KT) transition is analyzed by finite-size scaling of the helicity modulus at lattices of size L=32 through 128, giving
TKT\ensuremath{\approx}0.8108(1). The vorticity-vorticity correlation function is used to determine two different characteristic lengths, the
Z2 correlation length \ensuremath{\xi}, and the screening length \ensuremath{\lambda}, associated with the KT transition and free vortices. The temperature dependence of \ensuremath{\xi} is examined in order to determine
Tc and the correlation length exponent, \ensuremath{\nu}. The exponent is found to be consistent with the two-dimensional (2D) Ising value, \ensuremath{\nu}=1, and the obtained critical temperature is
Tc=0.8225(5). The determinations of both \ensuremath{\xi} and \ensuremath{\nu} are done carefully, first applying the techniques to the 2D Ising model, which serves as a convenient testing ground.
Using Monte Carlo and spin dynamics techniques we have investigated the critical behavior of the classical fully frustrated XY and fully frustrated Planar Rotor models in two dimensions in square lattices of size up to 256. We have found that the models have two transitions a Berezinskii-Kosterliz-Thouless transition at TBKT = 0:3655(5) and an Ising transition at TI = 0:3690(3) for the fully frustrated XY model and a Berezinskii-Kosterliz-Thouless transition at TBKT = 0:4410(5) and an Ising transition at TI = 0:4505(5) for the fully frustrated Planar Rotor model. After this _rst step we have simulated the Z2 (Ising) transition in both fully frustrated Planar Rotator and fully frustrated XY models. We have shown that they can be understood as a percolation problem if we de_ne properly the connectivity of the system. The critical exponent of the mean cluster size, is found as = 2:2(2) close to the two dimensional Ising case, = 2:389. The critical Fisher's exponent is obtained as _ = 1:8(1), close to the two dimensional Ising value _ = 2:055. Also, weobtain the critical Z2 temperatures for both models using the percolation criterion which agree very well with other calculations.Some preliminary calculations was done in the XY model in order to obtain the neutron scattering correlation function S(q; !) and the dispersion relation.
We study the phase diagram and critical behaviour of the two-dimensional square-lattice fully frustrated XY model (FFXY) and of two related models, a lattice discretization of the Landau–Ginzburg–Wilson Hamiltonian for the critical modes of the FFXY model, and a coupled Ising XY model. We present a finite-size-scaling analysis of the results of high-precision Monte Carlo simulations on L × L square lattices, up to L = O (10 3 ). In the FFXY model and in the other models, when the transitions are continuous, there are two very close but separate transitions. There is an Ising chiral transition characterized by the onset of chiral long-range order while spins remain paramagnetic. Then, as temperature decreases, the systems undergo a Kosterlitz–Thouless spin transition to a phase with quasi-long-range order. The FFXY model and the other models, in a rather large parameter region, show a crossover behaviour at the chiral and spin transitions that is universal to some extent. We conjecture that this universal behaviour is due to a multicritical point. The numerical data suggest that the relevant multicritical point is a zero-temperature transition. A possible candidate is the O(4) point that controls the low-temperature behaviour of the 4-vector model.
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.