Oksigen vakansiyalı tantal pentoksid (Ta₂O₅), kristal qəfəsdəki bəzi oksigen atomlarının çıxarıldığı və bu qüsurlar səbəbindən elektron zolaq quruluşu dəyişə bilən geniş qadağan olunmuş zonalı oksid yarımkeçirici sistemidir. Bu tədqiqat γ₁, B və δ struktur fazalarında %10 oksigen vakansiyası yaradıldıqda elektron zolaqların necə dəyişdiyini sıxlıq funksionalı nəzəriyyəsi (Density Functional Theory, DFT) ilə araşdırır. Hesablamalarda müəyyən istiqamətlərdə Fermi enerjisinin yaxınlığında olduqca düz elektron zolaqlar yaranır və bu zolaqların bəzi sahələrində elektronların zolaq-effektiv kütləsi sərbəst elektron kütləsinin \(10^3\) mərtəbəsini aşaraq bəzi nöqtələrdə \(10^4\) mərtəbəsinə çatır. Müəlliflər bu xüsusiyyəti ehtiyatla “heavy-fermion-like”, yəni ağır fermiona bənzər davranış adlandırırlar.

Tədqiqatın mühüm cəhəti ənənəvi ağır fermion sistemlərindən fərqli bir material sinfini araşdırmasıdır. Klassik ağır fermion birləşmələri əsasən nadir torpaq və ya aktinid elementləri ehtiva edən, güclü elektron korrelyasiyalarının və aşağı temperatur fizikasının həlledici olduğu metal və ya dar qadağan olunmuş zonalı sistemlərdir. Buradakı Ta₂O₅ tədqiqatı isə oksigen vakansiyalarının yaratdığı düz zolaqlardan törəyən çox böyük zolaq-effektiv kütlələri geniş qadağan olunmuş zonalı oksid sistemində araşdırır. Tədqiqat Kondo effektini, xüsusi istilik əmsalını, aşağı temperatur elektrik müqavimətini və ya birbaşa ağır fermion daşıyıcılarını eksperimental olaraq ölçmür; buna görə “ağır fermiona bənzər” ifadəsinin sərhədi qorunmalıdır.

Hesablamalarda VASP proqramı, PAW potensialları, PBE formasında GGA mübadilə-korrelyasiya yanaşması və 600 eV müstəvi-dalğa enerji kəsimi istifadə edilmişdir. γ₁ fazası üçün 1×1×2 superhüceyrə və 10×10×4 Monkhorst–Pack k-nöqtə şəbəkəsi; B fazası üçün 2×4×4, δ fazası üçün 6×6×12 k-nöqtə şəbəkəsi istifadə edilmişdir. γ₁ fazasında zolaq dispersiyasını yoxlamaq üçün HSE06 hibrid funksionalı ilə əlavə hesablamalar aparılmışdır.

γ₁ fazasında %10 tip-2 oksigen vakansiyası şəraitində GGA qadağan olunmuş zona eni 1,1 eV tapılmışdır. Γ–X, Γ–Y və Γ–Z istiqamətlərində düz zolaq enləri GGA ilə müvafiq olaraq 96,3, 41,3 və 23,4 meV; HSE ilə 110,0, 40,0 və 20,0 meV olaraq bildirilmişdir. B fazasında düz zolaq enləri Γ–X və Γ–Y üçün 29,8 və 76,6 meV; δ fazasında 76,2 və 62,1 meV-dir. B fazasının hesablanmış qadağan olunmuş zona eni 0,36 eV olduğu halda, %10 vakansiyalı δ fazası metal xassə göstərir.

γ₁ fazasındakı yüksək vakansiyalı quruluşun yalnız statik həndəsi optimallaşdırma ilə məhdud qalmaması üçün əlavə olaraq 30 K-də, NVT ansamblında, 2×2×1 superhüceyrə ilə 5 ps ab initio molekulyar dinamika (AIMD) simulyasiyası aparılmışdır. Ümumi enerji nəzərəçarpan sürüşmə göstərməmiş, başlanğıc və AIMD-dən sonrakı XRD pik mövqeləri böyük ölçüdə qorunmuş və Fermi səviyyəsi yaxınlığındakı düz zolaqlar simulyasiyanın sonunda da müşahidə edilmişdir. Bununla belə, bu yoxlama yalnız γ₁ fazasının göstərilən aşağı temperatur və qısa simulyasiya müddəti altındakı struktur/ədədi sabitliyinə aiddir; otaq temperaturunda uzunmüddətli eksperimental sabitliyin sübutu deyil.

Müəlliflər nəzəri proqnozun eksperimental yoxlanması üçün güclü maqnit sahəsində elektron siklotron rezonansını təklif edirlər. Məsələn, 100 GHz rezonans tezliyi qəbul edilərsə, \(m^\*\sim10^2\)–\(10^4m_0\) diapazonu təxminən 0,36–35,7 T maqnit sahəsi tələb edir. Lakin məqalədə bu təcrübə aparılmamışdır; verilən qiymətlər DFT-dən hesablanan effektiv kütlələrə əsaslanan eksperiment dizaynı təxminləridir.

