Critical properties in the two-dimensional XY model with four-fold and eight-fold crystal fields

We numerically investigate the critical properties of the two-dimensional XY model with four-fold and eight-fold crystal fields. Using extensive Monte Carlo simulations, we compute the specific heat, magnetization, orientational magnetization, and correlation length ratio to examine how these anisotropies influence the nature of phase transitions. The critical temperature and the critical exponent are systematically mapped across both compatible and competing field regimes, spanning from weak to strong field strengths. The numerical results reveal that the critical exponent exhibits a linear dependence on the inverse of the four-fold field strength in the strong field regime, while remaining finite and non-divergent in the weak field regime. Furthermore, the orientational magnetization effectively identifies and characterizes Ising-type phase transitions in the competing field regime.