Understanding the morphotropic phase boundary of perovskite solid solutions as a frustrated state

Ying Shi Teh,Jiangyu Li,Kaushik Bhattacharya

DOI: https://doi.org/10.1103/PhysRevB.103.144201

2020-11-02

Abstract:Perovskite solid solutions that have a chemical composition A(C$_x$D$_{1-x})$O$_3$ with transition metals C and D substitutionally occupying the B site of a perovskite lattice are attractive in various applications for their dielectric, piezoelectric and other properties. A remarkable feature of these solid solutions is the \emph{morphotropic phase boundary} (MPB), the composition across which the crystal symmetry changes. Critically, it has long been observed that the dielectric and piezoelectric as well as the ability to pole a ceramic increases dramatically at the MPB. While this has motivated much study of perovskite MPBs, a number of important questions about the role of disorder remain unanswered. We address these questions using a new approach based on the random-field Ising model with long-range interactions that incorporates the basic elements of the physics at the meso-scale. We show that the MPB emerges naturally in this approach as a frustrated state where stability is exchanged between two well-defined phases. Specifically, long-range interactions suppress the disorder at compositions away from MPB but are unable to do so when there is an exchange of stability. Further, the approach also predicts a number of experimentally observed features like the fragmented domain patterns and superior ability to pole at the MPB. The insights from this model also suggest the possibility of entirely new materials with strong ferroelectric-ferromagnetic coupling using an MPB.

Mesoscale and Nanoscale Physics,Applied Physics

What problem does this paper attempt to address?

The problems that this paper attempts to solve are related to the understanding of the nature of the morphotropic phase boundary (MPB) in perovskite solid solutions and its related phenomena. Specifically, the article mainly focuses on the following issues:

1. **Effect of compositional disorder**: Why does the effect of compositional disorder emerge near the MPB, while this disorder is suppressed far from the MPB, thus forming a well - defined crystal structure?
2. **Crystal structure at MPB**: Does a definite crystal structure exist at the MPB?
3. **Change in domain structure**: Why does the domain structure near the MPB become highly fragmented?
4. **Polarization ability**: Is it easier to achieve the polarization of ceramics at the MPB due to compositional disorder?
5. **Possibility of new phenomena**: Can the MPB be used to create new material phenomena, such as strong ferroelectric - ferromagnetic coupling?

To solve these problems, the author introduced a new method based on the random - field Ising model (RFIM) combined with long - range interactions. This method can describe the physical processes in perovskite solid solutions at the mesoscopic scale and explain the above problems. Through simulation and theoretical analysis, the author showed that the MPB naturally emerges as a "frustrated state", in which two different phases exchange in stability. In addition, this model also predicted some experimentally observed phenomena, such as the fragmented domain structure and the superior polarization ability at the MPB.

### Formula summary

1. **Total Hamiltonian**:
\[
W_{\text{tot}}(\{p_i\};\{c_i\})=\sum_{i = 1}^N h_{\text{loc}}(p_i; c_i)-\frac{J_e}{2}\sum_{<i,j>}p_i\cdot p_j+D_e W_{\text{dip}}(\{p_i\})-E_{\text{ext}}\cdot\sum_{i}p_i
\]
where:
- \( h_{\text{loc}}(p_i; c_i)\) represents the local preference energy.
- \( J_e\) is the exchange energy.
- \( D_e W_{\text{dip}}(\{p_i\})\) is the long - range electric dipole - electric dipole interaction.
- \( E_{\text{ext}}\) is the external electric field.

2. **Electric dipole - electric dipole interaction**:
\[
W_{\text{dip}}(\{p_i\})=\frac{1}{(d - 1)}\sum_{i,j = 1}^{N}\sum_{R}\left[\frac{p_i\cdot p_j - d(p_i\cdot r_{ij})(p_j\cdot r_{ij})/r_{ij}^2}{r_{ij}^d}\right]+\frac{2\pi}{d}\sum_{i = 1}^{N}|p_i|^2
\]

3. **Transition probability**:
\[
P_s=\frac{\exp(-\beta W_{\text{tot}}^{(s)})}{\sum_{r = 1}^{N_{\text{states}}}\exp(-\beta W_{\text{tot}}^{(r)})}
\]
where \(\beta\) is the inverse temperature and \(N_{\text{states}}\) is the number of possible states.

Through these formulas and models, the author successfully explained various phenomena near the MPB and provided a theoretical basis for further exploration of new multiferroic materials.