Towards more accurate natural orbital functional approximations: including 4-index cumulant contributions

कल्पना कीजिए कि आप एक अणु का एक आदर्श डिजिटल जुड़वां (digital twin) बनाने की कोशिश कर रहे हैं, जो इलेक्ट्रॉनों से बनी एक छोटी लेगो संरचना की तरह है। दशकों से, वैज्ञानिक 'रिड्यूस्ड डेंसिटी मैट्रिक्स फंक्शनल थ्योरी' (RDMFT) नामक एक उपकरण का उपयोग कर रहे हैं। RDMFT को एक सुपर-स्मार्ट रेसिपी बुक की तरह समझें जो आपको बताती है कि इलेक्ट्रॉन कैसे व्यवहार करते हैं, बिना हर एक इलेक्ट्रॉन को व्यक्तिगत रूप से ट्रैक किए—जो कि हवा के झोंकों के बीच समुद्र तट पर रेत के हर एक कण को गिनने की कोशिश करने जैसा होगा। लक्ष्य यह अनुमान लगाना है कि अणु कैसे आपस में जुड़े रहते हैं और, सबसे महत्वपूर्ण बात यह है कि जब वे टूटते हैं तो क्या होता है।

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

यह शोध पत्र इस विशिष्ट समस्या को हल करता है और PNOF5 नामक एक लोकप्रिय रेसिपी को ठीक करता है। लेखकों ने महसूस किया कि हालांकि PNOF5 कई चीजों में अच्छा है, लेकिन यह एक महत्वपूर्ण विवरण को छोड़ देता है कि इलेक्ट्रॉन जोड़ों में कैसे परस्पर क्रिया करते हैं जब एक अणु टूटने की ओर खिंच रहा होता है। यह ऐसा है जैसे रेसिपी मान लेती है कि इलेक्ट्रॉन केवल अपने सीधे साथियों के साथ ही संवाद करते हैं, और इस बात को नजरअंदाज कर देती है कि वे पास के अन्य जोड़ों से भी फुसफुसाहट (whisper) करते हैं। लेखकों ने एक नई विधि पेश की ताकि यह सुनिश्चित किया जा सके कि रेसिपी "टूटने के नियमों" का सम्मान करे: यह सुनिश्चित करना कि जब एक अणु विभाजित होता है, तो टुकड़े वास्तव में अलग हो जाएं और उनके स्पिन सही ढंग से स्थिर हो जाएं। इन सख्त नियमों को जोड़कर और फिर परिणाम को पॉलिश करके ताकि यह भौतिकी के नियमों का पालन करे, उन्होंने एक नया, बेहतर संस्करण बनाया। उनके सिमुलेशन दिखाते हैं कि नाइट्रोजन (N2) और ऑक्सीजन (O2) जैसे अणुओं के लिए, यह नया तरीका हमारे पास मौजूद सबसे उन्नत, महंगे सिमुलेशन के बिल्कुल समान व्यवहार करता है, और अंततः बॉन्ड के "स्नैप" (टूटना) को सही ढंग से पकड़ लेता है।

टूटे हुए रबर बैंड की कहानी

आइए इलेक्ट्रॉनों की दुनिया में उतरें, जहाँ चीजें थोड़ी अजीब हो जाती हैं। कल्पना कीजिए कि एक अणु एक डांस फ्लोर है जहाँ इलेक्ट्रॉन डांसर हैं। एक स्थिर अणु में, वे कसकर, समन्वित जोड़ों में नृत्य करते हैं। लेकिन क्या होता है जब आप अणु को अलग करने के लिए खींचते हैं, उसे तब तक खींचते हैं जब तक कि वह टूट न जाए? इसे "बॉन्ड डिसोसिएशन" (bond dissociation) कहा जाता है। वास्तविक दुनिया में, बॉन्ड टूटने के बाद, दोनों हिस्से स्वतंत्र, खुश परमाणुओं के रूप में दूर चले जाने चाहिए। वे अब हाथ नहीं पकड़े होने चाहिए, और उनके आंतरिक "स्पिन" (सोचिए कि वे छोटे घूमते हुए टॉप की तरह हैं) एक शांत, अनुमानित लय में स्थिर हो जाने चाहिए।

