🔬 materials science

Melting Behavior and Phase Stability of CaO from Neural Network Potentials: a Molecular Dynamics Study

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

मूल लेखक: Francesca Menescardi, Stefano de Gironcoli

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

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

CC BY 4.0

मूल लेखक: Francesca Menescardi, Stefano de Gironcoli

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

कल्पना कीजिए कि आप यह पता लगाने की कोशिश कर रहे हैं कि कैल्शियम ऑक्साइड (CaO) नामक पत्थर जैसे सख्त चीज़ का चट्टान जैसा ठोस रूप ठीक कब एक चिपचिपे तरल में बदल जाता है। वैज्ञानिक लंबे समय से इस तापमान को लेकर बहस कर रहे हैं। कुछ कहते हैं कि यह 2,800 डिग्री के आसपास है, अन्य कहते हैं कि यह 3,200 से ऊपर है। समस्या यह है कि CaO इतना गर्म और प्रतिक्रियाशील है कि यह ऐसा है जैसे आप किसी धातु के टुकड़े को पिघलाने की कोशिश कर रहे हों जबकि वह उसी बर्तन को खाने की भी कोशिश कर रहा हो जिसमें वह रखा है। एक वास्तविक लैब में सटीक माप लेना कठिन है।

इस समस्या को हल करने के लिए, शोधकर्ताओं ने CaO का एक डिजिटल ट्विन (Digital Twin) बनाया। असली चट्टानों को पिघलाने के बजाय, उन्होंने एक "स्मार्ट कंप्यूटर ब्रेन" (जिसे मशीन लर्निंग इंटरएटॉमिक पोटेंशियल कहा जाता है) बनाया जो जानता है कि CaO में हर एक परमाणु बिल्कुल कैसे व्यवहार करता है। इस 'ब्रेन' को एक अत्यंत सटीक नियम पुस्तिका की तरह समझें जो कंप्यूटर को बताती है कि परमाणु एक-दूसरे पर कैसे दबाव डालते हैं और कैसे खिंचाव महसूस करते हैं, लेकिन यह पुराने, धीमे भौतिकी सिमुलेशन की तुलना में दस लाख गुना तेज़ चलता है।

उन्होंने इस डिजिटल ब्रेन का उपयोग उत्तर खोजने के लिए इस प्रकार किया:

1. डिजिटल चट्टान को पिघलाने के दो तरीके

पिघलने के सटीक बिंदु को खोजने के लिए, उन्होंने अपने सिमुलेशन में दो अलग-अलग "खेल" आजमाए:

"दीवार में छेद" विधि (Void-Nucleated Melting):
कल्पना कीजिए कि एक आदर्श ईंट की दीवार है। यदि आप इसे गर्म करते हैं, तो यह अपने पिघलने के बिंदु के बहुत आगे तक ठोस रह सकती है क्योंकि वहां ढहने की शुरुआत करने के लिए कोई दरार नहीं है। इसे ठीक करने के लिए, शोधकर्ताओं ने अपनी डिजिटल दीवार के बीच में एक छेद किया। यह छेद एक कमजोर बिंदु के रूप में कार्य करता है। जैसे-जैसे उन्होंने दीवार को गर्म किया, छेद के चारों ओर तरल बनने लगा। छेद को बड़ा और बड़ा करके, उन्होंने वह तापमान खोज निकाला जहाँ दीवार हमेशा ढह जाती है।
- परिणाम: उन्होंने पाया कि पिघलने का बिंदु 3,055 केल्विन (लगभग 2,782°C) है। यह हाल के सबसे अच्छे वास्तविक-दुनिया के प्रयोगों से मेल खाता है।
"आधा-आधा" विधि (Two-Phase Coexistence):
कल्पना कीजिए कि एक लंबी ट्रेन का डिब्बा है जिसका अगला आधा हिस्सा जमी हुई बर्फ है और पिछला आधा हिस्सा उबलते पानी जैसा है। उन्होंने इस ट्रेन के डिब्बे को सिमुलेशन में डाला और बर्फ और पानी के बीच की सीमा को देखा। यदि बर्फ पिघलती है, तो पूरी चीज़ बहुत गर्म है। यदि पानी जमता है, तो यह बहुत ठंडा है। उन्होंने तापमान को तब तक समायोजित किया जब तक कि बर्फ और पानी पूरी तरह से संतुलित न हो गए।
- परिणाम: इस विधि ने एक कम संख्या दी, 2,847 केल्विन। पेपर में उल्लेख है कि यह विधि कभी-कभी तापमान को कम आंक देती है, लेकिन फिर भी यह एक उपयोगी जांच है।

