Next Article in Journal
Conjugated Polymer Nanoparticles for Bioimaging
Next Article in Special Issue
GaN-Based Laser Wireless Power Transfer System
Previous Article in Journal
LaAlO3:Mn4+ as Near-Infrared Emitting Persistent Luminescence Phosphor for Medical Imaging: A Charge Compensation Study
Previous Article in Special Issue
Diffusion-Driven Charge Transport in Light Emitting Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanical, Thermodynamic and Electronic Properties of Wurtzite and Zinc-Blende GaN Crystals

1
School of Mechanical and Electronic Engineering, Guilin University of Electronic Technology, Guilin 541004, China
2
EEMCS Faculty, Delft University of Technology, 2628 Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Materials 2017, 10(12), 1419; https://doi.org/10.3390/ma10121419
Submission received: 30 October 2017 / Revised: 7 December 2017 / Accepted: 8 December 2017 / Published: 12 December 2017
(This article belongs to the Special Issue Light Emitting Diodes and Laser Diodes: Materials and Devices)

Abstract

:
For the limitation of experimental methods in crystal characterization, in this study, the mechanical, thermodynamic and electronic properties of wurtzite and zinc-blende GaN crystals were investigated by first-principles calculations based on density functional theory. Firstly, bulk moduli, shear moduli, elastic moduli and Poisson’s ratios of the two GaN polycrystals were calculated using Voigt and Hill approximations, and the results show wurtzite GaN has larger shear and elastic moduli and exhibits more obvious brittleness. Moreover, both wurtzite and zinc-blende GaN monocrystals present obvious mechanical anisotropic behavior. For wurtzite GaN monocrystal, the maximum and minimum elastic moduli are located at orientations [001] and <111>, respectively, while they are in the orientations <111> and <100> for zinc-blende GaN monocrystal, respectively. Compared to the elastic modulus, the shear moduli of the two GaN monocrystals have completely opposite direction dependences. However, different from elastic and shear moduli, the bulk moduli of the two monocrystals are nearly isotropic, especially for the zinc-blende GaN. Besides, in the wurtzite GaN, Poisson’s ratios at the planes containing [001] axis are anisotropic, and the maximum value is 0.31 which is located at the directions vertical to [001] axis. For zinc-blende GaN, Poisson’s ratios at planes (100) and (111) are isotropic, while the Poisson’s ratio at plane (110) exhibits dramatically anisotropic phenomenon. Additionally, the calculated Debye temperatures of wurtzite and zinc-blende GaN are 641.8 and 620.2 K, respectively. At 300 K, the calculated heat capacities of wurtzite and zinc-blende are 33.6 and 33.5 J mol−1 K−1, respectively. Finally, the band gap is located at the G point for the two crystals, and the band gaps of wurtzite and zinc-blende GaN are 3.62 eV and 3.06 eV, respectively. At the G point, the lowest energy of conduction band in the wurtzite GaN is larger, resulting in a wider band gap. Densities of states in the orbital hybridization between Ga and N atoms of wurtzite GaN are much higher, indicating more electrons participate in forming Ga-N ionic bonds in the wurtzite GaN.

1. Introduction

Group-III nitride semiconductors have attracted considerable attention in recent years because of their great potential for technological applications such as short-wavelength light-emitting diodes (LEDs), photocatalysts, optoelectronic nanodevices as well as high-temperature, high-power and high-frequency electronic devices [1,2,3,4]. Among the Group-III nitride semiconductors, GaN is known as one of the most promising materials for third-generation semiconductors because of its excellent properties such as a wide direct band gap, low dielectric constant, corrosion resistance, radiation resistance, high luminous efficiency, high temperature resistance and high strength [5,6,7,8,9]. In addition, GaN can be made magnetic by introducing Mn magnetic dopants and is considered to be suitable for future spintronic devices [10,11,12,13].
Under ambient conditions, GaN presents a hexagonal wurtzite structure, while thin films of single-crystal zinc-blende GaN have been reported to be grown on β-SiC (100) substrates in a modified molecular-beam epitaxy system [14]. Zinc-blende is hoped to be more amenable to doping than the wurtzite GaN [15,16]. Most of the previous studies focused on the effect of surface adsorption [17,18,19,20], defects [21,22,23] and doping [24,25,26,27] on the chemical adsorption performance and their effect on electronic and optoelectronic properties of wurtzite GaN. Li et al. [28] pointed out that the density of states around the Fermi level for the wurtzite structure are much lower than that for the zinc-blende structure, which results in a wider band gap. Luo et al. [29] investigated the electronic structure and the elastic and thermal properties of zinc-blende GaN by first-principles calculations. Thus far, for the limitation of experimental methods, systematic and in-depth studies on the mechanical and thermodynamic of wurtzite and zinc-blende GaN crystals remain very limited, especially for the anisotropic behavior of mechanical properties. Fortunately, with the rapid development of numerical methods, calculations and predictions on many mechanical and physical properties of compounds and alloys based on density functional theory (DFT) have been successfully performed and their validity has been verified [30,31,32,33].
To help understand the physical properties of the wurtzite and zinc-blende GaN crystals, establish finite element models, analyze microelectronics reliability, and design components and devices, in this study, first-principles calculations based on DFT have been applied to investigate the elastic constants, elastic properties, elastic anisotropies, thermodynamic properties and electronic properties of wurtzite and zinc-blende GaN crystals.

2. Computational Method and Details

Under ambient conditions, GaN crystallizes in a hexagonal wurtzite structure with the space group P63 mc having lattice constants a = b, c and lattice angles α = β = 90°, γ = 120° [28]. The zinc-blende GaN can be described as a cubic structure with a = b = c and α = β = γ = 90°, having the space group F-43m [28]. The two structures for GaN crystals are shown in Figure 1. In this study, first principle calculations based on DFT calculations were performed using the Cambridge Sequential Total Energy Package (CASTEP) program [34], which has already been employed successfully to calculate the physical properties of numerous alloys and compounds. The generalized gradient approximation (GGA) [35] of the revised Perdew-Burke-Ernzerhof formalism [36] and the local density approximation (LDA) proposed by Ceperley and Alder and parameterized by Perdew and Zunger (CA-PZ) were performed to calculate the exchange-correlation potential. Meanwhile, the Broyden Fletcher Goldfarb Shanno algorithm [27] and Vanderbilt ultra-soft pseudopotentials [37] was used to optimize the crystal models. The electronic wave functions were expanded with an energy cut-off of 520 eV for all calculations. For the Brillouin-zone integration, a 12 × 12 × 12 Monkhorst-Pack mesh [38] was selected for the calculations based on the variation in total energy of less than 1 meV atom−1 [30]. The self-consistent field tolerance was set at 5.0 × 10−7 eV atom−1. It was considered that the system would reach the ground state when the convergence precision of the energy, maximum force, maximum displacement and maximum stress were less than 5.0 × 10−6 eV, 0.01 eV Å−1, 5.0 × 10−4 Å and 0.02 GPa, respectively. Moreover, the valence-electron configurations of Ga and N were 3d104s24p1 and 2s22p3, respectively, For the limitation of both GGA and LDA in calculating band gaps and structures, in this study, Heyd-Scuseria-Ernzerhof exchange-correlation function (HSE06 [39]) was applied to address the band gap problem of the two GaN crystals. The stiffness matrix Cij was obtained by linear fitting of four small strains (±0.001 and ±0.003) under nine deformation conditions [40], and elastic compliance matrix Sij was calculated as the inverse of Cij, i.e., Sij = [Cij]−1.

