First-principles study of bulk stacking, $J_{\rm eff}$ picture, magnetic Hamiltonian, $g$ factors, and structural distortions of $α$-RuCl$_3$

तकनीकी सारांश: $\alpha$-RuCl$_3$ के बल्क स्टैकिंग, $J_{\text{eff}}$ चित्र, चुंबकीय हैमिल्टनियन, $g$ कारक और संरचनात्मक विरूपण का प्रथम-सिद्धांत (First-Principles) अध्ययन

समस्या विवरण
$\alpha$-RuCl$_3$ किटाएव क्वांटम स्पिन लिक्विड अवस्था को साकार करने के लिए एक प्राथमिक उम्मीदवार है, फिर भी इसके मौलिक चुंबकीय गुण बहस का विषय बने हुए हैं। प्रमुख अनसुलझे मुद्दों में शामिल हैं:

1. बल्क संरचना: जबकि हालिया प्रयोगों से संकेत मिलता है कि निम्न-तापमान चरण $R\bar{3}$ स्पेस ग्रुप से संबंधित है (जिसमें $\sim$150 K पर $C2/m$ में संक्रमण होता है), इससे पहले कोई सैद्धांतिक अध्ययन नहीं हुआ था जो $R\bar{3}$ संरचना को निम्न-तापमान चरण के रूप में ऊर्जावान रूप से अनुकूल सिद्ध कर सके।
2. इलेक्ट्रॉनिक संरचना: $J_{\text{eff}} = 1/2$ चित्र की वैधता के लिए सत्यापन की आवश्यकता है, विशेष रूप से नेल (Néel) वेक्टर के सापेक्ष कोणीय संवेग क्वांटाइजेशन अक्ष के ओरिएंटेशन के संबंध में।
3. चुंबकीय हैमिल्टनियन: पिछले अध्ययन चुंबकीय विनिमय मापदंडों (जैसे, निकटतम पड़ोसी किटाएव विनिमय $K_1$ का चिह्न) में महत्वपूर्ण विसंगतियों की रिपोर्ट करते हैं और अक्सर द्वितीय-निकटतम पड़ोसी इंटरैक्शन की उपेक्षा करते हैं।
4. $g$-कारक अनिसोट्रॉपी (Anisotropy): $g$-कारक की अनिसोट्रॉपी ($g_{XY}$ बनाम $g_Z$) पर प्रयोगात्मक रिपोर्ट विरोधाभासी हैं, जो कम से लेकर उच्च अनिसोट्रॉपी तक विस्तृत हैं। इसके अलावा, इन कारकों की गणना करने वाली सैद्धांतिक विधियाँ (विशेष रूप से परमाणु कक्षीय प्रोजेक्शन) अपर्याप्त हो सकती हैं।
5. संरचनात्मक विरूपण: RuCl$_6$ अष्टफलक (octahedron) में विशिष्ट संरचनात्मक विरूपणों (त्रिकोणीय बनाम Cl त्रिकोणों का सापेक्ष घुमाव/twist) का चुंबकीय अनिसोट्रॉपी पर प्रभाव व्यवस्थित रूप से अलग नहीं किया गया है।

कार्यप्रणाली (Methodology)
लेखकों ने Quantum ESPRESSO फ्रेमवर्क के भीतर Local Density Approximation (LDA) और रोटेशनली इनवेरिएंट LDA+U विधियों का उपयोग करके कन्स्ट्रेंड डेंसिटी फंक्शनल थ्योरी (cDFT) का उपयोग किया।

• संरचनात्मक विश्राम (Structural Relaxation): विभिन्न हबार्ड $U$ मानों और विश्राम प्रकारों के तहत प्रयोगात्मक, फिक्स्ड-सेल और वेरिएबल-सेल विश्राम स्थितियों के लिए $C2/m$ और $R\bar{3}$ स्टैकिंग कॉन्फ़िगरेशन के लिए गणनाएँ की गईं ताकि ऊर्जा स्थिरता निर्धारित की जा सके।
• $J_{\text{eff}}$ विश्लेषण: तरंग फलनों (wavefunctions) को $J_{\text{eff}}$ आधार अवस्थाओं पर प्रोजेक्ट करके इलेक्ट्रॉनिक संरचना का विश्लेषण किया गया। महत्वपूर्ण रूप से, कोणीय संवेग क्वांटाइजेशन अक्ष को नेल वेक्टर के समानांतर संरेखित किया गया ताकि $J_{\text{eff}}$ चित्र का परीक्षण किया जा सके।
• चुंबकीय हैमिल्टनियन निर्माण: लेखकों ने विभिन्न चुंबकीय कॉन्फ़िगरेशन (जिगज़ैग, नेल, स्ट्रिपी, फेरोमैग्नेटिक, और एक विशिष्ट 1-बाय-4 कॉन्फ़िगरेशन) से कुल ऊर्जा को तीन मॉडल हैमिल्टोनियन पर फिट किया:
  1. पारंपरिक किटाएव-हाइजेनबर्ग (KH) मॉडल।
  2. स्पिन-आइस (SI) मॉडल।
  3. एक नया अनिसोट्रोपिक $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) मॉडल, जिसमें अनिसोट्रोपिक द्वितीय- और तृतीय-निकटतम पड़ोसी इंटरैक्शन शामिल हैं।
