Abstract
The formation of zigzag domain boundary in ferroelectric-ferroelastic Gd2(MoO4)3, which cannot be interpreted by the simple symmetry and energy considerations, was numerically studied on the basis of the time-dependent Ginzburg-Landau model. We assumed the uniformity of lamellar domain structure in the initial stage and applied electric fields at two edges of the system along the polar c direction but with opposite sense. One of the major factors to form the zigzag-like boundary is the lamellar structure caused by the minimization of the anisotropic elastic energy. Another is bending deformation induced by the non-uniform electric field through piezoelectricity.