🔬 materials science

Correlation-consistent Gaussian basis sets for copper solids from material-constrained atomic optimization

यह शोध पत्र कॉपर सॉलिड्स के लिए संख्यात्मक रूप से स्थिर, कोरिलेशन-कंसिस्टेंट गॉसियन बेस सेट्स उत्पन्न करने के लिए एक मटेरियल-कंस्ट्रेंड एटोमिक ऑप्टिमाइज़ेशन (MCAO) फ्रेमवर्क प्रस्तुत करता है, जो मानक एटोमली ऑप्टिमाइज़्ड सेट्स में निहित लीनियर डिपेंडेंस की समस्याओं को दूर करते हुए बल्क प्रॉपर्टीज और CO सोड्र्प्शन के लिए विश्वसनीय कम्प्लीट-बेसिस-सेट बेंचमार्क सक्षम बनाता है।

मूल लेखक: Jincheng Yu, Xiaoyu Zhang, Min-Ye Zhang, Yu Cao, Qiming Sun, Hong-Zhou Ye

प्रकाशित 2026-07-14

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

CC BY 4.0

मूल लेखक: Jincheng Yu, Xiaoyu Zhang, Min-Ye Zhang, Yu Cao, Qiming Sun, Hong-Zhou Ye

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

कल्पना कीजिए कि आप तांबे के तार का एक सटीक डिजिटल मॉडल बनाने की कोशिश कर रहे हैं। इसे करने के लिए, वैज्ञानिक "गौसियन बेस सेट" (Gaussian basis sets) नामक छोटे गणितीय निर्माण खंडों का उपयोग करते हैं। इन ब्लॉक्स को लेगो ब्रिक्स (Lego bricks) की तरह समझें। अंतरिक्ष में तैरते हुए एक अकेले तांबे के परमाणु के लिए, उसके इलेक्ट्रॉनों का सटीक वर्णन करने हेतु आपको बहुत बड़े, फूले हुए, विसरित (diffive) ब्रिक्स की आवश्यकता होगी। ये फूले हुए ब्रिक्स परमाणु के लिए तो बेहतरीन हैं, लेकिन धातु के एक ठोस ब्लॉक के लिए बहुत खराब हैं।

क्यों? क्योंकि जब आप एक ठोस धातु बनाने के लिए लाखों परमाणुओं को एक साथ पैक करते हैं, तो वे फूले हुए ब्रिक्स एक-दूसरे के ऊपर इतने अधिक ओवरलैप होने लगते हैं कि वे एक उलझे हुए, असंभव मलबे में बदल जाते हैं। यह एक मजबूत घर बनाने के लिए केवल विशाल, फूले हुए बादलों का उपयोग करने जैसा है; संरचना ढह जाती है क्योंकि हिस्से बहुत अधिक समान और नरम होते हैं। कंप्यूटर सिमुलेशन की दुनिया में, इस "उलझे हुए मलबे" को लीनियर डिपेंडेंस (linear dependence) कहा जाता है, और यह गणित को क्रैश कर देता है, जिससे तांबे के वास्तविक जीवन में व्यवहार (जैसे कि एक धातु के तार या सतह जहाँ गैस चिपकती है) का अध्ययन करना असंभव हो जाता है।

लंबे समय तक, वैज्ञानिकों को सटीकता (accuracy) और स्थिरता (stability) के बीच किसी एक को चुनना पड़ता था। या तो वे फूले हुए ब्रिक्स का उपयोग कर सकते थे (जो परमाणुओं के लिए सटीक थे, लेकिन सिमुलेशन क्रैश हो जाता था) या फूले हुए हिस्सों को काट देते थे (जो ठोसों के लिए स्थिर थे, लेकिन परिणाम अक्सर गलत होते थे)।

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

पुराना तरीका यह था कि कुत्ते की पूंछ काट दी जाए या उसे बांध दिया जाए (फूले हुए ब्रिक्स को काटना) ताकि वह लोगों से न टकराए। इससे पार्क सुरक्षित तो रहता था, लेकिन कुत्ता अब एक असली कुत्ते की तरह व्यवहार नहीं कर पाता था।

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

नए ब्रिक्स का परीक्षण
टीम ने इस नए तरीके को कॉपर (Cu) पर टेस्ट किया, जो एक ऐसी धातु है जिसे सिम्युलेट करना बहुत कठिन है क्योंकि इसमें "d" इलेक्ट्रॉन होते हैं (इन्हें धातु का "सीक्रेट सॉस" मान लें जो इसकी बिजली चालकता को नियंत्रित करता है)। उन्होंने नए ब्रिक्स बनाए, जिन्हें विभिन्न स्तरों के विवरण (डबल-जेटा, ट्रिपल-जेटा और क्वाड्रुपल-जेटा) के लिए MCAO-cc-pVXZ नाम दिया गया।