Ağır fermion anlayışının bu məqalədəki konteksti

Ənənəvi ağır fermion fizikasında yükdaşıyıcıların effektiv kütləsi sərbəst elektron kütləsinin yüzlərlə, hətta minlərlə qatına çata bilər. Mənbə giriş bölməsində CeAl₃ kimi klassik sistemlərdə aşağı temperatur Fermi-mayesi davranışını nümunə göstərir və elektrik müqaviməti üçün \(\Delta\rho(T)\propto T^2\), istilik tutumu üçün \(C(T)=\gamma T\) münasibətlərini xatırladır. Lakin Ta₂O₅ tədqiqatı bu termodinamik və ya nəqliyyat kəmiyyətlərini hesablamır və ölçmür.

Bu tədqiqatda ağır fermiona bənzər şərhin birbaşa əsası oksigen vakansiyalarının Fermi enerjisi yaxınlığında yaratdığı düz zolaqlar və həmin zolaqların çox kiçik əyrilik göstərən sahələrindən hesablanan çox böyük \(|m^\*|\) qiymətləridir. Buna görə tədqiqat klassik nadir-torpaq ağır fermion birləşmələrinin bütün çoxcisimli fizikasının Ta₂O₅-də göstərildiyini iddia etmir.

Düz Zolaq Elektronun Effektiv Kütləsini Niyə Böyüdür?

Düz elektron zolağında enerji \(E\), kristal impulsu \(k\) dəyişsə də çox az dəyişir; buna görə zolaq əyriliyinin böyüklüyü azalır. Zolaq nəzəriyyəsində effektiv kütlə əyriliyin tərsi ilə əlaqəli olduğundan kiçik \(|d^2E/dk^2|\) qiyməti çox böyük \(|m^\*|\) yaradır. Buna görə tədqiqatdakı qeyri-adi effektiv kütlə pikləri oksigen vakansiyalarının yaratdığı olduqca düz zolaq sahələrində meydana çıxır.

Öyrədici fon qeydi: Məqalə effektiv kütlə münasibətini ayrıca nömrələnmiş tənlik kimi yazmır, lakin \(E(k)\)-nın birinci və ikinci törəmələrindən sürət və effektiv kütləni hesabladığını açıq şəkildə bildirir. Bir əsas kristal istiqaməti boyunca istifadə olunan standart zolaq-nəzəriyyəsi münasibəti:

\[ v(k)=\frac{1}{\hbar}\frac{dE}{dk}, \]

\[ m^\*(k)= \frac{\hbar^2}{ \dfrac{d^2E}{dk^2} }. \]

şəklindədir. Birinci törəmə zolağın meylinin qrup sürətini necə müəyyən etdiyini; ikinci törəmə isə zolaq əyriliyinin effektiv kütləyə necə çevrildiyini göstərir. Bu iki tənlik Verianla-nın metodu izah etmək üçün verdiyi standart zolaq fizikası fonudur; məqalənin yeni tədqiqat nəticəsi kimi təqdim edilmir.

Formulda \(\hbar\) azaldılmış Plank sabitidir. \(dE/dk\) enerji zolağının meylini, \(d^2E/dk^2\) isə əyriliyini göstərir. Əyrilik sıfıra yaxınlaşdıqda effektiv kütlənin mütləq qiyməti çox böyüyür. Əyriliyin işarəsi dəyişdikdə \(m^\*\)-ın işarəsi də dəyişə bilər; məqalə buna görə ağır-kütlə sahələrini qiymətləndirərkən \(|m^\*|\)-dan istifadə edir. Mənfi zolaq-effektiv kütləsi hissəciyin həqiqi kütləsinin fiziki olaraq mənfi olması demək deyil; zolaq əyriliyinin işarəsini ifadə edir.

Oksigen vakansiyası növləri necə yaradıldı?

Ta₂O₅-dəki oksigen atomları Ta atomları ilə rabitə sayına görə iki koordinasiya sinfinə ayrılmışdır. Mənbədə tip 1 oksigen vakansiyası, iki Ta atomu ilə koordinasiyalı bir oksigenin çıxarılması; tip 2 oksigen vakansiyası isə üç Ta atomu ilə koordinasiyalı bir oksigenin çıxarılması kimi tərif edilir.

• γ₁ fazası: 1×1×2 superhüceyrə yaradılmış və bir üç-koordinasiyalı oksigen çıxarılaraq tip-2 vakansiya yaradılmışdır.
• B fazası: primitiv hüceyrədən iki iki-koordinasiyalı oksigen çıxarılmışdır; vakansiya tipi 1-dir.
• δ fazası: primitiv hüceyrədən bir iki-koordinasiyalı oksigen çıxarılmışdır; vakansiya tipi 1-dir.

