Представлено результати моделювання та експериментальні дослідження структурних, магнітних, електрокінетичних та енергетичних властивостей термометричного матеріалу Hf1-xNbxNiSn, а також функцій перетворення чутливих елементів термоелектричного термометра на його основі за температур 4.2–1000 К. Для випадку упорядкованого варіанту кристалічної структури термометричного матеріалу встановлено одночасне генерування донорних та акцепторних станів у забороненій зоні εg напівпровідника. Виявлено залежність між просторовим розташуванням атомів у вузлах елементарної комірки Hf1-xNbxNiSn та механізмами електропровідності, що дозволяє визначити умови синтезу матеріалів з максимальною ефективністю перетворення теплової енергії в електричну. Показано, що за концентрацій Hf1-xNbxNiSn, 0–0.02, переважно має місце заміщенням у позиції 4а атомів Hf (5d26s2) на атоми Nb (4d45s1), що генерує донорні стани. У районі концентрацій Hf1-xNbxNiSn, х=0.02–0.05, переважно відбувається заміщення атомів Ni (3d84s2) у позиції 4а на атоми Nb, що генерує акцепторні стани, а за концентрацій 0.05<х заміщенням атомів Hf на атоми Ni генерує додаткові донорні стани. Досліджений термометричний матеріал Hf1-xNbxNiSn є перспективним для виготовлення чутливих елементів термоелектричних термометрів. Проведено моделювання функцій перетворення термоелектричної пари Hf0.99Nb0.01NiSn-ПР 13, термоелектроди якої виготовлено з дослідженого термометричного матеріалу та платинороїду (ПР13) (додатна вітка).
Термометричні матеріали Hf1-xNbxNiSn, х=0.01–0.10, для виготовлення чутливих елементів перетворювачів температури отримано сплавленням шихти компонентів в електродуговій печі з вольфрамовим електродом (катод) в атмосфері очищеного аргону під тиском 0,1 кПа на мідному охолоджуваному поді (анод). Як гетер використовували попередньо сплавлений губчатий титан. Термічне оброблення сплавів Hf1-xNbxNiSn полягало у гомогенізуючому відпалюванні за температури 1073 К. Відпал зразків проводили на протязі 720 год. у вакуумованих кварцевих ампулах (до 1,0 Па) у муфельних електропечах з регулюванням температури з точністю ±10 ℃. Масиви дифракційних даних отримано на порошковому рентгенівському дифрактометрі STOE STADI-P (Cu Kα1-випромінювання). За допомогою пакета програм Fullprof розраховано структурні характеристики зразків Hf1-xNbxNiSn. Хімічний та фазовий склади зразків контролювали за допомогою металографічного аналізу (скануючий електронний мікроскоп Tescan Vega 3 LMU).
Для оптимізації параметрів кристалічної та електронної структур, енергетичних та кінетичних властивостей термометричних матеріалів Hf1-xNbxNiSn проведено розрахунки у межах теорії функціоналу густини (Density Functional Theory, далі DFT) за допомогою пакета програм Vienna Ab initio Simulation Package VASP v. 5.4.4 з потенціалами типу PAW. Розрахунок електронних кінетичних коефіцієнтів проводили з використанням коду Exciting (метод FLAPW – Full Potential Linearized Augumented Plane Waves) шляхом розв’язання лiнеаризованого рiвняння Больцмана у наближеннi сталого часу релаксацiї. Моделювання розподілу густини електронних станів (DOS) проводили методом Korringa-Kohn-Rostoker (KKR). Моделювання функцій перетворення термоелектричного термометра за температур 4.2–1000 К проведено методом FLAPW, пакет програм Elk, а реперними токами слугували результати експериментальних вимірювань.
Рентгеноструктурні дослідження кристалічної структури зразків термометричного матеріалу Hf1-xNbxNiSn, х=0–0.10, встановили, що дифрактограми індексуються у структурному типі MgAgAs і на них відсутні відбиття інших фаз. На основі отриманих дифракційних картин було розраховано зміну періоду елементарної комірки а(х) Hf1-xNbxNiSn, х=0–0.10. Очікували, що період комірки а(х) буде зменшуватися при збільшенні концентрації атомів Nb (rNb=0.146 нм), оскільки його атомний радіус є меншим, ніж атома Hf (rHf=0.158 нм). Однак, зміна періоду комірки а(х) Hf1-xNbxNiSn, х=0–0.10, є далекою від очікуваної і носить складний характер. Так, за концентрацій атомів Nb, х=0–0.02, значення періоду комірки а(х), як і прогнозували, стрімко зменшуються. Однак за концентрацій х=0.02–0.05 спостерігаємо такий самий стрімкий ріст залежності а(х), яка в районі концентрації х≈0.05 проходить через максимум і потім знову спадає. Отримані результати зміни періоду комірки а(х) Hf1-xNbxNiSn, х=0–0.10, вказують на складні структурні перетворення, які є наслідком одночасних змін у кількох кристалографічних позиціях напівпровідникового термометричного матеріалу, що зумовить перерозподіл густини електронних станів.
Моделювання розподілу густини електронних станів (DOS) для упорядкованого варіанту кристалічної структури Hf1-xNbxNiSn, х=0–0.10, показує положення рівня Фермі εF та ширину забороненої зони εg. Так, якщо в n-HfNiSn рівень Фермі εF знаходиться у забороненій зоні εg біля краю зони провідності εС, то вже за найменшої концентрації атомів Nb, х=0.005, рівень Фермі εF перетне край зони провідності εС: відбудеться перехід провідності діелектрик-метал.
Дослідження питомої магнітної сприйнятливості χ(х) термометричного матеріалу Hf1-xNbxNiSn, х=0–0.10, показали, що базовий напівпровідник n-HfNiSn є слабким діамагнетиком, на що вказують від’ємні значення питомої магнітної сприйнятливості χ за кімнатної температури. Легування n-HfNiSn атомами Nb робить напівпровідник Hf1-xNbxNiSn парамагнетиком Паулі, в якого питома магнітна сприйнятливість визначається виключно електронним газом і є пропорційною густині станів на рівні Фермі g(εF)(х). Дослідження встановили, що на ділянці концентрацій х=0–0.02 відбувається стрімке збільшення значень χ(х) Hf1-xNbxNiSn, викликане генеруванням донорних станів та збільшенням концентрації вільних електронів при заміщенні атомів Hf на атоми Nb. За більших концентрацій атомів Nb швидкість зміни магнітної сприйнятливості χ(х) Hf1-xNbxNiSn зменшується через появу та збільшення концентрації акцепторних станів, які захоплюють вільні електрони, зменшуючи їхню концентрацію. Отримані результати показують, що зовнішнє магнітне поле з напруженістю Н≤10 кГс не впливає на термометричні характеристики термоперетворювачів з даного матеріалу.
Характер зміни значень питомого електроопору ρ та коефіцієнта термо-ерс α Hf1-xNbxNiSn як зі зміною температури, так і концентрації атомів Nb узгоджується з висновками, зробленими на основі структурних досліджень та моделювання електронної структури. Оскільки електроопір r(Т,х) Hf1-xNbxNiSn практично лінійно збільшуються з ростом температури, що є результатом металізації електропровідності, даний матеріал є непридатним для отримання чутливих елементів термометрів опору.
Температурні залежності коефіцієнта термо-ерс a(Т,х) Hf1-xNbxNiSn показують, що основними носіями електричного струму термометричного матеріалу за всіх температур є вільні електрони. На це вказують від’ємні значення коефіцієнта термо-ерс a за всіх концентрацій та досліджених температур. Показано, що з ростом концентрації атомів Nb зменшується абсолютна зміна значень коефіцієнта термо-ерс a(Т,х) Hf1-xNbxNiSn. Причиною цього є зменшення ширини забороненої зони εg термометричного матеріалу. За характером поведінки коефіцієнта термо-ерс a(х,Т) Hf1-xNbxNiSn можна встановити концентрацію атомів Nb, за якої зміна значень коефіцієнта термо-ерс a буде найбільшою. У даному випадку ця умова досягається за концентрації напівпровідникового термометричного матеріалу Hf0.99Nb0.01NiSn, який став основою для формування електроду (від'ємної вітки) термоелектричного термометра.
- R. Marazza, R. Ferro, G. Rambaldi, Some phases in ternary alloys of titanium, zirconium, and hafnium, with a MgAgAs or AlCu2Mn type structure. J. Less-Common Met. Vol. 39,
p. 341–348, 1975. DOI: https://doi.org/10.1016/0022-5088(75)90207-6. - V.A. Romaka, Yu. Stadnyk, V. Krayovskyy, L. Romaka, O. Guk, V.V. Romaka, M. Mykyuchuk, A. Horyn. The latest heat-sensitive materials and temperature transducers. Lviv Polytechnic Publishing House, Lviv, 2020. DOI: https://opac.lpnu.ua/bib/1131184. [in Ukrainian].
- V.V. Romaka, P. Rogl, L. Romaka, Yu. Stadnyk, A. Grytsiv, O. Lakh, V. Krayovskii, Peculiarities of structural disorder in Zr- and Hf-containing Heusler and half-Heusler stannides, Intermetallics 35, 45 (2013). DOI: http://dx.doi.org/10.1016/j.intermet.2012.11.022
- B.I. Shklovskii, A.L. Efros. Electronic Properti es of Doped Semiconductors. Springer-Verlag, NY, 1984. DOI: http://doi10.1007/978-3-662-02403-4.
- T. Roisnel, J. Rodriguez-Carvajal. WinPLOTR: a Windows Tool for Powder Diffraction Patterns analysis, Mater. Sci. Forum, Proc. EPDIC7, vol.378–381, p.118–123, 2001. DOI:https://doi.org/10.4028/www.scientific.net/MSF.378-381.118.
- G. Kresse, J. Hafner. Ab initio molecular dynamics for liquid metals. Phys. Rev., B Vol. 47, p. 558–561, 1993. DOI: https://doi.org/10.1103/PhysRevB.47.558
- H.J. Monkhorst, J.K. Pack, Special points for Brillouin- zone integrations, Phys. Rev. B. Vol. 13, p. 5188–5192, 1976. DOI: https://doi.org/10.1103/PhysRevB.13.5188.
- K. Okhotnikov, T. Charpentier, S. Cadars, Supercell program: a combinatorial structure-generation approach for the local-level modeling of atomic substitutions and partial occupancies in crystals, J. Cheminform. Vol. 8(17), 1–15, 2016. DOI: 10.1186/s13321-016-0129-3.
- A. Gulans, S. Kontur, C. Meisenbichler, D. Nabok, P. Pavone, S. Rigamonti, S. Sagmeister, U. Werner, C. Draxl, Exciting – a full-potential all-electron package implemen- ting density-functional theory and many-body perturbation theory, J. Phys.: Condens Matter. Vol. 26, p. 363202, 1–24, 2014. DOI: 10.1088/0953-8984/26/36/363202.
- T.J. Scheidemantel, C. Ambrosch-Draxl, T. Thonhauser, H.V. Badding, and J.O. Sofo, Transport coefficients from first- principles calculations, Phys. Rev. B, Vol. 68, p. 125210, 2003. DOI: https://doi.org/10.1103/PhysRevB.68.125210.
- All-electron full-potential linearised augmented-plane wave (FP-LAPW) code – http://elk.sourceforge.net.
- M. Schruter, H. Ebert, H. Akai, P. Entel, E. Hoffmann, G.G. Reddy. Phys. Rev. B, vol.52, p.188–209, 1995. https://doi.org/10.1103/PhysRevB.52.188
- V.Moruzzi, J. Janak, A.Williams. Calculated Electronic Properties of Metals. NY, Pergamon Press, 1978. DOI: https://doi.org/10.1016/B978-0-08-022705-4.50002-8.