Moreau-Yosida-based Kohn-Sham Inversion for Periodic Systems

तकनीकी सारांश: आवधिक प्रणालियों के लिए मोरो–योसिडा-आधारित कोहन–शाम (Kohn–Sham) व्युत्क्रमण

समस्या विवरण
घनत्व-फलक सिद्धांत (density-functional theory - DFT) इलेक्ट्रॉनिक संरचना गणनाओं का एक आधार स्तंभ है, फिर भी सटीक विनिमय-सहसंबंध (exchange-correlation - xc) फलकों (functionals) का विकास एक महत्वपूर्ण चुनौती बना हुआ है। एक प्रमुख कठिनाई सटीक सार्वभौमिक घनत्व फलक और व्यापक रूप से उपयोग किए जाने वाले सन्निकटन (approximations) के बीच के संबंध में सीमित गणितीय अंतर्दृष्टि है। इन सन्निकटनों को सुधारने के लिए एक महत्वपूर्ण चरण "विपरीत कोहन–शाम (inverse Kohn–Sham - KS) समस्या" है: एक ज्ञात ग्राउंड-स्टेट घनत्व से प्रभावी KS विभव (विशेष रूप से xc विभव) का पुनर्निर्माण करना। जबकि कई व्युत्क्रमण योजनाएं मौजूद हैं, आवधिक प्रणालियों (ठोस-अवस्था सामग्रियों) के लिए सुदृढ़ और कुशल संख्यात्मक विधियां एक खुला प्रश्न बनी हुई हैं। मौजूदा दृष्टिकोण अक्सर एक कठोर कार्यात्मक-विश्लेषणात्मक आधार की कमी रखते हैं या आवधिक परिवेश में स्थिरता और अभिसरण (convergence) के साथ संघर्ष करते हैं।

कार्यप्रणाली
लेखक मोरो–योसिडा (Moreau–Yosida - MY) नियमितीकरण पर आधारित आवधिक प्रणालियों के लिए घनत्व-विभव व्युत्क्रमण हेतु एक गणितीय रूप से कठोर ढांचा विकसित करते हैं। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

1. कार्यात्मक-विश्लेषणात्मक ढांचा: सिद्धांत को आवधिक समरूप सोबोलेव स्थान (periodic homogeneous Sobolev space), $H^{-1}_{\text{per,hom}}(\Omega)$ के द्वैत (dual) के भीतर, और इसके द्वैत, $H^{1}_{\text{per,hom}}(\Omega)$ के लिए, घनत्वों के रूप में तैयार किया गया है। इन स्थानों का चयन महत्वपूर्ण है क्योंकि ये स्वाभाविक रूप से माध्य-क्षेत्र हार्ट्री ऊर्जा (mean-field Hartree energy) और विभवों की गेज स्वतंत्रता (gauge freedom) को संभालते हैं (जो एक स्थिरांक तक परिभाषित है)।
2. मोरो–योसिडा नियमितीकरण: सार्वभौमिक घनत्व फलक (विशेष रूप से गतिज-ऊर्जा बाधित-खोज फलक) को MY आवरण (envelope) का उपयोग करके नियमित किया जाता है। इसमें $F(\sigma) + \frac{1}{2\varepsilon}\|\sigma - \rho_{\text{gs}}\|^2_X$ के रूप में एक फलक का न्यूनीकरण शामिल है, जहाँ $\rho_{\text{gs}}$ लक्षित ग्राउंड-स्टेट घनत्व है और $\varepsilon > 0$ एक नियमितीकरण पैरामीटर है।
3. समीपस्थ मानचित्रण (Proximal Mapping): नियमित किए गए फलक का न्यूनाधिकतम (minimiser) "समीपस्थ घनत्व" $\rho_\varepsilon$ है। प्रभावी KS विभव को एक सीमा प्रक्रिया के माध्यम से पुनः प्राप्त किया जाता है:
   $$v_{\text{KS}} = \lim_{\varepsilon \to 0^+} \frac{1}{\varepsilon} J(\rho_\varepsilon - \rho_{\text{gs}})$$
   जहाँ $J$ घनत्व और विभव स्थानों के बीच द्वैत मानचित्रण (duality mapping) है।
4. सैद्धांतिक गारंटी: एक प्रमुख विश्लेषणात्मक योगदान यह प्रमाण है कि गतिज-ऊर्जा बाधित-खोज फलक चुनी गई टोपोलॉजी में निम्न-अर्ध-सतत (lower semicontinuous - l.s.c.) है। यह समीपस्थ बिंदु (proximal point) के अस्तित्व और व्युत्क्रमण योजना की वैधता को सुनिश्चित करता है।
5. संख्यात्मक कार्यान्वयन: समीपस्थ घनत्व को तीन अलग-अलग दृष्टिकोणों का उपयोग करके परिकलित किया जाता है:
   • प्रत्यक्ष न्यूनीकरण (Direct Minimisation): ऑर्बिटल्स पर ऊर्जा फलक का सीधे न्यूनीकरण करना (BFGS या न्यूटन विधियों का उपयोग करके)।
   • स्व-सुसंगत क्षेत्र (Self-Consistent Field - SCF): पुनरावृत्ति के माध्यम से KS-समान ऑयलर-लैग्रेंज समीकरणों को हल करना।
   • समीपस्थ बिंदु एल्गोरिदम (Proximal Point Algorithm): गैर-सुचारू उत्तल अनुकूलन (non-smooth convex optimisation) के लिए विशेष रूप से डिज़ाइन किया गया एक पुनरावृत्ति तंत्र।

प्रमुख योगदान
लेखक तीन प्राथमिक योगदान देते हैं:

1. सामान्य व्युत्क्रमण योजना: लेखक रिफ्लेक्सिव बेनाच स्पेस (reflexive Banach spaces) में घनत्वों के लिए एक सामान्य MY-आधारित व्युत्क्रमण एल्गोरिदम तैयार करते हैं। वे सिद्ध करते हैं कि यह योजना केवल फलक की निम्न-अर्ध-सततता और समीपस्थ मानचित्रण की गैर-संकुचनशीलता (non-expansiveness) पर निर्भर करती है, जो इनपुट घनत्व में त्रुटियों के प्रसार को सीमित करती है।
2. आवधिक प्रणालियों के लिए अनुप्रयोग: इस ढांचे को आवधिक समरूप सोबोलेव स्थानों के लिए विशिष्ट बनाया गया है। लेखक इस विशिष्ट टोपोलॉजी में गतिज-ऊर्जा फलक की निम्न-अर्ध-सततता को सिद्ध करते हैं, जिससे ज़ाओ-मोरिसन-पार (Zhao-Morrison-Parr - ZMP) विधि को पहली बार आवधिक परिवेश में एक ठोस कार्यात्मक-विश्लेषणात्मक आधार मिलता है। यह उप-विभेदक (subdifferential) में न्यूनतम मानदंड वाले तत्व के रूप में सटीक xc विभव के निर्धारण की अनुमति देता है।
3. संख्यात्मक विश्लेषण: यह शोध पत्र आवधिक एक-आयामी ग्रॉस-पिटाएव्स्की (Gross–Pitaevskii) समीकरण और तीन-आयामी बल्क सामग्रियों (सिलिकॉन, गैलियम आर्सेनाइड, सोडियम क्लोराइड, पोटेशियम क्लोराइड) का उपयोग करके व्युत्क्रमण योजना का व्यापक संख्यात्मक अन्वेषण प्रदान करता है।

परिणाम

• अभिसरण (Convergence): संख्यात्मक प्रयोग दर्शाते हैं कि समीपस्थ घनत्व $\rho_\varepsilon$, $\varepsilon \to 0$ होने पर, ग्राउंड-स्टेट घनत्व $\rho_{\text{gs}}$ की ओर दृढ़ता से (strongly) अभिसरित होता है। पुनर्निर्मित विभव संदर्भ xc विभव की ओर अभिसरित होता है।
• एल्गोरिदम प्रदर्शन: ग्रॉस–पिटाएव्स्की समीकरण के लिए, पूर्ण न्यूटन विधि (DFTK में कार्यान्वित) का उपयोग करके प्रत्यक्ष न्यूनीकरण ने SCF और क्वाजी-न्यूटन दृष्टिकोणों की तुलना में बेहतर अभिसरण दिखाया, विशेष रूप से छोटे $\varepsilon$ के लिए। समीपस्थ बिंदु एल्गोरिदम ने PDE सॉल्वर के समान अभिसरण गुण प्रदर्शित किए लेकिन इसमें अधिक कम्प्यूटेशनल ओवरहेड लगा।
• बल्क सामग्रियां: इस पद्धति ने विभिन्न बल्क सामग्रियों के लिए xc विभव को सफलतापूर्वक पुनर्निर्मित किया। परिणाम दर्शाते हैं कि अभिसरण दर लक्षित विभव की सुचारूता (smoothness) पर निर्भर करती; तीक्ष्ण विशेषताओं (जैसे परमाणु नाभिक के पास) को हल करने के लिए $\varepsilon$ के छोटे मानों की आवश्यकता होती है।
• मार्गदर्शक फलक (Guiding Functionals): अध्ययन ने मार्गदर्शक फलक की संरचना (हार्ट्री ऊर्जा, बाहरी विभव) के प्रभाव की जांच की। मार्गदर्शक फलक में बाहरी विभव को शामिल करने से पुनर्निर्मित विभव की सटीकता में महत्वपूर्ण सुधार हुआ।
• स्थिरता: व्युत्क्रमण योजना का परीक्षण विचलित इनपुट घनत्वों (फूरियर ट्रंकेशन द्वारा सिम्युलेटेड) के साथ किया गया। परिणामों ने दिखाया कि पुनर्निर्मित विभव की सटीकता इनपुट त्रुटि के परिमाण द्वारा सीमित है, जो नियमित किए गए दृष्टिकोण की स्थिरता की पुष्टि करती है।

महत्व
लेखक का दावा है कि इसकी प्राथमिक महत्ता आवधिक प्रणालियों के लिए KS व्युत्क्रमण हेतु एक गणितीय रूप से कठोर और संख्यात्मक रूप से स्थिर ढांचा प्रदान करने में निहित है। आवधिक समरूप सोबोलेव स्थानों में गतिज-ऊकी ऊर्जा फलक की निम्न-अर्ध-सततता स्थापित करके, लेखक सटीक विनिमय-सहसंबंध विभवों को पुनर्प्राप्त करने के लिए मोरो–योसिडा नियमितीकरण के उपयोग को मान्य करते हैं। यह कार्य अमूर्त उत्तल विश्लेषण (convex analysis) और व्यावहारिक DFT अनुप्रयोगों के बीच के अंतर को पाटता है, जो घनत्व फलकों के विश्लेषण और सुधार के लिए एक सुदृढ़ उपकरण प्रदान करता है। इसके अतिरिक्त, संख्यात्मक परिणाम सुझाव देते हैं कि MY-आधारित दृष्टिकोण मौजूदा व्युत्क्रमण योजनाओं के लिए एक व्यवहार्य विकल्प है, विशेष रूप से ठोस-अवस्था प्रणालियों के लिए जहाँ संख्यात्मक शोर और आवधिक सीमा स्थितियों के विरुद्ध मजबूती आवश्यक है।