Bu üç quruluş tədqiqat boyu %10 oksigen vakansiyası konsentrasiyası kimi götürülmüş və oksigen çıxarıldıqdan sonra qüsurlu strukturların ion mövqeləri tam relaksasiya edilmişdir.

DFT hesablama şərtləri

γ₁ fazasında zolaq quruluşu

%10 tip-2 oksigen vakansiyalı γ₁-Ta₂O₅ üçün GGA hesablaması 1,1 eV qadağan olunmuş zona eni verir. Fermi enerjisinin yaxınlığında düz zolaq yaranır. Tədqiqat bu düz zolağı Γ–X, Γ–Y və Γ–Z yüksək simmetriya istiqamətlərində ayrıca araşdırmışdır.

Daha kiçik zolaq eni araşdırılan impuls yolu boyunca enerjinin daha az dəyişdiyini və zolağın daha düz olduğunu ifadə edir. γ₁ fazasında bu üç istiqamət içində bildirilən ən dar zolaq Γ–Z istiqamətindədir.

E(k) əyrisi niyə hissə-hissə fit edildi?

Effektiv kütlə \(E(k)\) əyrisinin ikinci törəməsindən asılı olduğuna görə, fitləmədəki kiçik dalğalanmalar belə törəmələrdə böyüyə bilər. Məqalə buna görə \(E(k)\) verilənlərini tək bir qlobal polinomla fit etmək əvəzinə üst-üstə düşən hissələrə bölərək lokal fitləmə istifadə etmişdir.

Əlavə A-dakı nümunədə istifadə olunan polinom:

\[ E(k) = a k^{10} + b k^8 + c k^6 + d k^4 + e k^2 + f. \]

Burada \(a,b,c,d,e,f\) müvafiq lokal məlumat hissəsinə uyğunlaşdırılan əmsallardır. Yalnız cüt \(k\) qüvvələrinin istifadə edildiyi bu ifadə müvafiq zolaq hissəsinin hamar əyri ilə təmsil olunmasını təmin edir. Fit edilmiş \(E(k)\) daha sonra törəmələnərək elektron sürəti və zolaq əyriliyi hesablanmışdır.

Mənbə nümunə γ₁ Γ–Z verilənlərini 0–0.02, 0.01–0.03 və 0.02–0.04 olmaqla üç üst-üstə düşən seqmentə ayırdığını bildirir. Lakin eyni paraqrafda ümumi \(k\) intervalını 0–0.03 kimi yazır. Son seqmentin 0.04-ə qədər uzanması bu cümlə ilə uyğun deyil; Verianla hansı qiymətin yazı xətası olduğunu fərz etmir.

Effektiv kütlənin Brillouin zonasındakı paylanması

Mənbə şəkillərində sarı sahələr:

\[ |m^\*|\geq 10m_0 \]

şərtini; boz sahələr isə:

\[ |m^\*|\geq10^2m_0 \]

şərtini təmsil edir. Burada \(m_0\) sərbəst elektron kütləsidir.

γ₁ fazasında Γ–X, Γ–Y və Γ–Z istiqamətlərinin hamısında bəzi impuls sahələrində \(|m^\*|>10^3m_0\) qiymətləri tapılmışdır. Bəzi dar sahələrdə qiymət \(10^4m_0\) mərtəbəsinə çatır.

\(\log_{10}|m^\*/m_0|=4\) qiyməti effektiv kütlənin sərbəst elektron kütləsinin təxminən \(10^4\) qatı mərtəbəsində olduğunu bildirir. Buna görə cədvəldəki 4.131, 4.064 və 4.107 kimi qiymətlər birbaşa \(10^4\) mərtəbəsində ağır effektiv kütlə sahələrini göstərir.

B və δ fazalarında nə tapıldı?

%10 tip-1 vakansiyalı B-Ta₂O₅ üçün hesablanan qadağan olunmuş zona eni 0,36 eV-dir. %10 tip-1 vakansiyalı δ-Ta₂O₅ isə hesablamada aydın qadağan olunmuş zona göstərməmiş və metal xassə qazanmışdır. Hər iki fazada Fermi səviyyəsi yaxınlığında düz zolaqlar müşahidə edilmişdir.

B və δ fazalarının düz zolaqlarının müəyyən hissələrində də \(|m^\*|>10^3m_0\) qiymətləri və təxminən \(10^4m_0\) mərtəbəsində maksimumlar bildirilir.

Mənbənin hesabat qeydi: B və δ nəticələrini izah edən əsas mətn GGA hesablamalarından bəhs edir və görünən qrafiklərdə yalnız GGA/Fit əyriləri var. Buna qarşılıq Şəkil 4 və Şəkil 5 altyazılarında “GGA and HSE methods” ifadəsi keçir. Mənbə bu fərqi izah etmədiyinə görə Verianla B və δ fazaları üçün müstəqil HSE ədədi nəticələri varmış kimi yeni qiymət yaratmır.

