GEM: An implementation of the ghost-Gutzwiller approximation for simulating interacting quantum systems

ठोस पदार्थों में इलेक्ट्रॉन आमतौर पर एक सुव्यवस्थित भीड़ की तरह व्यवहार करते हैं, जो बिजली या चुंबकत्व बनाने के लिए सामग्रियों के माध्यम से सुचारू रूप से बहते हैं। लेकिन कुछ विशेष परिस्थितियों में, वे जिद्दी और अराजक हो सकते हैं, एक-दूसरे के आगे बढ़ने से इनकार कर सकते हैं और पदार्थ की ऐसी अवस्था बना सकते हैं जहाँ उनकी व्यक्तिगत अंतःक्रियाएं पूरे तंत्र पर हावी हो जाती हैं। यह दृढ़ता से सहसंबद्ध इलेक्ट्रॉनों (strongly correlated electrons) का क्षेत्र है, जो भौतिकी का वह कोना है जो बताता है कि क्यों कुछ सामग्रियां कुचालक होती हैं जबकि अन्य अतिचालक (superconductors) होती हैं, और क्यों कुछ धातुएं विशिष्ट तापमानों पर अपना चुंबकत्व खो देती हैं। दशकों से, वैज्ञानिक इन व्यवहारों को सिम्युलेट करने के लिए 'डायनामिकल मीन-फील्ड थ्योरी' (dynamical mean-field theory) नामक एक शक्तिशाली लेकिन गणनात्मक रूप से महंगी तकनीक पर निर्भर रहे हैं। यह विधि प्रत्येक इलेक्ट्रॉन के साथ इस तरह व्यवहार करती है जैसे वह एक लगातार बदलते परिवेश के साथ अंतःक्रिया कर रहा हो, जिसके लिए कणों के जटिल, समय-निर्भर नृत्य को ट्रैक करने के लिए अत्यधिक कंप्यूटिंग शक्ति की आवश्यकता होती है। हालांकि यह सटीक है, लेकिन यह दृष्टिकोण अक्सर बड़े सिस्टम, कम तापमान, या जब शोधकर्ताओं को किसी पैटर्न को खोजने के लिए कई अलग-अलग स्थितियों को स्कैन करने की आवश्यकता होती है, तो व्यावहारिक रूप से बहुत धीमा हो जाता है।

'गटज़विलर सन्निकटन' (Gutzwiller approximation) के रूप में जानी जाने वाली एक हल्की और तेज़ विधि ने इलेक्ट्रॉनों के औसत व्यवहार पर ध्यान केंद्रित करके एक शॉर्टकट प्रदान किया है, जो उच्च-ऊर्जा उतार-चढ़ाव को अनदेखा करती है। हालाँकि, इस सरलीकरण की एक कीमत चुकानी पड़ती है: यह ऊर्जा स्पेक्ट्रम के उन महत्वपूर्ण विवरणों को छोड़ देता है जो वास्तव में सामग्रियों के व्यवहार को परिभाषित करते हैं। इस अंतर को पाटने के लिए एक नया सैद्धांतिक ढांचा, 'घोस्ट-गटज़विलर सन्निकटन' (ghost-Gutzwiller approximation) विकसित किया गया। यह "घोस्ट" ऑर्बिटल्स (ghost orbitals) पेश करता है—काल्पनिक सहायक अवस्थाएं जो भारी-भरकम विधि की पूरी लागत के बिना छूटे हुए उच्च-ऊर्जा विवरणों को पुनः प्राप्त करने के लिए एक बफर के रूप में कार्य करती हैं। अब तक, यह आशाजनक मध्य मार्ग केवल सिद्धांत और बिखरे हुए कोड स्निपेट्स में मौजूद था, जिसमें शोधकर्ताओं के लिए वास्तविक दुनिया के मॉडलों पर इसे परीक्षण करने के लिए एक एकीकृत, ओपन-सोर्स सॉफ्टवेयर पैकेज की कमी थी।