3. Results and Discussion

3.1. Elastic Constants

The calculated lattice constants and equilibrium volumes in the ground state are listed in Table 1, where the experimental results of wurtzite and zinc-blende GaN crystals obtained from previous studies are also given for comparison. Obviously, the calculation results show good consistency with the experimental data, thus demonstrating the effectiveness of the calculations in this study. Moreover, the lattice parameters of wurtzite and zinc-blende structures are commonly overestimated by GGA while underestimated by LDA. The calculation results of the lattice constants obtained by LDA are closer to the experimental results as the maximum deviation of the lattice constant (herein, a, b and c) is about 1%, which is slightly less than that obtained by the GGA method.
The elastic constant is an essential property that can provide a link between the mechanical and dynamic behavior of crystals and give important information concerning the nature of the chemical bonding operating in solids [30]. In the stiffness matrix Cij of hexagonal crystals, there are five independent elastic constants, C11, C12, C13, C33 and C44, and the elastic constants for the hexagonal wurtzite GaN crystal should satisfy the following stability criteria [42]:
C 11 > 0 , ( C 11 C 12 ) > 0 , C 44 > 0 , ( C 11 + C 12 ) C 33 > 2 C 13 2
In the stiffness of cubic zinc-blende GaN, there are three independent elastic constants, C11, C44 and C12; restrictions on the elastic constants are [42]:
( C 11 + 2 C 12 ) > 0 , C 44 > 0 , ( C 11 C 12 ) > 0
The calculated elastic constants Cij and Sij for the two GaN crystals are listed in Table 2 and Table 3, respectively, where the previous experimental data of Cij are consistent with the results obtained in this work, thus proving the correctness of the simulation carried out in this paper. Clearly, according to the result shown in Table 2, the elastic constants of the two GaN crystals obey the above stability criteria.

3.2. Elastic Properties

The bulk modulus (B) is a parameter to measure how incompressible/resistant to compressibility the materials are. B is defined as the ratio of the increase in infinite small pressure to the resulting relative volume change and can be predicted from the calculated elastic constants within the Voigt–Reuss scheme [42]:
B V = 1 9 [ C 11 + C 22 + C 33 + 2 ( C 12 + C 13 + C 23 ) ]
B R = [ S 11 + S 22 + S 33 + 2 ( S 12 + S 13 + S 23 ) ] 1
G V = 1 15 [ ( C 11 + C 22 + C 33 ) ( C 12 + C 13 + C 23 ) + 3 ( C 44 + C 55 + C 66 ) ]
G R = 15 [ 4 ( S 11 + S 22 + S 33 ) 4 ( S 12 + S 13 + S 23 ) + 3 ( S 44 + S 55 + S 66 ) ] 1
where BV, BR and GV, GR stand for the upper and lower bounds of the bulk moduli and shear moduli of the polycrystalline aggregate, respectively. Moreover, the effective bulk modulus (B) and effective shear modulus (G) can be calculated using the Voigt-Reuss-Hill approximation [46], which are considered arithmetic means:
B = ( B V + B R ) / 2
G = ( G V + G R ) / 2
The effective elastic modulus (E) and Poisson’s ratio (ν) of the GaN polycrystalline aggregate can be expressed as follows [47]:
E = 9 B G 3 B + G
ν = 3 B 2 G 6 B + 2 G
The calculated B, G, E, ν and B/G for the two GaN polycrystals are listed in Table 4. Generally, the two GaN polycrystals have many common mechanical properties. The buck moduli of the two crystals are similar, while the other modulus values (i.e., G and E) for wurtzite GaN are a little larger, indicating that, compared to zinc-blende GaN, the stiffness of wurtzite GaN is slightly higher. According to the Pugh criterion, if B/G > 1.75, then ductile behavior may occur, otherwise the material behaves in a brittle mode [48]. Obviously, the ratio of B to G for hexagonal wurtzite GaN is 1.31 (GGA) or 1.61 (LDA), which is much less than 1.75, as given in Table 4, implying that wurtzite GaN may act like a brittle material. For zinc-blende GaN, the ratio is 1.59 (GGA) or 1.76 (LDA), i.e., it is close to 1.75. This implies that zinc-blende GaN might also be a brittle material while its brittleness is less than that of wurtzite GaN. Moreover, Poisson’s ratio (ν) is widely applied to evaluate the stability of the crystals under shear deformation, and a larger ν indicates that the materials would present a better plasticity. Zinc-blende GaN exhibits a better plastic behavior because of the larger ν (0.24 for GGA and 0.26 for LDA), while wurtzite GaN is more stable against shear stress for the smaller ν (0.20 for GGA and 0.24 for LDA). Furthermore, the value of ν implies the degree of directionality of the covalent bond. Its value is small (about 0.1) for covalent materials, while, for ionic materials, the typical value is 0.25 [49]. The calculated value of ν ranges from 0.20 to 0.27 for the two structures, indicating they are both ionic materials.

3.3. Elastic Anisotropy

