Enabling Ab-Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions

तकनीकी सारांश: बायस (Bias) के तहत एब-इनिटियो मॉलिक्यूलर डायनेमिक्स को CP2K+SMEAGOL इंटरफेस के माध्यम से सक्षम करना

समस्या विवरण
यद्यपि बाहरी बायस पोटेंशियल के तहत क्वांटम ट्रांसपोर्ट को मॉडल करने के लिए डेंसिटी फंक्शनल थ्योरी (DFT) को नॉन-इक्विलिब्रियम ग्रीन'स फंक्शन्स (NEGF) के साथ जोड़ना एक मानक दृष्टिकोण है, पारंपरिक पद्धतियां स्थिर प्रणालियों तक ही सीमित हैं। मौजूदा कार्यान्वयन आमतौर पर एक कैनोनिकल फ्रेमवर्क के भीतर कार्य करते हैं, जो निश्चित परमाणु स्थितियों पर ट्रांसपोर्ट गुणों की गणना करते हैं। यह सैद्धांतिक समझ और प्रयोगों के मैक्रोस्कोपिक परिणामों के बीच एक महत्वपूर्ण अंतर पैदा करता है, विशेष रूप से इलेक्ट्रोकेमिकल (EC) प्रणालियों के लिए जहाँ रूपांतरण स्वाभाविक रूप से ग्रैंड कैनोनिकल होते हैं। EC रूपांतरणों के प्राकृतिक वातावरण का सटीक वर्णन करने के लिए, सिमुलेशन को बायस के तहत मॉलिक्यूलर डायनेमिक्स (MD) और इलेक्ट्रॉनों के एक स्पष्ट ओपन-बाउंड्री विवरण की आवश्यकता होती है जो कम्प्यूटेशनल सेल में प्रवेश कर सकें और बाहर निकल सकें। इसके अलावा, मौजूदा DFT-NEGF कार्यान्वयन अक्सर ज्यामिति अनुकूलन (geometry optimization) या डायनेमिक्स के लिए आवश्यक करंट-प्रेरित बलों (current-induced forces) की गणना करने में अक्षम होते हैं, या वे विशिष्ट DFT कोड (जैसे SIESTA) तक सीमित होते हैं जिनमें कंडेंस्ड फेज सिमुलेशन के लिए अनुकूलित अत्याधुनिक पैकेज में मौजूद विशेषताएं नहीं होती हैं।

कार्यप्रणाली
लेखक लोकप्रिय मिक्स्ड गॉसियन/प्लेन-वेव इलेक्ट्रॉनिक स्ट्रक्चर पैकेज CP2K और NEGF कोड SMEAGOL के बीच एक नया इंटरफेस प्रस्तुत करते हैं। इस इंटरफेस को CP2K+SMEAGOL कहा जाता है, जिसे एक स्टैंडअलोन लाइब्रेरी के रूप में डिज़ाइन किया गया है जो किसी भी DFT सॉफ्टवेयर के साथ लिंक किया जा सकता है, और यह मूल SMEAGOL कोड का यथासंभव पुन: उपयोग करता है।

प्रमुख तकनीकी कार्यान्वयन शामिल हैं:

• करंट-प्रेरित बल (Current-Induced Forces): इंटरफेस बायस के तहत आयनों पर कार्य करने वाले नॉन-कंजर्वेटिव बलों की गणना को लागू करता है। ये बल परमाणु स्थितियों के सापेक्ष पूर्ण हैमिल्टनियन के ग्रेडिएंट्स का उपयोग करते हुए, परमाणु संवेग (atomic momenta) के समय व्युत्पन्न (time derivative) से प्राप्त होते हैं। महत्वपूर्ण रूप से, कार्यान्वयन में डेंसिटी मैट्रिक्स ($\Omega_{ij}$) के व्युत्पन्न से उत्पन्न होने वाला बल पद शामिल है, जो "इलेक्ट्रॉन विंड" घटक (आयनों से इलेक्ट्रॉनों का संवेग स्थानांतरण) की व्याख्या करता है। इस पद को अन्य कार्यान्वयनों में अक्सर उपेक्षित किया जाता है, लेकिन इसे गुणात्मक शुद्धता के लिए आवश्यक दिखाया गया है।
• डेंसिटी मैट्रिक्स निर्माण: बायस विंडो के भीतर बाउंड स्टेट्स को सही ढंग से भरने के लिए, इंटरफेस डेंसिटी मैट्रिक्स गणना को अपडेट करता है। सिंगल कॉन्टूर इवैल्यूएशन या लेफ्ट/राइट योगदान के सरल औसत के बजाय, यह एक वेटेड डबल-कॉन्टूर स्कीम (Brandbyge आदि के अनुसार) का उपयोग करता है जो नॉन-इक्विलिब्रियम योगदान द्वारा निर्धारित भार के साथ दोनों लीड्स के इंटीग्रल्स को जोड़ता है।
• सिमुलेशन सेटअप: यह विधि बायस के तहत बोर्न-ओपनहाइमर मॉलिक्यूलर डायमिक्स की अनुमति देती है। सिस्टम को सेमी-इनफिनिट इलेक्ट्रोड्स से जुड़े एक केंद्रीय "एक्सटेंडेड मॉलिक्यूल" क्षेत्र के साथ मॉडल किया गया है। लीड्स के सेल्फ-एनर्जी को स्टैटिक और केंद्रीय क्षेत्र के चार्ज डेंसिटी से स्वतंत्र माना जाता है, जो एक ऐसी धारणा है जो तब मान्य होती है जब केंद्रीय क्षेत्र चार्ज परिवर्तनों को स्क्रीन करने के लिए पर्याप्त बड़ा हो।

प्रमुख परिणाम और सत्यापन
लेखकों ने कई प्रणालियों में SIESTA+SMEAGOL के साथ किए गए संदर्भ गणनाओं के विरुद्ध CP2K+SMEAGOL इंटरफेस को मान्य किया:

1. जीरो-बायस फोर्सेस (अनंत Au वायर): इंटरफेस जीरो बायस पर परमाणु बलों को सफलतापूर्वक पुनरुत्पादित करता है, जो CP2K-ओनली गणनाओं से मेल खाता है। यह यह भी प्रदर्शित करता है कि एनर्जी डेंसिटी मैट्रिक्स ($\Omega$) पद की उपेक्षा करने से गुणात्मक रूप से गलत बल प्राप्त होते हैं।
2. पैरेलल-प्लेट कैपेसिटर: इंटरफेस कैपेसिटर के माध्यम से अपेक्षित पोटेंशियल ड्रॉप को पुनरुत्पादित करता है, जो लीड्स में एक स्थिर पोटेंशियल और वैक्यूम में एक लीनियर ड्रॉप दिखाता है। यह पुष्टि करता है कि हाई बायस पर सिमेट्रिक चार्ज डेंसिटी और सही हार्टी पोटेंशियलल उत्पन्न करने के लिए वेटेड डबल-कॉन्टूर स्कीम आवश्यक है, जबकि सिंगल-कॉन्टूर इवैल्यूएशन एसिमेट्रिक और गलत परिणाम देते हैं।
3. Au-H$_2$-Au जंक्शन ज्योमेट्री ऑप्टिमाइजेशन: लागू बायस (1.5 V तक) के तहत, इंटरफेस H-H बॉन्ड के विस्तार और परमाणु विस्थापन की सही भविष्यवाणी करता है, जो पिछले SIESTA+SMEAGOL परिणामों के अनुरूप है। इन परिवर्तनों को इलेक्ट्रिक फील्ड के बजाय इलेक्ट्रिक करंट के कारण माना गया है।
4. एक सॉल्वेटेड Au वायर का मॉलिक्यूलर डायनेमिक्स: लेखकों ने एक बड़े पैमाने के कंडेंस्ड फेज सिस्टम (166 जल अणुओं द्वारा सॉल्वेटेड एक मोनोएटमिक Au वायर) पर पहले DFT-NEGF मॉलिक्यूलर डायनेमिक्स सिमुलेशन किए।
   • स्थिरता: ओपन-सिस्टम की प्रकृति के बावजूद, डायनेमिक्स न्यूनतम लॉन्ग-टर्म एनर्जी ड्रिफ्ट के साथ स्थिर रहे।
   • इलेक्ट्रॉन माइग्रेशन: 1 V बायस पर 200 fs ट्रैजेक्टरी में, उन्होंने इलेक्ट्रॉन माइग्रेशन प्रभावों को देखा, जहाँ स्थानीय इलेक्ट्रॉन डेंसिटी मिनिममम केंद्रीय वायर परमाणु से हटकर दाएंतम परमाणु की ओर स्थानांतरित हो गया, जो दाएं इलेक्ट्रोड के पॉजिटिव चार्ज को स्क्रीन करने के अनुरूप है।
   • ट्रांसमिशन: सॉल्वेटेड वायर के लिए ट्रांसमिशन स्पेक्ट्रा एक आदर्श अनंत वायर के गुणात्मक रूप से समान थे, जिसमें अंतर जल के अणुओं द्वारा खोले गए अतिरिक्त टनलिंग चैनलों के कारण थे।

