Classical Monte Carlo algorithm for simulation of a pseudospin model for cuprates

A classical Monte Carlo algorithm based on the quasi-classical approximation is applied to the pseudospin Hamiltonian of the model cuprate. The model takes into account both local and non-local correlations, Heisenberg spin-exchange interaction, single-particle and correlated two-particle transfer. We define the state selection rule that gives both the uniform distribution of states in the phase space and the doped charge conservation. The simulation results show a qualitative agreement of a phase diagrams with the experimental ones.