The elastic anisotropy of materials has a significant effect on their physical properties, such as anisotropic deformation, crack propagation and elastic instability. In this section, the elastic anisotropy factor (A) is applied to explore the anisotropy degree of the two GaN monocrystals [56], as given in Table 5. If the value of anisotropy factor A is close to 1, then the mechanical behavior of the plane is prone to isotropy. The calculated results show that, in wurtzite GaN, A of planes containing the [001] axis (i.e., [0001] axis or orientation) is far from 1, indicating wurtzite GaN displays obvious anisotropic characteristics in these planes. In addition, zinc-blende GaN presents strong anisotropy in its crystal planes {100} and {110}.
To further investigate the anisotropic characteristics of monocrystal GaN crystals, three-dimensional surfaces revealing the elastic anisotropy were analyzed. The direction dependence of the elastic modulus of general monocrystals is given as follows [42]:
1 / E = S 11 l 1 4 + 2 S 12 ( l 1 l 2 ) 2 + 2 S 13 ( l 1 l 3 ) 2 + 2 S 14 ( l 1 2 l 2 l 3 ) + 2 S 15 ( l 3 l 1 3 ) + 2 S 16 ( l 1 3 l 2 ) + S 22 l 2 4 + 2 S 23 ( l 2 l 3 ) 2 + 2 S 24 ( l 2 3 l 3 ) + 2 S 25 ( l 1 l 2 2 l 3 ) + 2 S 26 ( l 1 l 2 3 ) + S 33 l 3 4 + 2 S 34 ( l 2 l 3 3 ) + 2 S 35 ( l 1 l 3 3 ) + 2 S 36 ( l 1 l 2 l 3 2 ) + S 44 ( l 2 l 3 ) 2 + 2 S 45 ( l 1 l 2 l 3 2 ) + 2 S 46 ( l 1 l 3 l 2 2 ) + S 55 ( l 1 l 3 ) 2 + 2 S 56 ( l 2 l 3 l 1 2 ) + S 66 ( l 1 l 2 ) 2
where l1, l2 and l3 are the direction cosines concerning the a, b and c directions of the crystal lattices, respectively, and Sij (i, j = 1, 2, 3) are the compliance coefficients, which are listed in Table 3. According to Equation (11) and the crystal features of hexagonal and cubic monocrystals, the direction dependences of the elastic moduli of the wurtzite and zinc-blende can be expressed as in Equations (12) and (13), respectively.
1 E = ( 1 l 3 2 ) 2 S 11 + l 3 4 S 33 + l 3 2 ( 1 l 3 2 ) ( 2 S 13 + S 44 )
1 E = S 11 + 2 ( S 11 S 12 0.5 S 44 ) ( l 1 2 l 2 2 + l 2 2 l 3 2 + l 3 2 l 1 2 )
The simulation results of direction dependences of the elastic moduli calculated by GGA and LDA are similar, and LDA method was applied for simulating the direction dependences of the elastic moduli for the higher accuracy in calculating structural parameters and elastic constants. Figure 2 shows the direction dependences of the elastic moduli and their projections on the main planes of GaN crystals. For wurtzite GaN, the shape of elastic modulus at plane (001) is a circle and the value is 324.8 GPa, while it is far from a circle at planes containing [001] axis (e.g., planes ( 01 1 ¯ 0 ) and ( 2 1 ¯ 1 ¯ 0 )) (Figure 2a,b); therefore, the elastic modulus of wurtzite GaN may exhibit an isotropic characteristic at the plane (001) while it shows significant anisotropy at the planes containing [001] axis, where the maximum and minimum values are 415.8 GPa and 267.8 GPa, respectively. For zinc-blende GaN, the elastic moduli along the [100] axis (i.e., Ex), [010] axis (i.e., Ey) and the [001] axis (i.e., Ez) are all the minimum value 175.0 GPa, and the maximum elastic moduli are located at the orientations <111> and the value is 388.5 GPa, as shown in Figure 2c. In addition, in the zinc-blende GaN, there are obvious anisotropic characteristics for all small miller indices crystallographic planes, and Figure 2d shows that the maximum and minimum elastic moduli in the projection of planes {100} are 297.7 and 175.0 GPa, respectively.
Moreover, shear moduli of hexagonal wurtzite GaN and cubic zinc-blende GaN can be expressed as Equation (14) [57] and Equations (15)–(17), respectively,
1 G = S 55 + ( S 11 S 12 S 44 2 ) ( 1 l 3 2 ) + 2 ( S 11 + S 33 2 S 13 S 44 ) ( 1 l 3 2 ) l 3 2
1 G = ( S 44 + 4 S 0 J )
S 0 = S 11 S 12 1 2 S 44
J = sin 2 θ cos 2 θ + 0.125 sin 4 θ ( 1 cos 4 ϕ )
where θ and φ are the Euler angles in the measurement systems. Figure 3a,b,d shows the direction dependences of the shear moduli and their projections on the main planes of GaN crystals. For wurtzite GaN, it is found that the shape of directional shear modulus at the (001) plane is a circle, while it is far from the circle at the planes containing the [001] axis. Therefore, the shear modulus of wurtzite GaN may exhibit an isotropic characteristic at the plane (001) and the isotropic value is 110.0 GPa, while it shows significant anisotropy at planes containing [001] axis, and the maximum and minimum values are 130.8 GPa and 98.9 GPa, respectively (see Figure 3a,b). For zinc-blende GaN, the shear moduli along the [100] axis (i.e., Ex), [010] axis (i.e., Ey) and the [001] axis (i.e., Ez) are all the maximum value 165.0 GPa, and the minimum shear moduli are located at the orientations <111> and the value is 81.0 GPa (see Figure 3d). In addition, in the zinc-blende GaN, there are also obvious anisotropic characteristics for all small miller indices crystallographic planes, and Figure 3e shows that the maximum and minimum shear moduli in the projection of planes {100} are 165.0 and 92.8 GPa, respectively. Obviously, to some extent, the shear modulus of GaN monocrystals has an inverse direction dependence compared to the elastic modulus. In summary, the elastic moduli and shear moduli of GaN monocrystals exhibit obvious anisotropy behavior, while the following calculation shows that the bulk moduli of the two monocrystals are nearly isotropic. Directional dependences of bulk moduli of hexagonal and cubic GaN monocrystal can be calculated using Equations (18) and (19), respectively,
1 B = ( S 11 + S 12 + S 13 ) ( S 11 + S 12 S 13 S 33 ) l 3 2
1 B = ( S 11 + 2 S 12 ) ( l 1 2 + l 2 2 + l 3 2 )
For wurtzite GaN, the maximum bulk modulus is located at crystal orientations which are vertical to [001] ([0001]) axis, e.g., [100] ( 2 1 ¯ 1 ¯ 0 ), [010] ( [ 01 1 ¯ 0 ] ), etc., and its value is 614.0 GPa. Besides, the minimum bulk modulus is in the direction of [100] axis and the value is 576.8 GPa, as shown in Figure 3c. Obviously, there is very slight anisotropy in the wurtzite GaN. For zinc-blende GaN, there is no anisotropy of bulk modulus, and the value is 602.4 GPa in all orientations (see Figure 3f).
Further, the curves of Poisson’s ratios at the planes containing [001] axis v(θ) in the wurtzite GaN can be expressed as Equation (20) [58], and Poisson’s ratios at the plane (hkl) of zinc-blende GaN can be described as Equation (21) [59],
v ( θ ) = s 12 sin 2 θ + s 13 cos 2 θ s 11 sin 4 θ + s 33 cos 4 θ + ( s 44 + 2 s 13 ) sin 2 θ cos 2 θ
ν ( h k l , θ ) = { S 12 + S 0 h 2 + k 2 + l 2 [ ( h 2 l h 2 + k 2 h 2 + k 2 + l 2 cos θ h k h 2 + k 2 sin θ ) 2 + ( k 2 l h 2 + k 2 h 2 + k 2 + l 2 × cos θ + h k h 2 + k 2 sin θ ) 2 + ( l h 2 + k 2 h 2 + k 2 + l 2 cos θ ) 2 ] } / [ S 11 + 2 S 0 ( h k ) 2 + ( h l ) 2 + ( l k ) 2 ( h 2 + k 2 + l 2 ) 2 ]
In the wurtzite GaN, Poisson’s ratios at the planes containing [001] axis are illustrated in Figure 4a; the maximum value is 0.31, and it is located at the directions vertical to [001] axis. Meanwhile, the minimum ν is 0.158, and the angle between the direction and [001] axis is 19.8°. In addition, ν along [001] axis of wurtzite GaN is 0.162. Additionally, for zinc-blende GaN, Poisson’s ratios at planes (100) and (111) are isotropic, as shown in Figure 4b,c, and the values are 0.35 and 0.18, respectively. However, ν at plane (110) exhibits obvious anisotropic behavior (see Figure 4d), where the maximum value 0.60 is located at [001] and [ 00 1 ¯ ] directions, while the minimum value is 0.10 (actually, −0.10) and it is in the orientations [ 1 ¯ 10 ] or [ 1 1 ¯ 0 ] .

