This thesis focuses on studying the interaction and evolution of the needle domains in barium titanate single crystal through the simulation and experimental methods. The results are expected to assist in the interpretation of the arrangement and evolution of the needle domains and improve the design of future material. The existence of a needle domain gives rise to an internal stress field and electric field around the needle tip, which influences the polarization switching process and crack propagation in a single crystal. A model is established to study the interaction and evolution of the needle domains in barium titanate single crystals using the theory of dislocations. Considering the electrical and mechanical incompatibility at the needle tip, the fields produced by a needle domain are represented using the fields due to the equivalent edge dislocations and line charges distributed over the needle tip. Accordingly, the dislocation fields derived by Barnett and Lothe for anisotropic piezoelectric media are used to analyse the stress and electric fields around the needle domains. The calculation of the modified Peach-Koehler force and the total energy due to the needle domains is used to study the interaction among the needle domains and the stability of the needle pattern. Through experimental observation, we further know about the microstructure and common pattern of needle domains. X-ray diffraction is used to map the shape of needle tip and detect the strain field around the needle tip. The interaction of pairs of needle domains in an infinite piezoelectric body is studied. It is found that the needle tip interactions tend to be dominated by the electrostatic terms. Additionally, the stability of groups of needle domains is investigated. Stable configurations of needle domains in a herringbone pattern are identified, consistent with experimental evidence. However, comb-like arrays of needles are found to be unstable. This is contrary to experimental observations, where needle domains are observed in ...