Domain structures play a key role in determining the physical properties of ferroelectric materials. The formation of these ferroelectric domains and domain walls are determined by the intrinsic nonlinearity and the nonlocal coupling of the polarization. Analogous to soliton excitations, domain walls can have high mobility when the domain wall energy is high. The domain wall can be describes by a continuum theory owning to the long range nature of the dipole-dipole interactions in ferroelectrics. The simplest form for the Landau energy is the so called ϕ model which can be used to describe a second order phase transition from a cubic prototype,where Pi (i =1, 2, 3) are the components of polarization vector, α's are the linear and nonlinear dielectric constants. In order to take into account the nonlocal coupling, a gradient energy should be included, for cubic symmetry the gradient energy is given by,