3.4.Thermodynamic Properties

The calculation results of phonon spectra can be used to compute lattice heat capacity (CV) as functions of temperature [60]. The temperature dependence of the energy can be calculated as follows:
E ( T ) = E t o t + E z p + ω exp ( ω k T ) 1 F ( ω ) d ω
where Etot is the total static energy at 0 K, which can be calculated by first principles, Ezp is the zero-point vibrational energy, ω is the frequency, k is the Boltzmann’s constant, T is the Kelvin temperature, ħ is Planck’s constant and F(ω) is the phonon density of states. Ezp can be expressed as follows:
E z p = 1 2 F ( ω ) ω d ω
The lattice contribution to the heat capacity, CV, is calculated as follows:
C V ( t ) = k ( ω k T ) 2 exp ( ω k T ) [ exp ( ω k T ) 1 ] 2 F ( ω ) d ω
The heat capacity calculated by the LDA method is illustrated in Figure 5. At low temperatures, power-law temperature dependence is clearly observed for the heat capacity, and the predicted power-law exponent is approximately 3 for the two GaN crystals according to the calculation result, which is consistent with the Debye heat capacity theory. Moreover, at 300 K, the calculated heat capacities of wurtzite and zinc-blende are 34.4 and 34.3 J mol−1 K−1, respectively, which are close to the experimental value of wurtzite GaN at 300 K (34.1 J mol−1 K−1) [61]. At high temperatures, the heat capacity of the two GaN crystals tends to converge at 49.4 mol−1 K−1, which is the Dulong–Petit limit. Obviously, heat capacities of the two GaN crystals are quite close as the temperature is increased from 0 to 1000 K, and the effect of crystal structures on the heat capacity of GaN is very limited.
Further, Debye temperature is one of important basic physical parameters, which has connections with many mechanical and physical properties such as elastic constants, heat capacity, bond strength, and so on [62]. By employing elastic constants and average sound velocity νm, Debye temperature ΘD can be approximately calculated by
Θ D = k ( 3 n N A ρ 4 π M ) 1 / 3 v m
v m = [ 1 3 ( 1 v l 3 + 2 v t 3 ) ] 1 / 3
where n is the number of atoms in the cell, NA is Avogadro constant, M is the molecular weight, ρ is the density, and νt and νl are longitudinal and transverse sound velocities, respectively, which can be obtained from Equations (27) and (28), respectively.
v t = ( G ρ ) 1 / 2
v l = ( 3 B + 4 G 3 ρ ) 1 / 2
The thermodynamic properties of GaN crystals calculated by LDA method are listed in Table 6. Obviously, the densities of the two structure GaN crystals are extremely close, and wurtzite GaN has higher νt, νl and νm. Besides, the measured Debye temperatures of GaN are scattered, and previous study reported that the experimental measured Debye temperatures of wurtzite GaN range from 586 K to 898 K [63]. In this study, the calculated Debye temperatures of wurtzite and zinc-blende GaN are 641.8 and 620.2 K, respectively. Compared with zinc-blende GaN, wurtzite GaN has the slightly higher Debye temperature, indicating stronger Ga–N ionic bonds in the wurtzite GaN. Therefore, wurtzite GaN exhibits slightly higher stiffness and melt point to a small degree.

3.5. Electronic Properties

Band gap is the energy required for a valence electron to become a conduction electron, which moves freely and serves as a charge carrier. According to previous experimental results, the band gaps of wurtzite and zinc-blende GaN crystals are 3.5 eV [55] and 3.1 eV [64], respectively. However, Band gaps of wurtzite and zinc-blende GaN crystals calculated by LDA are 1.66 eV and 1.48 eV, respectively, and they are 2.01 and 1.91 eV by GGA. Here, HSE06 scheme [39] was applied for band gap calculations, and the calculated band gaps of wurtzite and zinc-blende GaN are 3.62 eV and 3.06 eV, respectively, which are close to the experimental values. Clearly, the results also show that the band gap of wurtzite GaN is larger than that of zinc-blende GaN. In previous studies, it is reported that the band gap of wurtzite GaN is located at the Γ (G) point [65,66]. In the present study, the band structures of wurtzite and zinc-blende GaN crystals are shown in Figure 6. It is also found that the maximum value of valence band and the minimum value of conduction band are located at the G point for the two crystals, indicating that all of them are direct gap semiconductors. At the G point, the lowest energy of conduction band in the wurtzite GaN is larger than that in the zinc-blende GaN, which results in a wider band gap for the wurtzite GaN. Figure 7 shows DOSs of the GaN crystals around band gaps. According to the distribution characteristics of DOSs of GaN crystals, the two crystals exhibit similar atomic bonding and hybrid behavior. In the energy band with the energy ranging from −10 to 0 eV, DOSs of the two GaN crystals are primarily contributed by p orbits of N atoms (N-p). However, for the energy band in the conduction band, which is close to the band gap, e.g., the energies ranging from 3.62 eV to 20 eV for wurtzite GaN, DOSs of GaN crystals are mainly due to the hybridization between N-p and s and p orbits of Ga (Ga-s and Ga-p). Meanwhile, although the two crystals exhibit similar atomic bonding and hybridization behavior, the DOSs of Ga atoms, N atoms and wurtzite GaN are much larger than those of zinc-blende GaN, as observed in Figure 7. For example, the values of A1 (peak of DOS of GaN), B1 (peak of Ga-p) and C1 (peak of N-p) at −2.1 eV are 2.98, 0.48 and 2.36 electrons/eV, respectively, while A2, B2 and C2 at −2.1 eV for zinc-blende GaN are 1.15, 0.22 and 0.91 electrons/eV, respectively. Obviously, the values of A1, B1 and C1 are much larger than those of A2, B2 and C2, respectively. Therefore, it can be inferred that more electrons in corresponding atomic orbitals participate in forming Ga-N ionic bonds in the wurtzite GaN, which is why the wurtzite GaN has larger elastic modulus, higher Debye temperature, smaller Poisson’s ratio and shows more evident characteristic of brittleness.

4. Conclusions