• $g$-कारक गणना: दो विधियों की तुलना की गई:
  1. परमाणु कक्षीय प्रोजेक्शन (Atomic Orbital Projection): हबार्ड मैनिफोल्ड पर मोमेंट्स को प्रोजेक्ट करना।
  2. मोमेंट कैंटिंग विद वैनियर इंटरपोलेशन (Moment Canting with Wannier Interpolation): $J_{\text{eff}}$ मोमेंट को $1^\circ$ झुकाकर ट्रांसलेशनलली इक्विवेरिएंट वैनियर इंटरपोलेशन (WannierBerri) का उपयोग करके कुल चुंबकीय मोमेंट में परिवर्तन की गणना करना।
• विरूपण विश्लेषण (Distortion Analysis): एक पूर्ण अष्टफलक से दो विशिष्ट विरूपणों को अलग किया गया: त्रिकोणीय विरूपण (Ru-Cl-Z कोण को बदलना) और सापेक्ष घुमाव (XY प्ले में ऊपरी और निचले Cl त्रिकोणों को घुमाना)।

प्रमुख योगदान और परिणाम

1. $R\bar{3}$ स्टैकिंग की पुष्टि:
   अध्ययन प्रदर्शित करता है कि विभिन्न हबार्ड $U$ मानों और विश्राम प्रकारों के तहत $R\bar{3}$ संरचना $C2/m$ संरचना की तुलना में लगभग 5–6 meV/Ru कम ऊर्जा वाली है। यह प्रयोगात्मक रूप से देखे गए निम्न-तापमान $R\beta{3}$ चरण के लिए पहला सैद्धांतिक समर्थन प्रदान करता है।

2. परिष्कृत $J_{\text{eff}}$ चित्र:
   कंडक्शन बैंड मिनिमम अवस्थाओं को लगभग शुद्ध $J_{\text{eff}} = 1/2$ और $m_{\text{eff}} = -1/2$ चरित्र के रूप में दिखाया गया है, लेकिन केवल तभी जब कोणीय संमरक्षण क्वांटाइजेशन अक्ष को नेल वेक्टर के समानांतर सेट किया जाता है। यह परिप्रेक्ष्य, जिसे पहले उपेक्षित किया गया था, अवस्थाओं के कक्षीय चरित्र को सही ढंग से पहचानने के लिए आवश्यक है। हाइब्रिडाइजेशन और विरूपण के कारण $J_{\text{eff}}$ मोमेंट का परिमाण आदर्श $1/2$ से थोड़ा कम होने के बावजूद, कॉन्फ़िगरेशन में मजबूत पाया गया।

3. चुंबकीय हैमिल्टनियन और विनिमय पैरामीटर:

• पारंपरिक KH और स्पिन-आइस मॉडल DFT कुल ऊर्जा (विशेष रूप से नेल कॉन्फ़िगरेशन के लिए) को सटीक रूप से पुनरुत्पादित करने में विफल रहते हैं।
• प्रस्तावित A-$J_1J_2J_3$ मॉडल (और एक सरलीकृत A-$J_1J_2$-I-$J_3$ संस्करण) DFT ऊर्जाओं को सटीक रूप से पुनरुत्पादित करता है।
• द्वितीय-निकटतम पड़ोसी इंटरैक्शन महत्वपूर्ण हैं। द्वितीय-निकटतम पड़ोसी विनिमय के अनिसोट्रोपिक घटक ($K_2, \Gamma_2, \Gamma'_2$) नगण्य नहीं हैं और परिमाण में, उन अनिसोट्रोपिक निकटतम पड़ोसी पदों ($\Gamma'_1$) से बड़े हैं जिन्हें अक्सर पिछले मॉडलों में माना जाता है। इन पदों को शामिल करने से पारंपरिक मॉडल की तुलना में कुल ऊर्जा में रूट-मीन-स्क्वायर (RMS) त्रुटि में 3.6 गुना कमी आई।

