Abstract

The structural, mechanical, and electronic properties of the orthorhombic GaN (Pnma-GaN) are investigated at ambient pressure by using first-principles calculations method with the ultrasoft pseudopotential scheme. The elastic constants and phonon calculations reveal Pnma-GaN is mechanically and dynamically stable at ambient pressure. The calculated Young modulus of Pnma-GaN is 170 GPa, which is the three-fifths of wurtzite-GaN. Electronic structure study shows that Pnma-GaN is a direct semiconductor with band gap of 1.847 eV. The anisotropic calculation shows that wurtzite-GaN has a smaller elastic anisotropy than that of Pnma-GaN in Young’s modulus. In addition, when the composition of aluminum increases from 0 to 0.063 in the alloy, the band gap decreases initially and increases afterward for Pnma-Ga1−xAlxN, while, for wurtzite-Ga1−xAlxN, the band gap increases with the increasing composition x. Due to the structural porous feature, Pnma-GaN can also be expected to be a good hydrogen storage material.

1. Introductions

A complete acquisition of the physical properties of materials from the macro to the atomic scale is compulsory for device applications. Consequently, theoretical and experimental efforts are beneficial to the study of major properties of main materials and to the search for new materials. Theoretical studies are built on systematic models and computer simulations. One of the most precise microscopic theories in computational materials science is density functional theory (DFT) [1], which sufficiently defines the ground state physical properties of electronic system. It is well known that the GaN single crystal can be considered as the best candidate for wide-gap UV detectors and photodiodes [2, 3]. One of the best ways to improve their technology parameters in our opinion consists of the synthesis of an appropriate solid alloy to continuously vary their band energy gaps creating a continuous set of solid alloy. There have been many published articles of electronic and optical properties in GaN single crystal. The first band energy structure for wurtzite Ga1−xAlxN (x = 0.25, 0.50, and 0.75) has been calculated by Rubio et al. [4]. Then, Malachowski et al. [5] have calculated the band energy structure of Ga1−xAlxN (x = 0.15 and 0.35) solid alloys. The band structure of the Ga1−xAlxN (x = 0.15, 0.52, and 0.65) solid crystalline alloys was calculated by Kityk [6] first using another approach of calculations. This approach includes the use of the orthogonalisation procedure for the plane-wave basis set with respect to the core-like states and the applications of the PA screening procedure; and they found that this approach can be extended to different binary wide-gap semiconductors solid alloys during the estimation of the band energy structure. In [7], a good agreement was predicted with experimental results for the minimum of band gaps in the wurtzite-like Ga1−xAlxN structure. For the zinc-blend structure, it was evaluated that GaN will have a direct band gap of 3.1 eV at Γ point in good agreement with recent reported experiments on cubic GaN crystals (3.35 eV) [8]. Miwa and Fukumoto [9] calculated the GaN with wurtzite-like structure and zinc-blend structure using the first-principle pseudopotential calculations. Similar results were concluded in [10, 11].

Recently, we have found that Pnma-GaN is mechanically and dynamically stable at ambient pressure via elastic constants and phonon calculations. Previous work has not systematically investigated the structural and mechanical properties of the new Pnma phase of GaN and their Ga1−xAlxN alloys. In this work, we have first detailed the equilibrium crystal structure, mechanical properties, and electronic properties of GaN based on DFT.

2. Computational Methods

The self-consistent calculations were performed using density functional theory (DFT) [1, 12] within the Perdew-Burke-Ernzerhof (PBE) exchange-correlation and employing the generalized gradient approximation (GGA) [13] and the local density approximation (LDA) [14, 15] as implemented in the Cambridge Serial Total Energy Package (CASTEP) code [16]. A plane-wave basis set with energy cut-off of 380 eV was used for Pnma-GaN and their Ga1−xAlxN alloys in these calculations. For wurtzite-GaN, we used a plane-wave basis set with energy cut-off 400 eV, which is used with GGA and LDA. The Monkhorst-Pack [17] k-point separation with a grid of 0.025 Å−1 along the three axes in the reciprocal space was used in all calculations. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) [18] minimization scheme was used in geometry optimization. The total energy convergence tests showed that convergence to less than 1 meV/atom was achieved with the above calculation parameters. All energy on atoms converged to less than 1 meV/atom for the geometry relaxation; and all the components of the Hellmann-Feynman stress tensor were reduced to the order of 0.01 GPa at the given pressure. The elastic constants were determined by applying a set of given homogeneous strains with a finite value (0.1%) to the equilibrium lattice and calculating the resulting stress according to optimizing the structure under the strain. The bulk modulus, shear modulus, and Young’s modulus were estimated by employing Voigt-Reuss-Hill approximation [1921]. The phonons spectra of Pnma-GaN are using the linear response approach, or the so-called density functional perturbation theory (DFPT) is one of the most popular methods of ab initio calculation of lattice dynamics [22].