The mechanical, thermodynamic and electronic properties of wurtzite and zinc-Blende GaN crystals were investigated by first-principles calculations based on density functional theory. From the present study, the following conclusions can be drawn:
  • Generally, the two GaN polycrystals have common or similar elastic properties. The bulk moduli of the wurtzite and zinc-blende GaN polycrystals are very close, and polycrystalline wurtzite GaN has slightly larger shear modulus and elastic modulus but smaller Poisson’s ratio and shows more evident characteristic of brittleness.
  • Except for the bulk modulus, the anisotropic behavior of the two GaN monocrystal is quite different, or even opposite. For wurtzite GaN monocrystal, the maximum and minimum elastic moduli are located at orientations [001] and <111>, respectively, while they are in the orientations <111> and <100> in the monocrystal zinc-blende GaN, respectively. Compared to their respective elastic moduli, the shear moduli of the two GaN monocrystals have completely opposite direction dependences. Thus, the differences in their anisotropic behavior should draw enough attention.
  • For wurtzite GaN, Poisson’s ratios at the planes containing [001] axis are anisotropic. The maximum value is 0.31 and it is located at the directions vertical to [001] axis. For zinc-blende GaN, Poisson’s ratios at planes (100) and (111) are isotropic, while the Poisson’s ratio at plane (110) exhibits dramatically anisotropic phenomenon.
  • The Debye temperatures calculated based on elastic constants and average sound velocities of wurtzite and zinc-blende GaN are 641.8 and 620.2 K, respectively, and wurtzite GaN exhibits slightly higher stiffness and melt point for the higher Debye temperature. At the same temperature, the two GaN crystals have the same heat capacity. Zinc-blende GaN is hoped to be more amenable to doping than the wurtzite GaN [15,16], thus it can be a first decent approach of wurtzite GaN due to many common or similar properties.
  • The exchange-correlation functional theory HSE06 is recommended to calculate the band gaps. Band gaps are located at the G point for the two crystals, and the band gaps of wurtzite and zinc-blende GaN are 3.62 and 3.06 eV, respectively. At the G point, the lowest energy of conduction band in the wurtzite GaN is larger than that in the zinc-blende GaN, resulting in a wider band gap.
  • Densities of states in the orbital hybridization between Ga and N atoms of wurtzite GaN are much higher, indicating more electrons in corresponding atomic orbitals participate in forming Ga-N ionic bonds, which is why the wurtzite GaN has larger elastic modulus, higher Debye temperature, smaller Poisson’s ratio and shows more evident characteristic of brittleness.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (NSFC) under grant No. 51505095 and 51503049; Guangxi Natural Science Foundation under grant. No. 2016GXNSFBA380114; the Research Fund for Guangxi Key Laboratory of Manufacturing System & Advanced Manufacturing Technology under grant No. 15-140-30-003Z; and Innovation Project of Guet Graduate Education under grant No. 2018YJCX04.

Author Contributions

