Critical energy distribution function of the Baxter-Wu model

In this work we examine the critical finite-size scaling behavior of the energy probability distribution function and its corresponding Binder cumulant at critical point. Based on the results of Monte Carlo simulations at zero external magnetic field using the recently developed triangle-cluster algorithm, we calculate the energy distribution function's exponents and we derive the scaling relations for the energy Binder cumulant. Finally, we predict the exact form of the scaled energy distribution function. Most of our conclusions seem applicable not only to the Baxter-Wu model but also to other Ising-like models.