High-resolution Monte Carlo study of the multicritical point in the three-dimensional XXZ Heisenberg antiferromagnet

We use Monte Carlo simulations to study the XXZ Heisenberg antiferromagnet in a field in order to clearly determine the nature of the multicritical point. We use a hybrid sampling method with Metropolis and Wolff-cluster algorithms, along with histogram reweighting techniques. Staggered magnetization susceptibilities, Binder cumulants, and finite-size scaling are considered in an effort to detect a possible biconical phase. An analysis of the probability distribution of the magnetization allowed us to conclude that the multicritical point is bicritical and it is in the three-dimensional Heisenberg universality class.