4. $g$-कारक अनिसोट्रॉपी और गणना विधि:

• मोमेंट-कैंटिंग विधि को वैनियर इंटरपोलेशन के साथ मिलाने पर $g_{XY} \approx 2.22$ और $g_Z \approx 1.84$ के $g$-कारक प्राप्त होते हैं। ये मान एक छोटी अनिसोट्रॉपी का संकेत देते हैं, जो हालिया प्रयोगात्मक रिपोर्टों (जैसे, लोइडल एट अल.) के अनुरूप है और उच्च अनिसोट्रॉपी के पिछले दावों का खंडन करती है।
• परमाणु कक्षीय प्रोजेक्शन विधि $g$-कारक के परिमाण को काफी कम आंकती है ($g_{XY} \approx 1.54$), जिससे कुल चुंबकीय मोमेंट का एक गैर-नगण्य हिस्सा छूट जाता है। लेखक निष्कर्ष निकालते हैं कि सटीक कक्षीय चुंबकत्व गणना के लिए वैनियर इंटरपोलेशन दृष्टिकोण आवश्यक है।

5. संरचनात्मक विरूपण की भूमिका:
   विशिष्ट विरूपणों को अलग करके, अध्ययन पाता है कि ऊपरी और निचले Cl त्रिकोणों का सापेक्ष घुमाव (XY प्ले के भीतर) चुंबकीय अनिसोट्रॉपी और सही ग्राउंड-स्टेट मोमेंट दिशा ($\theta \sim 60^\circ$) को निर्धारित करने वाला प्रमुख कारक है। केवल त्रिकोणीय विरूपण प्रयोगात्मक ग्राउंड स्टेट को पुनरुत्पादित करने में विफल रहता है। यह ट्विस्ट विरूपण, जो बाहरी विक्षोभ के बिना बल्क सामग्री में मौजूद है, पिछले सैद्धांतिक मॉडलों में उपेक्षित रहा है।

महत्व
यह पेपर $\alpha$-RuCl$_3$ के संरचनात्मक और चुंबकीय विवरणों के लंबे समय से चले आ रहे विरोधाभासों को हल करके एक व्यापक प्रथम-सिद्धांत समझ प्रदान करने का दावा करता है। यह:

• $R\bar{3}$ निम्न-तापमान संरचना को सैद्धांतिक रूप से मान्य करता है।
• स्थापित करता है कि $J_{\text{eff}}$ चित्र को सही ढंग से व्याख्या करने के लिए क्वांटाइजेशन अक्ष को नेल वेक्टर के साथ संरेखित करना आवश्यक है।
• यह प्रदर्शित करता है कि सटीक चुंबकीय मॉडलिंग के लिए अनिसोट्रोपिक द्वितीय-निकटतम पड़ोसी विनिमय इंटरैक्शन को शामिल करना आवश्यक है।
• $g$-कारकों की गणना के लिए कार्यप्रणाली को सुधारता है, यह दिखाते हुए कि परमाणु कक्षीय प्रोजेक्शन अपर्याप्त है और वास्तविक अनिसोट्रॉपी छोटी है।
• Cl त्रिकोणों के सापेक्ष घुमाव को एक महत्वपूर्ण, पहले से उपेक्षित संरचनात्मक पैरामीटर के रूप में पहचानता है, जो चुंबकीय अनिसोट्रॉपी को नियंत्रित करता है, जिससे किटाएव सामग्रियों में चुंबकत्व को इंजीनियर करने का एक नया मार्ग प्रशस्त होता है।

लेखक निष्कर्ष निकालते हैं कि उनका कार्य $\alpha$-RuCl$_3$ के चुंबकत्व के प्रति एक अलग दृष्टिकोण प्रदान करता है और सुझाव देता है कि भविष्य के सैद्धांतिक मॉडलों और उच्च-दबाव प्रयोगों को सटीक व्याख्या के लिए इन विशिष्ट संरचनात्मक विरूपणों और विनिमय पदों को शामिल करना चाहिए।