2. "हीट बिल" (Enthalpy) की जांच करना

पिघलना केवल तापमान के बारे में नहीं है; यह इस बारे में भी है कि ठोस संरचना को तोड़ने के लिए आपको सिस्टम में कितनी ऊर्जा डालनी पड़ती है। शोधकर्ताओं ने इस "ऊर्जा बिल" (एन्थैल्पी ऑफ फ्यूजन) की गणना की।

उन्होंने पाया कि उनके डिजिटल ब्रेन ने लगभग 73 kJ/mol की ऊर्जा लागत का अनुमान लगाया।
यह संख्या वास्तविक दुनिया की केमिस्ट्री टेबल्स और अन्य उच्च-स्तरीय भौतिकी गणनाओं के सर्वोत्तम अनुमानों के साथ पूरी तरह मेल खाती है। इसने साबित कर दिया कि उनका डिजिटल ब्रेन सच बोल रहा है।

3. "दबाव" परीक्षण (High Pressure)

अंत में, उन्होंने पूछा: "क्या होता है अगर हम इस चट्टान को कुचल दें?" उन्होंने अपने डिजिटल CaO को 20 गीगापास्कल तक दबाया (यह समुद्र के तल के दबाव जैसा है, लेकिन हजार गुना अधिक)।

पुरानी धारणा: वैज्ञानिक पहले यह मानते थे कि जैसे-जैसे आप किसी सामग्री को दबाते हैं, "ओवरहीटिंग" (एक आदर्श क्रिस्टल को पिघलने के लिए आवश्यक अतिरिक्त गर्मी) का प्रतिशत समान रहता है।
नई खोज: शोधकर्ताओं ने पाया कि यह धारणा गलत थी। जैसे-जैसे उन्होंने CaO को अधिक दबाया, "ओवरहीटिंग" का अंतर वास्तव में बढ़ गया। सामान्य दबाव पर, एक आदर्श क्रिस्टल को पिघलने के लिए लगभग 17% अतिरिक्त गर्मी की आवश्यकता होती है। उच्च दबाव पर, उसे 24% अतिरिक्त गर्मी की आवश्यकता होती है।
क्यों? इसे एक भीड़ भरे डांस फ्लोर की तरह समझें। जब कम दबाव (कम भीड़) होता है, तो कुछ डांसर्स का हिलना-डुलना शुरू करना आसान होता है (पिघलना)। लेकिन जब बहुत अधिक दबाव (भारी भीड़) होता, तो भीड़ को अपना फॉर्मेशन तोड़ने और नाचना शुरू करने के लिए भारी मात्रा में ऊर्जा की आवश्यकता होती है, खासकर यदि वहां मदद के लिए कोई "कमजोर बिंदु" (दोष) न हो।

निष्कर्ष

इस पेपर ने केवल कैल्शियम ऑक्साइड के पिघलने के बिंदु का अनुमान नहीं लगाया; उन्होंने इसे साबित करने के लिए एक अत्यधिक सटीक, तेज़ कंप्यूटर मॉडल बनाया। उन्होंने पुष्टि की कि सामान्य दबाव पर CaO लगभग 3,055 K पर पिघलता है और दिखाया कि जब आप इसे दबाते हैं तो इसके पिघलने के नियम बदल जाते हैं। उनका नया "डिजिटल ब्रेन" अब वैज्ञानिकों के लिए अन्य चरम सामग्रियों का अध्ययन करने के लिए एक विश्वसनीय उपकरण है, बिना उन्हें वास्तविक लैब में पिघलाए।

तकनीकी सारांश: न्यूरल नेटवर्क पोटेंशियल से CaO का पिघलने का व्यवहार और चरण स्थिरता

समस्या विवरण
कैल्शियम ऑक्साइड (CaO) के सटीक गलनांक ( $T_m$ ) और उच्च-दबाव चरण स्थिरता को निर्धारित करना एक महत्वपूर्ण चुनौती बनी हुई है, क्योंकि उच्च वाष्प दबाव और प्रतिक्रियाशीलता जैसी प्रयोगात्मक सीमाएँ मौजूद हैं। $T_m$ के प्रयोगात्मक मान व्यापक रूप से भिन्न होते हैं (2843 K से 3223 K तक), और उच्च-दबाव गलनांक वक्र (melting curve) काफी हद तक अनछुआ है। पिछले कम्प्यूटेशनल अध्ययन एक समझौते (trade-off) का सामना कर रहे थे: अनुभवजन्य पोटेंशियल (जैसे, बोर्न-मेयर-हगिंस) में चरम स्थितियों के लिए आवश्यक क्वांटम सटीकता और स्थानांतरणीयता (transferability) का अभाव था, जबकि एब इनिटियो मॉलिक्यूलर डायनेमिक्स (AIMD) परिशुद्धता तो प्रदान करता था, लेकिन बड़े सिस्टम साइज और लंबी सिमुलेशन अवधि के लिए इसकी कम्प्यूटेशनल लागत अत्यधिक थी, जो गलनांक वक्रों को सटीक रूप से मैप करने और सुपरहीटिंग आर्टिफैक्ट्स से बचने के लिए आवश्यक थी।

कार्यप्रणाली
इस अंतर को पाटने के लिए, लेखकों ने PANNA 2.0 फ्रेमवर्क और LATTE डिस्क्रिप्टर का उपयोग करके CaO के लिए एक मशीन लर्निंग इंटरएटोमिक पोटेंशियल (MLIP) विकसित किया।

पोटेंशियल आर्किटेक्चर: मॉडल एक मल्टीलेयर परसेप्ट्रॉन (MLP) का उपयोग करता है जिसकी 256:124:1 हिडन लेयर संरचना है। यह अनुवाद, रोटेशन और परम्यूटेशन के प्रति अपरिवर्तनीय स्थानीय परमाणु परिवेश को एनकोड करने के लिए कार्टेशियन टेंसर कॉन्ट्रैक्शंस पर आधारित LATTE डिस्क्रिप्टर का उपयोग करता है। इस डिस्क्रिप्टर में 2-बॉडी, 3-बॉडी और 4-बॉडी टर्म शामिल हैं, जिनमें कुल 900 फीचर्स हैं।
प्रशिक्षण डेटासेट: पोटेंशियल को लगभग 12,000 कॉन्फ़िगरेशन पर प्रशिक्षित किया गया था जो AIMD सिमुलेशन (PBEsol फंक्शनल और Quantum ESPRESSO का उपयोग करके) से प्राप्त किए गए थे। डेटासेट विविध है, जिसमें शामिल हैं:
- ठोस और तरल बल्क संरचनाएं (300–6000 K)।
- ठोस-तरल इंटरफेस संरचनाएं (दो-चरण सह-अस्तित्व के लिए)।
- "वॉयडेड" (voided) ठोस संरचनाएं जिनमें केंद्रीय छिद्र हैं (1–40% आयतन के साथ) जो वॉयड-न्यूक्लेटेड मेल्टिंग के लिए हैं।
- संकुचित और विस्तारित संरचनाएं (20 GPa तक)।
  डेटा को प्रशिक्षण (75%), सत्यापन (15%), और परीक्षण (10%) सेटों में विभाजित किया गया था।
प्रशिक्षण प्रदर्शन: अंतिम पोटेंशियल ने ऊर्जा के लिए माध्य पूर्ण त्रुटि (MAE) <5 meV/atom और बलों के लिए <45 meV/Å प्राप्त की।
सिमुलेशन तकनीकें: प्रशिक्षित MLIP के साथ LAMDASS में आइसोबारिक (NPT) एन्सेम्बल का उपयोग करके बड़े पैमाने पर मॉलिक्यूलर डायनेमिक्स (MD) सिमुलेशन किए गए। $T_m$ $T_{m}$ निर्धारित करने के लिए दो विधियों का उपयोग किया गया:
1. वॉयड-न्यूक्लेटेड मेल्टिंग (VNM): विषम न्यूक्लिएशन को प्रेरित करने और सुपरहीटिंग को समाप्त करने के लिए एक 11 $\times$ 11 $\times$ 11 सुपरसेल में एक केंद्रीय वॉयड (void) पेश करना।
2. टू-फेज को-एग्ज़िस्टेंस (TPC): एक 6 $\times$ 6 $\times$ 60 सुपरसेल में सह-अस्तित्व वाले ठोस और तरल चरणों वाले सिस्टम का अनुकरण करना ताकि संतुलन तापमान पाया जा सके।

मुख्य परिणाम