Bu Tədqiqatdakı “Ağır Fermiona Bənzər” İfadəsi Nə Deməkdir?

Tədqiqatda “ağır fermiona bənzər” oksigen vakansiyalarının yaratdığı düz elektron zolaqlarda zolaq əyriliyindən hesablanan elektron effektiv kütləsinin bəzi impuls sahələrində \(10^3\)–\(10^4m_0\) mərtəbəsinə yüksəlməsini ifadə edir. Məqalə bunu klassik ağır fermion birləşmələrindəki bütün çoxcisimli mexanizmlərin Ta₂O₅-də eksperimental olaraq təsdiqləndiyi mənasında istifadə etmir; birbaşa nəqliyyat, xüsusi istilik, Kondo fizikası və ya siklotron rezonansı ölçümü aparılmamışdır.

Bu fərq xüsusilə vacibdir, çünki ənənəvi ağır fermion sistemlərində böyük effektiv kütlə adətən güclü elektron korrelyasiyalarının nəticəsi kimi müzakirə edilir. Bu məqalədə isə birbaşa hesablanan kəmiyyət oksigen vakansiyalı Ta₂O₅-in tək-hissəcik elektron zolaq dispersiyası və ondan törəyən zolaq-effektiv kütləsidir.

Müəlliflər nəticə bölməsində həmçinin elektron–fonon əlaqələnməsi və polaron yaranmasının Ta₂O₅-də effektiv kütləni daha da renormallaşdıra biləcəyini bildirirlər. Lakin bu mexanizmlər mövcud tədqiqatda ayrıca hesablanmamış və ya ölçülməmişdir; gələcək tədqiqat təklifi xarakteri daşıyır.

γ₁ fazasının dinamik sabitlik yoxlaması

%10 tip-2 oksigen vakansiyalı γ₁ fazasının dinamik sabitliyi üçün ab initio molekulyar dinamika tətbiq edilmişdir:

NVT, hissəcik sayının \(N\), həcmin \(V\) və temperaturun \(T\) sabit saxlandığı kanonik ansambldır. 5 ps ərzində ümumi enerji sabit bir orta qiymət ətrafında tərəddüd etmiş və nəzərəçarpan sistematik sürüşmə göstərməmişdir.

AIMD-dən sonra hesablanan XRD pik mövqelərinin başlanğıc quruluşla müqayisədə böyük ölçüdə qorunması kristal skeletin bu simulyasiya şərtləri altında ciddi yenidən qurulmaya məruz qalmadığını göstərir. AIMD-nin son kristal konfiqurasiyasından hesablanan zolaq quruluşunda Fermi səviyyəsi yaxınlığındakı düz zolaqlar da qorunmuş və mənbə dörd düz zolağın aydın şəkildə görünə bildiyini bildirmişdir.

Bu nəticə 30 K və 5 ps AIMD şərtlərindəki γ₁ quruluşu üçün ədədi dinamik sabitlik sübutudur. B və δ fazaları üçün eyni AIMD testi bildirilməmişdir. Bundan əlavə, 5 ps/30 K nəticəsi otaq temperaturunda uzunmüddətli eksperimental sabitliyin və ya sintez oluna bilmənin birbaşa sübutu deyil.

Oksigen vakansiyasının yaranma enerjisi

Əlavə C-də γ₁ fazasındakı %10 oksigen vakansiyası üçün yaranma enerjisi:

\[ E_{\mathrm{vf}} = E_{\mathrm{defect}} + 0.5E_{\mathrm{O_2}} - E_{\mathrm{perfect}} - T\Delta S. \tag{A1} \]

münasibəti ilə hesablanmışdır.

Burada:

• \(E_{\mathrm{vf}}\): oksigen vakansiyasının yaranma enerjisidir.
• \(E_{\mathrm{defect}}\): bir oksigen atomu çıxarılmış qüsurlu superhüceyrənin ümumi enerjisidir.
• \(E_{\mathrm{perfect}}\): qüsursuz superhüceyrənin ümumi enerjisidir.
• \(E_{\mathrm{O_2}}\): təcrid olunmuş O₂ molekulunun ümumi enerjisidir.
• \(T\Delta S\): temperatur və entropiya dəyişməsinin yaratdığı termodinamik düzəlişdir.

Mənbədə istifadə olunan statik DFT qiymətləri:

\[ E_{\mathrm{defect}} = -126.31649376\ \mathrm{eV}, \]

\[ E_{\mathrm{perfect}} = -136.63040186\ \mathrm{eV}, \]

\[ E_{\mathrm{O_2}} = -4.93232\ \mathrm{eV}. \]

Bu qiymətlərlə 0 K yaranma enerjisi:

\[ E_{\mathrm{vf}}(0\ \mathrm K) = 7.85\ \mathrm{eV} \]

olaraq tapılmışdır. 300 K və ətraf mühit təzyiqi üçün:

\[ T\Delta S = 0.6376\ \mathrm{eV} \]

termini daxil edildikdə:

\[ E_{\mathrm{vf}}(300\ \mathrm K) = 7.21\ \mathrm{eV} \]

alınmışdır.

Mənbə bu qiymətin ədəbiyyatda LSR fazasındakı tək bir oksigen vakansiyası üçün bildirilən təxminən 5 eV-dən təxminən 2 eV yüksək olduğunu və bunu yüksək vakansiya konsentrasiyasının qüsur energetikasına təsiri ilə əlaqələndirir. Buradakı təxminən 5 eV mövcud tədqiqatda eyni quruluş üzərində yenidən hesablanmış nəzarət qiyməti deyil, mənbənin Referans 48-dən apardığı müqayisədir.

Siklotron Rezonansı Bu Nəzəri Proqnozu Necə Sınaqdan Keçirə Bilər?

Siklotron rezonansı maqnit sahəsindəki yükdaşıyıcıların dairəvi hərəkət tezliyini elektromaqnit şüalanması ilə rezonansa gətirərək effektiv kütləni ölçməyə yarayan eksperimental üsuldur. Tədqiqat oksigen vakansiyalı Ta₂O₅-də hesablanan çox böyük \(m^\*\) qiymətlərinin gələcəkdə aşağı temperatur və güclü maqnit sahəsi şərtlərində bu üsulla sınaqdan keçirilə biləcəyini təklif edir; məqalənin özündə belə ölçmə aparılmamışdır.

Yükdaşıyıcının siklotron bucaq tezliyi:

\[ \omega_c = \frac{eB}{m^\*}. \tag{1} \]

Burada \(e\) elementar yük, \(B\) maqnit sahəsi, \(m^\*\) isə yükdaşıyıcının effektiv kütləsidir. Tənlikdə effektiv kütlə artdıqca eyni maqnit sahəsində siklotron tezliyi azalır. Buna görə çox ağır effektiv kütlələr mikrodalğa və ya uzaq infraqırmızı tezliklərə uyğun gələ bilər.

Siklotron tezliyi:

\[ \nu_c = \frac{\omega_c}{2\pi} \]

şəklində yazılır. Mənbə \(\nu_c=100\ \mathrm{GHz}\) üçün tələb olunan sahəni:

\[ 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 \]

şəklində verir.

\[ m^\*=10^2m_0 \]

üçün təxminən:

\[ B\approx0.36\ \mathrm T \]

və:

\[ m^\*=10^4m_0 \]

üçün təxminən:

\[ B\approx35.7\ \mathrm T \]

tələb olunur. Mənbə bu diapazonun mövcud davamlı yüksək maqnit sahəsi qurğularının imkanları daxilində olduğunu bildirir.

Səpilmə müddəti niyə kritikdir?

Yükdaşıyıcının rezonans hərəkətinin səpilmə ilə yox olmadan müşahidə edilə bilməsi üçün:

\[ \omega_c\tau>1 \]

şərti tələb olunur. Burada \(\tau\), yükdaşıyıcının orta sərbəst səpilmə müddətidir.

Hədd vəziyyəti:

\[ \omega_c\tau_0=1 \]

kimi tərif edildikdə:

\[ \tau_0(k) = \frac{m^\*(k)}{eB}. \tag{2} \]

alınır. Effektiv kütlə artdıqca tələb olunan minimum səpilmə müddəti də artır. Başqa sözlə, çox ağır yükdaşıyıcının siklotron orbitini kifayət qədər yaxşı müəyyən edə bilməsi üçün səpilmədən daha uzun müddət hərəkət etməsi lazımdır.

Kritik nümunə uzunluğu

Mənbə tələb olunan minimum xarakterik uzunluğu:

\[ L_0(k) = \tau_0(k)|v(k)|. \tag{3} \]

kimi qiymətləndirir. \(v(k)\), DFT zolaq quruluşundan əldə olunan elektron sürətidir.

Γ–X istiqaməti üçün \(L_0\), real fəzadakı X kristal istiqamətində rezonansın formalaşa bilməsi üçün tələb olunan minimum nümunə qalınlığı kimi şərh edilir.

γ₁ fazası üçün təklif olunan eksperiment şərtlərinin kəmiyyətləri

Mənbə, \(B=10\ \mathrm T\) üçün γ₁ fazasının \(|m^\*|\geq10^2m_0\) düz zolaq sahələrində siklotron tezliyinin təxminən:

\[ 10^{-2}\text{–}10^1\ \mathrm{GHz} \]

mərtəbəsinə düşdüyünü bildirir. Bu tezliklər sub-meV foton enerjiləri və uzaq infraqırmızı–mikrodalğa sahəsi ilə əlaqələndirilir.

Eyni sahələrdə tələb olunan kritik səpilmə müddəti:

\[ \tau_0 \sim 0.01\text{–}10\ \mathrm{ns} \]

diapazonuna yüksəlir.

γ₁ fazasının Γ–X istiqamətindəki kritik uzunluq \(L_0\), düz zolaq sahələrində təxminən:

\[ 10^3\ \mu\mathrm m \]

mərtəbəsinə çatır. Bu qiymət təxminən millimetr miqyasındadır. Mənbə yüksək keyfiyyətli monokristallar üçün bu tələbin çətin, lakin həyata keçirilə bilən ola biləcəyini müdafiə edir.

Γ–Y və Γ–Z istiqamətlərində tələb olunan kritik uzunluq təxminən:

\[ 10^2\ \mu\mathrm m \]

mərtəbəsindədir və müəlliflər bu istiqamətləri eksperimental baxımdan daha əlçatan hesab edirlər.

γ₁ müqayisəsində HSE, GGA ilə müqayisədə bəzi bir qədər daha geniş zolaqlar və buna bağlı olaraq bir qədər daha kiçik effektiv kütlələr verir. Nəticədə HSE ilə hesablanan \(\omega_c\) bir qədər böyük, \(\tau_0\) və \(L_0\) isə bir qədər kiçik ola bilər.

B və δ fazaları üçün eksperiment proqnozu

B və δ fazalarının Γ–X və Γ–Y düz zolaq sahələrində hesablanan davranış keyfiyyətcə γ₁ fazasına bənzəyir. B fazasında \(|m^\*|\geq10^2m_0\) sahələrində:

• \(\omega_c\) təxminən \(10^{-1}\)–\(10^1\) GHz,
• \(\tau_0\) təxminən 0.01–10 ns,
• \(L_0\) pik qiymətləri təxminən \(10^2\)–\(10^3\) µm

mərtəbələrindədir. δ fazasında da oxşar kəmiyyətlər bildirilmişdir. Hər iki quruluşda Γ–Y istiqamətindəki \(L_0\) adətən Γ–X istiqamətindəkindən daha kiçikdir və mənbə bunu siklotron rezonansının bu istiqamətdə daha əlçatan ola biləcəyi kimi şərh edir.

Tədqiqatın Metodu və Nəticələri

Hesablama axını

1. γ₁, B və δ Ta₂O₅ kristal quruluşları seçildi.
2. İki- və ya üç-koordinasiyalı oksigen atomları çıxarılaraq mənbədə %10 kimi tərif edilən qüsurlu quruluşlar yaradıldı.
3. Qüsurdan sonra ion quruluşları tam relaksasiya edildi.
4. PBE-GGA/PAW əsaslı VASP hesablamaları ilə elektron zolaq quruluşları və DOS əldə edildi.
5. γ₁ fazasında HSE06 zolaq hesablamaları GGA nəticələri ilə müqayisə üçün istifadə edildi.
6. Fermi səviyyəsi yaxınlığındakı düz zolaqlar Γ–X, Γ–Y və Γ–Z kimi yüksək simmetriya istiqamətlərində araşdırıldı.
7. \(E(k)\) əyriləri hissə-hissə yüksək dərəcəli polinomlarla fit edildi.
8. Fitin törəmələrindən elektron sürəti və zolaq-effektiv kütləsi hesablandı.
9. γ₁ fazasındakı %10 tip-2 vakansiya quruluşu üçün 30 K-də 5 ps AIMD aparıldı.
10. AIMD-dən sonra XRD və elektron zolaq quruluşu yenidən qiymətləndirilərək struktur və elektron xassənin qorunması yoxlanıldı.
11. γ₁ fazası üçün oksigen vakansiyasının yaranma enerjisi hesablandı.
12. DFT-dən əldə olunan \(m^\*(k)\) və \(v(k)\) qiymətlərindən istifadə etməklə gələcək siklotron rezonansı eksperimentinin \(\omega_c\), \(\tau_0\) və \(L_0\) tələbləri qiymətləndirildi.

Fazlar arasında əsas nəticələrin müqayisəsi

AIMD və qüsur energetikasının nəticəsi

• γ₁ fazasında 30 K, 5 ps NVT AIMD zamanı ümumi enerji nəzərəçarpan sürüşmə göstərməmişdir.
• Zaman ortalaması −504.98576 eV-dir.
• AIMD-dən sonrakı quruluşun 0 K əsas hal enerjisi −505.00368 eV-dir.
• Aradakı fərq 0.01792 eV-dir.
• Başlanğıc və AIMD-dən sonrakı hesablanmış XRD pik mövqeləri böyük ölçüdə üst-üstə düşmüşdür.
• Fermi səviyyəsi yaxınlığındakı düz zolaqlar AIMD-dən sonrakı zolaq quruluşunda qorunmuşdur.
• %10 vakansiyalı γ₁ quruluşu üçün oksigen vakansiyasının yaranma enerjisi 0 K-də 7,85 eV, 300 K/ətraf mühit təzyiqi termodinamik düzəlişi ilə 7,21 eV olaraq hesablanmışdır.

