Simulation of thickness dependence in ferroelectric thin films

The thickness dependence of coercive field (EC) and remanent polarization (Pr) in ferroelectric thin films has been numerically simulated using a two-dimensional four-state Potts model. In this model, each of the dipoles in the film is assigned to one of the four states corresponding to the four different mutually perpendicular orientations. Neighboring dipoles with the same orientation are then grouped together to form a domain. Four different kinds of domains exist. In the presence of the surface layer near the electrode/film interface, the thickness dependence of both coercive field and remanent polarization are simulated.