समस्या यह है कि हमारे वर्तमान कंप्यूटर मॉडल, जो इस व्यवहार की भविष्यवाणी करने की कोशिश करते हैं, अक्सर टूटने के दृश्य को गलत कर देते हैं। वे एक खराब फिल्म की तरह व्यवहार करते हैं जहाँ पर्दे गिरने के बाद भी अभिनेता मंच छोड़ने से इनकार कर देते हैं। अणु के दो हिस्से अजीब तरह से जुड़े रहते हैं, और उनके स्पिन भ्रमित हो जाते हैं, जिससे यह गलत उत्तर मिलता है कि बॉन्ड तोड़ने में कितनी ऊर्जा लगी। यह उन वैज्ञानिकों के लिए एक बड़ा सिरदर्द है जो नए सामग्रियां डिजाइन करने या रासायनिक प्रतिक्रियाओं को समझने की कोशिश कर रहे हैं।

यहाँ वैलेरी चुइको, पॉल डब्ल्यू. एयर्स और एडुआर्ड माटीटो की टीम आती है। उन्होंने PNOF5 नामक एक विशिष्ट रेसिपी को ठीक करने का निर्णय लिया। PNOF5 को इलेक्ट्रॉनों के जोड़े बनाने के निर्देशों के एक सेट के रूप में समझें। यह एक चतुर रेसिपी है, लेकिन इसमें एक अंधा बिंदु (blind spot) है: यह मानती है कि इलेक्ट्रॉन केवल अपने डांस पार्टनर की परवाह करते हैं और बाकी सबको नजरअंदाज करते हैं। यह तब ठीक काम करता है जब अणु पूरा होता है, लेकिन जब आप इसे अलग करने के लिए खींचना शुरू करते हैं, तो यह धारणा ढह जाती है। रेसिपी यह गणना करने में विफल रहती है कि अलग-अलग इलेक्ट्रॉन जोड़ों के बीच होने वाली सूक्ष्म, लंबी दूरी की फुसफुसाहट क्या है, जो टूटने को वास्तविक दिखाने के लिए आवश्यक है।

समाधान: गायब नियमों को जोड़ना

लेखकों का समाधान यह था कि रेसिपी में नियमों का एक नया सेट जोड़ा जाए। उन्होंने महसूस किया कि जब एक अणु टूटता है, तो टुकड़ों के "स्पिन" के व्यवहार के बारे में कुछ विशिष्ट, गैर-परक्राम्य (non-negotiable) नियम होते हैं। यदि आपके पास एक अणु दो टुकड़ों में टूट रहा है, तो प्रत्येक टुकड़े का कुल स्पिन एक एकल, अलग परमाणु के समान होना चाहिए। यह ऐसा ही है जैसे यदि आप एक खिलौना कार को दो भागों में तोड़ते हैं, तो प्रत्येक भाग के पहिये अभी भी सही गति से घूमने चाहिए, न कि बेतरतीब ढंग से या अटके हुए।

PNOF5 को ठीक करने के लिए, उन्होंने केवल अनुमान नहीं लगाया; उन्होंने गणित का उपयोग करके इलेक्ट्रॉन मानचित्र का "निकटतम संभव" संस्करण खोजा जो इन नए स्पिन नियमों का पालन करता है। कल्पना कीजिए कि आपके पास एक टूटे हुए फूलदान की धुंधली फोटो है, और आप जानते हैं कि टुकड़े वास्तव में कैसे दिखने चाहिए। उन्होंने एक गणितीय एल्गोरिदम का उपयोग करके फोटो को स्पष्ट (sharpen) किया, जिससे यह टूटने के ज्ञात नियमों से मेल खा सके। लेकिन एक पेच था: कभी-कभी, जब आप किसी फोटो को एक निश्चित तरीके से दिखने के लिए मजबूर करते हैं, तो आप अनजाने में भौतिकी के नियमों को तोड़ देते हैं (जैसे कि फूलदान के टुकड़ों को हवा में तैरते हुए दिखाना)। इसलिए, उन्होंने एक दूसरा कदम जोड़ा: एक "शुद्धिकरण" (purification) प्रक्रिया। यह चरण गणित की जांच करता है कि क्या नया, स्पष्ट फोटो अभी भी ब्रह्मांड के मौलिक नियमों का पालन करता है (विशेष रूप से, जिसे "N-representability" कहा जाता है, जो बस एक फैंसी तरीका है यह कहने का कि "क्या यह इलेक्ट्रॉन व्यवस्था वास्तव में अस्तित्व में है?")।