यहाँ उन्हें क्या मिला:

स्थिरता (Stability): जब उन्होंने इन नए ब्रिक्स का उपयोग तांबे के ब्लॉक को सिम्युलेट करने के लिए किया, तो गणित क्रैश नहीं हुआ। वह "उलझा हुआ मलबा" गायब हो गया था।
सटीकता (Accuracy): भले ही उन्होंने ब्रिक्स को अधिक सघन बनाया, फिर भी उन्होंने सटीकता नहीं खोई। जब उन्होंने तांबे के परमाणुओं के एक छोटे जोड़े (डिमर) को सिम्युलेट किया, तो ऊर्जा के परिणाम पुराने, सटीक परमाणु ब्रिक्स के लगभग समान थे (0.02 eV/atom के भीतर)।
वास्तविक दुनिया के गुण (Real-World Properties): उन्होंने तांबे के एक टुकड़े का सिमुलेशन किया और पाया कि परमाणुओं के बीच की दूरी (लैटिस कॉन्स्टेंट) और दबाव झेलने की क्षमता (बल्क मॉड्यूलस), "गोल्ड स्टैंडर्ड" के कंप्यूटर सिमुलेशन (प्लेन-वेव मेथड्स) के क्रमशः 0.01 Å और 1 GPa के भीतर थे।

"CO अधिशोषण पहेली" (CO Adsorption Puzzle) को सुलझाना
रसायन विज्ञान की सबसे बड़ी पहेलियों में से एक है "CO अधिशोषण पहेली"। यदि आप तांबे की सतह पर कार्बन मोनोऑक्साइड (CO) अणु रखते हैं, तो इसे तांबे के परमाणु के बिल्कुल शीर्ष (top site) पर चिपकना चाहिए। हालांकि, अधिकांश मानक कंप्यूटर प्रोग्राम अनुमान लगाते हैं कि यह तीन परमाणुओं के बीच के गड्ढे (hollow site) में चिपकता है।

उनके नए, स्थिर ब्रिक्स का उपयोग करके, शोधकर्ताओं ने यह देखने के लिए उच्च-स्तरीय सिमुलेशन (जिसे रैंडम-फेज एप्रोक्सिमेशन या RPA कहा जाता है) चलाए कि वास्तव में कौन सा स्थान पसंदीदा है।

उन्होंने पाया कि जब उन्होंने नए MCAO ब्रिक्स का उपयोग किया, तो सिमुलेशन ने सही ढंग से भविष्यवाणी की कि CO अणु टॉप साइट (top site) को पसंद करता है, जो प्रयोगों में देखी जाने वाली वास्तविकता से मेल खाता है।
उन्होंने यह भी पता लगाया कि कुछ पुराने तरीके, विशेष रूप से एक "लार्ज-कोर" स्यूडोपोटेंशियल (एक शॉर्टकट जो आंतरिक इलेक्ट्रॉनों को नकली कोर से बदल देता है), पूरी तरह से गलत थे। उस शॉर्टकट ने तांबे के परमाणुओं को बहुत मजबूती से चिपका दिया और भविष्यवाणी की कि CO गलत जगह पर चिपकेगा। नए तरीके ने दिखाया कि यह विशिष्ट शॉर्टकट ही अपराधी था, न कि स्वयं गौसियन ब्रिक्स का उपयोग करने का विचार।

इसका क्या अर्थ है
यह शोध पत्र सुझाव देता है कि आपको इन ब्रिक्स को ठोसों के लिए काम करने योग्य बनाने के लिए सटीक "परमाणु" तरीके को छोड़ने की आवश्यकता नहीं है। एक सरल "क्राउड-कंट्रोल" नियम जोड़कर, आप दोनों दुनिया का सर्वश्रेष्ठ प्राप्त कर सकते हैं: ऐसे ब्रिक्स जो जटिल रसायन विज्ञान के लिए पर्याप्त सटीक हैं और सुपरकंप्यूटर पर बिना क्रैश हुए चलने के लिए पर्याप्त स्थिर भी हैं।

लेखक आश्वस्त हैं कि ये नए MCAO-cc-pVXZ सेट तांबे और समान धातुओं के अध्ययन के लिए एक विश्वसनीय तरीका प्रदान करते हैं, जो धातु की सतहों पर गैस के अणुओं के चिपकने जैसी कठिन पहेलियों को हल करने का एक स्पष्ट मार्ग प्रदान करते हैं। वे यह भी नोट करते हैं कि हालांकि यह कॉपर के लिए बहुत अच्छा काम करता है, भविष्य के काम में यह देखना होगा कि क्या यह "क्राउड-कंट्रोल" ट्रिक अन्य प्रकार की धातुओं और सामग्रियों के लिए भी काम करती है।

