The metastable phase diagram of the Blume-Emery-Griffiths model in addition to the equilibrium phase diagram using the pair approximation

As a continuation of our previously published work, the metastable phase diagram of the Blume-Emery-Griffiths model with the arbitrary bilinear (J), biquadratic (K) and crystal field interaction (D) is presented in addition to the equilibrium phase diagram in (T/K, J/K) and (T/K, D/K) plane by using the pair approximation of the cluster variation method on a body centered cubic lattice. We also calculate the phase transitions for the unstable branches of order parameters. The calculated first- and second-order phase boundaries of the unstable branches of the order parameters are superimposed on the equilibrium phase diagram and metastable phase diagram. It is found that the metastable phase diagram and the first- and second-order phase boundaries for unstable branches of order parameters always exist at low temperatures, which are consistent with the experimental and theoretical works.