Abstract 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.
A model for the critical relaxation in the one-dimensional Ising-type S=1 spin system is presented. This model is equivalent to the Blume-Emery-Griffiths model and exhibits two simple critical points and one tricritical point. The kinetic behavior is studied using the real-space renormalization-group approach. In the two critical points we find that the critical slowing down is described by the dynamic exponent z, z=2. In each point this exponent belongs to the critical order parameter, while the second order parameter relaxes faster, with z=1 or 0. At the critical point the two order parameters relax with the same z, z=1.
We investigated the dynamic phase transitions (DPTs) in the mixed spin (2, 5/2) Blume-Emery Griffiths model with repulsive biquadratic interaction in the presence of a time-varying magnetic field. We used the path probability method to obtain the set of the dynamic equations. We numerically solved these dynamic equations to characterize the nature of first- and second-order phase transitions and to find the DPT temperatures as well as obtain the phases in the system. We constructed the dynamic phase diagrams (DPDs) in reduced temperature and amplitude of oscillating magnetic field plane. We observed that the DPDs display richer, complex and more topological various type of phase diagrams. In particular, DPDs exhibit the disordered phase, antiquadrupolar or staggered phase, six different ferrimagnetic phases, three different nonmagnetic phases, and numerous mixed phases. DPDs also display two dynamic tricritical points for only smaller values of crystal-field interactions, multiple critical end and double critical end points, one zero-temperature critical point, one inverse critical end point, and a quadruple point depending on interaction parameters. The system always shows the reentrant behaviors for the higher values of magnetic field amplitude, but it does not exhibit the dynamic tricritical behavior for higher values of crystal-field parameter.
We study the phase diagram of the 3-matrix model. This model corresponds to the Blume-Emery-Griffiths (BEG) model on a random surface. The BEG model has parameters J and K which are the spin and the lattice-gas coupling respectively. If K = 0, it is the usual tri-critical Ising model. If K/J = 3, this model becomes the 3-state Potts model. These models are solved in the context of the Coulomb gas representation. The intermediate region K/J = 1 is equivalent to the 3-matrix chain model and we can solve it exactly using the orthogonal polynomial method. We find the 3rd order critical line on which the spin ordering phase transition occurs. This line ends at the 4th order critical point which corresponds to the tri-critical point of the BEG model. The phase diagram is quite similar to that of the BEG model on a regular lattice.
Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the classical three component XY model (i.e., the anisotropic limit of an easy-plane Heisenberg ferromagnet), on square lattices of size up to
1922, for several temperatures below, at, and above
TKT. The temporal evolution of spin configurations was determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by a hybrid Monte Carlo algorithm. The neutron-scattering function S(q,\ensuremath{\omega}) was calculated from the resultant space-time displaced spin-spin correlation function. Pronounced spin-wave peaks were found both in the in-plane and the out-of-plane scattering function over a wide range of temperatures. The in-plane scattering function
Sxx also has a large number of clear but weak additional peaks, which we interpret to come from two-spin-wave scattering. In addition, we observed a small central peak in
Sxx, even at temperatures well below the phase transition. We used dynamic finite-size scaling theory to extract the dynamic critical exponent z. We find z=1.00(4) for all T\ensuremath{\le}
TKT, in excellent agreement with theoretical predictions, although the shape of S(q,\ensuremath{\omega}) is not well described by current theory. \textcopyright{} 1996 The American Physical Society.
Using a renormalization-group transformation in position space based on the Migdal-Kadanoff recursion relations, we investigate the three-dimensional semi-infinite cubic spin-
23 Blume-Emery-Griffiths model with single-site uniaxial anisotropy and nearest-neighbor pair interactions, both bilinear and biquadratic. The fixed-point topologies and corresponding phase diagrams predicted for this model are presented. Attention is focused on the renormalization-group description of the surface critical behaviors, in comparison with the mean-field approximation results. Three main classes of system phase diagrams are determined, featuring a rich variety of phase transitions associated with the surface with type being according to the values of bulk and surface interaction parameters. Surface critical, tricritical, bicritical, and multicritical points are also found.