Frequency Dependence of the Ferroelectric-Hysteresis Phase-Diagram: Monte Carlo Investigation

Abstract

In this work, Monte Caro simulation was used to investigate the dynamic hysteresis behavior of 2 dimensional DIFFOUR model with varying field amplitude, field frequency and temperature over extensive ranges. Metropolis algorithm and Ising spins with variable magnitude were used in the simulation. The period average polarization was used as a dynamic order parameter in specifying dynamic ordered and disordered states. Fourth order cumulant of the period average polarization was then used to extract the dynamic critical temperature for constructing the dynamic phase diagram. It was found that the dynamic phase transition curves strongly depend on the field frequency. At higher frequency, the system requires higher level of field amplitude (electrical driving work) or higher environmental temperatures (thermal fluctuation) for a transition from dynamic ferroelectric to dynamic paraelectric phases. Consequently, the hysteresis shape at fixed temperature and field amplitude can be of either asymmetric or symmetric depending on the field frequency. This work therefore illustrates prominent and interesting behavior of the ferroelectric system under the perturbation of external field and temperature.