The Pair Approximation method is applied to the antiferromagnetic Heisenberg–Ising spin- 1/2 bilayer with a simple cubic crystalline structure. The method allows for self-consistent calculations of thermodynamic quantities, based on the determination of Gibbs free-energy. In the paper the...