Application of elastostatic Green function tensor technique to electrostriction in cubic, hexagonal and orthorhombic crystals

J. Hlinka,E. Klotins

DOI: https://doi.org/10.1088/0953-8984/15/33/309

2002-06-18

Abstract:The elastostatic Green function tensor approach, which was recently used to treat electrostriction in numerical simulation of domain structure formation in cubic ferroelectrics, is reviewed and extended to the crystals of hexagonal and orthorhombic symmetry. The tensorial kernels appearing in the expressions for effective nonlocal interaction of electrostrictive origin are derived explicitly and their physical meaning is illustrated on simple examples. It is argued that the bilinear coupling between the polarization gradients and elastic strain should be systematically included in the Ginzburg-Landau free energy expansion of electrostrictive materials.

Materials Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to apply the elastic static Green's function tensor technique to the electrostrictive effect in materials with hexagonal and rhombic crystal symmetries. Specifically, the author aims to generalize a method that has been successfully used in the numerical simulation of the domain structure formation of cubic ferroelectrics to cover a wider range of crystal symmetries. This includes deriving the tensor kernels that appear in the expressions of the effective non - local interactions and explaining their physical meanings.

### Main problem decomposition

1. **Background and motivation**:
- The electrostrictive effect (electrostriction) plays a crucial role in ferroelectric materials, especially in the process of domain structure formation.
- The standard method eliminates the elastic degrees of freedom through the mechanical equilibrium condition, but the resulting effective energy terms are nonlinear and non - local.
- Using the Fourier transform can simplify the calculations in some cases, but it is necessary to deal with complex elastic Green's functions.

2. **Limitations of existing methods**:
- Existing methods are mainly applicable to crystals with cubic symmetry, and their applications to other symmetries (such as hexagonal and rhombic) have not been fully explored.
- For crystals with low symmetry, direct calculation of the elastic Green's function is very complicated.

3. **Research objectives**:
- Generalize the elastic static Green's function tensor technique to crystals with hexagonal and rhombic symmetries.
- Derive the effective interaction terms of the electrostrictive effect in these crystals.
- Explain the physical meaning of the tensor kernels and show their applications in different crystal symmetries.

### Formula summary

- **Total strain energy functional**:
\[
F = F_0\{P_i, P_{i,j}\} + F_1\{P_i, P_{i,j}, u_{ij}\} + F_2\{u_{ij}\}
\]
where \(F_0\) is the basic Landau free - energy density and gradient - energy density, \(F_1\) is the electrostrictive energy, and \(F_2\) is the pure elastic energy.

- **Electrostrictive stress**:
\[
f_{es} = -g_{ij}u_{ij}
\]
where \(g_{ij}\) is given by the electrostrictive tensor \(q_{ijkl}\):
\[
g_{ij} = q_{ijkl}P_k P_l
\]

- **Non - local effective interaction term**:
\[
F_{het} = -\frac{V}{2} \sum_{k \neq 0} \hat{n}_i g_{ij}(k) \Omega_{jk}(\hat{n}) g_{kl}(-k) \hat{n}_l
\]
or in a compact form:
\[
F_{het} = -\frac{V}{2} \sum_{k \neq 0} g_\alpha(k) B_{\alpha\beta}(\hat{n}) g_\beta(-k)
\]

- **Green's function tensor**:
\[
\Omega_{ij}(\hat{n}) = \delta_{ij} \left( \frac{1}{d_j} \right) - \frac{v_i v_j}{d_i d_j (1 + \sum_k \frac{v_k^2}{d_k})}
\]

- **Tensor kernel \(B_{\alpha\beta}\)**:
\[
B_{\alpha\beta}(\hat{n}) = \beta_{\alpha\beta} - \theta_\alpha \theta_\beta \left( 1 + \sum_k \frac{v_k^2}{d_k} \right)^{-1}
\]

### Conclusion

By generalizing the elastic static Green's function tensor technique, the author is able to deal with the electrostrictive effect in crystals with different symmetries. This method not only simplifies the calculations but also provides a new perspective for understanding the domain structure formation. Especially for rhombic crystals, this method is particularly important because these crystals are far from isotropic and...