Multicritical behavior of the ferromagnetic Blume-Emery-Griffiths model with repulsive biquadratic couplings

The ferromagnetic (J>0) version of the Blume-Emery-Griffiths model in the region of repulsive biquadratic couplings (K<0) is considered on a Cayley tree of coordination z, reducing the statistical problem to the analysis of a two-dimensional nonlinear discrete map. In order to investigate the effect of the coordination z on the system multicritical behavior, we study the particular case K/J=-3.5 with the inclusion of crystal fields (D≠0), but vanishing external magnetic fields (H=0), for two distinct lattice coordinations (z=4 and z=6). The thermodynamic solutions on the Bethe lattice (the central region of a large Cayley tree) are associated with the attractors of the two-dimensional map. The phase diagrams display several thermodynamic phases (paramagnetic, ferromagnetic, ferrimagnetic, and staggered quadrupolar). In some cases, there are regions of numerical costability of two different attractors of the map, associated with discontinuous phase transitions between the corresponding phases. To verify the thermodynamic stability of the phases and to locate the first-order boundaries, the analytical expression of the Gibbs free energy was obtained by the method proposed by Gujrati [Phys. Rev. Lett. 74, 809 (1995)PRLTAO0031-900710.1103/PhysRevLett.74.809]. For lower coordinations (z=4) the transition between the ferrimagnetic and the staggered quadrupolar phases is always continuous, while the transition between the ferromagnetic and the ferrimagnetic phases is discontinuous at low temperatures, turning into continuous for temperatures above a tricritical point. On the other hand, for higher coordinations (z=6), the transition between the ferromagnetic and the ferrimagnetic phases is always continuous. However, the transition between the ferrimagnetic and the staggered quadrupolar phases is continuous for higher temperatures and discontinuous for temperatures below a tricritical point, in agreement with previous results obtained in the mean-field approximation (infinity-coordination limit). In both cases, the occurrence and the thermodynamic stability of the ferrimagnetic phase is confirmed.