Semilocal exchange functionals from the exact-exchange condition for the hydrogen atom: Hydrogenic exactness and recovery of Rydberg-like bound states

तकनीकी सारांश: हाइड्रोजन परमाणु के लिए सटीक-विनिमय (Exact-Exchange) स्थिति से प्राप्त अर्ध-स्थानीय (Semilocal) विनिमय फलन

समस्या विवरण
मानक अर्ध-स्थानीय घनत्व फलन (GGAs जैसे PBE और meta-GGAs जैसे SCAN) एक मौलिक संरचनात्मक दोष से ग्रस्त हैं: परिमित प्रणालियों के बाहरी क्षेत्रों में उनका विनिमय विभव (exchange potentials, $v_x$) घातांकीय रूप से (exponentially) क्षय होता है, जिससे वे सटीक कूलम्बिक $-1/r$ पूंछ (tail) को पुनरुत्पादित करने में विफल रहते हैं। यह कमी विफलताओं के एक सुप्रसिद्ध परिवार का कारण बनती है, जिसमें बाउंड रिडबर्ग अवस्थाओं (bound Rydberg states) की अनुपस्थिति, अस्थिर ऋणात्मक आयन, व्यवस्थित रूप से कम अनुमानित आयनीकरण ऊर्जाएं, और स्व-अंतःक्रिया त्रुटि (self-interaction error, SIE) के कारण अत्यधिक आवेश विस्थानीकरण (charge delocalization) शामिल हैं। जबकि मॉडल विभव (जैसे LB94, mBJ) और रेंज-सेपरेटेड हाइब्रिड इन समस्याओं को संबोधित करते हैं, पूर्ववर्ती (मॉडल विभव) में विचर्य संगतता (variational consistency) का अभाव होता है (कोई मूल ऊर्जा फलन नहीं), और उत्तरवर्ती (हाइब्रिड) में उच्च गणनात्मक लागत आती है। मौजूदा अर्ध-स्थानीय प्रयास, जैसे कि आर्मिएंटो-कुमेल (AK13) फलन, एसिम्प्टोटिक्स (asymptotics) में सुधार करते हैं लेकिन वे हाइड्रोजन परमाणु के लिए सटीक स्थितियों से व्युत्पन्न नहीं हैं।

यह कार्य गिल और पॉपली (GP93) द्वारा 1993 में प्रतिपादित साधारण अवकल समीकरण (ordinary differential equation) की पुनरावृत्ति करता है, जो यह परिभाषित करता है कि एक GGA को हाइड्रोजन परमाणु के सटीक विनिमय विभव को पुनरुत्पादित करने के लिए किस स्थिति को संतुष्ट करना चाहिए। लेखक इस समीकरण को एक 'इनवर्स प्रॉब्लम' के रूप में पूर्णतः हल करते हैं ताकि एक अद्वितीय एन्हांसमेंट फैक्टर (enhancement factor) प्राप्त किया जा सके, जिसे GP93 फैक्टर ($F_{GP93}$) कहा जाता है।

कार्यप्रणाली तीन परतों में आगे बढ़ती है:

