Charting the thermodynamic stability of hybrid perovskite alloys with machine learning

तकनीकी सारांश: मशीन लर्निंग के साथ हाइब्रिड पेरोव्स्काइट मिश्र धातुओं की ऊष्मप्रवैगिकी स्थिरता का मानचित्रण

समस्या विवरण
हाइब्रिड पेरोव्स्काइट सौर सेल (PSCs) ने सिलिकॉन के बराबर बिजली रूपांतरण दक्षता हासिल की है, फिर भी इनका व्यावसायीकरण दीर्घकालिक परिचालन स्थिरता संबंधी समस्याओं और लेड (सीसा) विषाक्तता की चिंताओं के कारण बाधित है। चतुर्धातुक (quaternary) मिश्र धातुओं, विशेष रूप से $(Cs/FA)Pb(Br/I)_3$ और लेड-मुक्त $(Cs/FA)Sn(Br/I)_3$ (जहाँ FA फॉर्ममिडिनियम है) के माध्यम से संरचनात्मक इंजीनियरिंग, गुणों को ट्यून करने और स्थिरता में सुधार करने का एक मार्ग प्रदान करती है। हालाँकि, इन बहु-घटक प्रणालियों के विशाल विन्यास स्थान (configurational space) के साथ-साथ उन्मुखी योग्य (orientable) कार्बनिक धनायनों की जटिलता ने सटीक ऊष्मप्रवैगिकी मॉडलिंग को बाधित किया है। पारंपरिक डेंसिटी फंक्शनल थ्योरी (DFT) आवश्यक विन्यासों के नमूने लेने के लिए गणनात्मक रूप से अत्यधिक कठिन है, जबकि क्लस्टर एक्सपेंशन जैसे मौजूदा सन्निकटन (approximations) आणविक धनायन के ओरिएंटेशन को पकड़ने में विफल रहते हैं, और स्पेशल क्वासिरैंडम स्ट्रक्चर्स (SQS) सीमित विन्यास नमूनाकरण के कारण पूर्वाग्रह (biases) उत्पन्न करते हैं।

कार्यप्रणाली
लेखक एक मशीन लर्निंग (ML) त्वरित परमाणु मॉडल (atomistic modeling) वर्कफ़्लो प्रस्तुत करते हैं, जिसे इन चतुर्धातुक मिश्र धातुओं के पूर्ण संरचनात्मक स्थान में मुक्त ऊर्जा परिदृश्य (free energy landscapes) की गणना करने के लिए डिज़ाइन किया गया है। यह दृष्टिकोण कम्प्यूटेशनल बाधाओं को दूर करने के लिए दो-स्तरीय ML रणनीति का उपयोग करता है:

1. MACE इंटरएटोमिक पोटेंशियल: लेखकों ने DFT डेटा पर MACE (हायर-ऑर्डर इक्विवेरिएंट मैसेज पासिंग न्यूरल नेटवर्क) पोटेंशियल को प्रशिक्षित किया। इन पोटेंशियल का उपयोग कुशल संरचना विश्राम (structure relaxations) के लिए किया गया, जिसने DFT की तुलना में गणनाओं को लगभग छह आदेशों की मात्रा (six orders of magnitude) से तेज कर दिया। प्रशिक्षण डेटा में चार अलग-अलग चरणों ($Pm\bar{3}m$, $P4/mbm$, $I4/mcm$, और$Pnma$) से$2 \times 2 \times 2$ सुपरसेल्स और FA आणविक ओरिएंटेशन शामिल थे।
2. डायरेक्ट रिलैक्सेशन मॉडल (KRR): मुक्त ऊर्जा गणनाओं के लिए आवश्यक नमूनाकरण को और अधिक तेज करने के लिए, माध्यमिक कर्नेल रिज रिग्रेशन (KRR) मॉडल को MACE-रिलैक्स्ड डेटा पर प्रशिक्षित किया गया। ये मॉडल अनरिलैक्स्ड परमाणु संरचनाओं से सीधे रिलैक्स्ड ऊर्जा की भविष्यवाणी करते हैं, जिससे नमूनाकरण के दौरान स्पष्ट विश्राम (explicit relaxation) की आवश्यकता समाप्त हो जाती है। इस चरण ने अतिरिक्त तीन आदेशों की मात्रा का त्वरण प्रदान किया।
3. एक्टिव लर्निंग: सटीकता सुनिश्चित करने के लिए, विशेष रूप से स्थिरता विश्लेषण के लिए महत्वपूर्ण निम्न-ऊर्जा शासन (low-energy regime) में, एक एक्टिव लर्निंग योजना का उपयोग किया गया। इसमें निम्न-ऊर्जा मिश्र धातु विन्यासों की पहचान और सुधार के लिए पुनरावृत्त ऊर्जा न्यूनीकरण मोंटे कार्लो (MC) सिमुलेशन शामिल थे, जिससे माध्य पूर्ण त्रुटियों (MAEs) को काफी कम किया गया।
4. मुक्त ऊर्जा गणना: वांग-लैंडौ (Wang-Landau) एल्गोरिदम का उपयोग स्टेट्स के घनत्व (density of states) को नमूना लेने और 300 K (और संश्लेषण स्थितियों के लिए 150 °C) पर हेल्महोल्ट्ज़ मुक्त ऊर्जा के मिश्रण की गणना करने के लिए किया गया था। इसने फेज़ सेपरेशन (phase separation) के विरुद्ध ऊष्मप्रवैगिकी स्थिरता का आकलन करने के लिए कॉनवेक्स हल (convex hulls) और वक्रता विश्लेषण का निर्माण करने की अनुमति दी।

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

• वर्कफ़्लो विकास: यह शोध पत्र एक कठोर बहु-चरणीय ML वर्कफ़्लो स्थापित करता है जो उन्मुखी योग्य कार्बनिक धनायनों के उपचार सहित चतुर्धातुक हाइब्रिड पेरोव्स्काइट्स की जटिलता को संभालने में सक्षम है, जिसे पिछले तरीकों ने प्रभावी ढंग से मॉडल करने में संघर्ष किया था।
• त्वरण (Speedup): संयुक्त दृष्टिकोण पारंपरिक DFT की तुलना में लगभग नौ आदेशों की मात्रा का कुल त्वरण प्राप्त करता है, जिससे व्यापक मुक्त ऊर्जा मानचित्रण के लिए आवश्यक दो बिलियन से अधिक ऊर्जा विन्यासों का मूल्यांकन संभव हुआ।
• तुलनात्मक विश्लेषण: यह अध्ययन कई चरणों में लेड-आधारित $(Cs/FA)Pb(Br/I)_3$ और लेड-मुक्त $(Cs/FA)Sn(Br/I)_3$ चतुर्धातुक मिश्र धातुओं के बीच पहला व्यापक ऊष्मप्रवैगिक स्थिरता तुलना प्रदान करता है।

परिणाम

• मॉडल सटीकता: MACE पोटेंशियल ने रिलैक्स्ड संरचनाओं पर सिंगल-पॉइंट DFT गणनाओं की तुलना में Pb- और Sn-आधारित प्रणालियों के लिए क्रमशः 2.81 meV/f.u. और 3.07 meV/f.u. की MAE प्राप्त की। अंतिम डायरेक्ट रिलैक्सेशन मॉडल (एक्टिव लर्निंग के बाद) ने न्यूनतम-ऊर्जा विन्यासों पर लगभग 3.2 meV/f.u. की MAE प्राप्त की।
• स्थिरता परिदृश्य (Stability Landscapes):
  • Pb-आधारित प्रणाली: व्यापक स्थिर संरचना क्षेत्र प्रदर्शित करती है। क्यूबिक ($Pm\bar{3}m$) चरण में, संरचना स्थान के केंद्र में एक निम्न-ऊर्जा क्षेत्र मौजूद है। कॉनवेक्स हल पर स्थिर क्षेत्र विभिन्न A-साइट प्रतिस्थापनों में और Cs/FA अनुपात के आधार पर 60% Br सांद्रता तक विस्तृत हैं।
  • Sn-आधारित प्रणाली: अपने Pb समकक्ष की तुलना में संकीर्ण स्थिर संरचना क्षेत्र दिखाती है। Pb-आधारित क्यूबिक चरण में देखा गया केंद्रीय निम्न-ऊर्जा वेल (well) अनुपस्थित है; निम्नतम ऊर्जाएँ Cs-मुक्त किनारे पर पाई जाती हैं। कॉनवेक्स हल इंगित करता है कि संरचनात्मक इंजीनियरिंग के विकल्प कम हैं, और स्थिरता मुख्य रूप से उच्च आयोडीन सामग्री तक सीमित है।
  • तापमान प्रभाव: उन्नत तापमान (150 °C) पर स्थिरता क्षेत्र विस्तारित होते हैं, जो एंट्रोपिक स्थिरीकरण (entropic stabilization) के अनुरूप है, लेकिन Pb और Sn प्रणालियों के बीच मौलिक अंतर बना रहता है।
• एंट्रॉपी और स्थिरता: विश्लेषण से पता चलता है कि एंट्रॉपी-संचालित स्थिरीकरण सीमित है। अधिकतम स्थिरता उच्च आयोडीन सामग्री पर होती है, और जहाँ मिश्रण एंट्रॉपी उच्चतम है, वहाँ मध्य चतुर्धातुक संरचनाओं के पास कोई महत्वपूर्ण स्थिरीकरण नहीं देखा गया है। यह सुझाव देता है कि मिश्रण अकेले फेज़ सेपरेशन के विरुद्ध इन केंद्रीय चतुर्धातुक संरचनाओं को स्थिर नहीं कर सकता है।
• प्रायोगिक सत्यापन: Pb-आधारित प्रणाली के परिणाम वांग एट अल द्वारा प्रदान किए गए उच्च-थ्रूपुट प्रायोगिक डेटा के साथ गुणात्मक सहमति दिखाते हैं, विशेष रूप से क्यूबिक चरण में कम-Br, कम-Cs क्षेत्रों की स्थिरता के संबंध में। ऑर्थोरोम्बिक चरण में विसंगतियाँ गैर-संतुलन संश्लेषण स्थितियों और समान क्रिस्टल संरचनाओं के बीच संभावित फेज़ पहचान चुनौतियों के कारण हैं।

महत्व और दावे
यह शोध पत्र दावा करता है कि इसके परिणाम इन मिश्र धातुओं के मौलिक ऊष्मप्रवैगिक सीमाओं का मानचित्रण करके स्थिर पेरोव्स्काइट्स के डिजाइन का मार्गदर्शन करते हैं। लेखक निष्कर्ष निकालते हैं कि टिन-आधारित हाइब्रिड पेरोव्स्काइट्स की कम ऊष्मप्रवैगिक स्थिरता उनके लेड-आधारित समकक्षों की तुलना में उनके मुक्त ऊर्जा परिदृश्य में मौलिक अंतरों के कारण है, विशेष रूप से केंद्रीय निम्न-ऊर्जा क्षेत्रों की अनुपस्थिति और संकीर्ण स्थिर संरचना रेंज के कारण। फलस्वरूप, यह अध्ययन सुझाव देता है कि $(Cs/FA)Sn(Br/I)_3$ को अनुकूलित करने के लिए भविष्य के प्रायोगिक प्रयास उन संरचनाओं पर ध्यान केंद्रित करने चाहिए जिन्हें पहले ही Pb-आधारित प्रणाली के लिए अनुकूल माना गया है, न कि उस केंद्रीय संरचना स्थान पर जहाँ स्थिरता ऊष्मप्रवैगिक रूप से अनुकूल नहीं है। यह कार्य इस बात पर प्रकाश डालता है कि जबकि ML खोज को तेज कर सकता है, कुछ चतुर्धातुक क्षेत्रों की अंतर्निहित ऊष्मप्रवैगिक अस्थिरता को केवल मिश्रण एंट्रॉपी द्वारा दूर नहीं किया जा सकता है।