JOURNAL OF PHYSICS-CONDENSED MATTER, vol.19, no.37, 2007 (SCI-Expanded)
A bilayer Ising model consisting of two Bethe lattices, each of which is coupled with crystal fields of different strengths and each with a branching ratio of q Ising spins with one of the layers having only spin-1 and the other having only spin-3/2, is laid over the top of the other and the two layers are tied together via an interaction between the vertically aligned spins. After obtaining the ground-state (GS) phase diagrams on different possible planes depending on the given system parameters, the changes in the order-parameters and the free energy are investigated by use of the exact recursion relations in a pairwise approach to calculate the phase diagrams of the model. The ferromagnetic ordering in each of the layers and ferromagnetic or antiferromagnetic ordering of the adjacent nearest-neighbor (NN) spins of the layers are considered. The system presents both second- and first-order phase transitions. The lines of the first-order phase transitions end on either the stable or unstable tricritical points or at the isolated critical points. The model also displays one or two compensation temperatures when the bilinear interaction of the upper layer with spin-1 can compete with that of the lower layers with spin-3/2.