Ab initio calculations for the tetragonal PbZr0.5Ti0.5O3

Ab initio studies of structural, elastic and electronic properties of the tetragonal perovskite-type PbZr0.5Ti0.5O3 are presented using the pseudo-potential plane wave method within the density functional theory in generalized gradient approximation. The calculated equilibrium lattice parameters remain in a good agreement with the available experimental data. The bulk modulus obtained from the Birch-Murnaghan equation of state is calculated as B0=170 GPa, and the gap energy Eg=2.1 eV-3.5 eV. The some differences between calculated and nominal charges exist for all atoms. The biggest ones are on the Pb ions. They are caused by hybridization of the Pb 6s and O 2p states. The influence of the strain on the averaged over directions Young modulus in the 0.1%-0.3% range was studied.

💡 Research Summary

This paper presents a comprehensive first‑principles investigation of the tetragonal phase of lead zirconate‑titanate with composition PbZr0.5Ti0.5O3 (PZT). Using density‑functional theory within the generalized‑gradient approximation (GGA‑PBE) as implemented in the SIESTA code, the authors model a 40‑atom supercell in which Zr and Ti atoms alternate on the B‑site to approximate the 50/50 solid solution. A plane‑wave cutoff of 400 Ry and a 2 × 2 × 2 Monkhorst‑Pack k‑point mesh are employed, together with norm‑conserving Troullier‑Martins pseudopotentials.

Structural optimization yields lattice parameters a = 4.120 Å and c = 4.122 Å, which are in excellent agreement with experimental values (4.042 Å and 4.127 Å, respectively). The close match validates the computational setup for this complex perovskite.

Elastic properties are derived from two complementary approaches. First, the total‑energy versus volume curve is fitted to a third‑order Birch‑Murnaghan equation of state, giving a bulk modulus B0 = 170 GPa. Second, the full stress tensor is computed for small imposed strains (±0.3 %). The nine independent elastic constants (σ11, σ22, σ33, σ44, σ55, σ66, σ12, σ23, σ13) reveal pronounced anisotropy: the c‑axis (σ33) is about 15 % stiffer than the a‑b plane (σ11, σ22), while shear constants are considerably smaller than uniaxial ones. The stress‑strain relationships are nonlinear, with noticeable kinks near 0.05 % and 0.20 % strain, suggesting incipient domain‑wall or phase‑transition activity. Young’s modulus, extracted from the elastic constants, shows the highest values along the c‑direction and the lowest in the a‑b plane. Hooke’s law holds only in a narrow strain window of 0.07 %–0.15 %; outside this range the response deviates markedly. The average Young’s modulus is approximately 108 GPa, comparable to experimental reports (80–120 GPa) but higher than measured values at room temperature, a discrepancy attributed to defects and thermal effects.

Electronic structure analysis yields a total density of states (DOS) with a direct band gap of 2.1 eV when the conduction band edge is taken at the first unoccupied state. If the occupancy threshold is set to zero, the gap expands to about 3.5 eV, aligning well with experimental optical gaps (~3.4 eV). Bader‑type charge integration over spherical regions around each atom shows substantial deviation from nominal ionic charges. Lead atoms carry an effective charge of +3.25 e (≈50 % larger than the nominal +2), a result of strong hybridization between Pb 6s and O 2p orbitals. Zirconium and titanium exhibit charges of +3.53 e and +3.43 e, respectively, while oxygen bears –1.24 e. This charge redistribution is consistent with observed large LO‑TO phonon splittings and with electron‑paramagnetic‑resonance evidence for Pb³⁺ species in PZT.

The authors discuss the systematic overestimation of elastic constants and bulk modulus relative to experiment. They attribute this to the idealized, defect‑free, 0 K crystal model, noting that even a modest 10 % defect concentration can halve the elastic constants. Temperature effects, not captured in the static DFT calculations, also soften the material. Moreover, GGA‑PBE is known to underestimate band gaps, which partly explains the lower calculated gap before occupancy correction.

In conclusion, the study provides a detailed set of structural, mechanical, and electronic parameters for tetragonal PbZr0.5Ti0.5O3, establishing a reliable first‑principles baseline for future work on defect engineering, temperature‑dependent behavior, and device‑level modeling of PZT‑based MEMS, sensors, and actuators.