Revisit the phase diagram and piezoelectricity of lead zirconate titanate from first principles

Lead zirconate titanate (PbZr1-xTixO3, PZT) exhibits excellent piezoelectric properties in the morphotropic phase boundary (MPB) region of its temperature-composition phase diagram. However, the microscopic origin of its high piezoelectric response remains controversial. Here, we develop a machine-learning-based deep potential (DP) model of PZT using the training dataset from first principles density functional theory calculations. Based on DP-assisted large-scale atomic simulations, we reproduce the temperature-composition phase diagram of PZT, in good agreement with the experiment except the absence of structural transition from R3c to R3m. We find that the rhombohedral phase maintains R3c symmetry with slight oxygen octahedral tilting as increase of temperature, instead of appearing R3m symmetry. This discrepancy can trace back to the lack of experimental measurements to identify such slight octahedral tilting. More importantly, we clarify the atomic-level feature of PZT at the MPB, exhibiting the competing coupling of ferroelectric nanodomains with various polarization orientations. The high piezoelectric response is driven by polarization rotation of nanodomains induced by an external electric field.