Crystal field effect on a bilayer Bethe lattice

The influence of the crystal field on the phase diagrams of the bilayer spin-1 Ising model on the Bethe lattice is studied in terms of the intralayer coupling constants J(1) and J(2) of the two layers and interlayer coupling constant J(3) between the layers for given values of the coordination number q by using the recursion relation scheme. The six distinct ground-state configurations of the model are obtained on the (J(2)/J(1),J(3)/qJ(1)) plane with J(1)> 0, the ferromagnetic coupling, for given values of the crystal field. Then, the phase diagram of the system is obtained on the (kT/J(1),J(3)/J(1)) plane for given values of the crystal field and alpha=J(2)/J(1) with q=4 corresponding to the square lattice in real lattice systems. It was found that the system presents both first- and second-order phase transitions, therefore, tricritical points. The paramagnetic phase was also divided into two phases, P+ and P-, by studying the thermal behavior of the quadrupolar moments of the two layers.