FEATURES OF SIMULATION OF CHARACTERISTICS OF THERMOMETRIC MATERIAL Lu1-xZrxNiSb

2021;
: pp. 12-17
1
Lviv Polytechnic National University
2
Lviv Polytechnic National University
3
Lviv Polytechnic National University
4
Lviv Polytechnic National University
5
Ivan Franko National University of Lviv
6
Ivan Franko National University of Lviv
7
Ivan Franko National University of Lviv
8
Leibniz Institute for Solid State Research

The results of modeling the thermometric characteristics of the semiconductor solid solution Lu1–xZrxNiSb, which is a promising thermometric material for the manufacture of sensitive elements of thermoelectric and electro resistive thermocouples, are presented. Modeling of the electronic structure of Lu1–xZrxNiSb was performed by the Korringa–Kohn–Rostoker (KKR) method in the approximation of coherent potential and local density and by the full-potential method of linearized plane waves (FLAPW). KKR simulations were performed using the AkaiKKR software package in the local density approximation for the exchange-correlation potential with parameterization Moruzzi, Janak, Williams in the semi-relativistic one taking into account the spin-orbit interaction. The implementation of the method in the Elk software package was used to perform FLAPW calculations. To check the limits of the existence of the thermometric material Lu1–xZrxNiSb, both methods were used to calculate the change in the values of the period of the unit cell a(x) in the range x = 0 – 1.0. It is shown that there is an agreement between the change in the values of a(x) Lu1–xZrxNiSb calculated by the FLAPW method and the results of experimental studies. The obtained result indicates higher accuracy of modeling of structural parameters Lu1–xZrxNiSb by the FLAPW method in comparison with the KKR method. To study the possibility of obtaining thermometric material Lu1–xZrxNiSb and to establish the limits of its existence in the form of a continuous solid solution, modeling of thermodynamic characteristics in the approximation of harmonic oscillations of atoms within the theory of DFT density functional for a hypothetical solid solution Lu1–xZrxNiSb, x = 0 – 1.0. The change in the values of the enthalpy of mixing ΔH and the total energy E Lu1–xZrxNiSb, x = 0–1.0, allows us to state that the thermometric material exists in the form of a solid substitution solution in the concentration range x = 0 – < 0.20. At higher concentrations of impurity Zr atoms, x > 0.20, stratification occurs (spinoidal phase decay) and thermometric material does not exist. To understand the mechanisms of electrical conductivity of the thermometric material Lu1–xZrxNiSb, the methods of entry of impurity Zr atoms into the matrix of the basic semiconductor p-LuNiSb and their occupation of different crystallographic positions, as well as the presence of vacancies in them, were investigated. For this purpose, its electronic structure was modeled for different variants of the spatial arrangement of atoms and the presence of vacancies in crystallographic positions. It is shown that the most acceptable results of experimental studies are the model of the electronic structure of p-LuNiSb, which assumes the presence of vacancies in the crystallographic positions of 4a Lu atoms (~0.005) and 4c Ni atoms (~0.04). In this model of the spatial arrangement of atoms and the presence of vacancies at positions 4a and 4c, the LuNiSb compound is a semiconductor of the hole-type conductivity, in which the Fermi level eF is located near the level of the valence band eV. The kinetic characteristics of the semiconductor thermometric material Lu1–xZrxNiSb, in particular, the temperature dependences of the resistivity ρ(T, x) and the thermopower coefficient α(T, x) are modeled. It is established that at the lowest concentrations of impurity atoms Zr the Fermi level eF Lu1–xZrxNiSb passes from the bandgap to the conduction band eС. This is indicated by the negative values of the thermopower coefficient α(T, x) and the metallic conductivity type Lu1–xZrxNiSb. This changes the type of main current carriers from holes to electrons.

[1] L. Romaka, Yu. Stadnyk, V.A. Romaka, P. Klyzub, V. Krayovskyy, V. Pashkevych, A. Horyn, P. Garanuk, “Synthesis and electrical transport properties of Er1-xScxNiSb semiconducting solid solution”,.Phys. Chem. Sol. St., no. 22(1), p.146–152, 2021.
https://doi.org/10.15330/pcss.22.1.146-152
[2] M. Schruter, H. Ebert, H. Akai, P. Entel, E. Hoffmann, G.G. Reddy, “First-principles investigations of atomic disorder effects on magnetic and structural instabilities in transition-metal alloys”, Phys. Rev. B, vol. 52, p.188–209, 1995.
https://doi.org/10.1103/PhysRevB.52.188
[3] V. Moruzzi, J. Janak, A. Williams. Calculated Electronic Properties of Metals. NY, Pergamon Press, 1978.
[4] All-electron full-potential linearised augmented-plane wave (FP-LAPW) code. [Online]. Available: http://elk. sourceforge.net
[5] K. Momma, F. Izumi, “VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data”, J. Appl. Crystallogr., vol. 41, p.653–658, 2008.
https://doi.org/10.1107/S0021889808012016 

[6] I. Wolanska, K. Synoradzki, K. Ciesielski, K. Zaleski, P. Skokowski, D. Kaczorowski, Mater. Chem. Phys., vol.227, p.29–36, 2019. doi.org/10.1016/j.matchemphys.2019.01.056

[7] K. Synoradzki, K. Ciesielski, I. Veremchuk, H. Borrmann, P. Skokowski, D. Szymanski, Y. Grin, D. Kaczorowski, Materials, vol.12, p.1723–1730, 2019. doi:10.3390/ma12101723

[8] I. Karla, J. Pierre, R.V. Skolozdra, J. Alloys Compd., vol.265, p.42–50, 1998. doi.org/10.1016/S0925-8388(97)00419-2

[9] B. Shklovskii, A. Efros, Electronic Properties of Doped Semiconductors. Springer-Verlag, 1984. dx.doi.org/ 10.1007/978-3-662-02403-4