तकनीकी सारांश: कॉपर सॉलिड्स के लिए मटेरियल-कंस्ट्रेंड एटोमिक ऑप्टिमाइजेशन से प्राप्त कोरिलेशन-कंसिस्टेंट गॉसियन बेस सेट्स

समस्या विवरण
कोरिलेशन-कंसिस्टेंट गॉसियन बेस सेट्स (जैसे डनिंग के cc-pVXZ सीरीज़) व्यवस्थित आणविक क्वांटम केमिस्ट्री के लिए मानक हैं, जो विश्वसनीय कम्प्लीट-बेसिस-सेट (CBS) एक्सट्रैपोलेशन को सक्षम करते हैं। हालाँकि, इनका आवधिक ठोस (periodic solids), विशेष रूप से धात्विक और धातु-युक्त प्रणालियों में सीधा अनुप्रयोग, लीनियर डिपेंडेंस (linear dependence) के कारण गंभीर रूप से सीमित है। यह समस्या तब उत्पन्न होती है जब परमाणु अनुकूलन (atomic optimization) अक्सर विसरित प्रिमिटिव गॉसियन फलनों (बहुत छोटे घातांकों, $\alpha \ll 0.1$ a.u. के साथ) को उत्पन्न करता है जो ठोस जाली (lattice) में परमाणुओं को एक साथ लाने पर लगभग रैखिक रूप से निर्भर हो जाते हैं। इससे ओवरलैप-मैट्रिक्स कंडीशन नंबर ( $\kappa$ ) में बड़ी वृद्धि होती है, जिससे संख्यात्मक अस्थिरता (numerical instabilities) जैसे कि कठिन सेल्फ-कंसिस्टेंट फील्ड (SCF) अभिसरण और गैर-सुचारू पोटेंशियल एनर्जी सर्फेस उत्पन्न होते हैं।

मौजूदा उपचार, जैसे कि छोटे घातांकों को ट्रंकेट करना या सिस्टम-विशिष्ट अनुकूलन करना, अक्सर बेस सेट्स की ट्रांसफरेबिलिटी (transferability) से समझौता करते हैं या कोरिलेशन-कंसिस्टेंट संरचना को नष्ट कर देते हैं जो उच्च-सटीक सह-संबंधित तरंग-फलन (correlated-wavefunction) गणनाओं के लिए आवश्यक है। हालांकि कुछ हालिया प्रयास (जैसे GTH-cc-pVXZ) संयोजकता-सेट (valence-set) के आकार को सीमित करके कोरिलेशन कंसिस्टेंसी को बनाए रखते हैं, लेकिन ऐसे डिस्क्रीट चयन प्रोटोकॉल की प्रकृति सटीकता और स्थिरता के बीच इष्टतम संतुलन पर प्रश्न उठाती है।

कार्यप्रणाली: मटेरियल-कंस्ट्रेंड एटोमिक ऑप्टिमाइजेशन (MCAO)
इन चुनौतियों को संबोधित करने के लिए, लेखक मटेरियल-कंस्ट्रेंड एटोमिक ऑप्टिमाइजेशन (MCAO) पेश करते हैं। यह ढांचा सामग्री-निर्भर स्थिरता बाधा (material-dependent stability constraint) के साथ मानक परमाणु बेस-सेट अनुकूलन को संवर्धित करता है। MCAO का मूल एक लॉस फंक्शन है जिसे परमाणु ऊर्जा को न्यूनतम करने और प्रतिनिधि ठोसों में बड़े कंडीशन नंबरों को दंडित (penalize) करने के लिए डिज़ाइन किया गया है:

$L(\alpha; \{M\}, \kappa_0, \varepsilon_0) = E_{\text{atom}}(\alpha) + \frac{\varepsilon_0}{\kappa_0} \times \max_{M \in \{M\}} \kappa(\alpha; M)$

जहाँ:

$E_{\text{atom}}(\alpha)$ मानक परमाणु ऊर्जा है (वैलेंस के लिए HF, पोलराइजेशन के लिए MP2)।
$\kappa(\alpha; M)$ एक संदर्भ ठोस $M$ के लिए ओवरलैप मैट्रिक्स का कंडीशन नंबर है (इस कार्य में, फेस-सेंटर्ड क्यूबिक Cu)।
$\kappa_0$ और $\varepsilon_0$ उपयोगकर्ता द्वारा चुने गए पैरामीटर हैं जो लक्षित कंडीशन नंबर और दंड के ऊर्जा पैमाने को नियंत्रित करते हैं।

एक्सपोनेंट ट्रंकेशन के विपरीत, MCAO डिस्क्रीट कटऑफ के बजाय निरंतर अनुकूलन के माध्यम से उपयुक्त कॉम्पैक्टनेस निर्धारित करता है। यह प्रक्रिया मानक परमाणु अनुकूलन प्रोटोकॉल को यथासंभव संरक्षित करती है, यह सुनिश्चित करती है कि परिणामी बेस सेट्स कोरिलेटेड गणनाओं के लिए आवश्यक कोरिलेशन-कंसिस्टेंट संरचना बनाए रखें।

मुख्य योगदान और कार्यान्वयन
एक प्रूफ-ऑफ-कांसेप्ट के रूप में, लेखकों ने कॉपर (Cu) के लिए डनिंग-शैली के कोरिलेशन-कंसिस्टेंट बेस सेट्स उत्पन्न करने के लिए MCAO को लागू किया, जिन्हें MCAO-cc-pVXZ ( $X = \text{D, T, Q}$ ) कहा जाता है। यह अध्ययन व्यापक रेंज के न्यूक्लियर-पोटेंशियल और रिलेटिविस्टिक उपचारों को कवर करता है:

ऑल-इलेक्ट्रॉन (AE): नॉन-रिलेटिविस्टिक (NR) और स्केलर-रिलेटिवस्टिक (SFX2C-1e)।
इफेक्टिव कोर पोटेंशियल्स (ECP): हार्ड और सॉफ्ट ccECPs।
स्यूडोपोटेंशियल: स्मॉल-कोर (sc) और लार्ज-कोर (lc) GTH-PBE स्यूडोपोटेंशियल।

अनुकूलन [Ar] $3d^{10}4s^1$ कॉन्फ़िगरेशन को लक्षित करता है। वैलेंस सेट को HF स्तर पर अनुकूलित किया जाता है, और पोलराइजेशन सेट को MP2 स्तर पर, जो 11 वैलेंस इलेक्ट्रॉनों को सह-संबंधित (correlate) करता है। कंट्रैक्शन स्कीम गोलाकार औसत HF एटोमिक नेचुरल ऑर्बिटल्स (ANOs) का उपयोग करती है, जिसमें DZ, TZ, और QZ सेट्स बनाने के लिए आंशिक डीकन्ट्रैक्शन (partial decontraction) किया गया है।

परिणाम
MCAO-cc-pVXZ बेस सेट्स संख्यात्मक स्थिरता और सटीकता के बीच एक सफल संतुलन प्रदर्शित करते हैं:

संख्यात्मक स्थिरता: मानक परमाणु रूप से अनुकूलित सेट्स (जैसे cc-pVXZ, def2-XZVP) की तुलना में fcc Cu के लिए कंडीशन नंबर ( $\kappa$ ) काफी कम हो जाता है, जो उन्हें $10^9$ की अस्थिरता सीमा से नीचे लाता है। ये बेस सेट्स बिना किसी एड-हॉक ट्रंकेशन के Cu सॉलिड्स और सतहों के लिए स्थिर रहते हैं।
आणविक सटीकता: Cu डाइमर के लिए, MCAO सेट्स सभी ज़ेटा स्तरों और CBS लिमिट में मानक cc-pVXZ सेट्स के 0.02 eV/atom के भीतर बाइंडिंग ऊर्जा को पुनरुत्पादित करते हैं। इसके विपरीत, सरल एक्सपोनेंट ट्रंकेशन (जैसे संशोधित def2-TZVP) द्वारा निर्मित बेस सेट्स बाइंडिंग ऊर्जा को बढ़ा-चढ़ाकर बताते हैं, जो एक कम संतुलित विवरण का संकेत देता है।
सॉलिड-स्टेट गुण: बल्क Cu के लिए, MCAO सेट्स लैट्टिस कॉन्स्टेंट ( $a_0$ ) और बल्क मॉड्यूलस ( $B_0$ ) को प्लेन-वेव (PW) संदर्भ मानों (क्रमशः 0.01 Å और 1 GPa के भीतर) तक तेजी से अभिसरित करते हैं। बैंड स्ट्रक्चर, जिसमें सह-संबंधित गणनाओं के लिए महत्वपूर्ण लो-एनर्जी कंडक्शन बैंड शामिल हैं, सटीक रूप से पुनरुत्पादित होते हैं।
बल्क Cu में त्रुटि मूल्यांकन: CBS-एक्सट्रैपोलेटेड रैंडम फेज एप्रोक्सिमेशन (RPA@PBE) गणनाओं का उपयोग करके, लेखकों ने विभिन्न सन्निकटन (approximations) से होने वाली त्रुटियों का आकलन किया। GTH-PBE-sc स्यूडोपोटेंशियल और सॉफ्ट ccECPs, SFX2C-1e ऑल-इलेक्ट्रॉन संदर्भ के साथ अच्छी तरह से मेल खाते हैं। हालाँकि, GTH-PBE-lc स्यूडोपोटेंशियल ने इंटरैक्शन स्ट्रेंथ को काफी अधिक बताया ( $a_0$ को ~0.1 Å कम करके और $B_0$ को ~70 GPa अधिक करके), जो सह-संबंधित विधियों के लिए स्यूडोपोटेंशियल के पैरामीटराइजेशन की चुनौतियों को उजागर करता है।
Cu(111) पर CO अधिशोषण (Adsorption): इन बेस सेट्स को "CO अधिशोषण पहेली" पर लागू किया गया था, जहाँ सेमीलोकल DFT गलत तरीके से प्रायोगिक टॉप साइट के ऊपर हॉलो साइट वरीयता की भविष्यवाणी करता है। CBS-एक्सट्रैपोलेटेड RPA@PBE गणनाओं ने टॉप-साइट वरीयता की पुष्टि की। अध्ययन ने उच्च-गुणवत्ता वाले स्केलर-रिलेटिवस्टिक ऑल-इलेक्ट्रॉन बेंचमार्क प्रदान किए, जिससे पता चला कि अधिकांश स्यूडोपोटेंशियल और ECPs ऑल-इलेक्ट्रॉन संदर्भ के 0.02 eV के भीतर साइट-वरीयता ऊर्जा को पुनरुत्पादित करते हैं, जबकि लार्ज-कोर GTH-PBE-lc काफी विफल रहा।

महत्व और दावे
लेख दावा करते हैं कि MCAO सॉलिड्स के लिए गॉसियन बेस सेट्स उत्पन्न करने का एक मजबूत, व्यवस्थित मार्ग प्रदान करता है जो परमाणु अनुकूलन की सटीकता और कोरिलेशन-कंसिस्टेंट संरचना को बनाए रखता है और साथ ही लीनियर डिपेंडेंस को नियंत्रित करता है। लेखक तर्क देते हैं कि:

बलिदान के बिना स्थिरता: यह संभव है कि कोरिलेशन-कंसिस्टेंट फ्रेमवर्क को छोड़े बिना या ऐसी सिस्टम-विशिष्ट अनुकूलन का सहारा लिए बिना जो ट्रांसफरेबिलिटी से समझौता करती है, ठोसों में संख्यात्मक स्थिरता प्राप्त की जाए।
बेंचमार्किंग क्षमता: इन बेस सेट्स द्वारा सक्षम विश्वसनीय CBS एक्सट्रैपोलेशन, सह-संबंधित तरंग-फलन गणनाओं में स्यूडोपोटेंशियल, ECPs और रिलेटिविस्टिक सन्निकटन से उत्पन्न त्रुटियों के नियंत्रित अलगाव और मूल्यांकन की अनुमति देता है।
CO पहेली का समाधान: यह कार्य Cu(111) पर CO अधिशोषण के लिए उच्च-गुणवत्ता वाले स्केलर-रिलेटिवस्टिक ऑल-इलेक्ट्रॉन बेंचमार्क प्रदान करता है, जो ट्रांजिशन-मेटल सरफेस केमिस्ट्री के लिए RPA@PBE की वैधता का समर्थन करता है।

लेखक सीमाओं को स्वीकार करते हुए कहते हैं कि वर्तमान फॉर्मूलेशन कंडीशन नंबर (एक एकल मीट्रिक) पर निर्भर करता है और केवल प्रिमिटिव एक्सपोनेंट्स को अनुकूलित करता है। वे सुझाव देते हैं कि भविष्य के कार्य कंट्रैक्शन कोफिशिएंट्स को अनुकूलित करने के लिए विस्तारित हो सकते हैं और कई इल-कंडीशन्ड मोड्स को संभालने के लिए रेगुलाइजेशन को रिफाइन कर सकते हैं। फिर भी, वे तर्क देते हैं कि MCAO रणनीति भविष्य के सॉलिड-स्टेट गॉसियन बेस-सेट अनुकूलन के लिए एक मौलिक दृष्टिकोण प्रदान करती है।

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