“…Nevertheless, due to the bulky implement of these methods, the simulating supercells are generally small and the simulating time is limited. Recently, phase-field methods capable of simulating large-size systems and long-time period have been applied to simulate the mechanical behaviors of ferroelectric domain patterns (Li et al, 2002a;Wang et al, 2004;Su and Landis, 2007;Kontsos and Landis, 2010;Chen et al, 2012aChen et al, , 2012bWang et al, 2013). Although many works have been paid attention to two-dimensional geometries, the three-dimensional phase-field simulations (Li et al, 2002a;Wang et al, 2013) can capture the whole domain structure and complicated evolution of the domain pattern in ferroelectric thin films, with properly considering the effects of temperature, surfaces, external electric field and mechanical loads, etc.…”