3. Results and Discussion

The calculated lattice parameters of Pnma-GaN are presented in Table 1. Similarly, it is recognized within DFT that the lattice parameters of Pnma-GaN and wurtzite-GaN calculations are commonly underestimated for both GGA and LDA, while the LDA values are close to the experimental values for wurtzite-GaN, so LDA method can obtain more accurate results. The calculated lattice constants of Pnma-GaN are a = 5.900 Å, b = 3.232 Å, and c = 5.374 Å with LDA. The crystal structure of Pnma-GaN is shown in Figure 1(a). There are eight atoms in conventional cell, including four Ga atoms and four N atoms. There is an equivalent position for Ga or N atoms; that is to say, the four Ga atoms or N atoms are equivalent. Hence, it is reasonable that Al atom can replace any Ga atom. The Ga atoms occupy the crystallographic 4c sites in a conventional cell, which is (−0.1673, 0.2500, and 0.4008); and the N atoms occupy the crystallographic 4c sites in a conventional cell, which is (0.1731, 0.2500, and 0.4016), respectively. The crystal structure of Ga1−xAlxN (x = 0.031, 0.063, 0.094, and 0.125) is shown in Figures 1(b), 1(c), 1(d), and 1(e), respectively. In Ga1−xAlxN alloys, aluminum atoms substituted position is the lowest energy substitution method. In addition, wurtzite-GaN results were calculated using the same method; together with previous experimental results obtained other calculated results are presented in Table 1. The calculated lattice constants wurtzite-GaN (space group: P63mc) are a = b = 3.156 Å (3.225 Å), c = 5.148 Å (5.257 Å), and internal parameter u = 0.376 (0.375) with LDA (GGA), respectively. In addition, the experimental results of wurtzite-GaN are a = b = 3.190 Å, c = 5.189 Å [23], and u = 0.376. The calculated lattice parameters of wurtzite-GaN are in good agreement with the reported calculated results of other theoretical results [2326], while larger deviations from experimental results are found in both a and c directions 1.10% and 3.11% with GGA, respectively, but smaller deviations are found in both a and c directions 1.07% and 0.79% with LDA, respectively. Due to the structural porous feature (see Figures 1(b)1(e)), Pnma-GaN can also be expected to be a good hydrogen storage material. At ambient pressure, the enthalpy of Pnma-GaN is slightly greater 0.264 eV per formula unit than that of wurtzite-GaN.

The elastic properties of a solid are important. They are not only closely related to various fundamental solid-state phenomena such as interatomic bonding, equations of state, and phonon spectra, but also linked thermodynamically with specific heat, thermal expansion, Debye temperature, and Grüneisen parameter. Most importantly, knowledge of the elastic constants is essential for many practical applications related to the mechanical properties of a solid. The criteria for mechanical stability of orthorhombic phase are given by [27]:Obviously, the elastic constants of Pnma-GaN satisfy the generalized elastic stability criteria for a stable orthorhombic symmetry. To ensure the stability of Pnma-GaN, the phonon spectra were calculated at ambient conditions (see Figure 2(a)). There is no imaginary frequency, which means Pnma-GaN can be stable at ambient conditions. Calculated single crystal elastic constants , bulk modulus B, shear modulus G, and Young’s modulus for Pnma and wurtzite phases of GaN are listed in Table 2. The calculated elastic constants of wurtzite-GaN are in excellent agreement with the reported calculated results of other theoretical results [28, 35]. The calculated related enthalpies of Pnma-GaN, wurtzite-GaN, zinc blende-GaN, and NaCl-GaN are shown in Figure 2(b). All of them are smaller than the enthalpy of α-Ga and α-N2. That is to say, they can be synthesized in the experiment.

The Debye temperature () is a fundamental physical property and correlates with many physical properties of solids, such as specific heat and the thermal coefficient [36]. It is used to distinguish between high and low temperature regions for a solid. If the temperature , we expect all modes to have the energy of , and if , one expects the high-frequency modes to be frozen; that is, the vibrational excitations arise solely from the acoustic vibrations. The standard formula to evaluate the Debye temperature is given by [37]:where is Planck’s constant, is Boltzmann’s constant, is the number of atoms per formula unit, and , , and represent Avogadro’s number, molar mass, and density, respectively. The average sound velocity is described by the following expression:The transverse and longitudinal elastic wave velocity, and , can be achieved from Navier’s equation as follows:where and are bulk modulus and shear modulus, respectively. The calculated Debye temperature  K is evaluated from (2). The Debye temperature of Pnma-GaN is slightly smaller than that of wurtzite-GaN (600 K).

