Enhancement of the electromechanical response in ferroelectric ceramics by design

It is demonstrated based on continuum mechanics modeling and simulation that it is possible to obtain polycrystalline ceramic ferroelectric materials which beggars single crystals in electromechanical properties. The local inhomogeneities at the ferroelectric domain-scale level due to spontaneous polarization and the underlying anisotropy are taken into consideration in the framework of mathematical homogenization of physical properties in ferroelectric materials. The intrinsic randomness of the spatial distribution of polarization is shown to be judiciously employed for the design of better polycrystalline ferroelectrics. The noncollinear rotation of the net polarization-vectors embedded in crystallites of the ceramic ferroelectrics is demonstrated to play the key role in the enhancement of physical properties.

💡 Research Summary

The paper presents a comprehensive theoretical and computational framework for designing ferroelectric ceramic polycrystals with electromechanical performance that surpasses that of single crystals. Using continuum mechanics and mathematical homogenization, the authors model a representative unit‑cell composed of many grains, each characterized by a net polarization vector P. The orientation of P in each grain is described by three Euler angles (φ, θ, ψ) and is treated statistically: an initially uniform distribution (unpoled state) evolves into a Gaussian distribution after electric poling, with mean m and standard deviation σ reflecting processing conditions such as poling field strength and annealing temperature.

The constitutive relations for a piezoelectric medium (stress‑strain‑electric field coupling) are expressed in tensor form and transformed from the crystal coordinate system to the global frame via Euler rotation tensors. Periodic boundary conditions are imposed on the unit‑cell, and a three‑dimensional finite‑element model (20 × 20 × 20 mesh, 64 000 integration points) is solved to obtain the homogenized stiffness, piezoelectric, and dielectric tensors. The homogenization procedure yields effective macroscopic coefficients that depend explicitly on the statistical distribution of grain orientations.

First, the authors validate the method on single‑crystal BaTiO₃ by rotating the crystal through the full Euler space. The computed shear coefficients d₁₅ and d₃₁ decrease monotonically with rotation, while the longitudinal coefficient d₃₃ exhibits a pronounced peak at θ ≈ 50°, reaching 223.7 pC/N—more than twice the value for a crystal poled along