Numerical simulations study of a spin-1 Blume-Emery-Griffiths model on a square lattice

L. M. Luque*, S. A. Grigera, E. Albano

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

4 Citations (Scopus)


We use the Monte Carlo simulation technique to study the critical behavior of a three-state spin model, with bilinear and biquadratic nearest-neighbor interactions, known as the Blume-Emery-Griffiths model (BEG), in a square lattice. In order to characterize this model, we study the phase diagram, in which we identify three different phases: ferromagnetic, paramagnetic and quadrupolar, the later with one sublattice filled with spins and the other with vacancies. We perform our studies by using two algorithms: Metropolis update (MU) and Wang-Landau (WL). The critical scaling behavior of the model is complementary studied by applying results obtained by using both algorithms, while tricritical points and the tricritical scaling behavior is analyzed by means of WL measuring the joint density of states and using the method of field mixing in conjunction with finite-size scaling. Furthermore, motivated by the decoupling between spins observed within the quadrupolar phase, we further generalize the BEG model in order to study the behavior of the system by adding a next nearest neighbour (NNN) interaction between spins. We found that by increasing the strength of the (ferromagnetic) NNN interaction, a new ferromagnetic phase takes over that contains both quadrupolar and ferromagnetic order.
Original languageEnglish
Article number033210
Number of pages17
JournalJournal of Statistical Mechanics: Theory and Experiment
Publication statusPublished - 18 Mar 2019


  • Classical Monte Carlo simulations
  • Classical phase transitions
  • Numerical simulations
  • Phase diagrams


Dive into the research topics of 'Numerical simulations study of a spin-1 Blume-Emery-Griffiths model on a square lattice'. Together they form a unique fingerprint.

Cite this