An Efficient Method for Gibbs Free Energy Evaluation under Volume Compression

तकनीकी सारांश: आयतन संपीड़न के तहत गिब्स मुक्त ऊर्जा मूल्यांकन के लिए एक कुशल विधि

समस्या विवरण
गिब्स मुक्त ऊर्जा, $G(T, P)$, का सटीक मूल्यांकन, दबाव-तापमान चरण आरेखों (phase diagrams) के निर्माण और परिमित-तापमान एवं परिमित-दबाव स्थितियों में सामग्रियों की ऊष्मागतिक स्थिरता के सत्यापन के लिए अत्यंत महत्वपूर्ण है। जबकि शून्य-तापमान ($0$ K) संरचनात्मक विश्राम (relaxation) और ऊर्जा स्क्रीनिंग को उच्च-थ्रूपुट विधियों और मशीन लर्निंग द्वारा त्वरित किया गया है, उसके बाद का स्थिरता सत्यापन एक बाधा बना हुआ है। पारंपरिक दृष्टिकोण अर्ध-हार्मोनिक सन्निकटन (Quasi-Harmonic Approximation - QHA) पर निर्भर करते हैं, जिसके लिए मुक्त ऊर्जा सतह के निर्माण हेतु आयतन के घने ग्रिड (आमतौर पर 10–20 बिंदु) पर फोनन स्पेक्ट्रा की गणना करने की आवश्यकता होती है। संरचनात्मक रूप से जटिल सामग्रियों के लिए, जैसे कि $Ta_2O_5$ के पॉलीमॉर्फ्स जिनमें 77 परमाणुओं तक के प्रिमिटिव सेल होते हैं, यह घना आयतन नमूनाकरण (volume sampling) गणनात्मक रूप से अत्यधिक कठिन हो जाता है, जिससे उम्मीदवार पीढ़ी (candidate generation) से कठोर चरण-स्थिरता सत्यापन की ओर संक्रमण बाधित होता है।

कार्यप्रणाली
लेखक महत्वपूर्ण आयतन संपीड़न और विस्तार के तहत गिब्स मुक्त ऊर्जा का मूल्यांकन करने के लिए एक स्पार्स-वॉल्यूम ग्रुनेसेन इंटरपोलेशन (GI) विधि प्रस्तावित करते हैं। इसका मूल सिद्धांत यह है कि मध्यम आयतन सीमाओं के भीतर, फोनन आवृत्तियाँ लघुगणकीय आयतन निर्देशांक (logarithmic volume coordinate) के साथ सुचारू रूप से परिवर्तित होती हैं। यह विधि इस प्रकार कार्य करती है:

1. स्पार्स सैंपलिंग (Sparse Sampling): घने ग्रिड पर फोनन की गणना करने के बजाय, केवल स्थिर आयतनों के एक छोटे सेट ($N_{sp}$, आमतौर पर 3 बिंदु) पर स्पष्ट एब इनिशियो (ab initio) फोनन गणना की जाती है। इन बिंदुओं को लक्षित दबाव-तापमान विंडो के लिए गिब्स मुक्त ऊर्जा न्यूनीकरण द्वारा तय किए गए साम्यावस्था आयतन पथ को घेरने (bracket) के लिए चुना जाता है।
2. स्टैटिक-ZPE ब्रांच पुनर्निर्माण: स्पार्स आयतनों पर ज़ीरो-पॉइंट एनर्जी (ZPE) अनुपातों से एक प्रभावी ग्रुनेसेन पैरामीटर ($\gamma_{ZPE}$) प्राप्त किया जाता है। यह पैरामीटर एक घातांकीय स्केलिंग संबंध का उपयोग करके स्टैटिक-ZPE शाखा का पुनर्निर्माण करता है:
   $$ZPE(V) = ZPE(V_0) \exp\left[\gamma_{ZPE} \ln\left(\frac{V_0}{V}\right)\right]$$
3. परिमित-तापमान कंपन शाखा पुनर्निर्माण: अति-सरलीकरण से बचने के लिए, विधि खंडवार मोड-रिज़ॉलव्ड ग्रुनेसेन ढाल (piecewise mode-resolved Grüneisen slopes) का उपयोग करती है। दो स्पार्स आयतनों $V_a$ और $V_b$ के बीच स्थित लक्षित आयतन $V_m$ के लिए, प्रत्येक फोनन मोड $j$ की आवृत्ति को स्थानीय ढाल $\gamma_{ab}^j$ का उपयोग करके पुनर्गठित किया जाता है:
   $$\omega_j^{GI}(V_m) = \omega_j(V_a) \exp\left\{ \gamma_{ab}^j \left[ \ln\left(\frac{V_0}{V_m}\right) - \ln\left(\frac{V_0}{V_a}\right) \right] \right\}$$
   जहाँ $\gamma_{ab}^j$ की गणना अंत बिंदुओं पर आवृत्ति अनुपात से की जाती है।
