Controlling HER activity and stability of $γ$- and 6,6,12-Graphyne through engineered B-N doping: DFT and Reactive MD simulations

तकनीकी सारांश: इंजीनियर बी-एन (B-N) डोपिंग के माध्यम से $\gamma$- और 6,6,12-ग्राफिन की HER गतिविधि और स्थिरता को नियंत्रित करना

समस्या विवरण
हाइड्रोजन इवोल्यूशन रिएक्शन (HER) सतत हाइड्रोजन उत्पादन के लिए महत्वपूर्ण है, फिर भी व्यावहारिक अनुप्रयोगों के लिए ऐसे उत्प्रेरकों की आवश्यकता है जो हाइड्रोजन-समृद्ध स्थितियों में उच्च गतिविधि और दीर्घकालिक स्थिरता के बीच संतुलन बना सकें। जबकि उत्कृष्ट धातुएं (noble metals) मानक हैं, धातु-मुक्त, पृथ्वी-प्रचुर विकल्पों की आवश्यकता है। द्वि-आयामी कार्बन एलोट्रोप्स (allotropes), विशेष रूप से ग्राफिन (sp/sp$^2$ संकरित कार्बन शीट), एक ट्यूनेबल इलेक्ट्रॉनिक संरचना प्रदान करते हैं। हालांकि, दो विशिष्ट चुनौतियां ग्राफिन-आधारित HER उत्प्रेरकों के लिए बनी हुई हैं: (1) यह समझना कि बोरॉन-नाइट्रोजन (B-N) सह-डोपेंट (meta, ortho, para) की सापेक्ष ज्यामिति कैसे दोष ऊष्मागतिकी (defect thermodynamics) और फर्मी स्तर ($E_F$) के पास इलेक्ट्रॉनिक अवस्थाओं को एक साथ नियंत्रित करती है; और (2) यह निर्धारित करना कि 0 K पर हाइड्रोजन अधिशोषण मुक्त ऊर्जा ($\Delta G_{ads}$) को अनुकूलित करने के लिए अनुमानित मोटिफ (motifs), परिमित तापमान पर हाइड्रोजनीकरण-प्रेरित बॉन्ड विखंडन (bond scission) के विरुद्ध गतिज रूप से कितने सुदृढ़ रहते हैं।

कार्यप्रणाली
यह अध्ययन प्रिमिटिव (pristine), B-डोप्ड, N-डोप्ड और B-N सह-डोप्ड $\gamma$-ग्राफिन और 6,6,12-ग्राफिन का मूल्यांकन करने के लिए डेंसिटी फंक्शनल थ्योरी (DFT) और रिएक्टिव मॉलिक्यूलर डायनेमिक्स (ReaxFF) सिमुलेशन के संयोजन का उपयोग करते हुए एक बहु-स्तरीय गणनात्मक दृष्टिकोण अपनाता है।

• DFT गणना: SIESTA कोड का उपयोग करके GGA-PBE फंक्शनल्स और DZP बेस सेट के साथ की गई। अध्ययन ने इलेक्ट्रॉनिक बैंड संरचना, दोष निर्माण ऊर्जा ($E_{form}$), और कम्प्यूटेशनल हाइड्रोजन इलेक्ट्रोड ढांचे के भीतर गिब्स मुक्त ऊर्जा ऑफ हाइड्रोजन अधिशोषण ($\Delta G_{ads}$) की गणना की। स्पिन-पोलराइज्ड गणनाओं ने पुष्टि की कि संरचनाएं गैर-चुंबकीय रहती हैं।
• रिएक्टिव MD सिमुलेशन: परमाणु हाइड्रोजन के ग्राफिन शीट के साथ संपर्क का 300 K पर मूल्यांकन करने के लिए LAMMPS के साथ ReaxFF फोर्स फील्ड का उपयोग करके किए गए। इन सिमुलेशन को परिमित-तापमान संरचनात्मक स्थिरता और हाइड्रोजन अपटेक गतिकी (kinetics) का आकलन करने के लिए 500 ps तक NVT एन्सेम्बल में चलाया गया, जिसमें बॉन्ड टूटने और बनने का स्पष्ट रूप से मॉडलिंग किया गया।

मुख्य योगदान और परिणाम

1. इलेक्ट्रॉनिक संरचना और डोपिंग प्रभाव:

   • प्रिमिटिव $\gamma$-ग्राफिन एक अर्धचालक (semiconductor) है, जबकि 6,6,12-ग्राफिन एक एनिसोट्रोपिक डिरैक-जैसे सेमी-मेटैलिक फैलाव (dispersion) को प्रदर्शित करता है।
   • एकल B या N प्रतिस्थापन, केवल रिजिड बैंड शिफ्ट के बजाय डोपेंट-$\pi$ संकरण (hybridization) के माध्यम से निकट-$E_F$ अवस्थाओं को पुनर्गठित करता है। B एक एक्सेप्टर (acceptor) के रूप में कार्य करता है, और N एक डोनर (donor) के रूप में।
   • B-N सह-डोपिंग एक ज्यामिति-निर्भर विक्षोभ (perturbation) पेश करती है। B-N ऑर्थो (ortho) विन्यास को डोनर-एक्सेप्टर चार्ज मुआवजे के कारण सबसे अधिक ऊष्मागतिक रूप से अनुकूल पाया गया, जिसके बाद मेटा (meta) और पैरा (para) का स्थान आता है।