1. हाइड्रोजेनिक सटीकता (Hydrogenic Exactness): GP93 फैक्टर को हाइड्रोजन 1s घनत्व पर सटीक बनाया गया है। यूनिफॉर्म कोऑर्डिनेट स्केलिंग द्वारा, यह बिना किसी अनुभवजन्य पैरामीटर के सभी हाइड्रोजेनिक 1s आयनों के लिए सटीक आइजनवैल्यू ($\epsilon_{1s} = -Z^2/2$) प्रदान करता है।
2. एसिम्प्टोटिक व्यवहार (Asymptotic Behavior): बड़े ग्रेडिएंट शासन (large-gradient regime) में, GP93 फैक्टर $O[s(\ln s)^{2/3}]$ के रूप में अनबाउंड विकास प्रदर्शित करता है, जो AK13 के $s \ln s$ विकास से भिन्न है। यह विशिष्ट विकास यह सुनिश्चित करता है कि विनिमय विभव हाइड्रोजन घनत्व पर सटीक $-1/r$ पूंछ को पुनरुत्पादित करे।
3. स्विचिंग तंत्र (Switching Mechanism): GP93 फैक्टर के डाइवर्जेंट व्यवहार को बहु-इलेक्ट्रॉन प्रणालियों को अस्थिर करने से रोकने के लिए, लेखकों ने इसे एक "स्विच्ड" मेटा-GGA रूप में समाहित किया है।
   • $\alpha$-स्विच: यह एक-इलेक्ट्रॉन जैसी क्षेत्रों का पता लगाने के लिए इंडिकेटर $\alpha = (\tau - \tau^W)/\tau_{unif}$ का उपयोग करता है। GP93 टर्म केवल वहीं सक्रिय होती है जहाँ $\alpha \approx 0$ (एकल-कक्षक क्षेत्र) हो।
   • $q$-स्विच (ट्रिपल स्विच): परमाणु पूंछों (जहाँ पूंछ की आवश्यकता होती है) और सहसंयोजक बंध केंद्रों (जहाँ पूंछ अस्थिरता का कारण बनती है) के बीच अंतर करने के लिए, एक न्यूनीकरण लैपलेसियन इंडिकेटर $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$ पेश किया गया है। परमाणु पूंछों में $q \to +\infty$ होता है, जबकि बंध केंद्रों में $q \lesssim 0$ होता है। एक स्विच फंक्शन $D(q)$ सहसंयोजक बंध केंद्रों में GP93 टर्म को दबा देता है जबकि परमाणु पूंछों में इसे सुरक्षित रखता है।
   • कक्षक स्वतंत्रता (Orbital Independence): गतिज-ऊर्जा घनत्व $\tau$ को पर्ड्यू-कॉन्स्टेंटिन (PC07) निर्माण पर आधारित एक ऑर्बिटल-फ्री गतिज-ऊर्जा-घनत्व फलन (KEDF) द्वारा प्रतिस्थापित किया गया है, जो $\tau^W \le \tau$ को लागू करता है। यह संपूर्ण फलन को एक शुद्ध घनत्व फलन बनाता है, जिससे सामान्यीकृत कोहन-श}(\text{generalized Kohn–Sham, gKS}) सॉल्वर की आवश्यकता समाप्त हो जाती है।

मुख्य परिणाम

• हाइड्रोजेनिक सटीकता: फलन हाइड्रोजेनिक आयनों ($Z=1$ से $5$) के लिए मिलीहार्ट्री स्तर पर सटीक आइजनवैल्यू$\epsilon_{1s} = -Z^2/2$ को पुनरुत्पादित करता है।
• रिडबर्ग-जैसी बाउंड अवस्थाएँ: उन प्रणालियों में जिनमें एक-इलेक्ट्रॉन जैसी $s$-शेल HOMO होती हैं (H, He, और Li/Na/K बड़े-कोर स्यूडोपोटेंशियल के साथ), यह फलन एक श्रृंखला को बाउंड वर्चुअल कोहन-श $(\text{virtual Kohn–Sham})$ अवस्थाओं (4–8 स्तर) को पुनरुत्पादित करता है। इसके विपरीत, मानक फलन (PBE, SCAN, TPSS, M06-L) और ग्लोबल हाइब्रिड (B3LY light, PBE0) इन प्रणालियों में शून्य रिडबर्ग अवस्थाओं को बाउंड करते हैं।
• आयनीकरण ऊर्जाएं: परीक्षण किए गए $s$-शेल सिस्टम के लिए, उच्चतम-अधिकतम आइजनवैल्यू ($-\epsilon_{HOMO}$) प्रायोगिक आयनीकरण विभवों के काफी करीब है (औसत पूर्ण त्रुटि $\approx 0.80$ eV), जो PBE, SCAN, या B3LYP की तुलना में बेहतर है।
• चयनात्मकता (Selectivity): फलन $p$-शेल परमाणुओं (Ne, Ar) के लिए GP93 टर्म को सही ढंग से दबा देता है जहाँ $\alpha > 0$, जिसके परिणामस्वरूप शून्य बाउंड रिडबर्ग अवस्थाएँ प्राप्त होती हैं, जो डिज़ाइन के अनुरूप है।
• आणविक स्थिरता: ट्रिपल स्विच सहसंयोजक बंध केंद्रों (जैसे $H_2, N_2$) में संख्यात्मक डाइवर्जेंस को सफलतापूर्वक दबा देता है जहाँ अनस्विच्ड GP93 फैक्टर विफल हो जाता, जबकि परमाणु प्रणालियों के लिए सही एसिम्प्टोटिक टेल बनाए रखता है।
• सहसंबंध अनुकूलता (Correlation Compatibility): फलन को LYP सहसंबंध फलन के साथ जोड़ा गया है, जो एक-इलेक्ट्रॉन प्रणालियों के लिए शून्य हो जाता है, जिससे विनिमय घटक की हाइड्रोजेनिक सटीकता बनी रहती है।

