Tuning the piezoelectric and mechanical properties of the AlN system via alloying with YN and BN
Sukriti Manna, Geoff L. Brennecka, Vladan Stevanovi?, Cristian V. Ciobanu
11 Tuning the piezoelectric and mechanical properties of the AlN system via alloying with YN and BN
Sukriti Manna , Geoff L. Brennecka , Vladan Stevanovi ć and Cristian V. Ciobanu Dept. of Mechanical Engineering, Colorado School of Mines, Golden, Colorado 80401 Dept. of Metallurgical and Materials Engineering, Colorado School of Mines, Golden, Colorado 80401
Abstract:
Recent advances in microelectromechanical systems often require multifunctional materials, which are designed so as to optimize more than one property. Using density functional theory calculations for alloyed nitride systems, we illustrate how co-alloying a piezoelectric material (AlN) with different nitrides helps tune both its piezoelectric and mechanical properties simultaneously. Wurtzite AlN-YN alloys display increased piezoelectric response with YN concentration, accompanied by mechanical softening along the crystallographic c direction. Both effects increase the electromechanical coupling coefficients relevant for transducers and actuators. Resonator applications, however, require superior stiffness, thus leading to the need to decouple the increased piezoelectric response from a softened lattice. We show that co-alloying of AlN with YN and BN results in improved elastic properties while retaining some of the piezoelectric enhancements from YN alloying. This finding may lead to new avenues for tuning the design properties of piezoelectrics through composition-property maps. Keywords: piezoelectricity, electromechanical coupling, density functional theory, co-alloying ________________________________________ * To whom correspondence may be addressed, electronic mail: [email protected]
INTRODUCTION
Recent technological advances in microelectromechanical systems require materials to be more efficiently designed, often leading to optimization of more than one property.
1, 2
An electromechanical (piezoelectric) material requires the best piezoelectric coefficient, while materials for damping high mechanical loads require the high stiffness moduli.
3, 4
The high stiffness value is also associated with thermal stability, enabling the use of these materials over a large temperature range. Pseudo-binary alloys, i.e. solid solutions between two compounds, have been extensively investigated for use in electromechanical applications.
Most of these alloys have large miscibility gaps, which makes them less appealing for applications where mechanical stability is an important figure of merit. The presence of miscibility gaps is not particularly restrictive, since current non-equilibrium processing techniques allow for fabricating pseudo-binary alloys in which one can seek a given functionality by, e.g., changing composition (often through the miscibility gaps) so as to optimize piezoelectric coefficients, stiffness, electromechanical coupling constants, or electronic properties. However, achieving optimal or desired values for more than one such property is challenging since optimizing one property often changes another towards undesirable values. In such cases, more degrees of freedom may help create a framework in which advanced materials can be optimized for multiple functionalities. In this article, we focus on a piezoelectric system to show that alloying with more than one compound enables us to tune two properties, specifically, the piezoelectric coefficient and the mechanical stiffness. With the discovery of the anomalously large increase in the piezoelectric moduli of wurtzite AlN when alloyed with ScN or YN,
9, 12
AlN-based alloys have received a lot of attention as materials for improved telecommunications, sensors, and surface acoustic wave devices.
13, 14
The origin of the enhanced piezoelectricity was revealed to be an intrinsic effect of alloying, as opposed to textural or microstructural effects, with the obvious implication that controlling the piezoelectric enhancement can be done mainly by increasing composition of YN or ScN. However, such alloying softens the material significantly, which actually makes the material less attractive for resonant applications for which k Q is a common figure of merit wherein k is an electromechanical coupling coefficient and Q is a quality factor that is inversely proportional to mechanical dissipation. It is therefore desirable to control the mechanical properties as well, which could in principle be achieved by alloying with other species. Herein, we focus on AlN-YN alloys and advance the premise that further alloying this system with BN can increase the elastic constants of the material, thus leading to new avenues for piezoelectric design based on precise control YN and BN compositions. The choice for YN is motivated by the scarcity of literature reports on AlN-YN alloys,
12, 16 and by the fact that Y may be a less expensive alternative for Sc. While the elastic stiffening of AlN-YN upon alloying with BN would be enabled by the presence of shorter and stiffer B − N bonds, the effect of alloying with BN on the piezoelectric properties is presently unknown. Using density functional theory (DFT), we have found that the presence of BN, by itself, indeed leads to increases in the relevant elastic constant C , but does not change the piezoelectric coefficient e significantly. Furthermore, we have uncovered the combined effect of alloying with YN and BN on the properties of the AlN system for total alloying concentrations up to 50 at. %. The addition of BN to AlN-YN alloys counters to some extent the piezoelectric enhancement obtained via YN alloying, but increases the elastic modulus. Overall, the alloying of AlN with BN and YN leads to a high degree of control over the piezoelectric and elastic properties of the resulting alloy, and therefore over the coupling coefficients and various figures of merit commonly used in piezoelectric device design. METHODS
To assess the changes in elastic and piezoelectric properties driven by different levels of YN and BN in Al x - y Y x B y N alloys, we first performed DFT relaxations of atomic positions and lattice constants using the VASP package, with projector augmented waves in the generalized gradient approximation using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. Out of caution more than necessity, we introduced an on-site Coulomb interaction (parameter
𝑈 = 3 eV) for the yttrium d state since this state is explicit in the pseudopotential. When compared with calculations done with plain PBE ( 𝑈 = 0 eV), the results for lattice constants are nearly identical, while the piezo and stiffness tensor component vary by less than 2%. The parent structure is wurtzite AlN, in which Al is substituted with Y and B in desired amounts. Supercells with 32 atoms are simulated using special quasirandom structures (SQS) (Fig. 1), whose pair correlation functions are the same as those of random alloys up to second-nearest Figure 1.
Atomic structure of a representative SQS Al Y B N alloy structure. neighbors. A 5 × ×
21, 22 k -point mesh was used to sample the Brillouin zone, resulting in 38 irreducible k -points. A plane wave cutoff energy of 500 eV was used; the energy tolerance stopping criteria were 10 -8 and 10 -7 for the electronic and ionic convergence, respectively. These criteria were obtained by progressive tightening of the convergence until the relaxed SQS supercells stopped changing. This resulted in smoother curves for lattice parameters and properties, in which the only source of noise remained the use of small number of atoms coupled with a limited number of individual SQS supercells at each concentration (refer to Fig. SM-1 in Supplementary Material for an assessment of the dispersion of individual SQS values). Due to the variability of atomic-scale environments in the SQS supercells, we performed five calculations at each composition and used them to average the properties, including the lattice constants. For the separate alloying only with YN or with BN, the lattice constants a and c are plotted Figure 2. (a,b) Lattice constants a and c for AlN-YN and AlN-BN alloys as a function of dopant (Y or B) cation concentration x . Panel (c) shows the 𝑐/𝑎 ratio for these alloys. in Fig. 2(a,b) as functions of the cation concentration. Alloying with YN increases the lattice constants, while alloying with BN alone decreases them: these trends are simply due to the larger (smaller) cation radius of Y (B) compared to that of Al in AlN. The 𝑐/𝑎 ratio is more sensitive to doping for the case of Y compared to BN [Fig. 2(c)]. Our DFT calculations of mixing enthalpies for alloy wurtzite and rocksalt phases of AlN-YN revealed that wurtzite is more stable than rocksalt for up to ~72% Y content, above which rocksalt becomes more stable [Fig. 3(a)]. The mixing enthalpies in Fig. 3(a) are positive, meaning that mixing is not spontaneous. From the mixing enthalpies and combinatorial expression for mixing enthalpy in this pseudo-binary alloy, we have obtained the temperature-composition phase diagram in a manner similar to that described by, e.g., Burton et al. for other wurtzite systems. This phase diagram is shown in Fig. 3b, in which the temperature axis is rather qualitative since we have not considered lattice vibrations in deriving it; it is similar to the one that would be obtained in the presence of vibrations. As would be expected, the diagram shows a miscibility gap – the white space between the metastable regions in Fig. 3b. This means that compositions in that gap cannot be achieved under thermodynamic equilibrium. However, the current growth techniques (e.g. reactive dc magnetron co-sputtering) rely on non-equilibrium processes to
Figure 3.
Calculated mixing enthalpy curves for AlN alloyed with YN in the wurtzite ( w ) and rocksalt ( rs ) structures. overcome the solubility limits and create alloys with better control over properties through control of composition. RESULTS AND DISCUSSION
Focusing now on describing the piezoelectric properties of the alloys, we note that the piezoelectric coefficient e is evaluated at each Y ( x ) and B ( y ) composition via 𝑒 )) = 𝑒 ))+,-./012345 + ∗ ) - :; (1) where 𝑒 ))+,-./012345 is the clamped ion contribution to e , e is the electron charge, Z * is the axial component of the dynamical Born effective charge tensor, a is the equilibrium lattice constant, and 𝑑𝑢/𝑑𝜀 Figure 4.
The design premise behind co-alloying AlN with YN and BN. (a) Higher YN content in Al x Y x N increases the piezoelectric coefficient e but decreases the elastic modulus C . (b) Increasing concentration of BN in Al y B y N has a small effect on the piezoelectric coefficient, but increases the elastic modulus. Thus, co-alloying with YN and BN could lead to increases in both piezoelectric coefficient and stiffness relative to pure AlN. is the sensitivity of the internal parameter u to strains along the crystallographic c -axis.
23, 24
The elastic constants were computed by finite differences between the energies of the lattices subjected to small applied strains. Figure 4 shows the effects of separate BN and YN alloying on the properties of AlN. As seen in Fig. 4(a) the piezoelectric coefficient e increases with Y concentration up to 50%, while the elastic constant C decreases; this is consistent with previous works on the similar AlN-ScN system.
6, 8
The physical origin of this effect arises from the larger radius of the Y ion, which changes the internal parameter u of the alloyed wurtzite structure and its sensitivity to strain, 𝑑𝑢/𝑑𝜀 [refer to Eq. (1)]. The value of e that can be achieved by alloying with 50% Y is doubled compared to that of pure AlN (Fig. 4). Figure 4(b) shows that alloying with BN alone increases the elastic modulus, which is consistent to another report. Interestingly, we also note that substituting B for Al does not alter the piezoelectric coefficient significantly up to concentrations of 50% B.
When AlN-YN films are deposited, there is usually strain in the lattice stemming from the mismatch with the substrate, changing film texture with alloy concentration, or from insufficient annealing. It is therefore relevant to assess what the effect of this strain is on the piezoelectric coefficient, so we have re-done the DFT calculations for the situations were strain is applied in plane (along the a and b directions) or out of plane (along the c direction). The results of these calculations are displayed in Fig. 5: e can increase up to about 3.75 C/m (i.e., 2.5 times that of pure AlN) either by alloying with 50% Y and applying 2% tension along the a axis [Fig. 5(a)], or by alloying with 50% Y and applying 2% compression along the c axis [Fig. 5(b)]. The opposite trends of e with strain along different directions (Fig. 5) suggest that, for homogeneous and isotropic strain in the AlN-YN lattice, the effects of strain along c and a directions may cancel each other at least partially and thus not affect the piezoelectric coefficient significantly. This is Figure 5.
Variation of piezoelectric coefficient in Al x Y x N with strained applied along the a and c axes. important when alloying with a secondary species, such as BN, which would amount to roughly isotropic strain when B is (nearly) homogenously doped throughout the AlN-YN lattice. Based on these results, we have set out to analyze how co-alloying with both YN and BN changes properties, aiming to find composition regimes where the piezoelectric coefficient increases while the elastic modulus does not decrease to the same extent as in the case of alloying with Y alone. Figure 6 shows Figure 6.
The key piezoelectric ( e , red, left vertical axis) and elastic ( C , blue, right axis) properties of Al-Y-B nitrides as functions of the Y concentration at constant boron levels ranging from (a) 0.00% B to (h) 43.75% B. the variation of e and C with the Y concentration x , at a fixed level of B doping, y (such that x + y ≤ e and C ) are shown as functions of x and y in Fig. SM-2. Our starting premise that the presence of boron can mitigate or control the decrease in elastic modulus C brought about by Y alloying – while still retaining some increase of the piezoelectric coefficient – is borne out by the results in Fig. 6. For example, going from 0% B in Al Y N to 18.75% B concentration in Al Y B N results in a ~50 MPa increase of the elastic constant C with ~33% decrease in the piezoelectric coefficient. For ease of reference, some of the results shown in Fig. 6 are also tabulated in Table I. BN alloying has a strong effect on mechanical properties (Table I): increasing BN content from 31.25% to 50% at the expense of YN concentration nearly doubles the elastic modulus (from 238 GPa to 470 GPa). Table I.
Comparison between the properties of some of the Al-Y-B nitrides computed in this work and those of three commercial piezoelectrics.
Material e (C/m ) C (GPa) d (pC/N) k PZT-7A, Ref Na NbO , Ref.
33, 34 this work:
AlN 1.46 353 5 0.06 Al Y N 3.02 158 38 0.37 Al Y B N 2.51 161 22 0.29 Al Y B N 1.70 180 13 0.14 Al Y B N 1.85 209 13 0.14 Al Y B N 1.58 223 13 0.11 Al Y B N 1.93 238 13 0.14 Al Y B N 1.54 253 9 0.12 Al Y B N 1.73 328 8 0.08 Al B N 1.34 470 3 0.04 Next, we discuss the electromechanical coupling coefficients resulting from co-alloying of AlN with YN and BN. These coefficients are key parameters in designing transducers for energy harvesting and sensing applications.
5, 35
Applications such as pressure sensors, accelerometers, and gyroscopes commonly require the piezoelectric medium to operate in the longitudinal length mode, for which the relevant coupling coefficient is 𝑘 ))B = 𝑒 ))B 𝜀 )) 𝐶 )) + 𝑒 ))B where ε is the 33 component of the dielectric tensor. The static dielectric constants vary by about 10 -15 % from their value at x =12.5% (separate calculation) for the entire range of equimolar Y-B doping (i.e., x ≤
25% in Al x Y x B x N), so a constant value of ε = 11.0 was assumed in computing the coupling coefficients. High k coupling factor corresponds to transducers with better axial resolution, broader bandwidth, and higher sensitivity. Actuators based on cantilevers operate in the transverse length mode (bending),
35, 38 for which the relevant coupling coefficient is 𝑘 )DB = 𝑒 )DB 𝜀 )) 𝐶 DD + 𝑒 )DB , Figure 7.
The coupling constants (a) k and (b) k as functions of the Y and B concentrations in Al x - y Y x B y N with at least 50% Al. where e and C are components of the piezoelectric and elastic tensor, respectively. The two coupling coefficients are plotted in Fig. 7 (a,b), for the range of Y and B concentrations considered. High coupling constants require large piezoelectric coefficients and low elastic moduli. Increasing the amount of Y in the system increases the piezoelectric and coupling coefficients but softens the system elastically, which is true for both coupling modes addressed here. The two coupling coefficients exhibit similar trends (Fig. 7), with larger values achieved for higher Y concentrations. However, variations of k are significantly smaller, and larger values can also be achieved using equiatomic Y:B alloys. With the piezoelectric and elastic tensor components obtained from DFT calculations, it is useful to compare some of the properties of the co-alloyed AlN system with those of known piezoelectrics. Table I, which also includes the longitudinal piezoelectric modulus d , shows the properties of three well-known piezoelectrics in comparison with a few selected alloy compositions computed here (for completeness, piezoelectric moduli are also plotted in Fig. SM-3 for the entire range of compositions used). Direct inspection shows that AlN co-alloyed with Y and B can achieve parameters comparable to those of commercial piezoelectrics. More importantly, the simultaneous enhancement of e and C in comparison to the pure AlN values can be seen for boron levels of up 25%, illustrating that co-alloying with different species represents a valuable route for simultaneously engineering the piezoelectric and the elastic properties of a parent wurtzite structure such as that of AlN. CONCLUDING REMARKS
In conclusion, we have shown that co-alloying of AlN with YN and BN leads to superior control of both piezoelectric and mechanical properties. The idea of introducing BN into the systems arose as a way to stiffen the material, which otherwise would have becomes increasingly softer with increased YN content. In the process of assessing the influence of BN and YN composition on properties, we have separately calculated the effect of straining the lattice in different directions, and found that enhancements of the piezoelectric coefficient can be made by applied strain as well: this strain effect could become important for certain substrates, or for electromechanical applications where strain develops in the structure during operation. The results presented here could carry over to other parent phases and alloying phases as well, provided that each alloying agent offers a clear advantage in terms of engineering one specific property. In our case, YN leads to increasing the piezoelectric coefficient and BN leads to increasing the stiffness. The use of both YN and BN offers intrinsic and simultaneous control over the mechanical and piezoelectric properties of the Al-Y-B-N system, even though the maximal value for each individual parameter can be achieved in the presence of only one dopant. Acknowledgment.
The authors gratefully acknowledge the support of the National Science Foundation through Grant No. DMREF-1534503. The calculations were performed using the high-performance computing facilities of the Golden Energy Computing Organization and the National Renewable Energy Laboratory.
References S.-H. Park, H. B. Lee, S. M. Yeon, J. Park, and N. K. Lee, ACS Applied Materials & Interfaces , 24773 (2016). C. Pang, G.-Y. Lee, T.-i. Kim, S. M. Kim, H. N. Kim, S.-H. Ahn, and K.-Y. Suh, Nature Materials , 795 (2012). J. A. Rogers, T. Someya, and Y. Huang, Science , 1603 (2010). K. Uchino, Acta Materialia , 3745 (1998). S. Trolier-McKinstry and P. Muralt, Journal of Electroceramics , 7 (2004). F. Tasnádi, B. Alling, C. Höglund, G. Wingqvist, J. Birch, L. Hultman, and I. A. Abrikosov, Physical Review Letters , 137601 (2010). S.-B. Kim, H. Park, S.-H. Kim, H. C. Wikle, J.-H. Park, and D.-J. Kim, Journal of Microelectromechanical Systems , 26 (2013). M. A. Caro, S. Zhang, T. Riekkinen, M. Ylilammi, M. A. Moram, O. Lopez-Acevedo, J. Molarius, and T. Laurila, Journal of Physics: Condensed Matter , 245901 (2015). P. Mayrhofer, H. Riedl, H. Euchner, M. Stöger-Pollach, P. Mayrhofer, A. Bittner, and U. Schmid, Acta Materialia , 81 (2015). C. Tholander, J. Birch, F. Tasnádi, L. Hultman, J. Palisaitis, P. Å. Persson, J. Jensen, P. Sandström, B. Alling, and A. Žukauskaitė, Acta Materialia , 199 (2016). M. Akiyama, T. Kamohara, K. Kano, A. Teshigahara, Y. Takeuchi, and N. Kawahara, Advanced Materials , 593 (2009). A. Zukauskaite, C. Tholander, J. Palisaitis, P. O. A. Persson, V. Darakchieva, N. Ben Sedrine, F. Tasnadi, B. Alling, J. Birch, and L. Hultman, Journal of Physics D-Applied Physics (2012). G. Piazza, V. Felmetsger, P. Muralt, R. H. Olsson III, and R. Ruby, MRS Bulletin , 1051 (2012). Y. Zhang, W. Zhu, D. Zhou, Y. Yang, and C. Yang, Journal of Materials Science: Materials in Electronics , 472 (2015). R. R. Syms, H. Zou, and P. Boyle, Journal of Microelectromechanical Systems , 529 (2005). C. Tholander, I. A. Abrikosov, L. Hultman, and F. Tasnadi, Physical Review B , 094107 (2013). G. Kresse and J. Furthmüller, Computational Materials Science , 15 (1996). J. P. Perdew, K. Burke, and M. Ernzerhof, Physical Review Letters , 3865 (1996). G. Kresse and D. Joubert, Physical Review B , 1758 (1999). A. Van de Walle, P. Tiwary, M. De Jong, D. Olmsted, M. Asta, A. Dick, D. Shin, Y. Wang, L.-Q. Chen, and Z.-K. Liu, Calphad , 13 (2013). H. J. Monkhorst and J. D. Pack, Physical Review B , 5188 (1976). D. Vanderbilt, Journal of Physics and Chemistry of Solids , 147 (2000). B. P. Burton, A. van de Walle, and U. Kattner, Journal of Applied Physics , 113528 (2006). J. B. Bunn, R. Z. Voepel, Z. Wang, E. P. Gatzke, J. A. Lauterbach, and J. R. Hattrick-Simpers, Industrial & Engineering Chemistry Research , 1236 (2016). 4 D. Wu, K. Talley, S. Manna, A. Zakutayev, C. V. Ciobanu, P. Constantine, G. L. Brennecka, and C. Packard, submitted for publication (2017). S. Manna, K. Talley, A. Zakutayev, J. Mangum, G. L. Brennecka, V. Stevanovic, and C. V. Ciobanu, https://arxiv.org/abs/1708.00490 (2017). C. Hums, J. Blasing, A. Dadgar, A. Diez, T. Hempel, J. Christen, and A. Krost, Applied Physics Letters , 022105 (2007). A. M. Holder, S. Siol, P. Ndione, H. Peng, A. M. Deml, B. E. Matthews, L. T. Schelhas, M. F. Toney, R. G. Gordon, W. Tumas, A. Zakutayev, and S. Lany, Science Advances , e1700270 (2017). G. Zhao, X. Xu, H. Li, H. Wei, D. Han, Z. Ji, Y. Meng, L. Wang, and S. Yang, Scientific Reports , 26600 (2016). B. Narayanan, I. E. Reimanis, E. R. Fuller Jr, and C. V. Ciobanu, Physical Review B , 104106 (2010). F. Tasnadi, I. A. Abrikosov, and I. Katardjiev, Applied Physics Letters , 151911 (2009). X. Yang, G. Zeng, and W. Cao,
Analysis of mechanical and electrical damages in piezoelectric ceramics (INTECH Open Access Publisher, 2010). L. Egerton and D. M. Dillon, Journal of the American Ceramic Society , 438 (1959). D. Yang, L. Wei, X. Chao, Z. Yang, and X. Zhou, Physical Chemistry Chemical Physics , 7702 (2016). Q.-M. Wang, X.-H. Du, B. Xu, and L. E. Cross, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control , 638 (1999). M. Kim, J. Kim, and W. Cao, Applied Physics Letters , 132901 (2005). G. Wingqvist, Microsystem Technologies , 1213 (2012). K. Uchino, Smart Materials and Structures , 273 (1998)., 273 (1998).