🔬 materials science

A Unified microscopic picture of cation and anion migration in MAPbI $_3$

न्यूरल नेटवर्क पोटेंशियल के साथ मॉलिक्यूलर डायनेमिक्स सिमुलेशन का उपयोग करते हुए, यह अध्ययन प्रकट करता है कि MAPbI $_3$ में तीव्र आयन प्रवास (ion migration) MA इंटरस्टिशियल और चार्ज-डिपेंडेंट I-संबंधित दोषों की समन्वित सामूहिक गति द्वारा संचालित होता है, जबकि MA रिक्तियां अचल रहती हैं, जिससे हाइब्रिड पेरोव्स्काइट्स में आयनिक परिवहन की पारंपरिक समझ में संशोधन होता है।

मूल लेखक: Viren Tyagi, Geert Brocks, Shuxia Tao

प्रकाशित 2026-05-05

📖 5 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Viren Tyagi, Geert Brocks, Shuxia Tao

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

कल्पना कीजिए कि एक सौर सेल (solar cell) बना है एक विशेष क्रिस्टल से जिसे MAPbI3 कहा जाता है। इस क्रिस्टल को पत्थर के एक कठोर ब्लॉक के रूप में नहीं, बल्कि छोटे-छोटे निर्माण खंडों (building blocks) से बने एक नरम, स्पंजीदार स्पंज के रूप में सोचें। इस स्पंज के अंदर दो मुख्य प्रकार के ब्लॉक हैं: भारी धातु के ब्लॉक (लेड और आयोडीन) और हल्के, जैविक "अणुओं" के ब्लॉक (जिसे MA कहा जाता है, जो छोटे मिथाइलअमोनियम अणुओं की तरह हैं)।

समस्या यह है कि यह स्पंज एकदम सही नहीं है। कभी-कभी ब्लॉक्स गायब हो जाते हैं (वैकेंसी/रिक्तियां), और कभी-ते अतिरिक्त ब्लॉक्स को जहाँ उनकी जगह नहीं है वहाँ जबरदस्ती ठूँस दिया जाता है (इंटरस्टिशियल/अंतरालीय)। जब ये "दोष" (defects) इधर-उधर हिलने-डुलने लगते हैं, तो वे समय के साथ सौर सेल को खराब कर सकते हैं।

लंबे समय से, वैज्ञानिक इस बात को समझने की कोशिश कर रहे थे कि ये दोष ठीक कैसे चलते हैं और कितनी तेजी से चलते हैं। उन्हें प्रयोगों से जो संख्याएँ मिलीं, वे बहुत ही अजीब थीं, जैसे कि लोगों का एक समूह किसी कार की गति का अनुमान लगा रहा हो और उनके जवाब "पैदल चलने की गति" से लेकर "सुपरसोनिक" तक बदल रहे हों।

यह शोध पत्र एक अति-स्मार्ट कंप्यूटर सिमुलेशन (जो आर्टिफिशियल इंटेलिजेंस द्वारा संचालित है) का उपयोग करता है ताकि इन दोषों की गति को वास्तविक समय में देखा जा सके, जैसे कि एक हाई-स्पीड कैमरा डांस फ्लोर पर होने वाले नृत्य को फिल्मा रहा हो। उन्होंने क्या पाया, इसे यहाँ सरल भाषा में समझाया गया है:

1. "भूत" बनाम "भारी ढोने वाला" (The "Ghost" vs. The "Heavy Hauler")

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

2. आश्चर्यजनक नर्तक (The Surprising Dancer - द MA मॉलिक्यूल)

इस शोध पत्र का सबसे बड़ा आश्चर्य MA अणुओं से जुड़ा है। ये आयोडीन परमाणुओं की तुलना में बहुत बड़े और भारी होते हैं। आप उम्मीद कर सकते हैं कि ये बड़े अणु सुस्त, भारी और हिलने में कठिन होंगे—जैसे किसी कमरे में पियानो को धकेलने की कोशिश करना।

पुरानी धारणा: वैज्ञानिक सोचते थे कि ये बड़े अणु फंसे हुए हैं या बहुत धीरे चलते हैं।
नई खोज: सिमुलेशन ने दिखाया कि MA इंटरस्टिशियल (अतिरिक्त अणु) वास्तवं आयोडीन के भूतों जितने ही तेज़ हैं!

