Introduction
After the discovery of half-metallicity in the half-Heusler NiMnSb [1], Heusler compounds having the XYZ phase have proved their efficiency in the most revolutionized fields of technology such as optoelectronics [2-5], spintronics [6-8], and thermoelectric applications [9]. These compounds are mainly inter-metallic materials that consist of transition metals (X and Y) and a main group (Z) element. They have a non-centrosymmetric structure (C1b) conforming to an F-43m space group (No. 216) [10] with Wyckoff positions 4a (0,0,0), 4b (1/2,1/2,1/2) and 4c (1/4,1/4,1/4). Half-Heusler compounds are formed in a crystal of the MgAgAs type structure as four interpenetrating face-centered cubic superlattices with type 1, type 2 or type 3 structural phases [11]. In this work, we have studied the effect of pressure on the band gaps, the mechanical stability at ambient conditions and under pressure in addition to the thermoelectric properties of four lithium based half-Heusler compounds using first principles calculations [12]. These Li-based systems have been determined to be semiconducting with indirect large and narrow band gaps. The effect of pressure on all elastic parameters have been investigated. LiAlX(X = C, Si, Ge, Sn) systems have proved to be stable at high values of pressures. They also possess a high figure of merit (ZT), which makes them very attractive candidates in the fields of photovoltaics [13] and thermoelectric device applications [14].
1.1 Computational details
Using the full-potential linearized augmented plane-wave (FP-LAPW) method [15-17] implemented in the WIEN2k package [18], we have calculated the electronic properties, magnetic structures, mechanical stability and thermoelectric properties of a series of LiAlX (X = C, Si, Ge, Sn) half-Heusler compounds. The electronic exchange-correlation energy was treated using the generalized gradient approximation parameterized by Perdew-Burke-Ernzerhof (GGA-PBE) [19,20]. In addition, the modified Becke-Johnson (mBJ) approximation [21] was employed to predict the band gaps more precisely. The total energy dependence on the cell volume was fitted to the modified Murnaghan equation of state [22].
The cut-off parameter RMT*Kmax in the basis set was chosen to be 8 [23], where Kmax is the plane wave cut-off and RMT is the smallest muffin-tin radius. 2000 k-points have been used in the calculations, where the cut-off energy separating the core and the valence states was set at -6.0 Ry. For the elastic constants, we have employed the Thomas Charpin method [24]. While the thermoelectric properties were calculated using the Boltzmann theory as implemented in BoltzTrap code [25].
2 Results and discussion
2.1 Structural properties
The ternary half-Heusler compounds crystallize in three different structural phases labelled as Type 1, Type 2 and Type 3. These phases differ according to the positions of atoms in the cubic lattice structure. The optimal structural parameters of the LiAlX (X = C, Si, Ge, Sn) compounds were calculated in the nonmagnetic and ferromagnetic states for each structural phase. The total energy as function of the unit cell volume are shown in Figs. 1 and 2. We present as well the equilibrium lattice parameters, calculated band gaps, bulk modulus and its pressure derivative in Table I.
Compound | Type | ao(Å) | Bo(GPa) | B’0 | Enm—Efm | Gap(eV) | |||
---|---|---|---|---|---|---|---|---|---|
NM | FM | NM | NM | NM | GGA | GGA-mBJ | |||
1 | -576.9193 | -576.9183 | 4.91 | 127.4 | 3.69 | 10-3 | 0.81 | 1.64 | |
LiAlC | 2 | -576.6792 | -576.6892 | 0.81[34 | 1.63[34 | ||||
3 | -576.7492 | -576.7325 | |||||||
1 | -1081.8077 | -1081.8057 | 5.93 | 63 | 3.7 | 2 × 10-4 | 0.11 | 0.49 | |
LiAlSi | 2 | -1081.7311 | -1081.7352 | 0.118[34 | 0.51[34 | ||||
3 | -1081.7977 | -1081.7975 | |||||||
1 | -4698.9224 | -4698.9223 | 6.01 | 55.8 | 3.91 | 1.3 × 10-4 | 0.055 | 0.52 | |
LiAlGe | 2 | -4698.8462 | -4698.8562 | 0.053 | 0.52[34 | ||||
3 | -4698.8930 | -4698.8931 | |||||||
1 | -12859.0031 | -12858.9995 | 6.47 | 42.88 | 3.56 | 3.5 × 10-3 | 0.16 | 0.55 | |
LiAlSn | 2 | -12858.9586 | -12858.9584 | ||||||
3 | -12858.9967 | -12858.9961 |
Both Figs. 1 and 2 show that Type 1 structure is the most stable among all the structural phases. Furthermore, the minimum energy corresponds to the nonmagnetic state, which indicates a paramagnetic behavior of the studied compounds. The remaining properties will be estimated for the Type 1 nonmagnetic most stable structure for all four Li-based systems.
The resulting lattice parameters are in good agreement with other published experimental data [26,27]. From Table I, we notice that the lattice parameter increases downward in the Mendeleev periodic table from LiAlC to LiAlSn, while the bulk modulus decreases.
2.2 Mechanical properties
The study of the mechanical stability and durability of a crystal against applied forces is very important to determine the strength and resistivity of any material. The three independent elastic constants C11, C12 and C44 of the cubic LiAlX (X = C, Si, Ge, Sn) structures were calculated at normal pressure and under hydrostatic pressure. The calculations of the stability criteria and the evolution of Cij, bulk B, shear G and Young modulus E were also illustrated as function of pressure in Fig. 3. From this figure it can be observed that with the rise in compression the value of all elastic constants increases and we notice that the stability criteria are verified for all LiAlX alloys with changing pressure.
The constant C11, related to the linear elasticity, increases with pressure, signaling that pressure enhances the tensile strength of all compounds. The C44 elastic constant, which is related to shape deformation, becomes less affected against rising pressures, indicating less bond strength enhancement.
The stability of cubic systems at zero pressure can be studied using the Born-Huang stability conditions [28]
However, Born’s elastic stability criteria are inapplicable when the pressure is different from zero. The calculated values of the elastic constants under hydrostatic pressure were compared with the generalized stability criteria [29] using the following relations:
Figure 4 shows the pressure dependence of the generalized stability criteria for LiAlX(X = C, Si, Ge, and Sn) compounds. Our calculations show that
From the elastic constants, we have derived the rest of the mechanical parameters. The bulk values of LiAlX alloys that were extracted from the Murnaghan equation of state are very close to the values calculated using the elastic constants at zero pressure for all alloys. It is well known that the shear modulus G represents the resistance of a material against plastic deformation, while Young’s modulus Y [30] is associated with the stiffness of a solid. From Fig. 3, we notice that, at zero pressure, LiAlC has a larger shear modulus than the rest of materials, which confirms its high stiffness. By increasing the pressure, the values of G and E increase as well. In fact, the higher the pressure, the higher the stiffness of these materials as well as their resistivity against deformation.
The ductility and brittleness of a material is also explained by Pugh’s ratio (k) defined by B/G
[31]. A material is brittle if B/G < 1.75,otherwise it behaves as a ductile system. From the previous tables, it is evident that the four systems are brittle at ambient pressure. Even when increasing pressure, Pugh’s ratio increases gradually but the brittle behavior is still dominant in all four alloys. The negative Cauchy pressure
2.2.1 Electronic properties
The band structures of the semiconducting LiAlX alloys are presented in Fig. 6, where an indirect band gap characterizes the four Li-based alloys, with the valance band maximum at Γ point and the conduction band minimum at X point. The resulting band gaps calculated for LiAlC, LiAlSi and LiAlGe with either GGA or mBJ-GGA are very close to other theoretical data [35]. For LiAlSn, no theoretical or experimental data were available to compare the structural parameters of this compound.
Along Γ—X direction, the highest band gap was obtained for LiAlC with a value of 0.81 eV and 1.64 eV using GGA and mBJ-GGA, respectively, at zero pressure. When replacing C by Si, Ge or Sn, we found small narrow gaps that almost touch both valence and conduction bands. These gaps were corrected using the modified Becke-Johnson approximation and the estimated values were in the proximity of 0.5 eV.
LiAlX compounds possess three types of transitions: two of which are indirect along Γ—X and Γ—L directions and the last one is direct along Γ— Γ direction. Unlike the lattice constant, we notice that the Γ— X band gaps decrease when going from LiAlC to LiAlSn alloy at zero pressure.
Studying the influence of external pressure on the gap of the crystal structure was also one of the main interests in this paper. For this reason, we have tabulated the variations of the band gaps within different values of pressures for the four LiAlX (X = C, Si, Ge, Sn) compounds in , Tables II, III, IV and V, respectively. From these Tables, we notice that these half-Heuslers are characterized by an indirect band gap that does not change for the pressure values considered. Similarly, to LiAlSi, the direct band gap of LiAlSn and LiAlGe compounds is the largest among the three possible transitions at ambient conditions. However, the Γ— Γ band gap decreases gradually in these alloys with an increase of pressure. For LiAlC, the values of the band gap along Γ — L direction are the widest at ambient conditions.
LiAlC | ||||||
---|---|---|---|---|---|---|
GGA | GGA-mBJ | |||||
Γ — X | Γ — Γ | Γ — L | Γ — X | Γ — Γ | Γ — L | |
0 | 0.81 | 3.33 | 4.76 | 1.64 | 4.19 | 5.61 |
10 | 0.80 | 3.79 | 5.23 | 1.39 | 4.36 | 5.82 |
20 | 0.78 | 4.29 | 5.75 | 1.37 | 4.85 | 6.33 |
30 | 0.77 | 4.72 | 6.19 | 1.64 | 5.60 | 7.05 |
40 | 0.76 | 5.09 | 6.33 | 1.63 | 5.98 | 7.20 |
50 | 0.74 | 5.42 | 6.45 | 1.62 | 6.31 | 7.33 |
60 | 0.73 | 5.72 | 5.56 | 1.33 | 6.27 | 7.17 |
70 | 0.72 | 5.99 | 6.66 | 1.32 | 6.54 | 7.28 |
80 | 0.70 | 6.24 | 6.76 | 1.60 | 7.14 | 7.67 |
90 | 0.69 | 7.47 | 6.84 | 1.30 | 7.02 | 7.47 |
LiAlSn | ||||||
---|---|---|---|---|---|---|
GGA | GGA-mBJ | |||||
Γ — X | Γ — Γ | Γ — L | Γ — X | Γ — Γ | Γ — L | |
0 | 0.16 | 0.080 | 0.35 | 0.55 | 1.32 | 0.77 |
10 | 0.0 | 1.50 | 0.93 | 0.30 | 1.78 | 1.30 |
20 | 0.04 | 1.62 | 1.26 | 0.09 | 1.79 | 1.67 |
30 | 0.03 | 1.47 | 1.54 | 0.03 | 1.84 | 1.99 |
The energy gap along Γ — X direction decreases for all four compounds with increasing pressure. As for Γ — Γ and Γ — L directions, the band gaps increase for LiAlC, LiAlGe and LiAlSn with increasing pressure (see , Table II, IV and V, respectively) and decrease for the LiAlSi alloy (see Table III).
2.2.2 Thermoelectric properties
The Lithium based half-Heusler ompounds are known to be very promising in converting heat into power due to their large figure of merit (ZT) [36], which can be expressed as the following:
For this particular reason, we have estimated the Seebeck coefficient S,f igure of merit as well as the power factor PF for these LiAlX (X = C, Si, Ge, Sn) compounds, as illustrated in Fig. 7 and 8, respectively.
In the considered range of chemical potential, the Seebeck curve shows almost equal values in both negative and positive regions, meaning that both types of charge carriers may be involved. We notice as well that the large values of thermopower are obtained in low ranges of chemical potential (Fig. 8). The greatest value of the Seebeck reached 2755 micro
When the temperature increases from 300 K to 900 K, the maximum value of thermopower in all four compounds decreases dramatically for both n- and p-type regions.
In order to predict the amount of power generated by the current Li-based materials, we have estimated the variation of the power factor within temperature as can be seen in Fig. 8. The curve of the power factor shows a linear increment with increase of temperature for all four compounds. At 1000 K, LiAlX (X = C, Si, Ge, Sn) compounds show the maximum values 1.18 × 1012 WK-2m-1s-1, 7.5 × 1011 WK-2m-1s-1, 7.2 × 1011 WK-2m-1s-1 and 6.2 × 1011 WK-2m-1s-1, respectively.
The effective thermal conductivity has a contribution from two parameters namely electronic thermal conductivity k 𝑒 and lattice thermal conductivity
where A is a physical constant (A ~ 3.1 × 10-6), γ is the Grüneisen parameter θD the Debye temperature, V is the volume per atom, n is the number of atoms in the primitive unit cell, and m is the average mass of the atoms in the crystal.
The thermoelectric efficiency of the studied materials was also measured within temperature variation as can be seen in Fig. 7. An exponential increasing behavior is observed in the ZT curve of all compounds. The figure of merit increases sharply with increase of temperature until 700 K, where it shows a slow growth. The maximum values of figure of merit were found to be 0.73 for LiAlC and LiAlGe, and 0.72 for LiAlSi and LiAlSn alloys, respectively. These values correspond to a temperature in the range of 800 to 1000 K.
3 Conclusion
Studying the behavior of the four Li-based systems gave us the opportunity to discover their fascinating properties that might be very useful in different optical and thermoelectric applications. We have first illustrated the band structure and clarified the semiconducting nature of these compounds. Then, we have studied their mechanical stability and calculated the corresponding elastic parameters. All systems remained stable in a moderate range of pressure changes. We have then estimated band gaps of the three possible transitions in each compound under applied forces. LiAlX alloys showed no change of the natural indirect Γ - X gap within pressure variations. In the last part of our investigations, we have been concerned by the thermoelectric efficiency of these Li-based alloys. The Seebeck coefficient of these semiconductors showed that a p-n type of conduction is possible. We found that LiAlX alloys exhibit high values of power factor and figure of merit, which proves their efficiency in the thermoelectric device applications.
Author contributions
Khetir Mohammed: Investigation; resources. Amina Maafa: writing-original draft. Boukabrine Fouzia: Methodology; software; supervision; writing-original draft. Habib Rozale: Supervision. Asma Bouabça: Methodology; software. Abbes Chahed: Methodology; software