Piezoelectric constants for ZnO calculated using classical polarizable core-shell potentials

We demonstrate the feasibility of using classical atomistic simulations, i.e. molecular dynamics and molecular statics, to study the piezoelectric properties of ZnO using core-shell interatomic potentials. We accomplish this by reporting the piezoelectric constants for ZnO as calculated using two different classical interatomic core-shell potentials: that originally proposed by Binks and Grimes (1994 Solid State Commun. 89 921-4), and that proposed by Nyberg et al (1996 J. Phys. Chem. 100 9054-63). We demonstrate that the classical core-shell potentials are able to qualitatively reproduce the piezoelectric constants as compared to benchmark ab initio calculations. We further demonstrate that while the presence of the shell is required to capture the electron polarization effects that control the clamped ion part of the piezoelectric constant, the major shortcoming of the classical potentials is a significant underprediction of the clamped ion term as compared to previous ab initio results. However, the present results suggest that overall, these classical core-shell potentials are sufficiently accurate to be utilized for large scale atomistic simulations of the piezoelectric response of ZnO nanostructures.