Daoguo Yang conceived and designed the study; Hongbo Qin, Xinghe Luan and Chuang Feng performed the first principles calculation; Hongbo Qin wrote the manuscript; and Daoguo Yang and Guoqi Zhang reviewed and edited the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van de Walle, C.G. Wide-Band-Gap Semiconductors; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  2. Kanoun, M.B.; Goumri-Said, S.; Merad, A.E.; Merad, G.; Cibert, J.; Aourag, H. Zinc-blende AlN and GaN under pressure: Structural, electronic, elastic and piezoelectric properties. Semicond. Sci. Technol. 2004, 19, 1220–1231. [Google Scholar] [CrossRef]
  3. Soltani, A.; Cordier, Y.; Gerbedoen, J.C.; Joblot, S.; Okada, E.; Chmielowska, M.; Ramdani, M.R.; De Jaeger, J.C. Assessment of transistors based on GaN on silicon substrate in view of integration with silicon technology. Semicond. Sci. Technol. 2013, 28, 094003. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Wu, Z.F.; Gao, P.F.; Fang, D.-Q.; Zhang, S.-L. Enhanced visible light absorption in ZnO/GaN heterostructured nanofilms. J. Alloys Compd. 2017, 704, 478–483. [Google Scholar] [CrossRef]
  5. Nakamura, S. The roles of structural imperfections in InGaN-based blue light-emitting diodes and laser diodes. Science 1998, 281, 956–961. [Google Scholar] [CrossRef]
  6. Wang, X.; Chang, B.; Ren, L.; Gao, P. Influence of the p-type doping concentration on reflection-mode GaN photocathode. Appl. Phys. Lett. 2011, 98, 082109. [Google Scholar] [CrossRef]
  7. Chen, J.; Xue, C.; Zhuang, H.; Yang, Z.; Qin, L.; Li, H.; Huang, Y. Catalytic synthesis and optical properties of large-scale GaN nanorods. J. Alloys Compd. 2009, 468, L1–L4. [Google Scholar] [CrossRef]
  8. Paskova, T.; Hanser, D.A.; Evans, K.R. GaN Substrates for III-Nitride Devices. Proc. IEEE 2010, 98, 1324–1338. [Google Scholar] [CrossRef]
  9. Zhuang, D.; Edgar, J.H. Wet etching of GaN, AlN, and SiC: A review. Mater. Sci. Eng. R Rep. 2005, 48, 1–46. [Google Scholar] [CrossRef]
  10. Gao, X.; Man, B.; Zhang, C.; Leng, J.; Xu, Y.; Wang, Q.; Zhang, M.; Meng, Y. The important role of Ga vacancies in the ferromagnetic GaN thin films. J. Alloys Compd. 2017, 699, 596–600. [Google Scholar] [CrossRef]
  11. Sawicki, M.; Devillers, T.; Gałęski, S.; Simserides, C.; Dobkowska, S.; Faina, B.; Grois, A.; Navarro-Quezada, A.; Trohidou, K.N.; Majewski, J.A.; et al. Origin of low-temperature magnetic ordering in Ga1−xMnxN. Phys. Rev. B 2012, 85, 205204. [Google Scholar] [CrossRef]
  12. Stefanowicz, S.; Kunert, G.; Simserides, C.; Majewski, J.A.; Stefanowicz, W.; Kruse, C.; Figge, S.; Li, T.; Jakieła, R.; Trohidou, K.N.; et al. Phase diagram and critical behavior of the random ferromagnet Ga1−xMnxN. Phys. Rev. B 2013, 88, 081201. [Google Scholar] [CrossRef]
  13. Simserides, C.; Majewski, J.A.; Trohidou, K.N.; Dietl, T. Theory of ferromagnetism driven by superexchange in dilute magnetic semi-conductors. EPJ Web Conf. 2014, 75, 01003. [Google Scholar] [CrossRef]
  14. Paisley, M.J.; Sitar, Z.; Posthill, J.B.; Davis, R.F. Growth of cubic phase gallium nitride by modified molecular-beam epitaxy. J. Vac. Sci. Technol. A Vac. Surf. Films 1989, 7, 701–705. [Google Scholar] [CrossRef]
  15. Lei, T.; Fanciulli, M.; Molnar, R.J.; Moustakas, T.D.; Graham, R.J.; Scanlon, J. Epitaxial growth of zinc blende and wurtzitic gallium nitride thin films on (001) silicon. Appl. Phys. Lett. 1991, 59, 944–946. [Google Scholar] [CrossRef]
  16. Pankove, J.I. Perspective on gallium nitride. MRS Online Proc. Library Arch. 1989, 162. [Google Scholar] [CrossRef]
  17. Zywietz, T.K.; Neugebauer, J.; Scheffler, M. The adsorption of oxygen at GaN surfaces. Appl. Phys. Lett. 1999, 74, 1695–1697. [Google Scholar] [CrossRef]
  18. Adelmann, C.; Brault, J.; Mula, G.; Daudin, B.; Lymperakis, L.; Neugebauer, J. Gallium adsorption on (0001) GaN surfaces. Phys. Rev. B 2003, 67, 165419. [Google Scholar] [CrossRef]
  19. Rosa, A.L.; Neugebauer, J.; Northrup, J.E.; Lee, C.D.; Feenstra, R.M. Adsorption and incorporation of silicon at GaN (0001) surfaces. Appl. Phys. Lett. 2002, 80, 2008–2010. [Google Scholar] [CrossRef]
  20. Lee, D.S.; Lee, J.H.; Lee, Y.H.; Lee, D.D. GaN thin films as gas sensors. Sens. Actuators B Chem. 2003, 89, 305–310. [Google Scholar] [CrossRef]
  21. Pearton, S.; Zolper, J.; Shul, R.; Ren, F. GaN: Processing, defects, and devices. J. Appl. Phys. 1999, 86, 1–78. [Google Scholar] [CrossRef]
  22. Reshchikov, M.A.; Usikov, D.O.D.A.; Helava, H.; Makarov, Yu. Carbon defects as sources of the green and yellow luminescence bands in undoped GaN. Phys. Rev. B 2014, 90, 235203. [Google Scholar] [CrossRef]
  23. Li, L.; Yu, J.; Hao, Z.; Wang, L.; Wang, J.; Han, Y.; Li, H.; Xiong, B.; Sun, C.; Luo, Y. Influence of point defects on optical properties of GaN-based materials by first principle study. Comput. Mater. Sci. 2017, 129, 49–54. [Google Scholar] [CrossRef]
  24. Li, E.; Wu, B.; Lv, S.; Cui, Z.; Ma, D.; Shi, W. Field emission properties of Ge-doped GaN nanowires. J. Alloys Compd. 2016, 681, 324–329. [Google Scholar] [CrossRef]
  25. Kamimura, J.; Bogdanoff, P.; Ramsteiner, M.; Corfdir, P.; Feix, F.; Geelhaar, L.; Riechert, H. P-Type Doping of GaN Nanowires Characterized by Photoelectrochemical Measurements. Nano Lett. 2017, 17, 1529–1537. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, Y.P.; Zhang, Z.H.; Tan, S.T.; Hernandezmartinez, P.L.; Zhu, B.B.; Lu, S.P.; Kang, X.J.; Sun, X.W.; Demir, H.V. Investigation of p-type depletion doping for InGaN/GaN-based light-emitting diodes. Appl. Phys. Lett. 2017, 110, 033506. [Google Scholar] [CrossRef]
  27. Yang, M.; Chang, B.; Hao, G.; Wang, H.; Wang, M. Optoelectronic properties of GaN, AlN, and GaAlN alloys. Opt. Int. J. Light Electron Opt. 2015, 126, 3357–3361. [Google Scholar] [CrossRef]
  28. Li, S.; Ouyang, C. First principles study of wurtzite and zinc blende GaN: A comparison of the electronic and optical properties. Phys. Lett. A 2005, 336, 145–151. [Google Scholar] [CrossRef]
  29. Luo, B.; Wu, X.; Li, G. Electronic structure, elastic and thermal properties of semiconductor GaX (X = N, P, As, Sb) with zinc blende from first-principles calculation. Int. J. Mod. Phys. B 2014, 28, 1450183. [Google Scholar] [CrossRef]
  30. Chen, W.H.; Cheng, H.C.; Yu, C.F. The mechanical, thermodynamic, and electronic properties of cubic Au4Al crystal via first-principles calculations. J. Alloys Compd. 2016, 689, 857–864. [Google Scholar] [CrossRef]
  31. Guechi, A.; Merabet, A.; Chegaar, M.; Bouhemadou, A.; Guechi, N. Pressure effect on the structural, elastic, electronic and optical properties of the Zintl phase KAsSn, first principles study. J. Alloys Compd. 2015, 623, 219–228. [Google Scholar] [CrossRef]
  32. Liu, X.K.; Zhou, W.; Zheng, Z.; Peng, S.M. The elastic and thermodynamic properties of ZrMo2 from first principles calculations. J. Alloys Compd. 2014, 615, 975–982. [Google Scholar] [CrossRef]
  33. Piskorska-Hommel, E.; Winiarski, M.J.; Kunert, G.; Hommel, D. Polarization-dependent XAFS and density functional theory investigations of the quality of the epitaxial GaMnN structure. J. Alloys Compd. 2017, 725, 632–638. [Google Scholar] [CrossRef]
  34. Segall, M.D.; Lindan, P.J.D.; Probert, M.J.; Pickard, C.J.; Hasnip, P.J.; Clark, S.J.; Payne, M.C. First-principles simulation: Ideas, illustrations and the CASTEP code. J. Phys. Condens. Matter 2002, 14, 2717–2744. [Google Scholar] [CrossRef]
  35. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  36. Hammer, B.; Hansen, L.B.; Nørskov, J.K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 1999, 59, 7413–7421. [Google Scholar] [CrossRef] [Green Version]
  37. Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B 1990, 41, 7892–7895. [Google Scholar] [CrossRef]
  38. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  39. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Erratum: “Hybrid functionals based on a screened Coulomb potential”. [J. Chem. Phys. 2003, 118, 8207 (2003)]. J. Chem. Phys. 2006, 124, 219906. [Google Scholar] [CrossRef]
  40. Chen, W.H.; Yu, C.F.; Cheng, H.C.; Lu, S.T. Crystal size and direction dependence of the elastic properties of Cu3Sn through molecular dynamics simulation and nanoindentation testing. Microelectron. Reliab. 2012, 52, 1699–1710. [Google Scholar] [CrossRef]
  41. Yoshikawa, A.; Ohshima, E.; Fukuda, T.; Tsuji, H.; Oshima, K. Crystal growth of GaN by ammonothermal method. J. Cryst. Growth 2004, 260, 67–72. [Google Scholar] [CrossRef]
  42. Nye, J.F. Lindsay, R.B. Physical Properties of Crystals: Their Representation by Tensors and Matrices; Clarendon Press: Oxford, UK, 1985. [Google Scholar]
  43. Polian, A.; Grimsditch, M.; Grzegory, I. Elastic constants of gallium nitride. J. Appl. Phys. 1996, 79, 3343–3344. [Google Scholar] [CrossRef]
  44. Shimada, K.; Sota, T.; Suzuki, K. First-principles study on electronic and elastic properties of BN, AlN, and GaN. J. Appl. Phys. 1998, 84, 4951–4958. [Google Scholar] [CrossRef]
  45. Wright, A.F. Elastic properties of zinc-blende and wurtzite AlN, GaN, and InN. J. Appl. Phys. 1997, 82, 2833–2839. [Google Scholar] [CrossRef]
  46. Hill, R. The Elastic Behaviour of a Crystalline Aggregate. Proc. Phys. Soc. Sect. A 1952, 65, 349–354. [Google Scholar] [CrossRef]
  47. Kanoun, M.B.; Goumri-Said, S.; Reshak, A.H. Theoretical study of mechanical, electronic, chemical bonding and optical properties of Ti2SnC, Zr2SnC, Hf2SnC and Nb2SnC. Comput. Mater. Sci. 2009, 47, 491–500. [Google Scholar] [CrossRef]
  48. Pugh, S. XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Lond. Rdinb. Dublin Philos. Mag. J. Sci. 1954, 45, 823–843. [Google Scholar] [CrossRef]
  49. Bannikov, V.V.; Shein, I.R.; Ivanovskii, A.L. Electronic structure, chemical bonding and elastic properties of the first thorium-containing nitride perovskite TaThN3. Phys. Status Solidi RRL 2007, 1, 89–91. [Google Scholar] [CrossRef]
  50. Chen, X.; Liang, J.; Xu, Y.; Xu, T.; Jiang, P.; Yu, Y.; Lu, K. Structure and Debye temperature of wurtzite GaN. Mod. Phys. Lett. B 1999, 13, 285–290. [Google Scholar] [CrossRef]
  51. Xia, H.; Xia, Q.; Ruoff, A.L. High-pressure structure of gallium nitride: Wurtzite-to-rocksalt phase transition. Phys. Rev. B 1993, 47, 12925–12928. [Google Scholar] [CrossRef]
  52. Yonenaga, I. Hardness, yield strength, and dislocation velocity in elemental and compound semiconductors. Mater. Trans. 2005, 46, 1979–1985. [Google Scholar] [CrossRef]
  53. Nowak, R.; Pessa, M.; Suganuma, M.; Leszczynski, M.; Grzegory, I.; Porowski, S.; Yoshida, F. Elastic and plastic properties of GaN determined by nano-indentation of bulk crystal. Appl. Phys. Lett. 1999, 75, 2070–2072. [Google Scholar] [CrossRef]
  54. Kisielowski, C.; Krüger, J.; Ruvimov, S.; Suski, T.; Ager, J., III; Jones, E.; Liliental-Weber, Z.; Rubin, M.; Weber, E.; Bremser, M. Strain-related phenomena in GaN thin films. Phys. Rev. B 1996, 54, 17745–17753. [Google Scholar] [CrossRef]
  55. Vurgaftman, I.; Meyer, J.R.; Ram-Mohan, L.R. Band parameters for III–V compound semiconductors and their alloys. J. Appl. Phys. 2001, 89, 5815–5875. [Google Scholar] [CrossRef]
  56. Kei Lau, A.K.M. Elastic anisotropy factors for orthorhombic, tetragonal, and hexagonal crystals. Phys. Rev. B 1998, 58, 8980. [Google Scholar]
  57. Panda, K.B.; Ravi Chandran, K.S. Determination of elastic constants of titanium diboride (TiB2) from first principles using FLAPW implementation of the density functional theory. Comput. Mater. Sci. 2006, 35, 134–150. [Google Scholar] [CrossRef]
  58. Zhang, J.M.; Zhang, Y.; Xu, K.W.; Ji, V. Anisotropic elasticity in hexagonal crystals. Thin Solid Films 2007, 515, 7020–7024. [Google Scholar] [CrossRef]
  59. Zhang, J.M.; Zhang, Y.; Xu, K.W.; Ji, V. Young’s modulus surface and Poisson’s ratio curve for cubic metals. J. Phys. Chem. Solids 2007, 68, 503–510. [Google Scholar] [CrossRef]
  60. Baroni, S.; De Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 2001, 73, 515–562. [Google Scholar] [CrossRef]
  61. Chen, X.L.; Lan, Y.C.; Liang, J.K.; Cheng, X.R.; Xu, Y.P.; Xu, T.; Jiang, P.Z.; Lu, K.Q. Structure and Heat Capacity of Wurtzite GaN from 113 to 1073 K. Chin. Phys. Lett. 1999, 16, 107–108. [Google Scholar] [CrossRef]
  62. Kushwaha, A.K. Vibrational, elastic properties and sound velocities of ZnGa2O4 spinel. Comput. Mater.Sci. 2014, 85, 259–263. [Google Scholar] [CrossRef]
  63. Roder, C.; Einfeldt, S.; Figge, S.; Hommel, D. Temperature dependence of the thermal expansion of GaN. Phys. Rev. B 2005, 72, 085218. [Google Scholar] [CrossRef]
  64. Rubio, A.; Corkill, J.L.; Cohen, M.L.; Shirley, E.L.; Louie, S.G. Quasiparticle band structure of AlN and GaN. Phys. Rev. B Condens. Matter 1993, 48, 11810. [Google Scholar] [CrossRef] [PubMed]
  65. Dingle, R.; Sell, D.D.; Stokowski, S E.; Ilegems, M. Absorption, Reflectance, and Luminescence of GaN Epitaxial Layers. Phys. Rev. B Condens. Matter 1971, 3, 1211–1218. [Google Scholar] [CrossRef]
  66. Ilegems, M.; Dingle, R.; Logan, R.A. Luminescence of Zn- and Cd-doped GaN. J. Appl. Phys. 1972, 43, 3797–3800. [Google Scholar] [CrossRef]
