SURFACE SCIENCE, vol.726, 2022 (SCI-Expanded)
We investigated the order-disorder phase transitions with the four spin interacting Ising model Hamiltonian, which was obtained using first principle calculation results for the Si(100) surface. This is the first proposed effective spin Hamiltonian for dimer configurations. However, the four spin interacting Hamiltonian has not been established to date. All simulations were performed using the cellular automaton algorithm adapted for this Hamiltonian. As a result of the simulations made for three different parameter sets, the order-disorder phase transition temperature on the clean Si(001) surface was obtained on average at 203 +/- 4 K, which is consistent with the experimental results.