Crossover and universality breaking in the dilute Baxter-Wu model

The critical behavior of the Baxter-Wu model belongs to the universality class of the four-state Potts model. While the introduction of annealed vacancies does not alter the criticality of the four-state Potts model, the dilute Baxter-Wu model has remained the subject of several competing scenarios. Here we investigate the phase diagram of the spin-$1$ Baxter-Wu model in the presence of a crystal field using transfer-matrix calculations and large-scale Monte Carlo simulations. Our results reveal a systematic evolution of the effective critical behavior with increasing crystal field, accompanied by increasingly strong finite-size corrections near the crossover to the first-order regime. Along the line of continuous transitions, the central charge remains close to $c=1$, while the scaling dimensions systematically deviate from the spin-$1/2$ limit as the crystal field increases, consistent with either continuously varying effective critical exponents or a slow crossover between competing critical behaviors. The first-order regime is independently characterized through multicanonical simulations, which confirm the expected finite-size scaling and interfacial behavior. Taken together, our results provide a unified picture of the dilute spin-$1$ Baxter-Wu model, substantially narrowing the range of possible scenarios for the crossover between continuous and first-order phase transitions.