It is well known that the anisotropy of elasticity is an important implication in engineering science and crystal physics. The directional dependence of Young’s modulus and the projection in , , and planes for Pnma-GaN are illustrated in Figures 3(a) and 3(b). The maximum value of = 221 GPa and the minimal of = 106 GPa. The average value of all directions is 176 GPa, which is very close to CASTEP value (170 GPa; see Table 2). The ratio / = 2.09 shows Pnma-GaN exhibits a large elastic anisotropy. The directional dependence of Young’s modulus and the projection in , , and planes for wurtzite-GaN are shown in Figures 3(c) and 3(d). The maximal value of = 301 GPa and the minimal value of = 235 GPa. Its CASTEP value is 257 GPa, and this value is smaller than the average value of all directions by 1.15% (260 GPa; see Table 2). The ratio of / is 1.28; it shows that wurtzite-GaN exhibits a smaller elastic anisotropy than that of Pnma-GaN.

The electronic band structure of Pnma-GaN is shown in Figure 4(a). The dashed line represents the Fermi level (). The coordinate of high symmetry points in Brillouin zone for Pnma-GaN is , , , , , , , and . From the electronic band structure, it is shown that Pnma-GaN is a direct semiconductor and its band gap is 1.567 eV. It is known that the calculated band gap with DFT is usually underestimated by 30–50%, or even more; the true band gap must be larger than the calculated results. The valence band maximum and conduction band minimum are both at point. The total density of states and the partial wave density of states projected onto different atoms at ambient pressure are also shown in Figure 4. The main features of Pnma-GaN can be illustrated as follows: (i) the part present in the energy part (2.5 to 12.5 eV) of the DOS is mainly due to contributions from the electrons and electrons of Ga atom, and the contribution of the N atom orbital is the smallest; (ii) the states from −7.5 eV to Fermi energy mainly originate from electrons of N atom and with slight contributions from electrons of N atom; (iii) the peak present in the lower energy part (−17.5 to −10 eV) of the DOS is mainly due to contributions from the electrons of N, and the contribution from N- is very small.

The calculated band structures of Ga1−xAlxN (x = 0.031, 0.063, 0.094, and 0.125) for Pnma and wurtzite phases are plotted in Figures 5(a)5(d) and 5(e)5(h), respectively. The coordinate of high symmetry points in Brillouin zone for Pnma-Ga1−xAlxN is , , , and . The dashed line also represents the Fermi level (). Calculating the band structures of Ga1−xAlxN, aluminum atoms substituted position is the lowest energy substitution method. The values of band gap as a function of concentration for Pnma-Ga1−xAlxN and wurtzite-Ga1−xAlxN are listed in Table 3. Similarly, it is recognized within DFT that the energy band gap of electronic structure calculations is commonly underestimated for both GGA and LDA, while the LDA values are close to the experimental values for wurtzite-GaN, so we show the band structures of Pnma-Ga1−xAlxN and wurtzite-Ga1−xAlxN in Figure 5 with LDA. This is generally due to the circumstance that they have sample forms that are not adequately flexible to precisely reproduce both the exchange-correlation energy and equivalent potential [31]. The electronic energy band structure of Ga0.969Al0.031N, Ga0.937Al0.063N, Ga0.906Al0.094N, and Ga0.875Al0.125N for Pnma phase shows that they are direct semiconductors and its band gap is 1.740, 1.640, 1.710, and 1.741 eV, respectively. Furthermore, the electronic energy band structure of Ga0.969Al0.031N, Ga0.937Al0.063N, Ga0.906Al0.094N, and Ga0.875Al0.125N for wurtzite phase shows that they are direct semiconductors and its band gap is 2.217, 2.327, 2.489, and 2.648 eV, respectively. With the same composition x, the band gap of Ga1−xAlxN alloy for Pnma phase is smaller than that of wurtzite phase. With the increasing composition x, the band gap decreases initially and increases afterward for Pnma-Ga1−xAlxN, while, for wurtzite-Ga1−xAlxN, the band gap increases with the composition x increasing.

4. Conclusions

This study reports a detailed investigation of the structural and electronic properties of Pnma-GaN/Ga1−xAlxN semiconductor alloy using the DFT within the ultrasoft pseudopotential scheme in the frame of LDA and GGA. The elastic constants and phonon calculations reveal Pnma-GaN is mechanically and dynamically stable at ambient pressure and its Debye temperature is 510.55 K. Pnma-GaN is a direct semiconductor and its band gap is 1.847 eV within LDA at ambient pressure. The elastic anisotropy properties show that Pnma-GaN exhibits a large anisotropy in its Young’s modulus. When the composition of aluminum increases from 0 to 0.125 in the alloys, the band gap decreases initially and increases afterwards for Pnma-Ga1−xAlxN, while, for wurtzite-Ga1−xAlxN, the band gap increases as the composition x increases. Due to the structural porous feature, Pnma-GaN can also be expected to be a good hydrogen storage material.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the Natural Science Foundation of China (no. 61474089) and Open Fund of Key Laboratory of Complex Electromagnetic Environment Science and Technology, China Academy of Engineering Physics (no. 2015-0214. XY.K).