Quantum Monte Carlo study of honeycomb antiferromagnets under a triaxial strain

The honeycomb antiferromagnet under a triaxial strain is studied using the quantum Monte Carlo simulation. The strain dimerizes the exchange couplings near the corners, thus destructs the antiferromagnetic order therein. The antiferromagnetic region is continuously reduced by the strain. For the same strain strength, the exact numerical results give a much smaller antiferromagnetic region than the linear spin-wave theory. We then study the strained XY antiferromagnet, where the magnon pseudomagnetic field behaves quite differently. The 0th Landau level appears in the middle of the spectrum, and the quantized energies above (below) it are proportional to n13(n23), which is in great contrast to the equally spaced ones in the Heisenberg case. Besides, we find the antiferromagnetic order of the XY model is much more robust to the dimerization than the Heisenberg one. The local susceptibility of the Heisenberg case is extracted by the numerical analytical continuation, and no sign of the pseudo-Landau levels is resolved. It is still not certain whether the result is due to the intrinsic problem of the numerical analytical continuation. Thus, the existence of the magnon pseudo-Landau levels in the spin-12 strained Heisenberg Hamiltonian remains an open question. Our results are closely related to the two-dimensional van der Waals quantum antiferromagnets, and may be realized experimentally.