Phase transitions in a two-dimensional weakly diluted 4-state Potts model on a hexagonal lattice

The Monte Carlo method is used to study phase transitions in a two-dimensional weakly diluted 4-state Potts model on a hexagonal lattice. Systems with linear dimensions $L$ = 21 $\div$ 336 at spin concentrations $p$ = 1.00, 0.90 are considered. The use of the fourth-order binder cumulants method and histogram data analysis showed that the introduction of nonmagnetic impurities has a stabilizing role in the implementation of the second-order phase transition in the Potts model under study on a hexagonal lattice.