Abstract
Monte Carlo simulation of polarization switching in ferroelectric thin films has been carried out by using a three-dimensional realistic kinetic Ising model with the next-nearest interactions. The realistic kinetic Ising model adopted in this paper is a single-spin-filp model. The key feature of this model is that, different from the previous kinetic Ising models such as the Glauber one, it is no longer arbitrary in the dynamic aspect. The effects of switching electric field on the switching time have been investigated, and agreement with experimental observations has been obtained.