Implementation and application of a DFT$+U$$+V$ approach within the all-electron FLAPW method

यहाँ "all-electron FLAPW विधि के भीतर एक DFT+U+V दृष्टिकोण का कार्यान्वयन और अनुप्रयोग" नामक शोध पत्र का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

मानक डेंसिटी फंक्शनल थ्योरी (DFT) अक्सर उन सामग्रियों का सटीक वर्णन करने में विफल रहती है जिनमें प्रबल इलेक्ट्रॉनिक सहसंबंध (strong electronic correlations) होते हैं, विशेष रूप से वे जिनमें स्थानीयकृत (localized) $d$- और $f$-इलेक्ट्रॉन शामिल होते हैं। जबकि DFT+U विधि स्व-अंतःक्रिया त्रुटियों (self-interaction errors) को सुधारने और इलेक्ट्रॉनों को स्थानीयकृत करने के लिए ऑनसाइट कूलम्ब प्रतिकर्षण ($U$) को संबोधित करती है, यह केवल एक सिंगल-साइट सुधार तक ही सीमित रहती है। यह इंटरसाइट कूलम्ब अंतःक्रियाओं ($V$) की उपेक्षा करती है, जो निम्नलिखित के वर्णन के लिए अत्यंत महत्वपूर्ण हैं:

• पड़ोसी साइटों के बीच आवेश पुनर्वितरण (charge redistribution) और कक्षीय संकरण (orbital hybridization)।
• आवेश विप्रकीर्णन (charge disproportionation) और आवेश व्यवस्था (charge ordering)।
• आवेश स्थानांतरण प्रक्रियाएं (जैसे, धातु और लिगैंड कक्षकों के बीच)।
• सहसंयोजक बंधनों वाले प्रणालियों (जैसे, ग्राफीन) और चार्ज-ट्रांसफर इंसुलेटरों (जैसे, NiO) में गैर-स्थानीय सहसंबंध प्रभाव।

DFT+U+V के मौजूदा कार्यान्वयन अक्सर स्यूडोपोटेंशियल (pseudopotentials) पर निर्भर करते हैं, जिनमें ऑल-इलेक्ट्रॉन विवरण के लिए आवश्यक सटीकता की कमी हो सकती है या जो जटिल आधार सेटों (basis sets) के साथ संघर्ष कर सकते हैं। एक कठोर, ऑल-इलेक्ट्रॉन कार्यान्वयन की आवश्यकता है जो विविध सामग्री वर्गों में ऑनसाइट और इंटरसाइट दोनों सहसंबंधों को सटीक रूप से संभाल सके।

2. कार्यप्रणाली (Methodology)

लेखक FLEUR कोड का उपयोग करते हुए, फुल-पोटेंशियल लीनियराइज्ड ऑगमेंटेड प्लेन-वेव (FLAPW) विधि के भीतर DFT+U+V औपचारिक पद्धति का प्रथम-सिद्धांत (first-principles) कार्यान्वयन प्रस्तुत करते हैं।

• औपचारिकता विस्तार (Formalism Extension):

  • मानक DFT+U ऊर्जा फंक्शनल को एक इंटरसाइट टर्म ($E_V$) को शामिल करने के लिए विस्तारित किया गया है।
  • इस सुधार में एक सामान्यीकृत ऑक्यूपेशन मैट्रिक्स ($n_{IJ}$) शामिल है जो दो अलग-अलग परमाणुओं के परमाणु आधार फलनों (atomic basis functions) पर कॉन न्यूमैन-कोहन-शाम (Kohn-Sham) तरंग फलनों को प्रोजेक्ट करता है।
  • कार्यान्वयन रेडियल फलनों ($u$ और $\dot{u}$) के रेखीयकरण (linearization) से उत्पन्न होने वाले मिश्रित पदों (A-B और B-A योगदान) को ध्यान में रखता है, जो सटीक इंटरसाइट कपलिंग के लिए आवश्यक हैं। इसके लिए मफिन-टिन फंक्शन्स (MTFs) को प्रोजेक्टर के रूप में उपयोग किया जाता है, जो FLAPW आधार सेट के अनुरूप है।
  • डबल-काउंटिंग सुधार को मीन-फील्ड सीमा (mean-field limit) में लागू किया गया है।