2. ऊष्मागतिक गतिविधि ($\Delta G_{ads}$):

   • प्रिमिटिव जाली (lattices) कमजोर हाइड्रोजन बाइंडिंग प्रदर्शित करती हैं, जो उन्हें HER के लिए उप-इष्टतम बनाती हैं।
   • डोपिंग हाइड्रोजन अधिशोषण के लिए "हॉट स्पॉट्स" बनाती है। महत्वपूर्ण रूप से, सबसे सक्रिय साइटें स्वयं डोपेंट परमाणु नहीं हैं, बल्कि आस-पास के कार्बन साइट्स हैं, विशेष रूप से वे जो sp-संकरित एसिटिलीनिक श्रृंखलाओं में हैं।
   • अकेले N-डोपिंग, N की उच्च विद्युत ऋणात्मकता के कारण सीधे H-बाइंडिंग के लिए ऊष्मागतिक रूप से प्रतिकूल है; इसके बजाय, यह आसपास के ढांचे को ध्रुवीकृत (polarize) करती है।
   • B-N सह-डोपिंग $\Delta G_{ads}$ को सबेटियर ऑप्टिमम ($\approx 0$ eV) की ओर ट्यून करती है। 6,6,12-ग्राफिन के लिए, ऑर्थो प्रतिस्थापन विशिष्ट sp-समीपवर्ती साइटों को सीधे थर्मोन्यूट्रल विंडो ($\Delta G_{ads} \approx 0.03$ eV) के भीतर रखता है।

3. गतिज स्थिरता और गतिविधि-स्थिरता ट्रेड-ऑफ:

   • 300 K पर रिएक्टिव MD एक महत्वपूर्ण ट्रेड-ऑफ को प्रकट करता है: वे विन्यास जो इलेक्ट्रॉनिक परिदृश्य को दृढ़ता से सक्रिय करते हैं, वे ओवर-हाइड्रोजनीकरण के माध्यम से विनाशकारी बॉन्ड विखंडन को भी ट्रिगर कर सकते हैं।
   • $\gamma$-ग्राफिन: B-N ऑर्थो विन्यास अनियंत्रित संरचनात्मक क्षरण के बिना नियंत्रित हाइड्रोजन अपटेक को अनूठे रूप से बनाए रखता है। इसके विपरीत, B-N मेटा और पैरा प्रतिस्थापन बॉन्ड क्लेवेज (bond cleavage) की ओर ले जाते हैं। एकल B या N डोपिंग संरचनात्मक रूप से बरकरार रहती है।
   • 6,6,12-ग्राफिन: यह जाली आम तौर पर ओवर-हाइड्रोजनीकरण और बॉन्ड टूटने के प्रति अधिक संवेदनशील है। यहाँ तक कि प्रिमिटिव शीट भी प्रारंभिक क्षरण (incipient degradation) दिखाती है, और B-N सह-डोपिंग, विशेष रूप से मेटा और पैरा विन्यास के लिए, गंभीर नेटवर्क विरूपण और बॉन्ड विखंडन की ओर ले जाती है।

महत्व और दावे
यह शोध स्थापित करता है कि B-N ज्यामिति (विशेष रूप से डोपेंट्स की सापेक्ष स्थिति) ग्राफिन-आधारित HER उत्प्रेरकों के लिए एक प्राथमिक डिजाइन वेरिएबल है, जो चार्ज पुनर्वितरण और संरचनात्मक स्थिरता को सह-नियंत्रित करती है।

लेखकों का दावा है कि जबकि एक निकट थर्मोन्यूट्रल $\Delta G_{ads}$ प्राप्त करना एक आवश्यक ऊष्मागतिक मानदंड है, यह अपने आप में पर्याप्त नहीं है। व्यावहारिक उत्प्रेरकों को हाइड्रोजनीकरण-प्रेरित गिरावट के विरुद्ध गतिज सुदृढ़ता (kinetic robustness) को भी संतुष्ट करना चाहिए। अध्ययन B-N ऑर्थो-डोप्ड $\gamma$-ग्राफिन को एक आशाजनक उम्मीदवार के रूप में पहचानता है जो अनुकूल अधिशोषण ऊर्जा के साथ परिमित-तापमान स्थिरता को संतुलित करता है। इसके विपरीत, परिणाम बताते हैं कि 6,6,12-ग्राफिन, विशेष रूप से सह-डोप्ड होने पर, हाइड्रोजन-समृद्ध स्थितियों के तहत आंतरिक रूप से कम स्थिर है। यह कार्य इस बात पर जोर देता है कि उत्प्रेरक रूप से प्रासंगिक मोटिफ एक "BN-ध्रुवीकृत कार्बन वातावरण" (विशेष रूप से आस-पास के sp-कार्बन) है, न कि स्वयं हेटरोएटम्स (heteroatoms); और इलेक्ट्रॉनिक सक्रियण एवं संरचनात्मक अखंडता के बीच के अंतर्संबंध को एक साथ अनुकूलित किया जाना चाहिए।