The ground state of the square lattice bilayer quantum antiferromagnet with nearest (J1) and next-nearest (J2) neighbour intralayer interaction is studied by means of the dimer expansion method up to the 6-th order in the interlayer exchange coupling J3. The phase boundary between the spin-gap phase and the magnetically ordered phase is determined from the poles of the biased Padé approximants for the susceptibility and the inverse energy gap assuming the universality class of the 3-dimensional classical Heisenberg model. For weak frustration, the critical interlayer coupling decreases linearly with α(= J2/J1). The spin-gap phase persists down to J3 = 0 (single layer limit) for 0.45 < ∼ α < ∼ 0.65. The crossover of the short range order within the disordered phase is also discussed.KEYWORDS: bilayer Heisenberg antiferromagnet, frustration, Padé approximant, dimer expansion method, spingap state §1. IntroductionThe spin-1/2 square lattice Heisenberg model is now widely believed to have an antiferromagnetic long range order in the ground state.1, 2, 3, 4) It is, however, expected that the strong quantum fluctuation in this system may lead to the destruction of the long range order with the help of some additional mechanism. In this context, the square lattice antiferromagnetic Heisenberg model with nearest and nextnearest exchange interaction (hereafter called J 1 -J 2 model) 5,6,3,4,7,8,9,10,11,12,13,14,15,16,17,18) and the bilayer Heisenberg model 19,20,21,22,23,24,25,26,27,28) have been studied extensively. Considering the difference of the nature of the mechanism leading to the spin-gap phase in these two models, it must be most interesting to study their interplay in the bilayer J 1 -J 2 model.
30, 29)In the bilayer model, if the interlayer antiferromagnetic coupling is strong enough, the spins on both layers form interlayer singlet pairs and the quantum fluctuation is enhanced leading to the quantum disordered state. The dimer expansion study of this model has been quite successful 21,24,25,26,27) and it is shown that the transition between the Néel phase and the spin-gap phase belongs to the universality class of 3-dimensional classical Heisenberg model. This result is also confirmed by quantum Monte Carlo simulation.
23)On the other hand, in the J 1 -J 2 model, the competition between the nearest neighbour interaction J 1 and the nearest neighbour interaction J 2 introduces the frustration in the spin configuration which enhances the quantum fluctuation. The conclusion about the presence of the quantum disordered state in this model is, however, still controversial even in the most frustrated * e-mail: hida@riron.ged.saitama-u.ac.jp regime.In order to apply the dimer expansion method to the single layer J 1 -J 2 model, it is inevitable to start with the dimer configurations which break the translational symmetry as an unperturbed ground state.3, 4) In the bilayer J 1 -J 2 model, the unperturbed ground state can be taken as the interlayer dimers and the translational symmetry of the original Hamiltonian is ...