Figure 1. Crystal structures of: Wurtzite GaN (a); and zinc-blende GaN (b); and the relationship between a, b and c axes and crystallographic directions (c).
Figure 1. Crystal structures of: Wurtzite GaN (a); and zinc-blende GaN (b); and the relationship between a, b and c axes and crystallographic directions (c).
Materials 10 01419 g001
Figure 2. Direction dependences of elastic moduli: (a) wurtzite GaN; (b) projections on main crystal planes of wurtzite GaN, where the solid line denotes planes containing [001] axis and the dash line is plane [001]; (c) zinc-blende GaN; and (d) the projection on planes {111} of zinc-blende GaN.
Figure 2. Direction dependences of elastic moduli: (a) wurtzite GaN; (b) projections on main crystal planes of wurtzite GaN, where the solid line denotes planes containing [001] axis and the dash line is plane [001]; (c) zinc-blende GaN; and (d) the projection on planes {111} of zinc-blende GaN.
Materials 10 01419 g002
Figure 3. Directional dependence of elastic properties: (a) shear modulus of wurtzite GaN; (b) the projection of shear modulus of planes containing [001] axis in the wurtzite GaN; (c) bulk modulus of wurtzite GaN; (d) shear modulus of zinc-blende GaN; (e) the projection of shear modulus of planes {100} in the zinc-blende GaN; and (f) bulk modulus of zinc-blende GaN.
Figure 3. Directional dependence of elastic properties: (a) shear modulus of wurtzite GaN; (b) the projection of shear modulus of planes containing [001] axis in the wurtzite GaN; (c) bulk modulus of wurtzite GaN; (d) shear modulus of zinc-blende GaN; (e) the projection of shear modulus of planes {100} in the zinc-blende GaN; and (f) bulk modulus of zinc-blende GaN.
Materials 10 01419 g003
Figure 4. Poisson’s ratios v for the planes containing [001] ([0001]) axis in the wurtzite GaN (a); and three low index planes (100) (b); (111) (c); and (110) (d) in the zinc-blende GaN.
Figure 4. Poisson’s ratios v for the planes containing [001] ([0001]) axis in the wurtzite GaN (a); and three low index planes (100) (b); (111) (c); and (110) (d) in the zinc-blende GaN.
Materials 10 01419 g004
Figure 5. Calculation results of temperature-dependent heat capacities for wurtzite and zinc-blende GaN. The dash line denotes experimental data of wurtzite GaN obtained from Ref. [62].
Figure 5. Calculation results of temperature-dependent heat capacities for wurtzite and zinc-blende GaN. The dash line denotes experimental data of wurtzite GaN obtained from Ref. [62].
Materials 10 01419 g005
Figure 6. Band structures of GaN crystals: (a) Wurtzite GaN; and (b) zinc-blende GaN. The Fermi level is set to zero (see the red line).
Figure 6. Band structures of GaN crystals: (a) Wurtzite GaN; and (b) zinc-blende GaN. The Fermi level is set to zero (see the red line).
Materials 10 01419 g006
Figure 7. DOSs for the two GaN crystals and their atoms: (a1) wurtzite GaN and Ga atoms; (a2) N atoms in wurtzite GaN; (b1) zinc-blende GaN and Ga atoms; and (b2) N atoms in zinc-blende GaN. The mark * indicates the peak of DOS at −2.1 eV
Figure 7. DOSs for the two GaN crystals and their atoms: (a1) wurtzite GaN and Ga atoms; (a2) N atoms in wurtzite GaN; (b1) zinc-blende GaN and Ga atoms; and (b2) N atoms in zinc-blende GaN. The mark * indicates the peak of DOS at −2.1 eV
Materials 10 01419 g007
Table 1. Lattice constants of GaN crystals.
Table 1. Lattice constants of GaN crystals.
StructureMethodLattice Constants
a0 = b0 (Å)c0 (Å)c0/a0V0
Wurtzite GaNGGA (present)3.2425.2801.62948.075
LDA (present)3.1565.1451.63144.373
Expt. [41]3.1895.1851.62645.671
Zinc-blende GaNGGA (present)4.5824.5821.00096.221
LDA (present)4.4614.4611.00088.800
Expt. [15]4.4904.4901.00090.519
Table 2. Calculated elastic constant Cij, GPa.
Table 2. Calculated elastic constant Cij, GPa.
StructureC11 = C22C12C13 = C23C33C44 = C55C66RMSE *Note
Wurtzite GaN353.9382.5951.79380.4498.49135.6737.90GGA
374.35126.5680.85441.9498.86123.8921.35LDA
390 ± 15145 ± 20106 ± 20398 ± 20105 ± 10123 ± 10-Expt. [43]
Zinc-blende GaN 242.41122.05122.05242.41146.82146.82-GGA
286.93152.77152.77286.93164.97164.97-LDA
285161161285149149-LDA [44]
293159159293155155-LDA [45]
* Root Mean Square Error.
Table 3. Calculated elastic constant Sij, × 10−3/GPa.
Table 3. Calculated elastic constant Sij, × 10−3/GPa.
StructureS11=S22S12S13= S23S33S44S55S66Note
Wurtzite GaN3.026−0.659−0.3222.71610.15310.1537.371GGA
3.079−0.957−0.3882.40510.11510.1158.071LDA
Zinc-blende GaN 6.224−2.084−2.0846.2246.8116.8116.811GGA
5.712−2.027−2.0275.7126.0616.0616.061LDA
Table 4. Bulk modulus B, shear modulus G, elastic modulus E, Poisson’s ratio ν and B/G of polycrystalline GaN.
Table 4. Bulk modulus B, shear modulus G, elastic modulus E, Poisson’s ratio ν and B/G of polycrystalline GaN.
StructureB, GPaG, GPaE, GPaνB/GNote
Wurtzite GaN162.3124.1296.70.201.31GGA
196.3121.7302.60.241.61LDA
170 [50] 188 [51]116 [52]295 ± 3 [53]0.23 ± 0.06 [54]-Expt.
Zinc-blende GaN 162.2102.2253.40.241.59GGA
199.1113.1285.30.261.76LDA
203.7110.71281.180.271.84Expt. [55]
Table 5. Anisotropy factor A [56].
Table 5. Anisotropy factor A [56].
StructureSymmetryAnisotropy Factor ACalculation Results
GGALDA
Wurtzite GaNPlanes containing the [001] axis C 44 ( C 11 + 2 C 13 + C 33 ) / ( C 11 C 33 C 13 2 ) 0.630.61
Zinc-blende GaN{100} 2 C 44 / ( C 11 C 12 ) 2.442.46
{110} C 44 ( C L + 2 C 12 + C 11 ) / ( C L C 11 C 12 2 ) *1.851.85
* For cubic crystals CL = C66 + (C11 + C12)/2.
Table 6. The calculated thermodynamic properties of GaN.
Table 6. The calculated thermodynamic properties of GaN.
Structureρ, g/cmνt, m/sνl, m/sνm, m/sθD, K
Wurtzite GaN6.2684406.3867563.5814888.006641.8
Zinc-blende GaN6.2644249.6857474.8024724.200620.2

Share and Cite

MDPI and ACS Style

Qin, H.; Luan, X.; Feng, C.; Yang, D.; Zhang, G. Mechanical, Thermodynamic and Electronic Properties of Wurtzite and Zinc-Blende GaN Crystals. Materials 2017, 10, 1419. https://doi.org/10.3390/ma10121419

AMA Style

Qin H, Luan X, Feng C, Yang D, Zhang G. Mechanical, Thermodynamic and Electronic Properties of Wurtzite and Zinc-Blende GaN Crystals. Materials. 2017; 10(12):1419. https://doi.org/10.3390/ma10121419

Chicago/Turabian Style

Qin, Hongbo, Xinghe Luan, Chuang Feng, Daoguo Yang, and Guoqi Zhang. 2017. "Mechanical, Thermodynamic and Electronic Properties of Wurtzite and Zinc-Blende GaN Crystals" Materials 10, no. 12: 1419. https://doi.org/10.3390/ma10121419

APA Style

Qin, H., Luan, X., Feng, C., Yang, D., & Zhang, G. (2017). Mechanical, Thermodynamic and Electronic Properties of Wurtzite and Zinc-Blende GaN Crystals. Materials, 10(12), 1419. https://doi.org/10.3390/ma10121419

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop