Density-functional theory calculation of hydrogen solubility in cubic silicon carbide at finite temperatures

तकनीकी सारांश: परिमित तापमान पर क्यूबिक सिलिकॉन कार्बाइड में हाइड्रोजन घुलनशीलता का घनत्व-कार्यात्मक सिद्धांत (Density-Functional Theory) गणना

समस्या विवरण
फ्यूजन ऊर्जा तकनीक की प्रगति के लिए, विशेष रूप से प्रथम-दीवार (first-wall) सामग्रियों के लिए, प्रभावी ट्रिटियम पारगम्यता अवरोध (TPBs) अत्यंत महत्वपूर्ण हैं। सिलिकॉन कार्बाइड (SiC) एक प्रमुख TPB उम्मीदवार है; हालाँकि, SiC के लिए प्रयोगात्मक पारगम्यता माप छह क्रमों (orders of magnitude) तक भिन्न होते हैं। ये विसंगतियां आदर्श एकल क्रिस्टल और वास्तविक, दोषयुक्त सामग्रियों के बीच अंतर, साथ ही माप तकनीकों और नमूना तैयारी (जैसे, CVD पैरामीटर, विकिरण-प्रेरित दोष) में भिन्नता के कारण होती हैं। जबकि हाइड्रोजन (H) पारगम्यता को व्यापक रूप से मापा गया है, SiC में हाइड्रोजन घुलनशीलता पर शोध कम हुआ है। मौजूदा शोध में इस बात पर व्यापक डेटा का अभाव है कि तापमान और हाइड्रोजन आंशिक दबाव के फलन (function) के रूप में बाइंडिंग ऊर्जा कैसे बदलती है, विशेष रूप से इंटरस्टिशियल, रिक्तियों (vacancies) और गैर-स्टोइकोमेट्रिक (अमोर्फस) क्षेत्रों जैसे सूक्ष्म संरचनात्मक दोषों के प्रभाव के संबंध में। यह अंतराल सटीक मॉडल विकसित करने में बाधा डालता है जो हाइड्रोजन पारगम्यता की भविष्यवाणी कर सकें, क्योंकि पारगम्यता स्वाभाविक रूप से घुलनशीलता पर निर्भर करती है।

कार्यप्रणाली
लेखकों ने शुद्ध और दोषयुक्त $\beta$-SiC (क्यूबिक SiC) में हाइड्रोजन घुलनशीलता की भविष्यवाणी करने के लिए 'एब इनिटियो' (ab initio) ढांचे का विकास किया।

• कंप्यूटेशनल सेटअप: गणना DMol3 कोड का उपयोग करके की गई थी जिसमें जेनेलाइज्ड ग्रेडिएंट एप्रोक्सिमेशन (GGA-PBE) का उपयोग किया गया था। एक 2$\times$2$\times$2$\times$2 सुपरसेल (64 परमाणु) के साथ 4$\times$4$\times$4 k-पॉइंट मेश का उपयोग किया गया।
• संरचनात्मक मॉडल: अध्ययन ने चार प्रकार के इंटरस्टिशियल साइटों, शुद्ध क्रिस्टल में सिलिकॉन रिक्तियों ($V_{Si}$), कार्बन रिक्तियों ($V_C$), और गैर-स्टोइकोमेट्रिक अमोर्फस संरचनाओं को मॉडल किया। अमोर्फस क्षेत्रों को (001) प्लेन के साथ परमाणुओं को हटाकर कार्बन-समृद्ध (C-rich) और सिलिकॉन-समृद्ध (Si-rich) विक्षेपित क्षेत्रों के रूप में बनाया गया था, जो सुपरसेल के लगभग 25 वॉल्यूम % भाग को घेरते हैं।
• ऊष्मागतिकी (Thermodynamics): परिमित तापमान और दबावों को ध्यान में रखने के लिए, लेखकों ने 'एब इनिटियो' ऊष्मागतिकी का उपयोग किया। इसमें DFT कुल ऊर्जाओं को कंपन योगदान (शून्य-बिंदु ऊर्जा और फोनन डिस्पर्शन से प्राप्त परिमित-तापमान एन्ट्रॉपी/एन्थैल्पी) के साथ जोड़कर निर्माण की गिब्स मुक्त ऊर्जा ($\Delta G_f$) की गणना की गई। H के रासायनिक विभव (chemical potentials) को 0 से 1000 K के तापमान रेंज और 0.02 Pa से 0.1 MPa की आंशिक दबाव सीमा के लिए JANAF थर्मोकेमिकल डेटा और आदर्श गैस नियम का उपयोग करके प्राप्त किया गया।
• घुलनशीलता गणना: घुलनशीलता ($\theta$) को गैसीय $H_2$ और घुले हुए परमाणु H के रासायनिक विभवों को समान रखकर निकाला गया, जिसमें उन विभाजन फलनों (partition functions) का उपयोग किया गया जो गणना की गई $\Delta G_f$ को शामिल करते हैं। यह दृष्टिकोण ली एट अल. (Lee et al.) की विधियों को संशोधित करता है, जिसमें इलेक्ट्रॉनिक अपभ्रष्टता (degeneracy) पदों को हटाया गया है (जो अर्धचालकों के लिए कम प्रासंगिक हैं) और परिमित-तापमान मुक्त ऊर्जाओं को एकीकृत किया गया है।

मुख्य परिणाम

1. बाइंडिंग ऊर्जा और स्थिरता:

   • शुद्ध $\beta$-SiC में, सबसे अनुकूल इंटरस्टिशियल साइट चार Si परमाणुओं द्वारा समन्वित टेट्राहेड्रल साइट ($T_{Si}$) है, जिसकी एंडोथर्मिक बाइंडिंग ऊर्जा 2.87 eV है।
   • रिक्ति निर्माण ऊर्जा दर्शाती है कि $V_C$ (4.78 eV) $V_{Si}$ (7.89 eV) की तुलना में अधिक आसानी से बनता है। 4H-SiC के विपरीत, $\beta$-SiC में $V_{Si}$ अपघटन के विरुद्ध स्थिर है।
   • $V_{Si} + H$ कॉम्प्लेक्स, $V_C + H$ कॉम्प्लेक्स की तुलना में 0.91 eV से अधिक ऊष्मागतिक रूप से अनुकूल है। $V_{Si}$ में H बाइंडिंग स्थानीय वातावरण पर अत्यधिक निर्भर है, जो -1.88 से -0.97 eV के बीच है।
   • अमोर्फस संरचनाएं: C-समृद्ध अमोर्फस संरचना ने सबसे अनुकूल बाइंडिंग ऊर्जा (-2.15 eV) प्रदर्शित की, मुख्य रूप से वहां जहां H एक Si–H–Si बॉन्ड बनाता है। इसके विपरीत, Si-समृद्ध अमोर्फस संरचनाओं ने कम अनुकूल बाइंडिंग ऊर्जा (-0.10 से 0.56 eV) दिखाई।

2. ऊष्मागतिक स्थिरता ($\Delta G_f$):

   • सभी तापमानों $\le$ 1000 K पर, C-समृद्ध अमोर्फस संरचनाओं और $V_{Si}$ दोषों में H के लिए नकारात्मक $\Delta G_f$ पाया गया, जो ऊष्मागतिक स्थिरता और ट्रैपिंग (trapping) को दर्शाता है।
   • इंटरस्टिशियल H और $V_C + H$ विन्यास पूरे तापमान रेंज में सकारात्मक $\Delta G_f$ दिखाते हैं, जो ऊष्मागतिक अस्थिरता और ट्रैपिंग की कमी को दर्शाता है।

3. घुलनशीलता:

   • घुलनशीलता, $\Delta G_f$ के साथ दृढ़ता से सह-संबंधित है। उच्चतम घुलनशीलता C-समृद्ध अमोर्फस क्षेत्रों (विशेष रूप से वहां जहां H, Si–H–Si बॉन्ड बनाता है) और $V_{Si}$ दोषों में पाई गई।
   • 600 K पर, C-समृद्ध अमोर्फस संरचनाओं में घुलनशीलता $7.73 \times 10^4$ mol H m$^{-3}$ तक पहुँच गई, और $V_{Si}$ में यह $1.64 \times 10^{-2}$ mol H m$^{-3}$ थी।
   • शुद्ध इंटरस्टिशियल साइटों और $V_C$ दोषों का घुलनशीलता में नगण्य योगदान था।
   • अध्ययन में देखा गया कि जिन विन्यासों में नकारात्मक बाइंडिंग ऊर्जा होती है, उनमें बढ़ते तापमान के साथ स्थिर या घटती हुई घुलनशीलता होती है, जबकि सकारात्मक बाइंडिंग ऊर्जा वाले विन्यासों में बढ़ती हुई घुलनशीलता देखी जाती है।

महत्व और दावे
यह पत्र दावा करता है कि यह $\beta$-SiC में हाइड्रोजन व्यवहार के संबंध में पहला मात्रात्मक ऊष्मागतिक अंतर्दृष्टि प्रदान करता है, जो फ्यूजन रिएक्टर संचालन के लिए प्रासंगिक स्थितियों के तहत है, विशेष रूप से दोषयुक्त प्रणालियों के लिए घुलनशीलता डेटा की कमी को संबोधित करता है।

• विसंगतियों का स्पष्टीकरण: परिणाम बताते हैं कि प्रयोगात्मक पारगम्यता मानों में व्यापक भिन्नता विशिष्ट दोषों की उपस्थिति के कारण हो सकती है, विशेष रूप से $V_{Si}$ और गैर-स्टोइकोमेट्रिक C-समृद्ध क्षेत्र, जो शुद्ध क्रिस्टल की तुलना में हाइड्रोजन घुलनशीलता को काफी बढ़ा देते हैं।
• सूक्ष्म संरचनात्मक इंजीनियरिंग (Microstructural Engineering): लेखक प्रस्ताव देते हैं कि $V_{Si}$ सांद्रता और गैर-स्टोइकोमेट्रिक C-समृद्ध क्षेत्रों को कम करने के लिए निक्षेपण मापदंडों (जैसे, CVD प्रक्रियाओं में) को नियंत्रित करना ट्रिटियम पारगम्यता अवरोध प्रदर्शन को बेहतर बनाने के लिए एक व्यवहार्य रणनीति है।
• मॉडलिंग इनपुट: यह अध्ययन बहु-स्तरीय (multiscale) ट्रिटियम परिवहन मॉडलों के लिए आवश्यक इनपुट प्रदान करता है, जो तापमान, आंशिक दबाव और विशिष्ट दोष घनत्व के फलन के रूप में घुलनशीलता को मात्रात्मक रूप से निर्धारित करता है।
• सीमाएं: लेखक विनम्रतापूर्वक नोट करते हैं कि उनका मॉडल निश्चित दोष सांद्रता मानता है और अमोर्फस क्षेत्रों को पूर्ण सांख्यिकीय परीक्षण के बजाय सीमित सैंपलिंग के माध्यम से मानता है। वे सुझाव देते कि भविष्य के कार्यों में दोष घनत्वों को मान्य करने और ग्रेन बाउंडरी प्रभावों को शामिल करने के लिए प्रयोगात्मक मापों (जैसे, पॉज़िट्रॉन एनिहिलेशन स्पेक्ट्रोस्कोपी) को एकीकृत किया जाना चाहिए।