Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

कल्पना कीजिए कि आप यह अनुमान लगाने की कोशिश कर रहे हैं कि एक विशाल संगीत कार्यक्रम (कॉन्सर्ट) में लोगों की एक विशाल भीड़ कैसे व्यवहार करेगी। यदि आप केवल एक व्यक्ति को देखते हैं, तो आप 'मॉस पिट' (mosh pit) की अराजकता को देखने से चूक जाते हैं। यदि आप 1,00,000 प्रशंसकों के स्टेडियम में हर एक व्यक्ति को एक साथ ट्रैक करने की कोशिश करते हैं, तो आपका मस्तिष्क (या आपका कंप्यूटर) फट जाएगा। यह उन वैज्ञानिकों का दैनिक संघर्ष है जो "क्वांटम सामग्रियों" (quantum materials) का अध्ययन करते हैं। वे इलेक्ट्रॉन्स—जो कि हर चीज़ को बनाने वाले नन्हे, फुर्तीले कण हैं—को ठोस क्रिस्टल में एक साथ नाचते हुए समझना चाहते हैं। समस्या यह है कि इलेक्ट्रॉन अविश्वसनीय रूप से सामाजिक होते हैं; वे केवल अकेले नहीं बैठते, बल्कि वे जटिल, "सहसंबंधित" (correlated) तरीकों से एक-दूसरे को प्रभावित करते हैं। सही उत्तर पाने के लिए, आप केवल कुछ इलेक्ट्रॉन्स को नहीं देख सकते; आपको पूरे अनंत क्रिस्टल का अनुकरण (simulate) करना होगा। लेकिन एक ऐसे क्रिस्टल के लिए यह गणित करना जो अनंत तक फैला हुआ है, समुद्र में ज्वार आने के दौरान रेत के हर एक कण को गिनने की कोशिश करने जैसा है। पुराने तरीकों के साथ यह गणनात्मक रूप से असंभव है, जो जैसे-जैसे सिस्टम बड़ा होता जाता है, धीमा होता जाता है और अक्सर वास्तविक उत्तर मिलने से पहले ही रुक जाता है।

यहीं पर जेनी (Junjie Yang) और गार्नेट चैन (Garnet Chan) के नेतृत्व वाली शोधकर्ताओं की एक नई टीम एक चतुर शॉर्टकट के साथ आती है। उन्होंने इन विशाल गणनाओं को करने का एक नया तरीका विकसित किया है जो केवल थोड़ा तेज़ नहीं करता है—बल्कि यह खेल को पूरी तरह बदल देता है। क्रिस्टल के मानचित्रण के लिए आवश्यक बिंदुओं की भारी संख्या से घबराने के बजाय, उनकी विधि "रैखिक रूप से स्केल" (scales linearly) करती है। इसे इस तरह सोचें: यदि पुराने तरीके को 10 रेत के कण गिनने में एक घंटा लगा, 20 के लिए दो घंटे, और 1,000 के लिए 100 घंटे लगे, तो उनकी नई विधि 10 के लिए एक घंटा, 20 के लिए दो घंटे, और 1,000 के लिए केवल 100 घंटे लेगी। यह एक सीधी रेखा है, न कि एक ऐसा वक्र (curve) जो चार्ट से बाहर निकल जाए। "इंटरपोलेटिव सेपरेबल डेंसिटी फिटिंग" (Interpolative Separable Density Fitting - ISDF) नामक तकनीक का उपयोग करके, जिसे क्रिस्टल की समरूपता (symmetry) के साथ जोड़ा गया है, वे बिना किसी कठिनाई के क्रिस्टल की संरचना में 1,000 अलग-अलग बिंदुओं का अनुकरण करने में सफल रहे। उन्होंने कठोर डायमंड से लेकर जटिल चुंबकीय धातुओं तक, सब पर इसका परीक्षण किया, जिससे यह सिद्ध हुआ कि हम "ऊष्मप्रवैगिकी सीमा" (thermodynamic limit)—पदार्थ की "अनंत" अवस्था—में झाँकने में सक्षम हैं, जिसके लिए हमें शहर के आकार के सुपरकंप्यूटर की आवश्यकता नहीं है।

"मैजिक" शॉर्टकट का जादू

यह समझने के लिए कि इन वैज्ञानिकों ने वास्तव में क्या किया, कल्पना कीजिए कि आप सूप के एक बड़े बर्तन का स्वाद बताने की कोशिश कर रहे हैं। करने का पुराना तरीका (जिसे "गौसियन डेंसिटी फिटिंग" या GDF कहा जाता है) सूप के औसत स्वाद को जानने के लिए बर्तन की हर एक बूंद को चखने जैसा है। जैसे-जैसे बर्तन बड़ा होता जाता है, आपको अधिक और अधिक बूंदों को चखना पड़ता है, और इसमें लगने वाला समय विस्फोटक रूप से बढ़ता जाता है। यदि आप बर्तन का आकार दोगुना करते हैं, तो लगने वाला समय चार गुना या उससे भी अधिक हो सकता है।