SciPost Physics Codebases में एक नए सबमिशन में, शोधकर्ताओं की एक टीम GEM पेश करती है, जो एक सॉफ्टवेयर पैकेज है जिसे घोस्ट-गटज़विलर सन्निकटन को सुलभ और व्यावहारिक बनाने के लिए डिज़ाइन किया गया है। टीम, जिसका नेतृत्व सामुएल गिउली और फ्लैटआयरन इंस्टीट्यूट तथा विभिन्न अंतरराष्ट्रीय संस्थानों के सहयोगियों ने किया है, ने एक ऐसा टूल बनाया है जो वैज्ञानिकों को एक ट्यून करने योग्य सटीकता स्तर के साथ अंतःक्रिया करने वाले क्वांटम सिस्टम को सिम्युलेट करने की अनुमति देता है। मूल विचार सरल लेकिन शक्तिशाली है: इन घोस्ट ऑर्बिटल्स की संख्या को समायोजित करके, एक उपयोगकर्ता तेज़, अनुमानित गटज़विलर विधि और धीमी, सटीक डायनामिकल मीन-फील्ड थ्योरी के बीच सहजता से स्लाइड कर सकता है। जब घोस्ट ऑर्बिटल्स की संख्या कम होती है, तो गणना तेज़ होती है और कुल ऊर्जा और चुंबकीय क्रम जैसी स्थिर गुणों का एक अच्छा अनुमान प्रदान करती है। जैसे-जैसे घोस्ट ऑर्बिटल्स की संख्या बढ़ती है, सिमुलेशन अधिक जटिल गतिकी को कैप्चर करता है, और अंततः बहुत अधिक महंगी विधि के परिणामों के अनुरूप हो जाता है।

यह सॉफ्टवेयर पायथन (Python) में लिखा गया है और मौजूदा वैज्ञानिक उपकरणों, विशेष रूप से TRIQS इकोसिस्टम के साथ सहजता से काम करने के लिए डिज़ाइन किया गया है, जो क्वांटम सिमुलेशन के लिए एक मानक ढांचा है। GEM समीकरणों को सेट करने और 'सेल्फ-कंसिस्टेंसी लूप' (self-consistency loop) को प्रबंधित करने का भारी काम संभालता है, जहाँ सिमुलेशन अपने स्वयं के अनुमानों को तब तक बार-बार अपडेट करता है जब तक कि परिणाम स्थिर न हो जाएं। इसमें एक अंतर्निहित सॉल्वर शामिल है जो 'एक्जैक्ट डायगोनलाइजेशन' (exact diagonalization) का उपयोग करता है, जो एक छोटे क्वांटम सिस्टम की सटीक ऊर्जा अवस्थाओं को खोजने की एक तकनीक है, लेकिन इसे भविष्य में अन्य सॉलवर्स को स्वीकार करने के लिए भी बनाया गया है। यह मॉड्यूलरिटी का अर्थ है कि शोधकर्ता अध्ययन कर रहे विशिष्ट प्रश्न के आधार पर विभिन्न संख्यात्मक विधियों को प्लग इन कर सकते हैं, चाहे उसमें शून्य-तापमान ग्राउंड स्टेट शामिल हो या परिमित-तापमान ऊष्मागतिकी।

GEM की क्षमताओं का प्रदर्शन करने के लिए, लेखकों ने कई बेंचमार्क सिमुलेशन चलाए जो भौतिक परिदृश्यों की एक विस्तृत श्रृंखला को कवर करते हैं। एक परीक्षण में, उन्होंने बेथे लैटिस (Bethe lattice) पर एक सिंगल-ऑर्बिटल सिस्टम का मॉडल बनाया, जो जटिल सामग्रियों को सरल बनाने के लिए उपयोग किया जाने वाला एक सैद्धांतिक ढांचा है। उन्होंने तापमान की एक सीमा में सिस्टम की कुल ऊर्जा और एंट्रॉपी की गणना की, जिससे पता चला कि सॉफ्टवेयर उच्च सटीकता के साथ अपेक्षित ऊष्मागतिक व्यवहार को पुनरुत्पादित कर सकता है। दूसरे उदाहरण में, उन्होंने एक स्क्वायर लैटिस का चुंबकीय चरण आरेख (magnetic phase diagram) तैयार किया, जो दो-आयामी सामग्रियों में चुंबकत्व को समझने के लिए एक क्लासिक मॉडल है। विभिन्न तापमानों और अंतःक्रिया की ताकत पर सिस्टम को सिम्युलेट करके, उन्होंने सफलतापूर्वक उस महत्वपूर्ण तापमान की पहचान की जहाँ सामग्री चुंबकीय से गैर-चुंबकीय अवस्था में परिवर्तित होती है, जो अधिक मांग वाली विधियों के परिणामों से मेल खाता है।

