Using High Multipolar Orders to Reconstruct the Sound Velocity in Piezoelectrics from Lattice Dynamics

Information over the phonon band structure is crucial to predicting many thermodynamic properties of materials, such as thermal transport coefficients. Highly accurate phonon dispersion curves can be, in principle, calculated in the framework of density-functional perturbation theory. However, well-established techniques can run into trouble (or even catastrophically fail) in the case of piezoelectric materials, where the acoustic branches hardly reproduce the physically correct sound velocity. Here we identify the culprit in the higher-order multipolar interactions between atoms and demonstrate an effective procedure that fixes the aforementioned issue. Our strategy drastically improves the predictive power of perturbative lattice-dynamical calculations in piezoelectric crystals and is directly implementable for high-throughput generation of materials databases.