महत्व और दावे
लेख दावा करता है कि उन्होंने एक गैर-अनुभवजन्य (nonempirical), अर्ध-स्थानीय विनिमय फलन निर्मित किया है जो एक-इलेक्ट्रॉन जैसी $s$-शेल HOMOs वाली प्रणालियों के लिए सटीक $-1/r$ पूंछ और बाउंड रिडबर्ग अवस्थाओं को पुनरुत्पादित करता है, जो SCAN सहित सभी मौजूदा अर्ध-स्थानीय फलनों की गुणात्मक विफलता है।

लेखक अपने दावों के दायरे को स्पष्ट रूप से सीमित करते हैं:

• सफलता का डोमेन: फलन को एक-इलेक्ट्रॉन गुणों (रिडबर्ग अवस्थाएं, एसिम्प्टोटिक्स, हाइड्रोजेनिक स्तर, SIE) के लिए डिज़ाइन किया गया है जो एक-इलेक्ट्रॉन जैसी $s$-शेल HOMOs वाली प्रणालियों में होते हैं।
• सीमाएं:
  • ऊष्मागतिकी (Thermochemistry): यह फलन परमाणुकरण ऊर्जाओं (atomization energies) को लक्षित नहीं करता है। हाइड्रोजेनिक सटीकता को लागू करने से H परमाणु की पूर्ण ऊर्जा उन त्रुटि-क्षतिपूरक मानों से दूर चली जाती है जो PBE में अनुकूलित किए गए थे, जिससे यह H से जुड़ी ऊष्मागतिकी के लिए अनुपयुक्त हो जाता है।
  • $p$-शेल प्रणालियाँ: एसिम्प्टोटिक सुधार $p$-शेल HOMOs (जैसे Ne, Ar) के लिए विस्तारित नहीं होता है क्योंकि $\alpha$-स्विच तंत्र वहां प्रभावी नहीं है।
  • मात्रात्मक IP: He के लिए सटीक होने के बावजूद, भारी प्रणालियों के लिए आयनीकरण ऊर्जाओं में व्यवस्थित त्रुटियां हो सकती हैं क्योंकि यह हाइड्रोजन-संदर्भित अंशांकन (calibration) पर आधारित है।
  • प्रतिक्रिया गुण (Response Properties): स्विच सीमाओं की संवेदनशीलता के कारण उच्च डेरिवेटिव (जैसे ध्रुवीकरण क्षमताएं) इसके दायरे से बाहर हैं।
• कार्यान्वयन: लेख फलन को निश्चित घनत्वों पर गैर-स्व-संगत (non-self-consistently) रूप से प्रदर्शित करता है। यह नोट करता है कि गाऊसी आधारों (Gaussian bases) में पूर्ण स्व-संगत कार्यान्वयन लैपलेसियन पद ($v_{lapl}$) के कारण अस्थिर है, लेकिन तर्क देता है कि यह विधि वास्तविक-स्थान स्यूडोपोटेंशियल कार्यान्वयन के लिए स्वाभाविक रूप से उपयुक्त है जहाँ लैपलेसियन को स्थिरता से संभाला जा सकता है।

यह कार्य स्थापित करता है कि एक एकल गैर-अनुभवजन्य इनपुट—GP93 हाइड्रोजेनिक-सटीकता स्थिति—एक साथ सही एसिम्प्टोटिक टेल और सटीक हाइड्रोजेनिक स्तर उत्पन्न कर सकती है, जो नॉनलोकल सटीक विनिमय की आवश्यकता के बिना रिडबर्ग भौतिकी के लिए एक पथ प्रदान करती है।