4. मुक्त ऊर्जा संयोजन (Free Energy Assembly): पुनर्गठित आवृत्तियों और स्टैटिक ऊर्जाओं का उपयोग हिल्ट्ज़ मुक्त ऊर्जा (Helmholtz free energy) की गणना के लिए किया जाता है और तत्पश्चात गिब्स मुक्त ऊर्जा प्राप्त की जाती है। धात्विक प्रणालियों (जैसे Al) के लिए, थर्मल एक्सपेंशन पथ के साथ अवस्था घनत्व (DOS) के आधार पर इलेक्ट्रॉनिक उत्तेजनाओं ($\Delta F_{ele}$) के लिए एक सुधार शामिल किया जाता है।
5. सत्यापन: विधि का सत्यापन सरल तत्वों (C, Al, Si, Ge), बाइनरी ऑक्साइड (रुटाइल $TiO_2$, $\beta$-$PtO_2$) और जटिल $Ta_2O_5$ पॉलीमॉर्फ्स के लिए घने-ग्रिड QHA बेंचमार्क के विरुद्ध किया गया है।

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

• गणना दक्षता: यह विधि सबसे कम लागत वाली कार्यान्वयन में फोनन वॉल्यूम बिंदुओं की संख्या को ~20–21 से घटाकर केवल 3 कर देती है।
• सटीकता: यह गिब्स मुक्त ऊर्जाओं के पुनर्निर्माण में उच्च सटीकता प्राप्त करती है, जहाँ सरल प्रणालियों के लिए QHA बेंचमार्क के सापेक्ष माध्य पूर्ण त्रुटि (MAE) 0.53 meV/atom से नीचे रहती है।
• चरण स्थिरता संरक्षण: जटिल पॉलीमॉर्फ्स में बड़ी पूर्ण त्रुटियों के बावजूद, विधि मुख्य चरण-स्थिरता टोपोलॉजी और चरण सीमाओं को सफलतापूर्वक पुनरुत्पादित करती है।
• ऊष्मागतिक निरंतरता: पुनर्गठित मुक्त ऊर्जा सतहें तापीय विस्तार गुणांकों जैसे ऊष्मागतिक गुणों को सटीक रूप से पुनरुत्पादित करती हैं, जिसमें सिलिकॉन में नेगेटिव-टू-पॉजिटिव क्रॉसओवर व्यवहार भी शामिल है।

परिणाम

• सरल बेंचमार्क सिस्टम: C, Al, Si, Ge, रुटाइल $TiO_2$, और $\beta$-$PtO_2$ के लिए, GI विधि छह-प्रणाली औसत MAE 0.148 meV/atom प्रदान करती है, जिसमें अधिकतम त्रुटि 0.522 meV/atom है। Si और Ge के लिए GI और QHA चरण आरेखों के बीच चरण-लेबल समानता अत्यंत उच्च (क्रमशः 99.50% और 99.98%) है।
• जटिल प्रणालियाँ ($Ta_2O_5$): नौ $Ta_2O_5$ पॉलीमॉर्फ्स के लिए, विधि प्रमुख मुक्त ऊर्जा रुझानों और निम्न-दबाव $\gamma$ चरण तथा उच्च-दबाव B चरण के बीच प्रतिस्पर्धा को पकड़ती है। हालांकि जटिल मोड क्रॉसिंग और सॉफ्ट मोड्स के कारण पूर्ण त्रुटियां अधिक हैं (~2 से 17 meV/atom), फिर भी चरण-सीमा टोपोलॉजी QHA के अनुरूप बनी रहती है।
• गति वृद्धि (Speedup): फोनन गणनाओं में कमी के परिणामस्वरूप महत्वपूर्ण गति वृद्धि देखी गई है, जो 5.911$\times$ (Al) से 9.023$\times$ (Si) तक है, जबकि $Ta_2O_5$ ने 8.103$\times$ की गति वृद्धि प्राप्त की।
• अनुप्रयोग सीमा: विधि मध्यम संपीड़न (लगभग $0.8V_0$ तक) के लिए विश्वसनीय है। इस सीमा के परे, सुचारू ग्रुनेसेन धारणा के टूटने, मोड क्रॉसिंग और आसन्न संरचनात्मक परिवर्तनों के कारण त्रुटियां तेजी से बढ़ती हैं।

महत्व
यह शोध पत्र दावा करता है कि यह GI कार्यप्रणाली प्रथम-सिद्धांत चरण-आरेख निर्माण के लिए एक व्यावहारिक रूप से कुशल उपकरण के रूप में कार्य करती है, विशेष रूप से उन सामग्रियों के लिए जहाँ पारंपरिक आयतन-निर्भर QHA गणनाएँ गणनात्मक रूप से कठिन होती हैं। यह AI-सहायता प्राप्त उम्मीदवार पीढ़ी और कठोर परिमित-तापमान स्थिरता सत्यापन के बीच के अंतर को पाटती है। लेखक इस बात पर जोर देते हैं कि यह विधि संरचना पीढ़ी या $0$ K विश्राम का स्थान नहीं लेती है, बल्कि उसके बाद के गिब्स मुक्त ऊर्जा चेक की लागत को काफी कम कर देती है। यह दृष्टिकोण उन गैर-चुंबकीय क्रिस्टलीय सामग्रियों के लिए सबसे उपयुक्त है जहाँ जालक कंपन मुक्त ऊर्जा (lattice vibrational free energies) प्रमुख परिमित-तापमान योगदान हैं और जहाँ फोनन आवृत्तियाँ आयतन के साथ सुचारू रूप से बदलती हैं। लेखक विनम्रतापूर्वक नोट करते हैं कि चुंबकीय एंट्रॉपी, मजबूत इलेक्ट्रॉनिक सहसंबंधों, या बड़ी विन्यासगत एंट्रॉपी (configurational entropy) द्वारा संचालित प्रणालियाँ वर्तमान दायरे से बाहर हैं।