Published: On

Structure and Elastic Properties of (In,Ga)N

In the design of substrates for the epitaxial growth of InxGa1-xN alloys it is useful to know the elastic properties of the semiconductor as a function of composition. This application note shows the use of MedeA® with VASP 5.2 and the mechanical-thermal (MT) module in computing these properties. Judging by the results for the binaries GaN and InN, the level of accuracy is comparable with that achieved in experiments.

Structure and Elastic Properties of (In,Ga)N

Experimental Facts

InxGa1-xN alloys are important optoelectronic materials due to their direct band gap in the blue region of the visible electromagnetic spectrum. The wavelength of emitted light can be tuned by the ratio between Ga and In. At ambient conditions both GaN and InN crystallize in the hexagonal wurtzite structure as shown above.

Computed Results

Using the PBEsol semi-local exchange- correlation potential [1] and the S500 parameter setting in MedeA-VASP, the computed lattice parameters are in remarkable agreement with experiment for both GaN and InN. For GaN the computed value for the lattice parameter a is 3.179 Å, which is just slightly below the experimental range of 3.181 – 3.197 Å. The computed lattice parameter c of 5.181 Å falls into the experimental range of 5.176 – 5.280 Å. In the case of InN the computed values are also within the experimental range, albeit more on the longer side. It should be noted that this comparison is not completely consistent since the computed values do not contain any temperature effect whereas the experimental data are typically obtained at room temperature. A model structure and computed lattice parameters of an In0.125Ga0.875N are shown below.

It is remarkable that the computed lattice parameters are very close to the values obtained from a simple linear interpolation between GaN and InN, i.e. by applying Vegard’s law. Of course, this is just one of many possible distributions of In and Ga atoms in a wurtzite lattice. A more detailed modeling of the alloy could be obtained by using a larger supercell and exploring a range of different configurations. In reality, the distribution of the cations in this alloy depends on the conditions in the deposition process and subsequent annealing. With these caveats in mind, let us investigate the elastic coefficients of these three systems, namely pure GaN, pure InN, and the alloy model. The computed elastic coefficients and moduli are listed below. The computed values agree rather well with experimental data, as has also been found for a number of different inorganic crystalline systems (see other Application Notes on elastic coefficients).

Elastic coefficients and moduli for hexagonal GaN, InN, and In0.125Ga0.875N.
computed (expt) computed (expt) computed linear interpolation
The indices i,j (1,…,6) denote the sequence xx,yy,zz,yz,xz,xy. All values in GPa.
Experimental values (in parentheses) are from
C₁₁ 351 (390±15) 221 (190±7) 316 335
C₁₂ 142 (145±20) 120 (104±3) 124 139
C₁₃ 106 (106±20) 95 (121±7) 89 105
C₂₂ 315
C₂₃ 91
C₃₃ 397 (398±20) 232 (182±6) 349 376
C₄₄ 93 (105±10) 47 (10±1) 81 87
C₅₅ 81
C₆₆ 95
Bulk modulus 201 (210±10) 144 (140) 176 194
Young’s modulus 271 141 241 255
Shear modulus 106 53 95 99

Again note that the computed values do not contain effects due to thermal expansion, which are present in the experimental data taken at room temperature. Inclusion of thermal expansion would reduce the computed values. As expected, GaN is stiffer than InN and the values of the alloy are between those of the binaries. Due to the lower symmetry of the supercell for the alloy there are additional elastic coefficients in the fourth column, which can be backfolded to the hexagonal cell, e.g. C₂₂ should equal C₁₁ and C₂₃ should equal C₁₃. Within the numerical uncertainties this is clearly the case. It is interesting that the computed elastic coefficients and moduli are systematically smaller than values obtained by linear interpolation, as can be seen by comparing the last two columns. This is different from the lattice parameters, where Vegard’s law applies well for this class of ternary III-V semiconductors.


First-principles calculations provide an efficient way to obtain elastic coefficients for semiconductor alloys and other doped systems. Such calculations are probably much faster and more cost effective than experiments.

MedeA® modules used for this application


  1. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces, J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett. 100, 136406 (2008)
Full Version (pdf):