Transição de fase e universalidade dos modelos Blume-Capel e Baxter-Wu bidimensionais através das distribuições de probabilidade da energia e magnetizações

Abstract

The Blume-Capel model on a square lattice and the Baxter-Wu model on a triangular lattice have been studied through extensive Monte Carlo simulations using the single spin update according to the Metropolis algorithm. Both models have been treated with spin $S=1$ and spin $S=3/2$ in the presence of crystal field interactions. The transition temperature and correlation length critical exponent have been computed, as a function of the crystal field interaction, by employing the recent developed method of the zeros of the energy probability distribution. The important question of the universality class along the second-order transition lines could be resolved by looking at the corresponding universal probability distributions of the energy and magnetizations. Mixing field arguments turned out to be quite important to decide the possible change in the universality class of the models as the crystal field varies. Single histograms techniques have also been utilized in order to obtain simulational data close to the second-order transitions and multicritical points. As expected, the Blume-Capel model has first- and second-order transitions, where a tricritical point occurs for spin $S=1$, while an isolated double-critical-endpoint is present for spin $S=3/2$. Although the topology of the phase diagram of the Baxte-Wu model for spin $S=1$ is similar to that of the Blume-Capel model, with the tricritical point being replaced by a pentacritical point, the Baxter-Wu model has critical exponents that change along the tetracritical line as the crystal field is increased. On the other hand, the Baxter-Wu model with spin $S=3/2$ has a phase diagram that is completely different from the one obtained for the Blume-Capel, with the presence of a pentacritical point and a tetracritical endpoint. The present data also show that, despite the change of the probability distributions as the crystal field varies, the models with different spin values belong to the same universality class.