तीसरा और शायद सबसे महत्वपूर्ण प्रदर्शन एक मल्टी-ऑर्बिटल सिस्टम से संबंधित था, जहाँ इलेक्ट्रॉन परमाणु के चारों ओर कई अलग-अलग ऑर्बिटल्स में स्थित हो सकते हैं। यह सेटअप 'हुंड्स फिजिक्स' (Hund's physics) को समझने के लिए महत्वपूर्ण है, एक ऐसी घटना जहाँ परमाणु के भीतर इलेक्ट्रॉन स्पिन का संरेखण सामग्री के गुणों को गहराई से प्रभावित करता है। टीम ने एक तीन-ऑर्बिटल सिस्टम का अनुकरण किया और ट्रैक किया कि कैसे "क्वासिपार्टिकल वेट" (quasiparticle weight)—जो एक माप है कि इलेक्ट्रॉन सामग्री के माध्यम से कितनी स्वतंत्र रूप से चल सकते हैं—अंतःक्रिया की ताकत बढ़ाने के साथ बदलता है। उन्होंने पाया कि मध्यम संख्या में घोस्ट ऑर्बिटल्स का उपयोग करने से सॉफ्टवेयर हुंड्स कपलिंग के कारण होने वाले जटिल रेनोर्मलाइजेशन प्रभावों को पकड़ने में सक्षम हुआ, जिससे ऐसे परिणाम प्राप्त हुए जो गोल्ड-स्टैंडर्ड डायनामिकल मीन-फील्ड थ्योरी के बहुत करीब थे, लेकिन बहुत कम कंप्यूटेशनल लागत पर।

पेपर इस बात पर जोर देता है कि GEM केवल एक सैद्धांतिक अभ्यास नहीं है बल्कि जटिल सामग्रियों के चरण आरेख (phase diagrams) की खोज के लिए एक व्यावहारिक उपकरण है। लेखक दिखाते हैं कि घोस्ट ऑर्बिटल्स की एक छोटी संख्या के साथ भी, विधि गति और सटीकता के बीच एक विश्वसनीय संतुलन प्रदान करती है, जिससे उन प्रणालियों का अध्ययन करना संभव हो जाता है जो पहले सिम्युलेट करने के लिए बहुत बड़ी या बहुत जटिल थीं। एक ओपन-सोर्स कार्यान्वयन प्रदान करके, टीम को उम्मीद है कि अन्य शोधकर्ता इन परिणामों को दोहरा सकेंगे, नए मॉडल का परीक्षण कर सकेंगे, और संभावित रूप से इसे सामग्रियों की खोज के लिए बड़े वर्कफ़्लो में एकीकृत कर सकेंगे। यह सॉफ्टवेयर विस्तार योग्य है, जिसे अधिक उन्नत सॉलवर्स के साथ इंटरफेस करने और अंततः प्रथम-सिद्धांत (first-principles) इलेक्ट्रॉनिक संरचना कोड के साथ जुड़ने के लिए डिज़ाइन किया गया है, जिससे परमाणु संरचना से वास्तविक सामग्रियों का सीधा सिमुलेशन संभव होगा।

अंततः, GEM उन्नत क्वांटम सिमुलेशन के लोकतंत्रीकरण की दिशा में एक कदम है। यह कंप्यूटेशनल लागत और भौतिक निष्ठा (physical fidelity) के बीच के समझौते को नेविगेट करने का एक तरीका प्रदान करता है, जिससे शोधकर्ताओं को दृढ़ता से सहसंबद्ध इलेक्ट्रॉनों के व्यवहार की जांच करने के लिए एक लचीला उपकरण मिलता है। चाहे इसका उपयोग एक नए चुंबकीय चरण की स्थिरता की जाँच करने के लिए किया जाए या एक जटिल मिश्र धातु की ऊष्मागतिकी को समझने के लिए, यह सॉफ्टवेयर उन परिणामों के लिए एक स्पष्ट मार्ग प्रदान करता है जो कभी कई लोगों की पहुँच से बाहर थे। घोस्ट-गटज़विलर सन्निकटन को एक मानक, पुन: प्रयोज्य उपकरण बनाकर, लेखकों ने एक व्यापक समुदाय के लिए क्वांटम सामग्रियों की समृद्ध और अक्सर प्रति-सहज (counterintuitive) दुनिया का पता लगाने के द्वार खोल दिए हैं।