• पैरामीटर गणना (cRPA):

  • अंतःक्रिया मापदंडों ($U$ और $V$) को SPEX कोड के माध्यम से कन्स्ट्रेंड रैंडम-फेज एप्रोक्सिमेशन (cRPA) का उपयोग करके प्रथम-सिद्धांतों से प्राप्त किया गया है।
  • स्क्रीन किया गया कूलम्ब अंतःक्रिया दो अलग-अलग स्थानीयकृत आधार निरूपणों पर प्रोजेक्ट की जाती है ताकि उनकी प्रभावकारिता की तुलना की जा सके:
    1. मैक्सिमली लोकलाइज्ड वानियर फंक्शन्स (MLWFs): यह लो-एनर्जी मॉडल हैमिल्टोनियन के साथ एक सेतु प्रदान करता है और डेलोकलाइजेशन/हाइब्रिडाइजेशन प्रभावों को पकड़ता है।
    2. मफिन-टिन फंक्शन्स (MTFs): यह FLAPW आधार सेट के साथ सख्त निरंतरता सुनिश्चित करता है लेकिन मफिन-टिन त्रिज्या पर निर्भर करता है।
  • अध्ययन दोनों आधारों (जिन्हें w-cRPA और u-cRPA कहा जाता है) से प्राप्त मापदंडों की तुलना करता है।

• परीक्षण प्रणालियाँ (Test Systems):
  विभिन्न बॉन्डिंग और सहसंबंध शासनों का प्रतिनिधित्व करने वाली विविध सामग्रियों पर इस कार्यप्रणाली को मान्य किया गया है:

  1. ग्राफीन (Graphene): एक 2D सहसंयोजक प्रणाली जिसमें डेलोकलाइज्ड $\pi$-इलेक्ट्रॉन होते हैं।
  2. सिलिकॉन (Si) और जर्मेनियम (Ge): मजबूत $sp^3$ सहसंयोजक बंधन वाले बल्क मौलिक अर्धचालक।
  3. निकेल ऑक्साइड (NiO): एक प्रोटोटाइपिकल मॉट-चार्ज-ट्रांसफर इंसुलेटर जिसमें मजबूत $d$-$p$ हाइब्रिडाइजेशन होता है।

3. मुख्य योगदान (Key Contributions)

• ऑल-इलेक्ट्रॉन कार्यान्वयन: FLAPW फ्रेमवर्क (FLEUR) में DFT+U+V का पहला सफल एकीकरण, जो स्यूडोपोटेंशियल सन्निकटन के बिना उच्च-परिशुद्धता गणना सक्षम बनाता है।
• आधार सेट तुलना: cRPA पैरामीटर निष्कर्षण के लिए MTF बनाम वानियर फंक्शन प्रोजेक्टर का व्यवस्थित तुलनात्मक अध्ययन। अध्ययन दर्शाता है कि जबकि MTF पैरामीटर अपनी सघन स्थानीयकरण के कारण आम तौर पर बड़े होते हैं, वे मफिन-टिन त्रिज्या बढ़ने के साथ वानियर-आधारित मानों की ओर अभिसरित होते हैं।
• एकीकृत ढांचा (Unified Framework): यह प्रदर्शित करता है कि एक एकल ढांचा कमजोर रूप से सहसंबंधित 2D सामग्रियों (ग्राफीन) से लेकर दृढ़ रूप से सहसंबंधित चार्ज-ट्रांसफर इंसुलेटरों (NiO) तक के तंत्रों का प्रभावी ढंग से वर्णन कर सकता है।
• बेंचमार्किंग: प्रयोगात्मक डेटा, पिछले DFT+U+V स्यूडोपोटेंशियल अध्ययनों और उच्च-स्तरीय GW (ग्रीन फंक्शन) गणनाओं के विरुद्ध व्यापक सत्यापन।

4. मुख्य परिणाम (Key Results)

A. ग्राफीन (2D सहसंयोजक प्रणाली)

• फर्मी वेलोसिटी ($v_F$): मानक DFT $v_F$ को काफी कम आंकता है। DFT+U इसमें थोड़ा सुधार करता है, लेकिन DFT+U+V $10.74–11.11 \times 10^5$ m/s के मान देता है, जो प्रयोग (11.00 $\times 10^5$ m/s) के साथ उत्कृष्ट सहमति में है।
• बैंड संरचना: $V$ का समावेश डिरैक कोन (Dirac cone) के ढलान को ठीक करता है। कच्चे $U$ और $V$ मानों में महत्वपूर्ण अंतर के बावजूद, दोनों u-cRPA और w-cRPA पैरामीटर समान परिणाम देते हैं, जो मजबूती (robustness) का संकेत देते हैं।
• अंतर्दृष्टि: 2D सहसंयोजक प्रणालियों में सही क्वैसीपार्टिकल डिस्पर्शन को पकड़ने के लिए गैर-स्थानीय अंतःक्रियाएं महत्वपूर्ण हैं।

B. बल्क Si और Ge (सहसंयोजक अर्धचालक)

