Three-dimensional spin models are important for understanding real magnetic materials. In this work, we study the phase transition of the spin-half ( \(S=1/2\) ) antiferromagnetic quantum XXZ model with in-plane ferromagnetic interaction in three dimensions. The newly developed worm algorithm allows us to simulate a three-dimensional lattice system up to \(64 \times 64 \times 64\) . We simulated the spin system in the parametric space of temperature (T) and z-directional spin–spin interaction strength ( \(\Delta \) ). A tentative phase diagram is obtained by the finite-size scaling of superfluid stiffness and the analysis of compressibility. We also obtained the correlation-length critical exponent ( \(\nu \) ) with finite-size scaling theory and found \(\nu =0.65(3)\) , which is consistent with the classical model.