Numerical Modeling of Solvent Diffusion through the Transition Metal Dichalcogenides based Nanomaterials

यहाँ गीतिका साहू के शोध पत्र "Numerical Modeling of Solvent Diffusion through the Transition Metal Dichalcogenides based Nanomaterials" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

ट्रांजिशन मेटल डाइचैल्कोजेनाइड्स (TMDCs) बहुमुखी 2D सामग्रियां हैं जिन्हें उनके इलेक्ट्रॉनिक और ऑप्टिकल गुणों को ट्यून करने के लिए शून्य-आयामी क्वांटम डॉट्स (QDs) तक छोटा किया जा सकता है। हालाँकि, TMDC-आधारित QDs के संश्लेषण में एक महत्वपूर्ण चुनौती सॉल्वोथर्मल विधियों (solvothermal methods) के माध्यम से कण आकार वितरण (particle size distribution) और एकरूपता पर सटीक नियंत्रण प्राप्त करना है।

संश्लेषण प्रक्रिया में TMDCs के कमजोर वैन डेर वाल्स (van der Waals) परतों के बीच विलायक (solvent) का अंतर्वेशन (intercalation) शामिल है, जिससे एक्सफोलिएशन (परत अलग होना) होता है। यह प्रक्रिया जटिल प्रयोगात्मक मापदंडों द्वारा नियंत्रित होती है, जिनमें शामिल हैं:

• विलायक का चयन (pH, आणविक आकार, चयनात्मकता)।
• विलायक की विसरणशीलता (Diffusivity)।
• अभिक्रिया समय (एक्सपोज़र की अवधि)।
• तापमान और अग्रदूत (precursor) सांद्रता।

वर्तमान प्रयोगात्मक दृष्टिकोण यह अनुमान लगाने में संघर्ष करते हैं कि ये पैरामीटर अंतिम नैनोकण के आकार को निर्धारित करने के लिए कैसे परस्पर क्रिया करते हैं। एक समान QD संश्लेषण के लिए इन मापदंडों को अनुकूलित करने हेतु एक भविष्य कहनेवाला संख्यात्मक ढांचे (predictive numerical framework) की आवश्यकता है।

2. कार्यप्रणाली (Methodology)

यह अध्ययन एक संख्यात्मक सिमुलेशन (numerical simulation) दृष्टिकोण का उपयोग करता है जो सॉल्वोथर्मल बाथ में एक नैनोकण प्रणाली के विकास को सिम्युलेट करने के लिए संशोधित फिक के प्रसार नियम (Modified Fick's Law of Diffusion) को डायनेमिक बॉन्ड परकोलेशन मॉडल (DBPM) के साथ जोड़ता है।

• सिस्टम सेटअप:

  • मॉडल TMDC नैनोकणों को स्टैक्ड परतों (सैंडविच संरचना) के रूप में दर्शाता है जहाँ परतों के भीतर मजबूत सहसंयोजक बंध (covalent bonds) होते हैं और परतों के बीच कमजोर वैन डेर वाल्स बल होते हैं।
  • प्रारंभिक अवस्था: 50 कॉन्फ़िगरेशन, प्रत्येक में 50 परतें (कुल मोटाई ~50 nm)। संरचनात्मक विकार (structural disorder) की नकल करने के लिए परत की चौड़ाई यादृच्छिक रूप से भिन्न होती है (1–15 nm)।
  • विलायक का प्रसार केवल अंतर-परत स्थानों (परतों के बीच) तक ही सीमित है, न कि एक एकल परत के भीतर।

• गणितीय ढांचा (Mathematical Framework):

  • विसरण समीकरण (Diffusion Equation): एक साइट $(i, j)$ पर विलायक स्तर ($P$) को संशोधित फिक के दूसरे नियम के परिमित अंतर सन्निकटन (finite difference approximation) का उपयोग करके अपडेट किया जाता है:
    $$P^{n+1}_{i,j} = P^n_{i,j} + s \sum (W \cdot P^n_{neighbors} - 4W \cdot P^n_{i,j})$$
    जहाँ $s$ एक विसरणशीलता चर (diffusivity variable) है (जो विलायक की शक्ति/अभिक्रिया स्थितियों का प्रतिनिधित्व करता है) और $W$ एक यादृच्छिक हॉपिंग प्रायिकता (hopping probability) है (जो सामग्री के दोषों/विषमता का प्रतिनिधित्व करता है)।
  • बंध टूटने की प्रक्रिया (Bond Breaking Mechanism): परतों के बीच एक बंध तब "टूटा हुआ" (एक्सफोलिएटेड) माना जाता है जब स्थानीय विलायक स्तर एक सीमा से अधिक हो जाता है ($P_{i,j} > 0.9$)। जब एक परत के सभी बंध टूट जाते हैं, तो कण छोटे टुकड़ों में विभाजित हो जाता है।

• सिमुलेशन पैरामीटर:

  • विसरणशीलता चर ($s$): विभिन्न विलायकों या अभिक्रिया तीव्रता को सिम्युलेट करने के लिए 0.1 से 0.9 तक बदला गया।
  • पुनरावृत्ति (Iterations - $I$): 100 पुनरावृत्तियों तक, जो अभिक्रिया समय के प्रॉक्सी के रूप में कार्य करती है।
  • विश्लेषण मेट्रिक्स: औसत कण आकार ($\bar{d}$), एवलांच सांख्यिकी (नए बंध टूटना), शैनन एंट्रॉपी ($E$) आकार वितरण की एकरूपता के लिए, और आकार में सापेक्ष परिवर्तन ($\Delta \bar{d} / \Delta I$)।

3. मुख्य योगदान (Key Contributions)

1. भविष्य कहनेवाला ढांचा (Predictive Framework): विलायक विसरणशीलता और अभिक्रिया समय (पुनरावृत्तियों) को सीधे नैनोकण आकार के विकास और एक्सफोलिएशन गतिकी से जोड़ने वाला एक संख्यात्मक मॉडल स्थापित किया।
2. महत्वपूर्ण थ्रेशोल्ड की पहचान: एक महत्वपूर्ण विसरणशीलता थ्रेशोल्ड ($s \approx 0.6$) की पहचान की जो यह निर्धारित करता है कि सिस्टम स्थिर रहेगा या तीव्र एक्सफोलिएशन से गुजरेगा।
3. एंट्रॉपी-आधारित एकरूपता विश्लेषण: आकार वितरण की गैर-एकरूपता से एकरूपता में संक्रमण को मापने के लिए शैनन एंट्रॉपी (Shannon Entropy) का उपयोग किया, और उन विशिष्ट पुनरावृत्ति बिंदुओं की पहचान की जहाँ सिस्टम संतृप्ति (saturation) तक पहुँच जाता है।
4. पावर लॉ सहसंबंध (Power Law Correlation): अधिकतम एंट्रॉपी (अधिकतम आकार उतार-चढ़ाव) के लिए आवश्यक पुनरावृत्ति और न्यूनतम आकार परिवर्तन दर (संतृप्ति की शुरुआत) के लिए आवश्यक पुनरावृत्ति के बीच एक पावर-लॉ सहसंबंध की खोज की।

4. मुख्य परिणाम (Key Results)

A. आकार विकास और विसरणशीलता ($s$)

• कम विसरणशीलता ($s \leq 0.3$): औसत कण आकार ($\bar{d}$) 100 पुनरावृत्तियों के बाद भी स्थिर (~18 nm) रहता है। विलायक की शक्ति अंतर-परत वैन डेर वाल्स बलों को दूर करने के लिए अपर्याप्त है।
• मध्यवर्ती विसरणशीलता ($0.4 \leq s \leq 0.6$): $\bar{d}$ धीरे-धीरे कम होता है। एक्सफोलिएशन होता है लेकिन विलंबित होता है, जिससे महत्वपूर्ण आकार कमी प्राप्त करने के लिए अधिक पुनरावृत्तियों की आवश्यकता होती है।
• उच्च विसरणशीलता ($s \geq 0.7$): तीव्र एक्सफोलिएशन होता है। $s=0.9$ के लिए, $\bar{d}$ 100 पुनरावृत्तियों के भीतर 18 nm से घटकर ~6.5 nm हो जाता है, जिसमें महत्वपूर्ण कमी पहले 25 पुनरावृत्तियों में ही हो जाती है।

B. एवलांच सांख्यिकी (बंध टूटना)

• उच्च $s$ मानों के लिए नए टूटे हुए बंधों की संख्या ($\bar{b}_n$) जल्दी चरम (peak) पर पहुँचती है, जो तीव्र प्रारंभिक विखंडन का संकेत देती है।
• उच्च $s$ के लिए, सिस्टम एक ऐसी स्थिति में पहुँच जाता है जहाँ टूटने के लिए कम बड़े कण शेष रह जाते हैं, जिससे $\bar{b}_n$ बाद के पुनरावृत्तियों में घट जाता है।
• औसत आकार और कुल टूटे हुए बंधों के बीच संबंध घातांकीय क्षय (exponential decay) का अनुसरण करता है: $\bar{d} = A e^{-m\bar{b}_t}$।

C. सर्वाधिक बार होने वाला आकार ($d^*$) और स्विचिंग व्यवहार

• सिस्टम या तो प्रारंभिक बड़े आकार (18 nm) या एक छोटे संतृप्त आकार (3 nm) द्वारा हावी रहता है।
• एक "स्विच" एक महत्वपूर्ण पुनरावृत्ति ($I_c$) पर होता है जहाँ प्रभावी आकार बड़े से छोटे की ओर बदल जाता है।
• महत्वपूर्ण थ्रेशोल्ड: यह स्विच व्यवहार केवल तभी शुरू होता है जब $s > 0.6$ हो। इससे नीचे, सिस्टम बड़े कणों द्वारा हावी रहता है।
• $I_c$ और $s$ के बीच का संबंध रैखिक है लेकिन $s=0.68$ पर ढलान बदल देता है, जो विसरण गतिकी के दो अलग-अलग शासन (regimes) को दर्शाता है।

D. एंट्रॉपी और संतृप्ति (Entropy and Saturation)

• शैनन एंट्रॉपी ($E$): आकार वितरण की अव्यवस्था (गैर-एकरूपता) को मापता है।
  • कम $s$: $E$ निम्न रहता है (कोई परिवर्तन नहीं)।
  • मध्यवर्ती $s$: $E$ बढ़ता है, जो मिश्रित आकारों (गैर-एकरूपता) को दर्शाता है।
  • उच्च $s$: $E$ पहले तेजी से बढ़ता है और फिर घटता है, जो यह दर्शाता है कि सिस्टम एक समान, छोटे आकार के वितरण की ओर बढ़ रहा है।
• संतृप्ति बिंदु: अध्ययन ने अधिकतम एंट्रॉपी ($I''$) के लिए पुनरावृत्ति और आकार परिवर्तन की न्यूनतम दर ($I'$) के लिए पुनरावृत्ति के बीच एक सहसंबंध स्थापित किया। यह सहसंबंध (~2.01 घातांक के साथ पावर लॉ) पुष्टि करता है कि एक विशिष्ट पुनरावृत्ति बिंदु मौजूद है जहाँ सिस्टम आकार के उतार-चढ़ाव से संतृप्त, एक समान अवस्था में परिवर्तित होता है।

5. महत्व और निष्कर्ष (Significance and Conclusion)

यह अध्ययन TMDC-आधारित क्वांटम डॉट्स के संश्लेषण को अनुकूलित करने के लिए एक मजबूत कम्प्यूटेशनल उपकरण प्रदान करता है।

• अनुकूलन (Optimization): यह प्रदर्शित करता है कि केवल अभिक्रिया समय बढ़ाना पर्याप्त नहीं है; प्रभावी एक्सफोलिएशन प्राप्त करने के लिए विलायक विसरणशीलता को एक महत्वपूर्ण थ्रेशोल्ड ($s > 0.6$) से ऊपर होना चाहिए।
• एकरूपता नियंत्रण (Uniformity Control): एंट्रॉपी और आकार परिवर्तन की दर की निगरानी करके, शोधकर्ता एक समान आकार वितरण प्राप्त करने के लिए आवश्यक सटीक "अभिक्रिया समय" (पुनरावृत्तियों) की भविष्यवाणी कर सकते हैं, जिससे ओवर-प्रोसेसिंग या अंडर-प्रोसेसिंग से बचा जा सकता है।
• भविष्य का दायरा: लेखक नोट करते हैं कि वास्तविक दुनिया के नैनोस्केल सिस्टम के लिए सटीकता को और बेहतर बनाने के लिए मॉडल में सतह ऊर्जा (surface energy) (विशेष रूप से डैंगलिंग बॉन्ड्स और सर्फेक्टेंट प्रभाव) को शामिल करने की आवश्यकता है।

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