Tədqiqatın dəstəklədiyi nəticələr

• Araşdırılan %10 oksigen vakansiyası modellərində γ₁, B və δ Ta₂O₅ fazalarında Fermi səviyyəsi yaxınlığında düz zolaqlar yaranır.
• Bu düz zolaqların bəzi sahələrində hesablanan zolaq-effektiv kütləsi sərbəst elektron kütləsinin \(10^3\) qatını aşır və ətraflı zolaq analizində təxminən \(10^4\) mərtəbəsinə çatır.
• γ₁ fazasında GGA və HSE zolaq nəticələri detallarda fərqlənsə də, düz zolaq və böyük effektiv kütlə davranışını keyfiyyətcə qoruyur.
• B və δ fazalarında da düz zolaq və böyük effektiv kütlə sahələri hesablanmışdır.
• %10 vakansiyalı δ fazası hesablamada metal hala keçir.
• γ₁ fazasının %10 tip-2 vakansiyalı konfiqurasiyası 30 K-dəki 5 ps AIMD testində struktur və elektron xassələrini qorumuşdur.
• Siklotron rezonansı nəzəri effektiv kütlə proqnozlarını gələcəkdə sınamaq üçün mənbənin təklif etdiyi eksperimental üsuldur.

Tədqiqatın dəstəkləmədiyi ümumiləşdirmələr

• Məqalə ağır fermiona bənzər davranışı eksperimental olaraq müşahidə etməmişdir.
• Siklotron rezonansı ölçməsi aparılmamışdır; yalnız eksperiment şərtləri hesablanmışdır.
• Tədqiqat otaq temperaturunda ağır fermion davranışının müşahidə edildiyini göstərmir.
• 30 K-dəki 5 ps AIMD uzunmüddətli və ya otaq temperaturundakı sabitliyi sübut etmir.
• Klassik ağır fermionlarda mühüm olan Kondo effekti və ya digər çoxcisimli korrelyasiya mexanizmləri birbaşa hesablanmamışdır.
• Elektron–fonon əlaqələnməsinin və polaron yaranmasının effektiv kütləyə töhfəsi mövcud hesablamanın nəticəsi deyil; müəlliflərin gələcək tədqiqat təklifidir.
• B və δ fazaları üçün γ₁ fazasındakı ilə eyni AIMD sabitlik testi bildirilməmişdir.
• Yüksək oksigen vakansiyası konsentrasiyasının eksperimental olaraq nə qədər asan yaradıla və ya saxlanıla biləcəyi bu tədqiqatda eksperimental şəkildə göstərilməmişdir.
• Tədqiqatın nəticə bölməsindəki \(10^3\)–\(10^5m_0\) ifadəsi əsas mətn və Cədvəl 1-dəki təxminən \(10^4m_0\) pik kəmiyyətləri ilə tam uyğun gəlmir; Verianla \(10^5m_0\) qiymətini müstəqil təsdiqlənmiş nəticə kimi istifadə etmir.

Metodoloji və hesabat məhdudiyyətləri

• Nəzəri tədqiqat: Əsas nəticələr ilk-prinsiplər elektron quruluş hesablamalarından əldə olunur.
• Vakansiya konsentrasiyası: Əsas nəticələr %10 oksigen vakansiyası modellərinə yönəlib.
• AIMD əhatəsi: Dinamik sabitlik yoxlaması yalnız γ₁ tip-2 quruluşu üçün, 30 K və 5 ps şərtlərində aparılmışdır.
• Effektiv kütlə fitləməsi: Çox böyük effektiv kütlə qiymətləri hissə-hissə fit edilmiş \(E(k)\) əyrisinin törəmələrindən əldə edilir; mənbə ayrıca xəta zolağı və ya fitləmə qeyri-müəyyənliyi bildirmir.
• Əlavə A interval uyğunsuzluğu: Mətn ümumi \(k\) intervalını 0–0.03 kimi yazarkən son fitləmə seqmentini 0.02–0.04 verir.
• Nəticə bölməsində kəmiyyət fərqi: Nəticədəki \(10^5m_0\) yuxarı mərtəbəsi ətraflı nəticə və cədvəldə görünən \(10^4m_0\) mərtəbəsindən fərqlidir.
• B/δ HSE hesabatı: bəzi şəkil altyazılarında HSE ifadəsi yer alsa da, əsas mətn və görünən əyrilər bu iki faza üçün əsasən GGA nəticələrini təqdim edir.
• Eksperimental tələblər: təklif olunan siklotron rezonansı ağır effektiv kütlə sahələrində uzun səpilmə müddətləri, güclü maqnit sahələri və bəzi istiqamətlərdə yüzlərlə mikrometr–millimetr miqyaslı keyfiyyətli nümunələr tələb edə bilər.

Mənbə və Metod Qeydi

