The antiferromagnetic Ising model for a bilayer Bethe lattice

The totally antiferromagnetic Ising model is analyzed on a bilayer Bethe lattice in detail by studying the order-parameters, response functions, i.e. susceptibility and specific heat, and free energy by using the recursion relations in a pairwise approach. The ground state phase diagrams of the model are also obtained on the (_J2/|_J1|,_J3/q|_J1|)$(J_2/|J_1|,J_3/q|J_1|)$ plane for given values of H/q|_J1|$H/q\left|J_1\right|$ and on the (H/q|_J1|,_J3/q|_J1|)$(H/q|J_1|,J_3/q|J_1|)$ plane for given _J2/|_J1|$J_2/\left|J_1\right|$. As a result, we have obtained the temperature-dependent phase diagrams for various values of the coordination number q   on the (_J3/|_J1|,kT/|_J1|)$(J_3/|J_1|,\mathit{kT}/|J_1|)$ and (H/|_J1|,kT/|_J1|)$(H/|J_1|,\mathit{kT}/|J_1|)$ planes for given values of the rest of the system parameters.