परिवेशी दबाव पर गलनांक तापमान:
- VNM विधि ने $T_m = 3055 \pm 11$ K दिया। यह मान हाल के प्रयोगात्मक डेटा और पिछले MLIP अध्ययनों के करीब है, जो सुझाव देता है कि यह परीक्षण की गई विधियों में सबसे सटीक अनुमान है।
- TPC विधि ने $T_m = 2847 \pm 15$ K दिया। लेखक नोट करते हैं कि यह कम अनुमान (underestimation) CaO और अन्य सामग्रियों के लिए TPC विधि की एक ज्ञात सीमा है।
ऊष्मागतिक गुण (Thermodynamic Properties):
- फ्यूजन की एन्थैल्पी ( $\Delta H_f$ ): गणना की गई $\sim$ 73 kJ/mol (3055 K पर) और $\sim$ 71.6 kJ/mol (2847 K पर) है। ये मान ऊष्मागतिक आकलन और AIMD गणनाओं के साथ सहमत हैं।
- फ्यूजन की एंट्रॉपी ( $\Delta S_f$ ): गणना किए गए मान ( $\sim$ 24.05–25.14 J·mol $^{-1}$ ·K $^{-1}$ ) NIST-JANAF तालिकाओं और MgO एवं BeO से प्राप्त अनुमानों के अनुरूप हैं।
- आयतन विस्तार (Volume Expansion): सिमुलेशन $T_m$ पर $\sim$ 29% आयतन वृद्धि की भविष्यवाणी करता है, जो घनत्व प्रोफाइल के अनुरूप है जो दर्शाता है कि तरल चरण ठोस चरण की तुलना में कम घना है।
- तापीय विस्तार (Thermal Expansion): MLIP तापीय विस्तार के रुझानों को पुनरुत्पादित करता है जो AIMD और प्रयोगात्मक डेटा के अनुरूप है, जो पिछले अनुभवजन्य BMH पोटेंशियल सिमुलेशन में देखे गए तीव्र ढलान को ठीक करता है।
उच्च-दबाव गलनांक वक्र (20 GPa तक):
- TPC विधि का उपयोग करते हुए, लेखकों ने 20 GPa तक गलनांक वक्र की गणना की।
- अध्ययन पुष्टि करता है कि ओवरहीटिंग अनुपात ( $\eta = T_s/T_m - 1$ ) दबाव के तहत स्थिर नहीं है। यह 0 GPa पर 17.3% से बढ़कर 20 GPa पर 24.0% हो जाता है।
- अनुभवजन्य फिट ( $T = aP^b + c$ ) दिखाते हैं कि जबकि $T_m$ और तापीय अस्थिरता तापमान ( $T_s$ ) समान स्केलिंग नियमों का पालन करते हैं, $T_s$ उच्च दबाव में पूर्ण क्रिस्टल में न्यूक्लिएशन साइटों की कमी के कारण दबाव के साथ अधिक तीव्रता से बढ़ता है।

महत्व और दावे
यह शोध पत्र MLIP-आधारित सिमुलेशन को आयनिक ऑक्साइड में चरण स्थिरता की जांच के लिए एक मजबूत और कुशल ढांचे के रूप में स्थापित करता है। अनुभवजन्य पोटेंशियल की कम्प्यूटेशनल लागत पर एब इनिटियो सटीकता प्राप्त करके, यह अध्ययन पिछले अनुभवजन्य मॉडलों (विशेष रूप से तापीय विस्तार ढलान और ओवरहीटिंग अनुपात के संबंध में) में पाई गई ऊष्मागतिक विसंगतियों को हल करता है। इसके परिणाम 20 GPa तक CaO के कुछ ही कम्प्यूटेशनल निर्धारणों में से एक के रूप में उच्च-दबाव गलनांक वक्र प्रदान करते हैं। महत्वपूर्ण रूप से, यह कार्य यह सत्यापित करता है कि दबाव के तहत $T_m$ का अनुमान लगाने के लिए एक निश्चित ओवरहीटिंग अनुपात मानना गलत है, क्योंकि तापीय अस्थिरता सीमा और वास्तविक गलनांक बिंदु के बीच का अनुपात दबाव के साथ महत्वपूर्ण रूप से बदलता है। अध्ययन पुष्टि करता है कि परिवेशी दबाव पर CaO के लिए गलनांक तापमान 3000 K से ऊपर है।

1. डिजिटल चट्टान को पिघलाने के दो तरीके

2. "हीट बिल" (Enthalpy) की जांच करना

3. "दबाव" परीक्षण (High Pressure)

निष्कर्ष

इसी तरह के और पेपर