Stability and low-energy orientations of interphase boundaries in multiaxial ferroelectrics: Phase-field simulations

Yang Zhang,Fei Xue,Bo Wang,Jia-Mian Hu,Shuai Dong,Jun-Ming Liu,Long-Qing Chen
DOI: https://doi.org/10.1103/PhysRevB.105.014108
2022-01-25
Abstract:The coexistence of different ferroelectric phases enables the tunability of the macroscopic properties and extensive applications from piezoelectric transducers to nonvolatile memories. Here we develop a thermodynamic model to predict the stability and low-energy orientations of boundaries between different phases in ferroelectrics. Taking lead zirconate titanate and bismuth ferrite as two examples, we demonstrate that the low-energy orientations of interphase boundaries are largely determined by minimizing the electrostatic and elastic energies. Phase-field simulations are employed to analyze the competition between the interfacial energy and the electrostatic and elastic energies. Our simulation results demonstrate that the lowering of crystal symmetry could occur due to the electrical and mechanical incompatibilities between the two phases, which can be used to explain the experimentally observed low-symmetry phases near morphotropic phase boundaries. Our work provides theoretical foundations for understanding and controlling the interphase boundaries in ferroelectric materials for multifunctional applications.
Materials Science,Computational Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is related to the stability and low - energy orientation of interphase boundaries (IBs) in ferroelectric materials. Specifically, the author aims to predict the stability of interphase boundaries between different ferroelectric phases and their low - energy orientations through a thermodynamic model, and to explain the formation mechanism of the low - symmetry phase observed in experiments near the interphase boundaries. ### Specific problem description 1. **Stability of interphase boundaries in multiaxial ferroelectric materials**: - Due to the coexistence of multiple phases, ferroelectric materials can achieve macroscopic property regulation and wide - ranging applications, such as piezoelectric transducers and non - volatile memories. - The focus of the research is on how to predict and understand the stability and low - energy orientation of these interphase boundaries in multiaxial ferroelectric materials. 2. **Determinants of low - energy orientation**: - The author proposes that the low - energy orientation of interphase boundaries is mainly determined by minimizing electrostatic energy and elastic energy. - Through phase - field simulation, the competitive relationship between interface energy, electrostatic energy, and elastic energy was analyzed, and the influence of these energy terms on the stability of interphase boundaries was explored. 3. **Change in crystal symmetry**: - It was found that in cases of electrical and mechanical incompatibility, crystal symmetry may be reduced, which can be used to explain the appearance of the low - symmetry phase near the interphase boundaries observed in experiments. ### Main methods - **Thermodynamic model**: A thermodynamic model including long - range electrostatic energy and elastic energy was established. Electrical and mechanical incompatibility factors \(Q(n)\) and \(B(n)\) were introduced, and the total incompatibility factor \(E(n)=Q(n)+B(n)\) was defined. - **Phase - field simulation**: Phase - field simulation was used to calculate the energy changes of different interphase boundaries under different conditions, verifying the predictions of the theoretical model. ### Key formulas - **Electrostatic energy density**: \[ f_{\text{elec}}=\frac{1}{2}\varepsilon_0\kappa_{ij}E_iE_j - P_iE_i \] where \(\varepsilon_0\) is the vacuum permittivity, \(\kappa_{ij}\) is the background dielectric tensor, \(E_i\) is the electric field component, and \(P_i\) is the spontaneous polarization component. - **Elastic energy density**: \[ f_{\text{elas}}=\frac{1}{2}s_{ijkl}\sigma_{ij}\sigma_{kl}-\sigma_{ij}\varepsilon_{ij}^0 \] where \(s_{ijkl}\) is the fourth - order compliance tensor, \(\sigma_{ij}\) is the stress component, and \(\varepsilon_{ij}^0\) is the self - strain component. - **Electrical incompatibility factor**: \[ Q(n)=\left(\frac{\Delta P\cdot n}{\varepsilon_0\kappa_{nn}}\right)^2 \] - **Mechanical incompatibility factor**: \[ B(n)=c_{ijkl}n_in_k(\varepsilon_{ij}^0-\varepsilon_{ij}^{0'})n_jn_l \] Through these formulas and methods, the author successfully predicted and explained the stability and low - energy orientation of interphase boundaries in ferroelectric materials, providing a theoretical basis for understanding and controlling interphase boundaries in ferroelectric materials.