Size-driven ferroelectrics in 2D pseudo-spin Ising model

We performed Monte Carlo simulation of the pseudo-spin Ising model in a two-dimensional ferroelectric system. The Metropolis method was employed by studying the size dependence of dielectric properties along the electric-field direction and the internal energy of the two-dimensional ferroelectric system. The critical phenomena in the vicinity of the phase transition temperature were considered by the scaling function obeying the values of critical exponents in view of the large-scale lattice.