Orijinal başlıq: Possible Heavy-Fermion-like Behavior in Wide-Bandgap Semiconductor Ta2O5

Müəlliflər: Hui-Min Tang, Yong Yang, Yoshiyuki Kawazoe.

Müəllif sırası: Hui-Min Tang → Yong Yang → Yoshiyuki Kawazoe.

Məsul müəllif: Yong Yang; mənbədə yyanglab@issp.ac.cn ünvanı ilə “Correspondence” kimi göstərilmişdir.

Bərabər töhfə: Mənbədə bərabər töhfə və ya bərabər birinci müəlliflik bəyanı yoxdur.

Qurumlar: 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.

Nəşriyyat: MDPI.

Cild / nömrə / məqalə nömrəsi: 8(3), 68.

DOI: 10.3390/physics8030068.

Göndərilmə: 30 iyun 2026.

Reviziya: 28 avqust 2026.

Qəbul: 3 sentyabr 2026.

Nəşr: 19 sentyabr 2026.

Mənbə növü: Resenziyalı tədqiqat məqaləsi; nəzəri/ədədi kondensə olunmuş maddə və hesablama materiallar fizikası tədqiqatı.

Jurnal bölməsi: Condensed Matter Physics.

Lisenziya: Creative Commons Attribution, CC BY 4.0.

Əsas metod: VASP əsaslı DFT; PAW; PBE-GGA; müqayisəli HSE06; zolaq quruluşu/DOS; hissə-hissə \(E(k)\) fitləməsi; törəmə əsaslı sürət və effektiv kütlə analizi; AIMD; hesablanmış XRD; qüsurun yaranma enerjisi; siklotron rezonansı eksperiment dizaynı.

Əsas sistem: %10 oksigen vakansiyalı γ₁, B və δ Ta₂O₅ fazaları.

Maliyyələşdirmə: National Natural Science Foundation of China, layihə nömrələri 12074382, 12464012, 11804062 və 11474285.

Məlumatların əlçatanlığı: Mənbə tədqiqatın orijinal töhfələrinin məqalə daxilində olduğunu və əlavə sualların məsul müəllifə yönəldilə biləcəyini bildirir. Ayrı açıq xam məlumat deposu bildirilməmişdir.

Maraqlar toqquşması: Müəlliflər maraqlar toqquşmasının olmadığını bəyan etmişlər.

Müəllif töhfələri: Konseptuallaşdırma Y.Y.; metod H.-M.T. və Y.Y.; proqram təminatı H.-M.T. və Y.Y.; yoxlama Y.Y.; formal analiz və tədqiqat H.-M.T. və Y.Y.; resurslar H.-M.T., Y.Y. və Y.K.; məlumatların kurasiyası H.-M.T. və Y.Y.; ilkin layihə H.-M.T. və Y.Y.; nəzərdən keçirmə/redaktə Y.Y. və Y.K.; vizuallaşdırma, nəzarət və layihə idarəçiliyi Y.Y.; maliyyələşdirmənin əldə edilməsi H.-M.T., Y.Y. və Y.K. olaraq bildirilmişdir.

Hesablama infrastrukturu: Tədqiqat Chinese Academy of Sciences Hefei Branch Supercomputing Center, Hefei Advanced Computing Center və Tohoku University Institute for Materials Research nəzdindəki hesablama resurslarından istifadə etmişdir; HPCI System Research Project hp200246 də qeyd edilmişdir.

Mənbədaxili uyğunsuzluq 1: Abstract/əsas nəticələr və Cədvəl 1 təxminən \(10^4m_0\) maksimum mərtəbəsini dəstəklədiyi halda, nəticə bölməsi \(10^3\)–\(10^5m_0\) ifadəsini istifadə edir. Bu fərq düzəldilmədən görünən şəkildə qorunmuşdur.

Mənbədaxili uyğunsuzluq 2: Əlavə A, \(k\) intervalını 0–0.03 kimi müəyyənləşdirərkən fitləmə seqmentlərini 0–0.02, 0.01–0.03 və 0.02–0.04 kimi verir.

Mənbədaxili hesabat qeyri-müəyyənliyi 3: B və δ fazalarının bəzi şəkil altyazılarında GGA və HSE birlikdə qeyd olunsa da, əsas nəticə mətni və görünən əyrilər əsasən GGA hesablamalarını təqdim edir. Verianla bu fazalar üçün mənbədə açıq şəkildə verilməyən HSE rəqəmləri yaratmamışdır.

Əsas şərh sərhədi: Tədqiqat çox böyük zolaq-effektiv kütlələrinə əsaslanaraq “possible heavy-fermion-like behavior” təklif edir. Bu nəticə ağır fermionun eksperimental müşahidəsi və ya klassik ağır fermion many-body fizikasının tam təsdiqi deyil. Siklotron rezonansı nəzəri proqnozu sınamaq üçün gələcəkdə tətbiqi təklif olunan eksperimentdir.