यह कैसे संभव है?
शोध पत्र बताता है कि ये बड़े अणु अकेले नहीं चलते। वे एक "ग्रुप हग" (समूह गलबंदी) में चलते हैं। कल्पना कीजिए कि डांस फ्लोर पर तीन लोग हैं। एक व्यक्ति के अकेले रास्ता बनाने की कोशिश करने के बजाय, वे सभी एक समन्वित (coordinated) गति में एक साथ घूमते और बदलते हैं। एक व्यक्ति आगे बढ़ता है, दूसरे स्थान बनाने के लिए घूमते हैं, और अचानक पूरा समूह खिसक जाता है। यह टीमवर्क भारी MA अणुओं को आयोडीन के नन्हे परमाणुओं की तरह तेज़ दौड़ने में मदद करता है।

3. वह जो अपनी जगह पर टिका रहता है

यहाँ एक अपवाद है: MA रिक्तियां (Vacancies) (जहाँ एक MA अणु गायब है)। सिमुलेशन ने दिखाया कि ये खाली स्थान लगभग स्थिर (immobile) हैं। भले ही सिमुलेशन में तापमान को बहुत ऊँचा कर दिया गया हो, ये खाली स्थान हिले भी नहीं। ऐसा लगता है जैसे वह खाली जगह फर्श से चिपकी हुई है। इससे पता चलता है कि यदि आप सौर सेल में MA को चलते हुए देखते हैं, तो वह संभवतः अतिरिक्त अणु हैं, न कि खाली स्थान।

4. संख्याएँ भ्रमित करने वाली क्यों थीं?

शोध पत्र सुझाव देता है कि पिछले प्रयोगों में अलग-अलग उत्तर मिलने का कारण यह था कि वे अलग-अलग चीजों को माप रहे थे।

तेज़ गति (0.15–0.20 eV ऊर्जा अवरोध) वह है जो क्रिस्टल के गहराई में होती है (बल्क डिफ्यूजन), जिस पर यह अध्ययन केंद्रित था।
अन्य अध्ययनों में रिपोर्ट की गई धीमी गति क्रिस्टल के कणों के किनारों पर या उनके बीच की सीमाओं पर हो सकती है, जहाँ चीजें फंस जाती हैं और अलग तरह से चलती हैं।

बड़ी तस्वीर (The Big Picture)

यह अध्ययन इन सामग्रियों के बारे में हमारी समझ के नियमकोश (rulebook) को फिर से लिखता है। यह हमें बताता है कि:

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

यह समझकर कि ये दोष इतने गतिशील हैं और विशिष्ट तरीकों से चलते हैं, वैज्ञानिक अब इन दोषों को "पैसिवेट" (बंद करने) या उन्हें हिलने से रोकने के तरीके बेहतर ढंग से डिजाइन कर सकते हैं, जिससे सौर सेल और लाइटों को लंबे समय तक चलने में मदद मिलेगी।

तकनीकी सारांश: MAPbI $_3$ में धनायन और ऋणायन प्रवास का एक एकीकृत सूक्ष्म चित्र

समस्या विवरण
हैलराइड पेरोव्स्काइट उपकरणों के जीवनकाल को बढ़ाने के लिए दोष निष्क्रियता (defect passivation) और दोष गतिशीलता (defect mobility) पर नियंत्रण अत्यंत महत्वपूर्ण है। यद्यपि प्रोटोटाइप सामग्री MAPbI $_3$ (मिथाइलअमोनियम लेड आयोडाइड) में आयन परिवहन के अध्ययन के लिए प्रयोगात्मक मापों का व्यापक रूप से उपयोग किया गया है, फिर भी रिपोर्ट किए गए माइग्रेशन बैरियर्स (migration barriers) के परिणामों में एक महत्वपूर्ण प्रसार (0.10 से 0.94 eV के बीच) दिखाई देता है। इसके अलावा, प्रवास करने वाले विशिष्ट दोषों की प्रकृति और संरचना—चाहे वे रिक्तियां (vacancies), अंतराकाशी (interstitials), या विस्तारित दोष (extended defects) हों—को प्रयोगात्मक रूप से पहचानना कठिन बना हुआ है। पारंपरिक कम्प्यूटेशनल दृष्टिकोण, जो नज़ल एस्टिक इलास्टिक बैंड (Nudged Elastic Band - NEB) विधि पर निर्भर करते हैं, इन 'सॉफ्ट' सामग्रियों में संघर्ष करते हैं क्योंकि यहाँ कई प्रतिस्पर्धी निम्न-ऊर्जा पथ मौजूद होते हैं और मिथाइलअमोनियम (MA) धनायन की समरूपता (symmetry) बहुत कम होती है। इसके अतिरिक्त, मानक मॉलिक्यूलर डायनेमिक्स (MD) सिमुलेशन अक्सर क्लासिकल फोर्स फील्ड्स की सीमाओं के कारण विश्वसनीय बैरियर निष्कर्षण के लिए आवश्यक सटीकता की कमी प्रदर्शित करते हैं।

