Ab Initio Piezoelectric Properties Of Wurtzite Zno-Based Alloys: Impact Of The C/A Cell Ratio

The piezoelectric (PE) and stiffness tensors of 32-atom supercells of ZnO-based alloys have been obtained from ab initio simulations using density functional perturbation theory in the local density approximation. Low concentration for substituents to Zn, O, or both were considered in unstrained and biaxially strained supercells. The d(33) coefficient for unstrained Zn15YO15N and Zn15LaO15N alloys are, respectively, 17.5 and 18 pC N-1, whereas e(33) is 1.7 C m(-2) for both alloys. These values are significantly improved compared to simulated values for pristine ZnO (d(33) = 11.4 pC N-1 and e(33) = 1.3 C m(-2)). Applying 2% tensile strain on Zn15YO15N results in an increase of the e(33) coefficient to 2.1 C m(-2), a 62% increase over the value calculated for pristine ZnO. We confirm for a variety of ternary and quaternary ZnO-based alloys that a linear relation is verified between the e(33) coefficient and the cell ratio c/a, described by a slope approximate to -9 C m(-2). Our results also indicate that the PE coefficients follow the same trends with respect to changes in c/a caused by variations in chemical composition or by applying biaxial strain. Based on this correlation, we propose a simple method to identify promising candidates among piezoelectric alloys in the wurtzite family, effectively reducing the intensive computational resources needed to obtain optimal PE performance for applications compatible with the many requirements of thin film growth and processing.