नई विधि, जिसे लेखक FFTISDF कहते हैं (एक लंबा नाम जिसका अर्थ है Fast Fourier Transform Interpolative Separable Density Fitting), एक मास्टर शेफ की तरह है जो जानता है कि पूरे बर्तन का स्वाद ठीक से जानने के लिए आपको केवल कुछ विशिष्ट, रणनीतिक चम्मचों को चखने की आवश्यकता है। वे सबसे महत्वपूर्ण "इंटरपोलेशन पॉइंट्स" (interpolation points)—यानी मुख्य चम्मचों—को चुनने के लिए एक गणितीय ट्रिक का उपयोग करते हैं जो पूरे सूप का सटीक प्रतिनिधित्व करते हैं। केवल इन प्रमुख बिंदुओं पर ध्यान केंद्रित करके और क्रिस्टल के दोहराव वाले पैटर्न (जैसे फर्श पर दोहराए जाने वाले टाइल्स) का उपयोग करके, वे अविश्वसनीय गति के साथ पूरे अनंत बर्तन के स्वाद का पुनर्निर्माण कर सकते हैं।

क्वांटम भौतिकी की दुनिया में, यह "सूप" इलेक्ट्रॉन्स के बीच की परस्पर क्रिया है। लेखकों ने इस विधि का उपयोग चार बहुत अलग सामग्रियों में इलेक्ट्रॉन्स की ऊर्जा की गणना करने के लिए किया:

1. डायमंड (हीरा): एक कठोर, वाइड-गैप सेमीकंडक्टर (एक बहुत ही मजबूत इंसुलेटर की तरह)।
2. कार्बन डाइऑक्साइड (CO2): एक आणविक क्रिस्टल जहाँ अणु कमजोर बलों द्वारा जुड़े होते हैं।
3. निकल मोनोऑक्साइड (NiO): एक जटिल चुंबकीय सामग्री जहाँ इलेक्ट्रॉन "स्ट्रॉन्गली कोरिलेटेड" हैं, जिसका अर्थ है कि वे इस बारे में बहुत चूजी हैं कि वे किसके साथ समय बिताएंगे।
4. CaCuO2: एक क्यूप्रेट सुपरकंडक्टर, एक ऐसी सामग्री जो बिना किसी प्रतिरोध के बिजली का संचालन कर सकती है, जो ऊर्जा अनुसंधान का "होली ग्रिल" है।

परिणाम: "शायद" से "लगभग पहुँच गए" तक

टीम ने केवल एक तेज़ कार नहीं बनाई; उन्होंने एक ऐसा हाईवे बनाया जहाँ पुराने कार चल भी नहीं सकते थे। उन्होंने 1,000 k-points का उपयोग करके इन सामग्रियों पर सिमुलेशन चलाए। क्रिस्टल भौतिकी की भाषा में, एक "k-point" इलेक्ट्रॉन के संवेग (momentum) के स्नैपशॉट की तरह है। एक अनंत क्रिस्टल की पूर्ण तस्वीर पाने के लिए, आपको अनंत संख्या में स्नैपशॉट की आवश्यकता होती है। चूंकि आप अनंत नहीं कर सकते, वैज्ञानिक आमतौर पर कुछ (जैसे 8 या 16) करते हैं और अनुमान लगाते हैं कि अनंत संस्करण कैसा दिखेगा।

उनकी नई विधि के साथ, उन्हें उतना अनुमान लगाने की आवश्यकता नहीं पड़ी। वे वास्तव में 1,000 स्नैपशॉट के साथ सिमुलेशन चला सके।

यहाँ उन्हें क्या मिला:

• गति और पैमाना (Speed and Scale): डायमंड क्रिस्टल के लिए, उन्होंने दिखाया कि उनकी विधि k-points की संख्या के साथ पूरी तरह से रैखिक रूप से स्केल करती है। जबकि अन्य विधियाँ (जैसे GDF) अधिक पॉइंट्स जोड़ने पर धीमी और महंगी होने लगीं, FFTISDF एक स्थिर, तेज़ गति से चलती रही। उन्होंने गणना की कि डायमंड के एक विशाल 1,000 k-point सिमुलेशन के लिए, उनकी विधि को केवल 2.11 GB मेमोरी की आवश्यकता थी, जबकि पुराने तरीके को चौंका देने वाले 881.09 GB की आवश्यकता होती। यह एक लैपटॉप और हार्ड ड्राइवों के एक गोदाम के बीच के अंतर जैसा है।
• सटीकता (Accuracy): उन्होंने अपने "शॉर्टकट" परिणामों की तुलना "हर बूंद चखने वाले" परिणामों से करके सटीकता का परीक्षण किया। उन्होंने पाया कि एक विशिष्ट पैरामीटर (जिसे $c_{IP}$ कहा जाता है, जो यह नियंत्रित करता है कि वे कितने "चम्मच" चखते हैं) को समायोजित करके, वे त्रुटि को प्रति परमाणु $10^{-4}$ परमाणु इकाइयों से कम कर सकते हैं। यह अविश्वसनीय रूप से सटीक है।
• ऊष्मप्रवैगिकी सीमा (The Thermodynamic Limit): अंतिम लक्ष्य पदार्थ की "ऊष्मप्रवैगिकी सीमा" (अनंत आकार) में ऊर्जा खोजना था। अधिक और अधिक k-points (1,000 तक) के साथ अपने सिमुलेशन चलाकर और फिर गणितीय रूप से परिणामों को अनंत तक खींचकर, वे अंततः CCSD(T) (एक बहुत ही उच्च-स्तरीय क्वांटम केमिस्ट्री विधि) जैसी जटिल सामग्रियों के लिए विश्वसनीय उत्तर प्राप्त करने में सक्षम थे।
  • डायमंड के लिए, उन्होंने सफलतापूर्वक ऊर्जा को ऊष्मप्रवैगिकी सीमा तक पहुँचाया (extrapolate किया), जिससे यह दिखाया कि उनकी विधि उस "एनर्जी-इन्क्रीजिंग रिजीम" को संभाल सकती है जो आमतौर पर अन्य कंप्यूटरों को रोक देता है।
  • CO2 के लिए, उन्होंने एक चतुर एक्सट्रपलेशन ट्रिक का उपयोग किया। उन्होंने क्रिस्टल के छोटे हिस्सों के लिए ऊर्जा की गणना की और फिर पूरे के लिए ऊर्जा की भविष्यवाणी करने के लिए एक गणितीय रेखा का उपयोग किया। उन्होंने पाया कि उनकी अनुमानित ऊर्जा वास्तविक वैश्विक मान के लगभग समान थी, जिसमें प्रति परमाणु त्रुटि केवल $6 \times 10^{-5}$ परमाणु इकाइयाँ थी।
  • चुंबकीय सामग्रियों NiO और CaCuO2 के लिए, उन्होंने गणना की कि इलेक्ट्रॉन चुंबकत्व बनाने के लिए कैसे परस्पर क्रिया करते हैं। उन्होंने पाया कि NiO के लिए, चुंबकीय संपर्क शक्ति ($J$) ऊष्मप्रवैगिकी सीमा में -14.05 meV थी। यह वास्तविक दुनिया के प्रयोगों (जो लगभग -17 से -19 meV के आसपास हैं) के बहुत करीब है, जबकि पुरानी विधियों ने केवल -7.59 meV का मान दिया था। CaCuO2 के लिए, उन्होंने -165.44 meV की चुंबकीय शक्ति पाई, जो प्रयोगात्मक डेटा के साथ बहुत अच्छी तरह मेल खाती है।

यह क्यों मायने रखता है

यह पेपर यह दावा नहीं करता है कि इसने भौतिकी के हर रहस्य को सुलझा लिया है। लेखक सावधानीपूर्वक नोट करते हैं कि उनकी विधि को सटीक होने के लिए अभी भी बहुत सारे "इंटरपोलेशन पॉइंट्स" (परमाणुओं की संख्या का लगभग 10 गुना) की आवश्यकता होती है, जो बहुत बड़े क्रिस्टल के लिए गणनात्मक रूप से भारी हो सकता है। वे यह भी उल्लेख करते हैं कि "रियल-स्पेस ग्रिड" (डिजिटल कैनवास जिस पर वे पेंट करते हैं) भारी धातुओं वाली सामग्रियों के लिए अभी भी काफी बड़ा है।

हालाँकि, मुख्य निष्कर्ष एक बड़ी छलांग है। उन्होंने सिद्ध किया है कि अनंत आवधिक प्रणालियों (infinite periodic systems) पर रैखिक स्केलिंग के साथ उच्च-स्तरीय, "एब इनीशियो" (प्रथम सिद्धांतों से) क्वांटम गणनाएं करना संभव है। इसका मतलब है कि जैसे-जैसे हम बड़ी और अधिक जटिल सामग्रियां बनाएंगे, हम उस दीवार से नहीं टकराएंगे जहाँ गणित असंभव हो जाता है। अब हम हजारों बिंदुओं के साथ क्रिस्टल में इलेक्ट्रॉन्स के व्यवहार का अनुकरण कर सकते हैं, जिससे हमें यह समझने के लिए अधिक स्पष्ट और सटीक चित्र मिलता है कि ये सामग्रियां वास्तव में कैसे काम करती हैं।

संक्षेप में, उन्होंने एक ऐसे गणितीय प्रश्न को जो "हल करना बहुत कठिन था" उसे "बिल्कुल सही" वाले प्रश्न में बदल दिया। उन्होंने केवल गिनने का तेज़ तरीका नहीं खोजा; उन्होंने केवल कुछ लहरों को देखकर पूरे समुद्र को देखने का तरीका खोजा, और उन्होंने इसे इतनी सटीकता के साथ किया कि पानी वास्तविक महसूस होता है।