कार्यप्रणाली (Methodology)
इन चुनौतियोंों को संबोधित करने के लिए, लेखकों ने डेंसिटी फंक्शनल थ्योरी (DFT) डेटा पर प्रशिक्षित मशीन लर्निंग फोर्स फील्ड्स (MLFFs) द्वारा संचालित मॉलिक्यूलर डायनेमिक्स (MD) सिमुलेशन का उपयोग किया।

फोर्स फील्ड विकास: अध्ययन में Allegro मॉडल का उपयोग किया गया, जो एक स्थानीय E(3)-इक्विवेरिएंट ग्राफ न्यूरल नेटवर्क पोटेंशियल (NNP) है, जिसे NequIP पैकेज के माध्यम से कार्यान्वित किया गया है।
प्रशिक्षण डेटा: NNPs को VASP पैकेज के भीतर SCAN+rVV10 फंक्शनल (वैन डेर वाल्स इंटरैक्शन सहित) का उपयोग करके उत्पन्न DFT डेटा पर प्रशिक्षित किया गया था। डेटासेट्स का निर्माण लघु-समय MD रन से "ऑन-द-फ्लाई" सैंपलिंग दृष्टिकोण का उपयोग करके किया गया था, जिससे विविध संरचनाओं का समावेश सुनिश्चित हो सके।
दोष मॉडलिंग: सिमुलेशन में 2 $\times$ 2 $\times$ 2 सुपरसेल्स (96 परमाणु) में पॉइंट डिफेक्ट्स और उत्पादन रन (production runs) के लिए 6 $\times$ 6 $\times$ 6 सुपरसेल्स (2592 परमाणु) का उपयोग किया गया। अध्ययन में अंतर्निहित स्थितियों के तहत सबसे प्रमुख चार्ज अवस्थाओं में आयोडाइड (I) और MA रिक्तियों और अंतराकाशी (interstitials) पर विचार किया गया: $I^-_I$ , $V^+_I$ , $I^0_{MA}$ , $I^+_{MA}$ , $V^0_{MA}$ , और $V^-_{MA}$ । न्यूट्रल और चार्ज्ड MA दोषों के लिए अलग-अलग NNPs को वैलेंस या कंडक्शन बैंड के पास फर्मी स्तरों वाले DFT डेटा का चयन करके प्रशिक्षित किया गया था।
सत्यापन (Validation): NNPs की सटीकता का सत्यापन CI-NEB के माध्यम से गणना किए गए माइग्रेशन बैरियर्स की DFT बेंचमार्क के साथ तुलना करके किया गया, जिसमें विचलन 0.11 eV से कम पाया गया। फोर्स तुलनाओं ने 0.99 का निर्धारण गुणांक ( $R^2$ ) दिखाया।
सिमुलेशन प्रोटोकॉल: 500 K से 600 K के बीच पांच तापमानों पर तीन स्वतंत्र 2-नैनोसेकंड MD सिमुलेशन किए गए। आइंस्टीन संबंध (Einstein relation) का उपयोग करके डिफ्यूजन गुणांक निकाले गए, जिसमें ट्रेसर डिफ्यूजन के बजाय विशिष्ट प्रजातियों के सभी परमाणुओं को शामिल किया गया ताकि 'किक-आउट' प्रक्रियाओं को समझा जा सके।

मुख्य परिणाम