उन्होंने इस दो-चरणीय प्रक्रिया को—पहले स्पिन नियमों को ठीक करें, फिर भौतिकी को शुद्ध करें—बार-बार चलाया जब तक कि संख्याएं बदलना बंद नहीं हो गईं। यह एक गिटार के तार को ट्यून करने जैसा है: आप इसे कसते हैं, पिच की जांच करते हैं, तनाव को समायोजित करते हैं, और तब तक दोहराते हैं जब तक कि यह सही सुर न निकाल दे।

परिणाम: स्नैप को सही पकड़ना

जब उन्होंने नाइट्रोजन ($N_2$), ऑक्सीजन ($O_2$) और नाइट्रिक ऑक्साइड ($NO^+$) जैसे अणुओं पर अपने सुधरे हुए तरीके का परीक्षण किया, तो परिणाम आश्चर्यजनक थे। "ब्रेकअप ज़ोन" में—जब परमाणुओं को दूर खींचा गया था, जैसे कि 3.5 Å या उससे अधिक की दूरी पर—उनके नए तरीके ने ऊर्जा मानों को कंप्यूटर केमिस्ट्री के स्वर्ण मानक (gold standard), जिसे CASSCF कहा जाता है, के लगभग पूरी तरह से मेल खाया।

उदाहरण के लिए, नाइट्रोजन ($N_2$) के सिमुलेशन में, पुराने PNOF5 रेसिपी ने -107.315 Hartrees (ऊर्जा की एक इकाई) की ऊर्जा दी, जबकि नए संशोधित संस्करण ने -107.438 Hartrees दी। "गोल्ड स्टैंडर्ड" CASSCF गणना ने भी -107.438 Hartrees ही दिया। ऑक्सीजन और नाइट्रिक ऑक्साइड के लिए भी यही हुआ। नए तरीके ने न केवल ऊर्जा को सही किया; इसने "डिलोकलाइजेशन इंडेक्स" (delocalization index) को भी ठीक किया, जो यह मापता है कि इलेक्ट्रॉन दो टुकड़ों के बीच कितने साझा किए जाते हैं। पुराने रेसिपी में, टुकड़े अभी भी इलेक्ट्रॉन साझा कर रहे थे जब उन्हें अलग हो जाना चाहिए था। नए रेसिपी में, वह साझाकरण गायब हो गया, ठीक वैसे ही जैसे वास्तविक दुनिया में होना चाहिए।

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

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

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

लेखक स्वीकार करते हैं कि अभी भी काम करना बाकी है। उनका तरीका कम्प्यूटेशनल रूप से भारी (computationally heavy) है, जिसका अर्थ है कि इसमें बहुत अधिक कंप्यूटर शक्ति की आवश्यकता होती है, विशेष रूप से बड़े अणुओं के लिए। वे यह भी नोट करते हैं कि कुछ बहुत जटिल स्थितियों के लिए, उनका वर्तमान दृष्टिकोण हर एक विवरण को कैप्चर करने में विफल हो सकता है। लेकिन उन्होंने जो रास्ता साफ किया है वह स्पष्ट है: यदि हम अणुओं के टूटने की भविष्यवाणी करना चाहते हैं, तो हमें यह सुनिश्चित करना होगा कि हमारे मॉडल टुकड़ों के भौतिकी का सम्मान करें। और इस नए तरीके के साथ, उन्होंने दिखाया है कि यह संभव है कि बॉन्ड के "स्नैप" को बिल्कुल सही तरीके से पकड़ा जाए।