प्रदर्शन
प्राथमिक कम्प्यूटेशनल बाधा ग्रीन'स फंक्शन के इंटीग्रल के रूप में इलेक्ट्रॉन डेंसिटी का मूल्यांकन करना है, जिसके लिए प्रति SCF स्टेप 96 एनर्जी पॉइंट कंप्यूटेशन (64 रियल एक्सिस, 32 कॉम्प्लेक्स एक्सिस) की आवश्यकता होती है, जो एक सिंगल हैमिल्टोनियन डायगोनालाइजेशन को प्रतिस्थापित करता है।

• पैरेलरालाइजेशन: कोड रियल-स्पेस कंप्यूटेशन और k-पॉइंट्स के लिए MPI, और DO लूप्स एवं LAPACK रूटीन के लिए OpenMP का उपयोग करता है। प्रदर्शन 32 OpenMP थ्रेड्स तक अच्छी तरह से स्केल करता है।
• लागत: एक सॉल्वेटेड Au वायर सिस्टम के लिए, CP2K+SMEAGOL गणना प्रति SCF स्टेप एक मानक CP2K गणना की तुलना में लगभग 3.2 गुना धीमी थी, जबकि इसमें 8 गुना अधिक कोर (2048 बनाम 256) का उपयोग किया गया था। लेखक नोट करते हैं कि यह लागत में एक ऑर्डर-ऑफ-मैग्निट्यूड वृद्धि को दर्शाता है, जिससे बिना और अधिक त्वरण (जैसे ScaLAPACK या मशीन लर्निंग फोर्स फील्ड) के लंबे MD ट्राजेक्टरीज चुनौतीपूर्ण हो जाते हैं।

महत्व और दावे
यह पेपर दावा करता है कि CP2K+SMEAGOL अब तक उपलब्ध सबसे फीचर-रिच और पूर्ण DFT-NEGF कार्यान्वयन है, जो नेटिव CP2K कार्यान्वयन से बेहतर है जो बिना फोर्सेस के केवल $\Gamma$-पॉइंट ट्रांसपोर्ट तक सीमित है।

इस कार्य का प्राथमिक महत्व यह प्रदर्शन करना है कि DFT-NEGF का उपयोग वास्तविक ऑपरेटिंग स्थितियों के तहत बड़े पैमाने के कंडेंस्ड फेज सिस्टम के लिए मॉलिक्यूलर डायनेमिक्स सिमुलेशन करने के लिए किया जा सकता है। लेखकों की जानकारी के अनुसार, यह ऐसे पहले गणनाओं में से एक है। यह इंटरफेस इलेक्ट्रो-कैटेलिटिक प्रतिक्रियाओं और इलेक्ट्रोकेमिकल प्रणालियों के अध्ययन को सक्षम बनाता है जहाँ लागू पोटेंशियल, करंट फ्लो और परमाणु डायनेमिक्स के बीच का परस्पर संबंध महत्वपूर्ण है, जो सूक्ष्म सैद्धांतिक मॉडलों और मैक्रोस्कोपिक प्रयोगात्मक स्थितियों के बीच के अंतर को पाटता है।