माइग्रेशन बैरियर्स और गतिशीलता: सिमुलेशन से पता चला कि आयोडाइड अंतराकाशी ( $I^-_I$ ) और रिक्तियां ( $V^+_I$ ) अत्यधिक गतिशील हैं, जिनकी सक्रियता ऊर्जा ( $E_a$ ) क्रमशः 0.18 eV और 0.15 eV है। उल्लेखनीय रूप से, MA अंतराकाशी ( $I^0_{MA}$ और $I^+_{MA}$ ) भी 0.18 eV और 0.20 eV की सक्रियता ऊर्जा के साथ तुलनीय गतिशीलता प्रदर्शित करते हैं।
MA रिक्ति की अचल प्रकृति: अंतराकाशी के विपरीत, MA रिक्तियों ( $V^0_{MA}$ और $V^-_{MA}$ ) को सिम्युलेटेड तापमान और समय सीमा (600 K तक) के भीतर अचल पाया गया। यह उनके सक्रियता ऊर्जा के लिए लगभग 0.4 eV के निचले स्तर (lower bound) का सुझाव देता है।
माइग्रेशन तंत्र:
- आयोडाइड दोष: प्रवास एक "किक-आउट" तंत्र के माध्यम से होता है जहाँ एक पड़ोसी आयोडाइड परमाणु रिक्ति को भरने के लिए घूमता है, या एक अंतराकाशी Pb-I-Pb ब्रिज के बीच हॉप करता है।
- MA अंतराकाशी: अपने आणविक स्वरूप के बावजूद, MA अंतराकाशी तीन पड़ोसी MA आयनों के शामिल एक समन्वित गति (concerted motion) के माध्यम से प्रवास करते हैं। एक Pb-I पिंजरे (cage) में दो MA आयन और तीसरा पड़ोसी अपने C-N अक्षों को संरेखित करने के लिए घूमते हैं, जिससे केंद्रीय आयन को बगल के पिंजरे में जाने की अनुमति मिलती है।
चार्ज स्टेट निर्भरता: आयोडाइड से संबंधित दोषों का डिफ्यूजन चार्ज अवस्था पर निर्भर करता है, जबकि MA से संबंधित दोषों का डिफ्यूजन चार्ज अवस्था से लगभग स्वतंत्र है। यह इस तथ्य के कारण है कि MA अंतराकाशी और रिक्तियां बहुत उथले (shallow) इम्प्योरिटी लेवल बनाती हैं जो कमरे के तापमान पर आयनित हो जाते हैं, जिससे विसरित चार्ज (delocalized charge) गति में बाधा नहीं डालता है।
प्रयोगों के साथ निरंतरता: गणना किए गए सक्रियता बैरियर्स (0.15–0.20 eV) प्रयोगात्मक रूप से रिपोर्ट किए गए विस्तृत रेंज के निचले स्तर के अनुरूप हैं।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह MAPbI $_3$ में आयन प्रवास का एक एकीकृत सूक्ष्म चित्र प्रदान करता है जो प्रयोगात्मक सक्रियता ऊर्जा के विसंगतियों को सुलझाता है। लेखक तर्क देते हैं कि प्रयोगों में देखे गए निम्न मान (0.15–0.20 eV) आयोडाइड दोषों और MA अंतराकाशी द्वारा संचालित अंतर्निहित बल्क डिफ्यूजन के अनुरूप हैं। अन्य अध्ययनों में रिपोर्ट किए गए उच्च मान संभवतः जटिल दोषों के निर्माण या इंटर-ग्रेन ट्रांसपोर्ट जैसे अतिरिक्त योगदानों को दर्शाते हैं, न कि अंतर्निहित बल्क माइग्रेशन बैरियर्स को।

एक प्राथमिक योगदान पारंपरिक समझ को संशोधित करना है कि MA अपनी आणविक संरचना के कारण स्वाभाविक रूप से धीमा होता है। अध्ययन यह प्रदर्शित करता है कि MA अंतराकाशी, एक सामूहिक, समन्वित माइग्रेशन तंत्र के कारण आयोडाइड दोषों के लगभग समान गतिशील हैं। इसके अलावा, यह कार्य इस बात पर प्रकाश डालता है कि A-कैटायन प्रजाति (MA) हैलाइड गति को महत्वपूर्ण रूप से प्रभावित नहीं करती है, क्योंकि MAPbI $_3$ में हैलाइड दोषों की गतिशीलता सभी-अकार्बनिक CsPb(I $_x$ Br $_{1-x}$ ) $_3$ पेरोव्स्केट्स के समान है। अंततः, यह अध्ययन हाइब्रिड पेरोव्स्केट्स में तेज़ आयनिक प्रवास को सक्षम करने के लिए सामूहिक आणविक गति के महत्व को रेखांकित करता है।