• लैटिस स्थिरांक और मॉड्युली: मानक DFT लैटिस स्थिरांक को थोड़ा अधिक आंकता है। केवल $U$ जोड़ने से लैटिस विस्तार और सॉफ्टनिंग (सहमति बिगड़ना) होती है। $V$ जोड़ने से यह प्रभाव संतुलित होता है, जिससे लैटिस स्थिरांक और बल्क मॉड्युली प्रयोगात्मक मानों के करीब आ जाते हैं।
• बैंड गैप: DFT बैंड गैप को बहुत कम आंकता है। DFT+U गैप को खोलता है लेकिन अक्सर संरचनात्मक सटीकता की कीमत पर। DFT+U+V संरचनात्मक मापदंडों में सुधार करते हुए सफलतापूर्वक गैप को खोलता है।
• पैरामीटर संवेदनशीलता: अध्ययन नोट करता है कि Si और Ge के लिए, GW परिणामों से मेल खाने के लिए $U$ के तुल्य $V$ मानों की आवश्यकता होती है, जो सुझाव देता है कि $sp^3$ प्रणालियों के लिए केवल निकटतम-पड़ोसी अंतःक्रियाएं पर्याप्त नहीं हो सकती हैं और अगले-निकटतम-पड़ोसी पदों की आवश्यकता हो सकती है।

C. NiO (चार्ज-ट्रांसफर इंसुलेटर)

• इलेक्ट्रॉनिक संरचना: मानक DFT वैलेंस बैंड के शीर्ष पर सही कक्षीय लक्षण (orbital character) को दोहराने में विफल रहता है (Ni $d$-प्रभुत्व के बजाय O $p$-प्रभुत्व की भविष्यवाणी करता है)। DFT+U गैप में सुधार करता है लेकिन स्पेक्ट्रल वेट वितरण के साथ संघर्ष करता है।
• DFT+U+V की सफलता: Ni और O के बीच इंटरसाइट $V$ को शामिल करने से वैलेंस बैंड मैक्सिमम पर O $p$-चरित्र और कंडक्शन बैंड मिनिमम पर Ni $d$-चरित्र सही ढंग से प्राप्त होता है। यह -7 eV पर सैटेलाइट पीक को भी कैप्चर करता है, जो प्रयोगात्मक XPS डेटा के अनुरूप है।
• चुंबकीय गुण: यह विधि चुंबकीय क्षण (magnetic moment) की सटीक भविष्यवाणी करती है और AFM-2 ग्राउंड स्टेट को स्थिर करती है। $V$ को बढ़ाने से सुपerexchange अंतःक्रियाएं मजबूत होती हैं, जिससे एंटीफेरोमैग्नेटिक चरण और अधिक स्थिर होता है, हालांकि हाइब्रिडाइजेशन बढ़ने के कारण चुंबकीय क्षणों में थोड़ी कमी आती है।
• तुलना: DFT+U+V संरचनात्मक और इलेक्ट्रॉनिक गुणों के लिए DFT+U से बेहतर प्रदर्शन करता है और विभिन्न हाइब्रिड फंक्शनल (PBE0, HSE06) की सटीकता के बराबर या उससे अधिक है।

5. महत्व और भविष्य की दिशा (Significance and Outlook)

• सटीकता बनाम लागत: DFT+U+V दृष्टिकोण GW गणनाओं की उच्च सटीकता और DFT की कम्प्यूटेशनल दक्षता के बीच एक बेहतर समझौता प्रदान करता है, जिससे यह जटिल सामग्रियों के लिए स्केलेबल बन जाता है।
• भौतिक अंतर्दृष्टि: यह कार्य पुष्टि करता है कि इंटरसाइट कूलम्ब अंतःक्रियाएं केवल मॉट इंसुलेटर के लिए एक सुधार नहीं हैं, बल्कि सहसंयोजक और मिश्रित-बॉन्डिंग प्रणालियों में चार्ज ट्रांसफर, हाइब्रिडाइजेशन और गैर-स्थानीय सहसंबंधों का वर्णन करने के लिए आवश्यक हैं।
• भविष्य की दिशाएं: लेखक हबार्ड मापदंडों के स्व-सुसंगत (self-consistent) मूल्यांकन को लागू करने और बड़े ऊर्जा रेंज में बेहतर विवरण के लिए आवृत्ति-निर्भर (frequency-dependent/dynamical) $U$ और $V$ का पता लगाने की योजना बना रहे हैं।

निष्कर्षतः, यह शोध पत्र एक सुदृढ़ ऑल-इलेक्ट्रॉन DFT+U+V ढांचे को स्थापित करता है जो विविध प्रकार की सहसंबंधित सामग्रियों के लिए इलेक्ट्रॉनिक संरचना गणनाओं की भविष्य कहने वाली शक्ति को महत्वपूर्ण रूप से बढ़ाता है, जो मानक DFT और उच्च-स्तरीय मेनी-बॉडी परटर्बेशन थ्योरी के बीच के अंतर को पाटता है।