Kislorod vakansiyali tantal pentoksid (Ta₂O₅), kristall panjaradagi ayrim kislorod atomlari chiqarib tashlangan va shu nuqsonlar tufayli elektron zona tuzilishi oʻzgarishi mumkin boʻlgan keng taqiqlangan zonali oksid yarimoʻtkazgich tizimidir. Ushbu tadqiqot γ₁, B va δ strukturaviy fazalarida %10 kislorod vakansiyasi hosil qilinganda elektron zonalar qanday oʻzgarishini zichlik funksionali nazariyasi (Density Functional Theory, DFT) yordamida oʻrganadi. Hisob-kitoblarda ayrim yoʻnalishlarda Fermi energiyasi yaqinida juda tekis elektron zonalar hosil boʻladi va bu zonalarning ayrim sohalarida elektronlarning zona-samarali massasi erkin elektron massasining \(10^3\) tartibidan oshib, ayrim nuqtalarda \(10^4\) tartibiga yetadi. Mualliflar bu xususiyatni ehtiyotkorlik bilan “heavy-fermion-like”, ya’ni ogʻir fermionga oʻxshash xatti-harakat deb ataydilar.
Tadqiqotning muhim jihati shundaki, u an’anaviy ogʻir fermion tizimlaridan farq qiladigan materiallar sinfini oʻrganadi. Klassik ogʻir fermion birikmalari koʻpincha noyob yer yoki aktinid elementlarini oʻz ichiga oladigan, kuchli elektron korrelyatsiyalari va past harorat fizikasi hal qiluvchi boʻlgan metall yoki tor taqiqlangan zonali tizimlardir. Bu yerdagi Ta₂O₅ tadqiqoti esa kislorod vakansiyalari hosil qilgan tekis zonalardan kelib chiqadigan juda katta zona-samarali massalarni keng taqiqlangan zonali oksid tizimida oʻrganadi. Tadqiqot Kondo effektini, solishtirma issiqlik koeffitsiyentini, past haroratdagi elektr qarshiligini yoki ogʻir fermion tashuvchilarini bevosita tajribada oʻlchamaydi; shuning uchun “ogʻir fermionga oʻxshash” iborasining chegarasi saqlanishi lozim.
Hisob-kitoblarda VASP dasturi, PAW potensiallari, PBE shaklidagi GGA almashinuv-korrelyatsiya yondashuvi va 600 eV tekis-toʻlqin energiya kesimi ishlatilgan. γ₁ fazasi uchun 1×1×2 superyacheyka va 10×10×4 Monkhorst–Pack k-nuqta toʻri; B fazasi uchun 2×4×4, δ fazasi uchun 6×6×12 k-nuqta toʻri qoʻllangan. γ₁ fazasidagi zona dispersiyasini tekshirish uchun HSE06 gibrid funksionali bilan qoʻshimcha hisob-kitoblar bajarilgan.
γ₁ fazasida %10 tip-2 kislorod vakansiyasi sharoitida GGA taqiqlangan zona kengligi 1,1 eV deb topilgan. Γ–X, Γ–Y va Γ–Z yoʻnalishlaridagi tekis zona kengliklari GGA bilan mos ravishda 96,3, 41,3 va 23,4 meV; HSE bilan 110,0, 40,0 va 20,0 meV deb xabar qilingan. B fazasida tekis zona kengliklari Γ–X va Γ–Y uchun 29,8 va 76,6 meV; δ fazasida 76,2 va 62,1 meV. B fazasining hisoblangan taqiqlangan zona kengligi 0,36 eV boʻlsa, %10 vakansiyali δ fazasi metall xususiyat koʻrsatadi.
γ₁ fazasidagi yuqori vakansiyali tuzilma faqat statik geometriya optimallashtirish bilan cheklanib qolmasligi uchun 30 K da, NVT ansamblida, 2×2×1 superyacheyka bilan 5 ps ab initio molekulyar dinamika (AIMD) simulyatsiyasi ham oʻtkazilgan. Umumiy energiyada sezilarli siljish kuzatilmagan, boshlangʻich va AIMD dan keyingi XRD choʻqqi holatlari katta darajada saqlangan va Fermi sathi yaqinidagi tekis zonalar simulyatsiya oxirida ham koʻringan. Biroq bu tekshiruv faqat γ₁ fazasining koʻrsatilgan past harorat va qisqa simulyatsiya vaqtidagi strukturaviy/sonli barqarorligiga taalluqlidir; u xona haroratidagi uzoq muddatli tajribaviy barqarorlikning dalili emas.
Mualliflar nazariy bashoratni tajribada sinash uchun kuchli magnit maydon ostida elektron siklotron rezonansini taklif qiladilar. Masalan, 100 GHz rezonans chastotasi faraz qilinsa, \(m^\*\sim10^2\)–\(10^4m_0\) oraligʻi taxminan 0,36–35,7 T magnit maydonni talab qiladi. Ammo maqolada bu tajriba bajarilmagan; keltirilgan qiymatlar DFT dan hisoblangan samarali massalarga asoslangan tajriba dizayni baholaridir.
Ogʻir fermion tushunchasining ushbu maqoladagi konteksti
An’anaviy ogʻir fermion fizikasida zaryad tashuvchilarning samarali massasi erkin elektron massasidan yuzlab, hatto minglab marta katta boʻlishi mumkin. Manba kirish qismida CeAl₃ kabi klassik tizimlarda past haroratdagi Fermi-suyuqlik xatti-harakatini misol qilib keltiradi va elektr qarshilik uchun \(\Delta\rho(T)\propto T^2\), issiqlik sigʻimi uchun \(C(T)=\gamma T\) munosabatlarini eslatadi. Biroq Ta₂O₅ tadqiqoti bu termodinamik yoki transport kattaliklarini hisoblamaydi ham, oʻlchamaydi ham.
Ushbu tadqiqotda ogʻir fermionga oʻxshash talqinning bevosita asosi kislorod vakansiyalari Fermi energiyasi yaqinida hosil qiladigan tekis zonalar va bu zonalarning juda kichik egrilik koʻrsatadigan sohalaridan hisoblangan juda katta \(|m^\*|\) qiymatlaridir. Shu sababli tadqiqot klassik noyob-yer ogʻir fermion birikmalarining barcha koʻp-jismli fizikasi Ta₂O₅ da namoyon boʻlganini da’vo qilmaydi.
Nega Tekis Zona Elektronning Samarali Massasini Oshiradi?
Tekis elektron zonada energiya \(E\), kristall impulsi \(k\) oʻzgarganiga qaramay juda kam oʻzgaradi; shu sababli zona egriligining kattaligi kamayadi. Zona nazariyasida samarali massa egrilikning teskarisi bilan bogʻliq boʻlganligi sababli kichik \(|d^2E/dk^2|\) qiymati juda katta \(|m^\*|\) ni hosil qiladi. Shu bois tadqiqotdagi favqulodda samarali massa choʻqqilari kislorod vakansiyalari hosil qilgan juda tekis zona sohalarida paydo boʻladi.
Tushuntiruvchi fon izohi: Maqola samarali massa munosabatini alohida raqamlangan tenglama sifatida yozmaydi, biroq \(E(k)\) ning birinchi va ikkinchi hosilalaridan tezlik va samarali massani hisoblaganini aniq bildiradi. Bitta asosiy kristall yoʻnalishi boʻylab ishlatiladigan standart zona-nazariyasi munosabati:
\[ v(k)=\frac{1}{\hbar}\frac{dE}{dk}, \]
\[ m^\*(k)= \frac{\hbar^2}{ \dfrac{d^2E}{dk^2} }. \]
koʻrinishidadir. Birinchi hosila zona qiyaligi guruh tezligini qanday belgilashini; ikkinchi hosila esa zona egriligi samarali massaga qanday aylanishini koʻrsatadi. Bu ikki tenglama Verianla usulni tushuntirish uchun keltirgan standart zona fizikasi fonidir; maqolaning yangi tadqiqot natijasi sifatida taqdim etilmaydi.
Formulada \(\hbar\) keltirilgan Plank doimiysidir. \(dE/dk\) energiya zonasining qiyaligini, \(d^2E/dk^2\) esa egriligini koʻrsatadi. Egrilik nolga yaqinlashganda samarali massaning mutlaq qiymati juda kattalashadi. Egrilik ishorasi oʻzgarganda \(m^\*\) ishorasi ham oʻzgarishi mumkin; shu sababli maqola ogʻir-massa sohalarini baholashda \(|m^\*|\) dan foydalanadi. Manfiy zona-samarali massa zarrachaning haqiqiy massasi fizik jihatdan manfiy degani emas; u zona egriligining ishorasini ifodalaydi.
Kislorod vakansiyasi turlari qanday hosil qilindi?
Ta₂O₅ ichidagi kislorod atomlari Ta atomlari bilan bogʻlanishlar soniga koʻra ikki koordinatsiya sinfiga ajratilgan. Manbada tip 1 kislorod vakansiyasi, ikki Ta atomi bilan koordinatsiyalangan kislorodning chiqarilishi; tip 2 kislorod vakansiyasi esa uch Ta atomi bilan koordinatsiyalangan kislorodning chiqarilishi sifatida ta’riflanadi.
- γ₁ fazasi: 1×1×2 superyacheyka yaratilgan va bitta uch-koordinatsiyalangan kislorod chiqarilib tip-2 vakansiya hosil qilingan.
- B fazasi: boshlangʻich yacheykadan ikkita ikki-koordinatsiyalangan kislorod chiqarilgan; vakansiya turi 1.
- δ fazasi: boshlangʻich yacheykadan bitta ikki-koordinatsiyalangan kislorod chiqarilgan; vakansiya turi 1.
Bu uch tuzilma tadqiqot davomida %10 kislorod vakansiyasi konsentratsiyasi sifatida koʻrib chiqilgan va kislorod chiqarilgandan soʻng nuqsonli tuzilmalarning ion holatlari toʻliq relaksatsiya qilingan.
DFT hisoblash shartlari
| Parametr | Manbada ishlatilgan qiymat / usul | Ilmiy roli |
|---|---|---|
| Dastur | VASP | Birinchi prinsiplardan elektron tuzilish hisob-kitoblari |
| Asosiy usul | Density Functional Theory (DFT) | Elektronlarning umumiy energiyasi va zona tuzilishi |
| Elektron–ion oʻzaro ta’siri | PAW | Ion yadrolari va valent elektronlarning ifodalanishi |
| Almashinuv/korrelyatsiya | PBE-GGA | Asosiy DFT hisoblash yondashuvi |
| Tekis-toʻlqin energiya kesimi | 600 eV | Bazis toʻplamining maksimal kinetik energiya chegarasi |
| B fazasi k-toʻri | 2×4×4 | Brillouin zonasini namunalash |
| δ fazasi k-toʻri | 6×6×12 | Brillouin zonasini namunalash |
| γ₁ fazasi | 1×1×2 superyacheyka; 10×10×4 k-toʻri | %10 tip-2 kislorod vakansiyasini modellashtirish |
| k-toʻrini hosil qilish | Monkhorst–Pack | Davriy teskari fazoni namunalash |
| Qoʻshimcha zona tekshiruvi | HSE06 | Ayniqsa γ₁ fazasidagi zona dispersiyasini GGA bilan taqqoslash |
γ₁ fazasidagi zona tuzilishi
%10 tip-2 kislorod vakansiyali γ₁-Ta₂O₅ uchun GGA hisob-kitobi 1,1 eV taqiqlangan zona kengligini beradi. Fermi energiyasi yaqinida tekis zona hosil boʻladi. Tadqiqot bu tekis zonani Γ–X, Γ–Y va Γ–Z yuqori simmetriya yoʻnalishlarida alohida tahlil qilgan.
| Yoʻnalish | GGA zona kengligi | HSE zona kengligi |
|---|---|---|
| Γ–X | 96,3 meV | 110,0 meV |
| Γ–Y | 41,3 meV | 40,0 meV |
| Γ–Z | 23,4 meV | 20,0 meV |
Kichikroq zona kengligi tekshirilayotgan impuls yoʻli boʻylab energiya kamroq oʻzgarishini va zonaning tekisroq ekanini bildiradi. γ₁ fazasida ushbu uch yoʻnalish ichida xabar qilingan eng tor zona Γ–Z yoʻnalishidadir.
Nega E(k) egri chizigʻi boʻlakma-boʻlak fit qilindi?
Samarali massa \(E(k)\) egri chizigʻining ikkinchi hosilasiga bogʻliq boʻlganligi sababli, fittingdagi kichik tebranishlar ham hosilalarda kuchayishi mumkin. Shuning uchun maqola \(E(k)\) ma’lumotlarini bitta global polinom bilan fit qilish oʻrniga bir-birini qoplaydigan boʻlimlarga ajratib, mahalliy fittingdan foydalangan.
A ilovadagi misolda ishlatilgan polinom:
\[ E(k) = a k^{10} + b k^8 + c k^6 + d k^4 + e k^2 + f. \]
Bu yerda \(a,b,c,d,e,f\) tegishli mahalliy ma’lumot boʻlimiga moslashtirilgan koeffitsiyentlardir. Faqat juft \(k\) darajalari ishlatilgan ushbu ifoda tegishli zona boʻlagini silliq egri chiziq bilan ifodalashni ta’minlaydi. Fit qilingan \(E(k)\) keyin hosila olinib, elektron tezligi va zona egriligi hisoblangan.
Manba namuna sifatidagi γ₁ Γ–Z ma’lumotlarini 0–0.02, 0.01–0.03 va 0.02–0.04 boʻlgan uchta bir-birini qoplaydigan segmentga ajratganini bildiradi. Biroq shu paragrafning oʻzida umumiy \(k\) oraligʻini 0–0.03 deb yozadi. Oxirgi segmentning 0.04 gacha uzayishi bu jumla bilan mos kelmaydi; Verianla qaysi qiymat yozuv xatosi ekanini taxmin qilmaydi.
Samarali massaning Brillouin zonasida taqsimlanishi
Manba tasvirlaridagi sariq sohalar:
\[ |m^\*|\geq 10m_0 \]
shartini; kulrang sohalar esa:
\[ |m^\*|\geq10^2m_0 \]
shartini ifodalaydi. Bu yerda \(m_0\) erkin elektron massasidir.
γ₁ fazasida Γ–X, Γ–Y va Γ–Z yoʻnalishlarining barchasida ayrim impuls sohalarida \(|m^\*|>10^3m_0\) qiymatlari topilgan. Ayrim tor sohalarda qiymat \(10^4m_0\) tartibiga yetadi.
| Yoʻnalish | k koordinatasi | \(\log_{10}|m^\*/m_0|\) |
|---|---|---|
| Γ–X | (0.045, 0, 0) | 4.131 |
| Γ–X | (0.102, 0, 0) | 4.064 |
| Γ–Y | (0, 0.015, 0) | 3.970 |
| Γ–Y | (0, 0.051, 0) | 3.684 |
| Γ–Y | (0, 0.099, 0) | 3.854 |
| Γ–Z | (0, 0, 0.020) | 4.107 |
\(\log_{10}|m^\*/m_0|=4\) qiymati samarali massaning erkin elektron massasidan taxminan \(10^4\) marta katta tartibda ekanini anglatadi. Shu sababli jadvaldagi 4.131, 4.064 va 4.107 kabi qiymatlar bevosita \(10^4\) tartibidagi ogʻir samarali massa sohalarini koʻrsatadi.
B va δ fazalarida nima topildi?
%10 tip-1 vakansiyali B-Ta₂O₅ uchun hisoblangan taqiqlangan zona kengligi 0,36 eV. %10 tip-1 vakansiyali δ-Ta₂O₅ esa hisoblashda aniq taqiqlangan zona koʻrsatmagan va metall xususiyatga ega boʻlgan. Har ikkala fazada ham Fermi sathi yaqinida tekis zonalar kuzatilgan.
| Faza | Γ–X | Γ–Y |
|---|---|---|
| B | 29,8 meV | 76,6 meV |
| δ | 76,2 meV | 62,1 meV |
B va δ fazalarining tekis zonalaridagi ayrim boʻlimlarda ham \(|m^\*|>10^3m_0\) qiymatlari va taxminan \(10^4m_0\) tartibidagi maksimumlar xabar qilingan.
Manba hisobot izohi: B va δ natijalarini bayon qiluvchi asosiy matn GGA hisob-kitoblari haqida soʻz yuritadi va koʻrinadigan grafiklarda faqat GGA/Fit egri chiziqlari mavjud. Bunga qarshi ravishda 4-rasm va 5-rasm tagyozuvlarida “GGA and HSE methods” iborasi uchraydi. Manba bu farqni izohlamaganligi sababli Verianla B va δ fazalari uchun mustaqil HSE sonli natijalari bordek yangi qiymatlar yaratmaydi.
Ushbu Tadqiqotdagi “Ogʻir Fermionga Oʻxshash” Iborasi Nimani Anglatadi?
Tadqiqotda “ogʻir fermionga oʻxshash” kislorod vakansiyalari hosil qilgan tekis elektron zonalarda zona egriligidan hisoblangan elektron samarali massasining ayrim impuls sohalarida \(10^3\)–\(10^4m_0\) tartibiga chiqishini anglatadi. Maqola buni klassik ogʻir fermion birikmalaridagi barcha koʻp-jismli mexanizmlar Ta₂O₅ da tajribaviy tasdiqlangan degan ma’noda ishlatmaydi; transport, solishtirma issiqlik, Kondo fizikasi yoki siklotron rezonansi bevosita oʻlchanmagan.
Bu farq ayniqsa muhim, chunki an’anaviy ogʻir fermion tizimlarida katta samarali massa koʻpincha kuchli elektron korrelyatsiyalarining natijasi sifatida muhokama qilinadi. Ushbu maqolada esa bevosita hisoblangan kattalik kislorod vakansiyali Ta₂O₅ ning bir-zarrachali elektron zona dispersiyasi va undan hosil qilingan zona-samarali massadir.
Mualliflar xulosa boʻlimida, shuningdek, elektron–fonon bogʻlanishi va polaron hosil boʻlishi Ta₂O₅ da samarali massani yanada qayta renormallashtirishi mumkinligini ta’kidlaydilar. Biroq bu mexanizmlar joriy tadqiqotda alohida hisoblanmagan yoki oʻlchanmagan; ular kelajak tadqiqot uchun taklif sifatidadir.
γ₁ fazasining dinamik barqarorligini tekshirish
%10 tip-2 kislorod vakansiyali γ₁ fazasining dinamik barqarorligi uchun ab initio molekulyar dinamika qoʻllangan:
| Shart | Manba qiymati |
|---|---|
| Faza | γ₁-Ta₂O₅ |
| Kislorod vakansiyasi | %10, tip 2 |
| Superyacheyka | 2×2×1 |
| Harorat | 30 K |
| Davomiylik | 5 ps |
| Ansambl | NVT |
| Oʻrtacha umumiy energiya | −504.98576 eV |
| AIMD dan keyingi tuzilmaning 0 K asosiy holat energiyasi | −505.00368 eV |
| Energiya farqi | 0.01792 eV |
NVT — zarrachalar soni \(N\), hajm \(V\) va harorat \(T\) doimiy saqlanadigan kanonik ansambldir. 5 ps davomida umumiy energiya barqaror oʻrtacha qiymat atrofida tebrangan va sezilarli tizimli siljish koʻrsatmagan.
AIMD dan keyin hisoblangan XRD choʻqqi holatlarining boshlangʻich tuzilmaga nisbatan katta darajada saqlanishi kristall karkas ushbu simulyatsiya sharoitlarida jiddiy qayta tuzilishga uchramaganini koʻrsatadi. AIMD ning yakuniy kristall konfiguratsiyasidan hisoblangan zona tuzilishida Fermi sathi yaqinidagi tekis zonalar ham saqlangan va manba toʻrtta tekis zona aniq koʻrinishini bildirgan.
Bu natija 30 K va 5 ps AIMD sharoitidagi γ₁ tuzilmasi uchun sonli dinamik barqarorlik dalilidir. B va δ fazalari uchun xuddi shunday AIMD sinovi xabar qilinmagan. Bundan tashqari, 5 ps/30 K natijasi xona haroratidagi uzoq muddatli tajribaviy barqarorlik yoki sintezlanish imkoniyatining bevosita dalili emas.
Kislorod vakansiyasining hosil boʻlish energiyasi
C ilovada γ₁ fazasidagi %10 kislorod vakansiyasi uchun hosil boʻlish energiyasi:
\[ E_{\mathrm{vf}} = E_{\mathrm{defect}} + 0.5E_{\mathrm{O_2}} - E_{\mathrm{perfect}} - T\Delta S. \tag{A1} \]
munosabati orqali hisoblangan.
Bu yerda:
- \(E_{\mathrm{vf}}\): kislorod vakansiyasining hosil boʻlish energiyasi.
- \(E_{\mathrm{defect}}\): bitta kislorod atomi chiqarilgan nuqsonli superyacheykaning umumiy energiyasi.
- \(E_{\mathrm{perfect}}\): nuqsonsiz superyacheykaning umumiy energiyasi.
- \(E_{\mathrm{O_2}}\): izolyatsiyalangan O₂ molekulasining umumiy energiyasi.
- \(T\Delta S\): harorat va entropiya oʻzgarishi hosil qiladigan termodinamik tuzatish.
Manbada ishlatilgan statik DFT qiymatlari:
\[ E_{\mathrm{defect}} = -126.31649376\ \mathrm{eV}, \]
\[ E_{\mathrm{perfect}} = -136.63040186\ \mathrm{eV}, \]
\[ E_{\mathrm{O_2}} = -4.93232\ \mathrm{eV}. \]
Ushbu qiymatlar bilan 0 K dagi hosil boʻlish energiyasi:
\[ E_{\mathrm{vf}}(0\ \mathrm K) = 7.85\ \mathrm{eV} \]
deb topilgan. 300 K va atrof-muhit bosimi uchun:
\[ T\Delta S = 0.6376\ \mathrm{eV} \]
hadi kiritilganda:
\[ E_{\mathrm{vf}}(300\ \mathrm K) = 7.21\ \mathrm{eV} \]
olingan.
Manba bu qiymat adabiyotda LSR fazasidagi bitta kislorod vakansiyasi uchun xabar qilingan taxminan 5 eV dan taxminan 2 eV yuqori ekanini va buni yuqori vakansiya konsentratsiyasining nuqson energetikasiga ta’siri bilan bogʻlaydi. Bu yerdagi taxminan 5 eV joriy tadqiqotda ayni tuzilma uchun qayta hisoblangan nazorat qiymati emas, balki manbadagi 48-Referensdan olingan taqqoslashdir.
Siklotron Rezonansi Ushbu Nazariy Bashoratni Qanday Sinashi Mumkin?
Siklotron rezonansi magnit maydondagi zaryad tashuvchilarning doiraviy harakat chastotasini elektromagnit nurlanish bilan rezonansga keltirish orqali samarali massani oʻlchashga xizmat qiladigan tajribaviy usuldir. Tadqiqot kislorod vakansiyali Ta₂O₅ da hisoblangan juda katta \(m^\*\) qiymatlarini kelajakda past harorat va kuchli magnit maydon sharoitlarida shu usul bilan sinash mumkinligini taklif qiladi; maqolaning oʻzida bunday oʻlchov bajarilmagan.
Tashuvchining siklotron burchak chastotasi:
\[ \omega_c = \frac{eB}{m^\*}. \tag{1} \]
Bu yerda \(e\) elementar zaryad, \(B\) magnit maydon, \(m^\*\) esa tashuvchining samarali massasidir. Tenglamada samarali massa ortgani sari ayni magnit maydondagi siklotron chastotasi kamayadi. Shu sababli juda ogʻir samarali massalar mikrotoʻlqin yoki uzoq infraqizil chastotalarga mos kelishi mumkin.
Siklotron chastotasi:
\[ \nu_c = \frac{\omega_c}{2\pi} \]
koʻrinishida yoziladi. Manba \(\nu_c=100\ \mathrm{GHz}\) uchun zarur maydonni:
\[ B = \frac{2\pi\nu_cm_0}{e} \left( \frac{m^\*}{m_0} \right) \approx 3.57\times10^{-3} \left( \frac{m^\*}{m_0} \right) \mathrm T \]
koʻrinishida beradi.
\[ m^\*=10^2m_0 \]
uchun taxminan:
\[ B\approx0.36\ \mathrm T \]
va:
\[ m^\*=10^4m_0 \]
uchun taxminan:
\[ B\approx35.7\ \mathrm T \]
talab qilinadi. Manba bu oraliq mavjud uzluksiz yuqori magnit maydon inshootlarining imkoniyatlari doirasida ekanini bildiradi.
Sochilish vaqti nega muhim?
Tashuvchining rezonans harakati sochilish tufayli yoʻqolmasdan kuzatilishi uchun:
\[ \omega_c\tau>1 \]
sharti talab qilinadi. Bu yerda \(\tau\) tashuvchining oʻrtacha erkin sochilish vaqtidir.
Chegaraviy holat:
\[ \omega_c\tau_0=1 \]
deb ta’riflanganda:
\[ \tau_0(k) = \frac{m^\*(k)}{eB}. \tag{2} \]
hosil boʻladi. Samarali massa ortgan sari zarur minimal sochilish vaqti ham ortadi. Boshqacha aytganda, juda ogʻir tashuvchi siklotron orbitasini yetarlicha yaxshi belgilashi uchun sochilmasdan uzoqroq vaqt harakat qilishi kerak.
Kritik namuna uzunligi
Manba zarur minimal xarakterli uzunlikni:
\[ L_0(k) = \tau_0(k)|v(k)|. \tag{3} \]
deb baholaydi. \(v(k)\) — DFT zona tuzilishidan olingan elektron tezligi.
Γ–X yoʻnalishi uchun \(L_0\) haqiqiy fazodagi X kristall yoʻnalishida rezonans rivojlanishi uchun zarur minimal namuna qalinligi sifatida talqin qilinadi.
γ₁ fazasi uchun taklif qilingan tajriba sharoitlarining kattaliklari
Manba \(B=10\ \mathrm T\) uchun γ₁ fazasining \(|m^\*|\geq10^2m_0\) tekis zona sohalarida siklotron chastotasi taxminan:
\[ 10^{-2}\text{–}10^1\ \mathrm{GHz} \]
tartibiga tushishini xabar qiladi. Bu chastotalar sub-meV foton energiyalari va uzoq infraqizil–mikrotoʻlqin sohasi bilan bogʻlanadi.
Ayni sohalarda zarur kritik sochilish vaqti:
\[ \tau_0 \sim 0.01\text{–}10\ \mathrm{ns} \]
oraligʻigacha koʻtariladi.
γ₁ fazasining Γ–X yoʻnalishidagi kritik uzunlik \(L_0\), tekis zona sohalarida taxminan:
\[ 10^3\ \mu\mathrm m \]
tartibiga yetadi. Bu qiymat taxminan millimetr miqyosidadir. Manba yuqori sifatli monokristallar uchun bu talab qiyin, ammo amalga oshirilishi mumkin boʻlishi mumkinligini ta’kidlaydi.
Γ–Y va Γ–Z yoʻnalishlarida zarur kritik uzunlik taxminan:
\[ 10^2\ \mu\mathrm m \]
tartibidadir va mualliflar bu yoʻnalishlarni tajribaviy jihatdan qulayroq deb hisoblaydilar.
γ₁ taqqoslashida HSE GGA ga nisbatan ayrim biroz kengroq zonalar va shunga mos ravishda biroz kichikroq samarali massalarni beradi. Natijada HSE bilan hisoblangan \(\omega_c\) biroz kattaroq, \(\tau_0\) va \(L_0\) esa biroz kichikroq boʻlishi mumkin.
B va δ fazalari uchun tajriba bahosi
B va δ fazalarining Γ–X va Γ–Y tekis zona sohalarida hisoblangan xatti-harakat sifat jihatdan γ₁ fazasiga oʻxshaydi. B fazasida \(|m^\*|\geq10^2m_0\) sohalarida:
- \(\omega_c\) taxminan \(10^{-1}\)–\(10^1\) GHz,
- \(\tau_0\) taxminan 0.01–10 ns,
- \(L_0\) choʻqqi qiymatlari taxminan \(10^2\)–\(10^3\) µm
tartibidadir. δ fazasida ham shunga oʻxshash kattaliklar xabar qilingan. Har ikkala tuzilmada Γ–Y yoʻnalishidagi \(L_0\) odatda Γ–X yoʻnalishidagidan kichikroq va manba buni siklotron rezonansi shu yoʻnalishda qulayroq boʻlishi mumkinligi sifatida talqin qiladi.
Tadqiqot Usuli va Natijalari
Hisoblash jarayoni
- γ₁, B va δ Ta₂O₅ kristall tuzilmalari tanlandi.
- Ikki- yoki uch-koordinatsiyalangan kislorod atomlari chiqarilib, manbada %10 deb ta’riflangan nuqsonli tuzilmalar hosil qilindi.
- Nuqson hosil qilingandan soʻng ion tuzilmalari toʻliq relaksatsiya qilindi.
- PBE-GGA/PAW asosidagi VASP hisob-kitoblari bilan elektron zona tuzilmalari va DOS olindi.
- γ₁ fazasida HSE06 zona hisob-kitoblari GGA natijalarini taqqoslash uchun ishlatildi.
- Fermi sathi yaqinidagi tekis zonalar Γ–X, Γ–Y va Γ–Z kabi yuqori simmetriya yoʻnalishlarida tahlil qilindi.
- \(E(k)\) egri chiziqlari boʻlakma-boʻlak yuqori darajali polinomlar bilan fit qilindi.
- Fit hosilalaridan elektron tezligi va zona-samarali massa hisoblandi.
- γ₁ fazasidagi %10 tip-2 vakansiyali tuzilma uchun 30 K da 5 ps AIMD bajarildi.
- AIMD dan keyin XRD va elektron zona tuzilishi qayta baholanib, strukturaviy va elektron xususiyatlarning saqlanishi tekshirildi.
- γ₁ fazasi uchun kislorod vakansiyasining hosil boʻlish energiyasi hisoblandi.
- DFT dan olingan \(m^\*(k)\) va \(v(k)\) qiymatlaridan foydalanib kelajakdagi siklotron rezonansi tajribasining \(\omega_c\), \(\tau_0\) va \(L_0\) talablari baholandi.
Fazalar orasidagi asosiy natijalarni taqqoslash
| Faza | Vakansiya turi | Elektron holati | Tekis zona kengligi | Ogʻir samarali massa natijasi |
|---|---|---|---|---|
| γ₁ | Tip 2 | GGA taqiqlangan zona kengligi 1,1 eV | Γ–X: 96,3 meV; Γ–Y: 41,3 meV; Γ–Z: 23,4 meV | Har uchala yoʻnalishda ham \(>10^3m_0\); ayrim sohalarda \(10^4m_0\) tartibi |
| B | Tip 1 | Taqiqlangan zona kengligi 0,36 eV | Γ–X: 29,8 meV; Γ–Y: 76,6 meV | Ayrim sohalarda \(>10^3m_0\); maksimum taxminan \(10^4m_0\) tartibida |
| δ | Tip 1 | %10 vakansiyada metall holat | Γ–X: 76,2 meV; Γ–Y: 62,1 meV | Ayrim sohalarda \(>10^3m_0\); maksimum taxminan \(10^4m_0\) tartibida |
AIMD va nuqson energetikasi natijasi
- γ₁ fazasida 30 K, 5 ps NVT AIMD davomida umumiy energiyada sezilarli siljish kuzatilmagan.
- Vaqt boʻyicha oʻrtacha qiymat −504.98576 eV.
- AIMD dan keyingi tuzilmaning 0 K asosiy holat energiyasi −505.00368 eV.
- Farq 0.01792 eV.
- Boshlangʻich va AIMD dan keyin hisoblangan XRD choʻqqi holatlari katta darajada ustma-ust tushgan.
- Fermi sathi yaqinidagi tekis zonalar AIMD dan keyingi zona tuzilishida saqlangan.
- %10 vakansiyali γ₁ tuzilmasi uchun kislorod vakansiyasining hosil boʻlish energiyasi 0 K da 7,85 eV, 300 K/atrof-muhit bosimi termodinamik tuzatishi bilan 7,21 eV deb hisoblangan.
Tadqiqot qoʻllab-quvvatlaydigan natijalar
- Tekshirilgan %10 kislorod vakansiyasi modellarida γ₁, B va δ Ta₂O₅ fazalarida Fermi sathi yaqinida tekis zonalar hosil boʻladi.
- Bu tekis zonalarning ayrim sohalarida hisoblangan zona-samarali massa erkin elektron massasining \(10^3\) baravaridan oshadi va batafsil zona tahlilida taxminan \(10^4\) tartibiga yetadi.
- γ₁ fazasida GGA va HSE zona natijalari tafsilotlarda farq qilsa ham, tekis zona va katta samarali massa xatti-harakatini sifat jihatdan saqlaydi.
- B va δ fazalarida ham tekis zona va katta samarali massa sohalari hisoblangan.
- %10 vakansiyali δ fazasi hisoblashda metall holatga oʻtadi.
- γ₁ fazasining %10 tip-2 vakansiyali konfiguratsiyasi 30 K dagi 5 ps AIMD sinovida strukturaviy va elektron xususiyatlarini saqlagan.
- Siklotron rezonansi nazariy samarali massa bashoratlarini kelajakda sinash uchun manba tomonidan taklif qilingan tajribaviy usuldir.
Tadqiqot qoʻllab-quvvatlamaydigan umumlashtirishlar
- Maqolada ogʻir fermionga oʻxshash xatti-harakat tajribada kuzatilmagan.
- Siklotron rezonansi oʻlchovi bajarilmagan; faqat tajriba sharoitlari hisoblangan.
- Tadqiqot xona haroratida ogʻir fermion xatti-harakati kuzatilganini koʻrsatmaydi.
- 30 K dagi 5 ps AIMD uzoq muddatli yoki xona haroratidagi barqarorlikni isbotlamaydi.
- Klassik ogʻir fermionlarda muhim boʻlgan Kondo effekti yoki boshqa koʻp-jismli korrelyatsiya mexanizmlari bevosita hisoblanmagan.
- Elektron–fonon bogʻlanishi va polaron hosil boʻlishining samarali massaga hissasi joriy hisoblash natijasi emas; bu mualliflarning kelajak tadqiqot boʻyicha taklifidir.
- B va δ fazalari uchun γ₁ fazasidagiga teng AIMD barqarorlik sinovi xabar qilinmagan.
- Yuqori kislorod vakansiyasi konsentratsiyasini tajribada qanchalik oson hosil qilish yoki saqlash mumkinligi ushbu tadqiqotda tajribaviy koʻrsatilmagan.
- Tadqiqotning xulosa boʻlimidagi \(10^3\)–\(10^5m_0\) iborasi asosiy matn va 1-jadvaldagi taxminan \(10^4m_0\) choʻqqi kattaliklari bilan toʻliq mos kelmaydi; Verianla \(10^5m_0\) qiymatini mustaqil tasdiqlangan natija sifatida ishlatmaydi.
Usuliy va hisobot cheklovlari
- Nazariy tadqiqot: Asosiy natijalar birinchi-prinsiplar elektron tuzilish hisob-kitoblaridan kelib chiqadi.
- Vakansiya konsentratsiyasi: Asosiy natijalar %10 kislorod vakansiyasi modellariga qaratilgan.
- AIMD qamrovi: Dinamik barqarorlik tekshiruvi faqat γ₁ tip-2 tuzilmasi uchun, 30 K va 5 ps sharoitida bajarilgan.
- Samarali massa fittingi: Juda katta samarali massa qiymatlari boʻlakma-boʻlak fit qilingan \(E(k)\) egri chizigʻining hosilalaridan olinadi; manba alohida xatolik chizigʻi yoki fitting noaniqligini xabar qilmaydi.
- A ilovadagi oraliq nomuvofiqligi: Matn umumiy \(k\) oraligʻini 0–0.03 deb yozgan holda, oxirgi fitting segmentini 0.02–0.04 deb beradi.
- Xulosa boʻlimidagi kattalik farqi: Xulosadagi \(10^5m_0\) yuqori tartibi batafsil natijalar va jadvalda koʻrinadigan \(10^4m_0\) tartibidan farq qiladi.
- B/δ HSE hisobot berilishi: ayrim rasm tagyozuvlarida HSE iborasi uchrasa-da, asosiy matn va koʻrinadigan egri chiziqlar bu ikki faza uchun asosan GGA natijalarini taqdim etadi.
- Tajribaviy talablar: taklif qilingan siklotron rezonansi ogʻir samarali massa sohalarida uzoq sochilish vaqtlarini, kuchli magnit maydonlarni va ayrim yoʻnalishlarda yuzlab mikrometr–millimetr miqyosidagi sifatli namunalarni talab qilishi mumkin.
Manba va Usul Haqida Izoh
Asl sarlavha: Possible Heavy-Fermion-like Behavior in Wide-Bandgap Semiconductor Ta2O5
Mualliflar: Hui-Min Tang, Yong Yang, Yoshiyuki Kawazoe.
Mualliflar tartibi: Hui-Min Tang → Yong Yang → Yoshiyuki Kawazoe.
Mas’ul muallif: Yong Yang; manbada yyanglab@issp.ac.cn manzili bilan “Correspondence” sifatida koʻrsatilgan.
Teng hissa: Manbada teng hissa yoki teng birinchi mualliflik bayonoti yoʻq.
Muassasalar: School of Physical Science and Technology, Guangxi Normal University; Key Laboratory of Photovoltaic and Energy Conservation Materials, Institute of Solid State Physics, Hefei Institutes of Physical Sciences, Chinese Academy of Sciences; Science Island Branch of Graduate School, University of Science and Technology of China; New Industry Creation Hatchery Center, Tohoku University; Department of Centre for Inter Disciplinary Research, SRM University-AP.
Jurnal: Physics.
Nashriyot: MDPI.
Jild / son / maqola raqami: 8(3), 68.
DOI: 10.3390/physics8030068.
Topshirilgan: 30 iyun 2026.
Tahrir: 28 avgust 2026.
Qabul qilingan: 3 sentabr 2026.
Nashr qilingan: 19 sentabr 2026.
Manba turi: Taqrizdan oʻtgan tadqiqot maqolasi; nazariy/sonli kondensirlangan modda va hisoblash materiallari fizikasi tadqiqoti.
Jurnal boʻlimi: Condensed Matter Physics.
Litsenziya: Creative Commons Attribution, CC BY 4.0.
Asosiy usul: VASP asosidagi DFT; PAW; PBE-GGA; qiyosiy HSE06; zona tuzilishi/DOS; boʻlakma-boʻlak \(E(k)\) fittingi; hosilaga asoslangan tezlik va samarali massa tahlili; AIMD; hisoblangan XRD; nuqson hosil boʻlish energiyasi; siklotron rezonansi tajriba dizayni.
Asosiy tizim: %10 kislorod vakansiyali γ₁, B va δ Ta₂O₅ fazalari.
Moliyalashtirish: National Natural Science Foundation of China, loyiha raqamlari 12074382, 12464012, 11804062 va 11474285.
Ma’lumotlarning mavjudligi: Manba tadqiqotning original hissalari maqola ichida ekanini va qoʻshimcha savollar mas’ul muallifga yoʻnaltirilishi mumkinligini bildiradi. Alohida ochiq xom ma’lumotlar ombori xabar qilinmagan.
Manfaatlar toʻqnashuvi: Mualliflar manfaatlar toʻqnashuvi yoʻqligini bayon qilganlar.
Mualliflar hissasi: Konseptuallashtirish Y.Y.; usul H.-M.T. va Y.Y.; dasturiy ta’minot H.-M.T. va Y.Y.; tekshirish Y.Y.; formal tahlil va tadqiqot H.-M.T. va Y.Y.; resurslar H.-M.T., Y.Y. va Y.K.; ma’lumotlarni kuratsiya qilish H.-M.T. va Y.Y.; dastlabki qoʻlyozma H.-M.T. va Y.Y.; koʻrib chiqish/tahrirlash Y.Y. va Y.K.; vizuallashtirish, nazorat va loyiha boshqaruvi Y.Y.; moliyalashtirishni jalb qilish H.-M.T., Y.Y. va Y.K. deb xabar qilingan.
Hisoblash infratuzilmasi: Tadqiqot Chinese Academy of Sciences Hefei Branch Supercomputing Center, Hefei Advanced Computing Center va Tohoku University Institute for Materials Research tarkibidagi hisoblash resurslaridan foydalangan; HPCI System Research Project hp200246 ham qayd etilgan.
Manba ichidagi nomuvofiqlik 1: Abstract/asosiy natijalar va 1-jadval taxminan \(10^4m_0\) maksimal tartibini qoʻllab-quvvatlasa, xulosa boʻlimi \(10^3\)–\(10^5m_0\) iborasini ishlatadi. Bu farq tuzatilmasdan koʻrinadigan tarzda saqlangan.
Manba ichidagi nomuvofiqlik 2: A ilova \(k\) oraligʻini 0–0.03 deb belgilagan holda, fitting segmentlarini 0–0.02, 0.01–0.03 va 0.02–0.04 deb beradi.
Manba ichidagi hisobot noaniqligi 3: B va δ fazalarining ayrim rasm tagyozuvlarida GGA va HSE birga tilga olinsa-da, asosiy natijalar matni va koʻrinadigan egri chiziqlar asosan GGA hisob-kitoblarini taqdim etadi. Verianla bu fazalar uchun manbada aniq berilmagan HSE sonlarini yaratmagan.
Asosiy talqin chegarasi: Tadqiqot juda katta zona-samarali massalardan kelib chiqib “possible heavy-fermion-like behavior” ni taklif qiladi. Bu natija ogʻir fermionning tajribaviy kuzatuvi yoki klassik ogʻir fermion many-body fizikasining toʻliq tasdigʻi emas. Siklotron rezonansi nazariy bashoratni sinash uchun kelajakda qoʻllash taklif qilingan tajribadir.
Izoh qoldiring
E-pochta manzilingiz chop etilmaydi. Majburiy maydonlar * bilan belgilangan